Changeset 15293

Timestamp:

2021-09-27T17:43:00+02:00 (2 years ago)

Author:

clem

Message:

try to pass debug mode with outputs activated. One issue was that the shear was not calculated at the top and bottom levels. Another issue is that the shear is not defined in the haloes but solving this one probably requires to wait for the new xios to be released

Location:

NEMO/trunk/src/OCE/ZDF

Files:

• zdfphy.F90 (modified) (1 diff)
• zdfsh2.F90 (modified) (1 diff)

NEMO/trunk/src/OCE/ZDF/zdfphy.F90

@@ -381 +381 @@
       IF( iom_use('avt_k') .OR. iom_use('avm_k') .OR. iom_use('eshear_k') .OR. iom_use('estrat_k') ) THEN
          IF( l_zdfsh2 ) THEN
-            DO_2D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1 )
-               zsh2(ji,jj,1  ) = 0._wp
-               zsh2(ji,jj,jpk) = 0._wp
-            END_2D
             CALL iom_put( 'avt_k'   ,   avt_k       * wmask )
             CALL iom_put( 'avm_k'   ,   avm_k       * wmask )

NEMO/trunk/src/OCE/ZDF/zdfsh2.F90

@@ -96 +96 @@
          END_2D
       END DO
+      DO_2D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1 ) ! set p_sh2 to 0 at the surface and bottom for output purpose
+         p_sh2(ji,jj,1)   = 0._wp
+         p_sh2(ji,jj,jpk) = 0._wp
+      END_2D
       !
    END SUBROUTINE zdf_sh2