Changeset 12918 for branches/UKMO/AMM15_v3_6_STABLE_package_collate_PS44/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

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/DIA/diawri.F90 (modified) (6 diffs)

branches/UKMO/AMM15_v3_6_STABLE_package_collate_PS44/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

@@ -138 +138 @@
       REAL(wp), POINTER, DIMENSION(:,:)   :: z2d      ! 2D workspace
       REAL(wp), POINTER, DIMENSION(:,:,:) :: z3d      ! 3D workspace
+      REAL(wp), DIMENSION(jpi,jpj)        :: zw2d 
+      REAL(wp) :: zmdi
       REAL(wp), POINTER, DIMENSION(:,:,:) :: zrhd , zrhop  ! 3D workspace
       !!----------------------------------------------------------------------
@@ -147 +149 @@
       CALL wrk_alloc( jpi , jpj, jpk , zrhd , zrhop )
       !
+      zmdi=1.e+20                               !  missing data indicator for masking
+
       ! Output the initial state and forcings
       IF( ninist == 1 ) THEN                       
@@ -160 +164 @@
       ENDIF
 
-      CALL iom_put( "ssh" , sshn )                 ! sea surface height
+      zw2d(:,:)=sshn(:,:)*tmask(:,:,1) + zmdi*(1.0-tmask(:,:,1))
+      CALL iom_put( "ssh" , zw2d )                 ! sea surface height
       if( iom_use('ssh2') )   CALL iom_put( "ssh2", sshn(:,:) * sshn(:,:) )   ! square of sea surface height
       
@@ -169 +174 @@
             DO ji = 1, jpi
                jkbot = mbkt(ji,jj)
-               z2d(ji,jj) = tsn(ji,jj,jkbot,jp_tem)
+               z2d(ji,jj) = tsn(ji,jj,jkbot,jp_tem)*tmask(ji,jj,jkbot) + zmdi*(1.0-tmask(ji,jj,jkbot))
             END DO
          END DO
@@ -204 +209 @@
          
       CALL iom_put( "uoce", un(:,:,:)         )    ! 3D i-current
-      CALL iom_put(  "ssu", un(:,:,1)         )    ! surface i-current
+      zw2d(:,:)=un(:,:,1)*umask(:,:,1) + zmdi*(1.0-umask(:,:,1))
+      CALL iom_put(  "ssu", zw2d(:,:)         )    ! surface i-current
       IF ( iom_use("sbu") ) THEN
          DO jj = 1, jpj
@@ -221 +227 @@
       
       CALL iom_put( "voce", vn(:,:,:)         )    ! 3D j-current
-      CALL iom_put(  "ssv", vn(:,:,1)         )    ! surface j-current
+      zw2d(:,:)=vn(:,:,1)*vmask(:,:,1) + zmdi*(1.0-vmask(:,:,1))
+      CALL iom_put(  "ssv", zw2d(:,:)         )    ! surface j-current
       IF ( iom_use("sbv") ) THEN
          DO jj = 1, jpj