Changeset 9707 for branches/UKMO/dev_r5518_GO6_diag_bitcomp/NEMOGCM/NEMO/OPA_SRC/IOM/iom.F90

Timestamp:

2018-05-31T11:48:34+02:00 (6 years ago)

Author:

frrh

Message:

Work out N-fold mask points for U and V grids based on
value flipping during N-fold lbc_lnk operations.
It's much easier this way.
Tidy up remaining code.

File:

• branches/UKMO/dev_r5518_GO6_diag_bitcomp/NEMOGCM/NEMO/OPA_SRC/IOM/iom.F90 (modified) (1 diff)

branches/UKMO/dev_r5518_GO6_diag_bitcomp/NEMOGCM/NEMO/OPA_SRC/IOM/iom.F90

@@ -1400 +1400 @@
       IF ( ln_mskland ) THEN
          ! mask land points, keep values on coast line -> specific mask for U, V and W points
+
+
          SELECT CASE ( cdgrd )
+         ! The masks applied here are specifically used to mask out duplicate 
+         ! data poits in wrap columns and N-fold rows in order to ensure bit 
+         ! comparison of diagnostics which have not undergone an explicit
+         ! lbc_lnk prior to writing and which are thus prone to junk in 
+         ! in these duplicate points since they are not usually explicitly 
+         ! involved in field computation. 
          CASE('T')   
-            zmask(:,:,:) = tmask_i_3d(:,:,:)
-         CASE('U')
-            zmask(:,:,:) = umask_i_3d(:,:,:) 
-         CASE('V') 
-            zmask(:,:,:) = vmask_i_3d(:,:,:) 
+            zmask(:,:,:) = tmask_i_diag(:,:,:)
+         CASE('U')      
+            zmask(:,:,:) = umask_i_diag(:,:,:)
+    CASE('V')   
+            zmask(:,:,:) = vmask_i_diag(:,:,:)
          CASE('W')   
-            zmask(:,:,2:jpk  ) = tmask_i_3d(:,:,1:jpkm1) + tmask_i_3d(:,:,2:jpk)   
-            zmask(:,:,1) = tmask_i_3d(:,:,1)
-         END SELECT
+            zmask(:,:,2:jpk  ) = tmask_i_diag(:,:,1:jpkm1) + tmask_i_diag(:,:,2:jpk)   
+            zmask(:,:,1) = tmask_i_diag(:,:,1)
+        END SELECT
          !
 #if ! defined key_xios2