Changeset 15260 for NEMO/branches/NERC/dev_r4.0.6_GEOMETRIC/tests

Timestamp:

2021-09-14T18:31:24+02:00 (3 years ago)

Author:

acc

Message:

Untested, initial port of GEOMETRIC changes. Includes a bug fix to eken diagnostic in diawri.F90 which, if confirmed, is also relevant to current 4.0-HEAD. #2722

File:

• NEMO/branches/NERC/dev_r4.0.6_GEOMETRIC/tests/CANAL/MY_SRC/diawri.F90 (modified) (1 diff)

NEMO/branches/NERC/dev_r4.0.6_GEOMETRIC/tests/CANAL/MY_SRC/diawri.F90

@@ -212 +212 @@
       ENDIF
 
-      CALL iom_put( "woce", wn )                   ! vertical velocity
+      !                                            ! vertical velocity
+      IF( ln_zad_Aimp ) THEN   ;   CALL iom_put( "woce", wn + wi )   ! explicit plus implicit parts
+      ELSE                     ;   CALL iom_put( "woce", wn )
+      ENDIF
+
       IF( iom_use('w_masstr') .OR. iom_use('w_masstr2') ) THEN   ! vertical mass transport & its square value
-         ! Caution: in the VVL case, it only correponds to the baroclinic mass transport.
-         z2d(:,:) = rau0 * e1e2t(:,:)
+         !                     ! Caution: in the VVL case, it only correponds to the baroclinic mass transport.
          DO jk = 1, jpk
-            z3d(:,:,jk) = wn(:,:,jk) * z2d(:,:)
+            IF( ln_zad_Aimp ) THEN
+               z3d(:,:,jk) = rau0 * e1e2t(:,:) * ( wn(:,:,jk) + wi(:,:,jk) )
+            ELSE
+               z3d(:,:,jk) = rau0 * e1e2t(:,:) * wn(:,:,jk)
+            ENDIF
          END DO
          CALL iom_put( "w_masstr" , z3d )  
-         IF( iom_use('w_masstr2') )   CALL iom_put( "w_masstr2", z3d(:,:,:) * z3d(:,:,:) )
+         IF( iom_use('w_masstr2') )   CALL iom_put( "w_masstr2", z3d * z3d )
       ENDIF