Changeset 11908 for NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfcpl.F90

Timestamp:

2019-11-14T20:20:07+01:00 (4 years ago)

Author:

mathiot

Message:

ENHANCE-02_ISF_nemo: fix issue initialisation after ice sheet coupling step

File:

• NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfcpl.F90 (modified) (3 diffs)

NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfcpl.F90

@@ -61 +61 @@
       ! check presence of variable needed for coupling
       ! iom_varid return 0 if not found
-      id = 0
+      id = 1
       id = id * iom_varid(numror, 'ssmask', ldstop = .false.)
       id = id * iom_varid(numror, 'tmask' , ldstop = .false.)
@@ -145 +145 @@
       !!----------------------------------------------------------------------
       !!
-      INTEGER :: ji, jj, jd          !! loop index
+      INTEGER :: ji, jj, jd, jk      !! loop index
       INTEGER :: jip1, jim1, jjp1, jjm1
       !!
@@ -195 +195 @@
       !
       ! recompute the vertical scale factor, depth and water thickness
+      IF(lwp) write(numout,*) 'isfcpl_ssh : recompute scale factor from sshn (new wet cell)'
+      IF(lwp) write(numout,*) '~~~~~~~~~~~'
+      DO jk = 1, jpk
+         e3t_n(:,:,jk) =  e3t_0(:,:,jk) * ( ht_0(:,:) + sshn(:,:) ) &
+             &                          / ( ht_0(:,:) + 1._wp - ssmask(:,:) ) * tmask(:,:,jk)   &
+             &          + e3t_0(:,:,jk)                               * (1._wp -tmask(:,:,jk))
+      END DO
+      e3t_b(:,:,:) = e3t_n(:,:,:)
       CALL dom_vvl_zgr()
       !