Changeset 10490

Timestamp:

2019-01-09T20:49:56+01:00 (4 years ago)

Author:

mathiot

Message:

change the deepest e3w_0 inside the ice shelf => correct depth when using e3_to_depth

File:

• utils/tools/DOMAINcfg/src/domzgr.f90 (modified) (1 diff)

utils/tools/DOMAINcfg/src/domzgr.f90

@@ -1827 +1827 @@
             !       ... on ik / ik-1 
                e3w_0  (ji,jj,ik  ) = e3t_0  (ji,jj,ik) !2._wp * (gdept_0(ji,jj,ik) - gdepw_0(ji,jj,ik)) 
+               gdept_0(ji,jj,ik-1) = gdept_0(ji,jj,ik) - e3w_0(ji,jj,ik)
                e3t_0  (ji,jj,ik-1) = gdepw_0(ji,jj,ik) - gdepw_1d(ik-1)
-! The next line isn't required and doesn't affect results - included for consistency with bathymetry code 
-               gdept_0(ji,jj,ik-1) = gdept_1d(ik-1)
+               e3w_0  (ji,jj,ik-1) = gdept_0(ji,jj,ik-1) - gdept_1d(ik-2)
+               gdepw_0(ji,jj,ik-1) = gdepw_0(ji,jj,ik) - e3t_0(ji,jj,ik-1)
             ENDIF 
          END DO