Changeset 12918 for branches/UKMO/AMM15_v3_6_STABLE_package_collate_PS44/NEMOGCM/NEMO/OPA_SRC/ZDF

Timestamp:

2020-05-13T14:39:24+02:00 (4 years ago)

Author:

petesykes

Message:

Land masking for operational diagnostics

File:

• branches/UKMO/AMM15_v3_6_STABLE_package_collate_PS44/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl.F90 (modified) (4 diffs)

branches/UKMO/AMM15_v3_6_STABLE_package_collate_PS44/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl.F90

@@ -501 +501 @@
 
       REAL(wp) :: zdt             ! timestep variable  
+      REAL(wp) :: zmdi
+      REAL(wp), DIMENSION(jpi,jpj    ) :: zw2d
 
       CHARACTER(len=1) :: cmld
@@ -510 +512 @@
 
       !!----------------------------------------------------------------------
-      
+      zmdi=1.e+20                               !  missing data indicator for masking
+
       IF( kt == nit000 ) THEN
          REWIND( numnam_ref )              ! Namelist namzdf_mldzint in reference namelist 
@@ -550 +553 @@
 
                IF( iom_use( "mldzint_"//cmld ) ) THEN
-                  CALL iom_put( "mldzint_"//cmld, hmld_zint(:,:) )
+                  zw2d(:,:) = hmld_zint(:,:)*tmask(:,:,1) + zmdi*(1.0-tmask(:,:,1))
+                  CALL iom_put( "mldzint_"//cmld, zw2d(:,:) )
                ENDIF
 
@@ -582 +586 @@
                         WRITE(numout,*) 'zdf_mxl_zint (25h) : Outputting 25h data at i_cnt_25h=',i_cnt_25h
                      ENDIF
-                     CALL iom_put( "mldzint25h_"//cmld , hmld_zint_25h(:,:,jn) / 25._wp   )
+                     zw2d(:,:) = (hmld_zint_25h(:,:,jn)/25._wp)*tmask(:,:,1) + zmdi*(1.0-tmask(:,:,1))
+                     CALL iom_put( "mldzint25h_"//cmld , zw2d(:,:) )
                      ! Reset array
                      hmld_zint_25h(:,:,jn) = hmld_zint(:,:)