Changeset 10164 for NEMO/branches/2018/dev_r10057_ENHANCE03_ZTILDE/src/OCE/DIA/diawri.F90

Timestamp:

2018-10-02T16:02:24+02:00 (6 years ago)

Author:

jchanut

Message:

ztilde update, #2126: Add higher order filtering + alternative scale factor decomposition

File:

• NEMO/branches/2018/dev_r10057_ENHANCE03_ZTILDE/src/OCE/DIA/diawri.F90 (modified) (3 diffs)

NEMO/branches/2018/dev_r10057_ENHANCE03_ZTILDE/src/OCE/DIA/diawri.F90

@@ -58 +58 @@
    USE diurnal_bulk    ! diurnal warm layer
    USE cool_skin       ! Cool skin
+   USE domvvl
 
    IMPLICIT NONE
@@ -124 +125 @@
 
       ! Output of initial vertical scale factor
-      CALL iom_put("e3t_0", e3t_0(:,:,:) )
+      CALL iom_put("e3t_0", e3t_0(:,:,:))
       CALL iom_put("e3u_0", e3u_0(:,:,:) )
       CALL iom_put("e3v_0", e3v_0(:,:,:) )
@@ -132 +133 @@
       CALL iom_put( "e3v" , e3v_n(:,:,:) )
       CALL iom_put( "e3w" , e3w_n(:,:,:) )
+      !
+      IF ( (ln_vvl_ztilde).OR.(ln_vvl_layer) ) THEN
+         IF( iom_use("e3t_star") ) THEN    ! Barotropic cell thickness anomaly
+            z3d(:,:,:) = (e3t_n(:,:,:)-tilde_e3t_n(:,:,:)-e3t_0(:,:,:))*tmask(:,:,:) 
+            CALL iom_put( "e3t_star" , z3d(:,:,:) )
+         ENDIF
+         IF( iom_use("e3t_tilde") )  THEN  ! Baroclinic cell thickness anomaly
+            CALL iom_put( "e3t_tilde" , tilde_e3t_n(:,:,:) )
+         ENDIF
+      ENDIF
+
       IF( iom_use("e3tdef") )   &
          CALL iom_put( "e3tdef"  , ( ( e3t_n(:,:,:) - e3t_0(:,:,:) ) / e3t_0(:,:,:) * 100 * tmask(:,:,:) ) ** 2 )