Changeset 9529 for branches/2018/dev_r8864_nemo_v3_6_ZTILDE/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

Timestamp:

2018-04-30T16:38:28+02:00 (6 years ago)

Author:

jchanut

Message:

ztilde stability update

File:

• branches/2018/dev_r8864_nemo_v3_6_ZTILDE/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90 (modified) (3 diffs)

branches/2018/dev_r8864_nemo_v3_6_ZTILDE/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

@@ -198 +198 @@
       !!----------------------------------------------------------------------
       INTEGER  ::   jk                     ! dummy loop indices
+      LOGICAL  ::   ll_e3_dep
       REAL(wp) ::   zt, zw                 ! temporary scalars
       REAL(wp) ::   zsur, za0, za1, zkth   ! Values set from parameters in
@@ -209 +210 @@
       ! Set variables from parameters
       ! ------------------------------
+       ll_e3_dep = .TRUE.
        zkth = ppkth       ;   zacr = ppacr
        zdzmin = ppdzmin   ;   zhmax = pphmax
@@ -312 +314 @@
       ENDIF
 
-      IF ( ln_isfcav ) THEN
+      IF ( ln_isfcav .OR. ll_e3_dep ) THEN
 ! need to be like this to compute the pressure gradient with ISF. If not, level beneath the ISF are not aligned (sum(e3t) /= depth)
 ! define e3t_0 and e3w_0 as the differences between gdept and gdepw respectively