Changeset 12614 for NEMO/branches/2020/dev_r12527_Gurvan_ShallowWater/src/SWE/domain.F90

Timestamp:

2020-03-26T15:59:52+01:00 (4 years ago)

Author:

gm

Message:

first Shallow Water Eq. update

Location:

NEMO/branches/2020/dev_r12527_Gurvan_ShallowWater/src/SWE

Files:

• . (added)
• domain.F90 (copied) (copied from NEMO/branches/2020/dev_r12527_Gurvan_ShallowWater/src/OCE/DOM/domain.F90) (3 diffs)

NEMO/branches/2020/dev_r12527_Gurvan_ShallowWater/src/SWE/domain.F90

@@ -143 +143 @@
 
       CALL dom_msk( ik_top, ik_bot )    ! Masks
-      !
-      ht_0(:,:) = 0._wp  ! Reference ocean thickness
+
+      ht_0(:,:) = 0._wp                 ! Reference ocean thickness
       hu_0(:,:) = 0._wp
       hv_0(:,:) = 0._wp
-      DO jk = 1, jpk
+      hf_0(:,:) = 0._wp
+      DO jk = 1, jpkm1
          ht_0(:,:) = ht_0(:,:) + e3t_0(:,:,jk) * tmask(:,:,jk)
          hu_0(:,:) = hu_0(:,:) + e3u_0(:,:,jk) * umask(:,:,jk)
          hv_0(:,:) = hv_0(:,:) + e3v_0(:,:,jk) * vmask(:,:,jk)
+         hf_0(:,:) = hf_0(:,:) + e3f_0(:,:,jk) * ssfmask(:,:)     ! CAUTION : only valid in SWE, not with bathymetry
       END DO
-      !
+      !                                 ! Inverse of reference ocean thickness
+      r1_ht_0(:,:) =  ssmask(:,:) / ( ht_0(:,:) + 1._wp -  ssmask(:,:) )
+      r1_hu_0(:,:) = ssumask(:,:) / ( hu_0(:,:) + 1._wp - ssumask(:,:) )
+      r1_hv_0(:,:) = ssvmask(:,:) / ( hv_0(:,:) + 1._wp - ssvmask(:,:) )
+      r1_hf_0(:,:) = ssfmask(:,:) / ( hf_0(:,:) + 1._wp - ssfmask(:,:) )
+      
       !           !==  time varying part of coordinate system  ==!
       !
@@ -190 +197 @@
       !
       IF( ln_meshmask    )   CALL dom_wri       ! Create a domain file
+
       IF( .NOT.ln_rstart )   CALL dom_ctl       ! Domain control
       !
@@ -200 +208 @@
          WRITE(numout,*) 
       ENDIF
+      
       !
    END SUBROUTINE dom_init