Changeset 9529 for branches/2018/dev_r8864_nemo_v3_6_ZTILDE/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.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/DIA/diawri.F90 (modified) (1 diff)

branches/2018/dev_r8864_nemo_v3_6_ZTILDE/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

@@ -157 +157 @@
          CALL iom_put( "e3tdef"  , ( ( fse3t_n(:,:,:) - e3t_0(:,:,:) ) / e3t_0(:,:,:) * 100 * tmask(:,:,:) ) ** 2 )
       CALL iom_put("tpt_dep", fsdept_n(:,:,:) )
-
-
+      IF( iom_use("gdepwt") )  THEN
+         z3d(:,:,1) = (gdepw_n(:,:,1)-sshn(:,:))*tmask(:,:,1)
+         DO jk=2,jpk
+            z3d(:,:,jk) = (gdepw_n(:,:,jk)-sshn(:,:))*tmask(:,:,jk-1)
+         END DO
+         CALL iom_put("gdepwt" , z3d(:,:,:) )
+      END IF
 
       CALL iom_put( "ssh" , sshn )                 ! sea surface height