Changeset 15556 for NEMO/trunk/src/OCE/DOM/domzgr.F90

Timestamp:

2021-11-29T16:23:06+01:00 (2 years ago)

Author:

jchanut

Message:

#2638: Add the possibility to read bottom levels at U/V/F points in the mesh file. Store fe3mask (i.e. fmask as it is prior updating it for lateral boundary conditions). All is this is only needed to ensure a correct update of parent grid variables with AgRIF. This also anticipates the possible use of coarsened meshes.

File:

• NEMO/trunk/src/OCE/DOM/domzgr.F90 (modified) (10 diffs)

NEMO/trunk/src/OCE/DOM/domzgr.F90

@@ -74 +74 @@
       INTEGER  ::   ikt, ikb            ! top/bot index
       INTEGER  ::   ioptio, ibat, ios   ! local integer
+      INTEGER  ::   is_mbkuvf           ! ==0 if mbku, mbkv, mbkf to be computed
       REAL(wp) ::   zrefdep             ! depth of the reference level (~10m)
       REAL(wp), DIMENSION(jpi,jpj  ) ::   zmsk
@@ -97 +98 @@
             &              e3t_0   , e3u_0   , e3v_0 , e3f_0    ,   &    ! vertical scale factors
             &              e3w_0   , e3uw_0  , e3vw_0           ,   &    ! vertical scale factors
-            &              k_top   , k_bot            )                  ! 1st & last ocean level
+            &              k_top   , k_bot                      ,   &    ! 1st & last ocean level
+            &              is_mbkuvf, mbku, mbkv, mbkf )                 ! U/V/F points bottom levels
             !
       ELSE                          !==  User defined configuration  ==!
          IF(lwp) WRITE(numout,*)
          IF(lwp) WRITE(numout,*) '          User defined vertical mesh (usr_def_zgr)'
+         is_mbkuvf = 0
          !
          CALL usr_def_zgr( ln_zco  , ln_zps  , ln_sco, ln_isfcav,   & 
@@ -177 +180 @@
 
       !                                ! top/bottom ocean level indices for t-, u- and v-points (f-point also for top)
-      CALL zgr_top_bot( k_top, k_bot )      ! with a minimum value set to 1
+      CALL zgr_top_bot( k_top, k_bot, is_mbkuvf )      ! with a minimum value set to 1
       !
       !                                ! ice shelf draft and bathymetry
@@ -220 +223 @@
       &                 pe3t  , pe3u  , pe3v   , pe3f ,            &   ! vertical scale factors
       &                 pe3w  , pe3uw , pe3vw         ,            &   !     -      -      -
-      &                 k_top  , k_bot    )                            ! top & bottom ocean level
+      &                 k_top  , k_bot  ,                          &   ! top & bottom ocean level
+      &                 k_mbkuvf  , k_bot_u  , k_bot_v  , k_bot_f  )   ! U/V/F points bottom levels
       !!---------------------------------------------------------------------
       !!              ***  ROUTINE zgr_read  ***
@@ -235 +239 @@
       REAL(wp), DIMENSION(:,:,:), INTENT(out) ::   pe3w , pe3uw, pe3vw         !    -       -      -
       INTEGER , DIMENSION(:,:)  , INTENT(out) ::   k_top , k_bot               ! first & last ocean level
+      INTEGER                   , INTENT(out) ::   k_mbkuvf                    ! ==1 if mbku, mbkv, mbkf are in file
+      INTEGER , DIMENSION(:,:)  , INTENT(out) ::   k_bot_u , k_bot_v, k_bot_f  ! bottom levels at U/V/F points
       !
       INTEGER  ::   ji,jj,jk     ! dummy loop index
@@ -322 +328 @@
       ENDIF
       !
+      IF( iom_varid( inum, 'mbku', ldstop = .FALSE. ) > 0 ) THEN
+         IF(lwp) WRITE(numout,*) '          mbku, mbkv & mbkf read in ', TRIM(cn_domcfg), ' file'
+         CALL iom_get( inum, jpdom_global, 'mbku', z2d )
+         k_bot_u(:,:) = NINT( z2d(:,:) )
+         CALL iom_get( inum, jpdom_global, 'mbkv', z2d )
+         k_bot_v(:,:) = NINT( z2d(:,:) )
+         CALL iom_get( inum, jpdom_global, 'mbkf', z2d )
+         k_bot_f(:,:) = NINT( z2d(:,:) )
+         k_mbkuvf = 1
+      ELSE
+         k_mbkuvf = 0
+      ENDIF
+      !
       ! reference depth for negative bathy (wetting and drying only)
       IF( ll_wd )  CALL iom_get( inum,  'rn_wd_ref_depth' , ssh_ref   )
@@ -330 +349 @@
 
 
-   SUBROUTINE zgr_top_bot( k_top, k_bot )
+   SUBROUTINE zgr_top_bot( k_top, k_bot, k_mbkuvf )
       !!----------------------------------------------------------------------
       !!                    ***  ROUTINE zgr_top_bot  ***
@@ -342 +361 @@
       !!                                     (min value = 1)
       !! ** Action  :   mbkt, mbku, mbkv :   vertical indices of the deeptest 
-      !!                                     ocean level at t-, u- & v-points
+      !!                mbkf                 ocean level at t-, u-, v- & f-points
       !!                                     (min value = 1 over land)
       !!----------------------------------------------------------------------
       INTEGER , DIMENSION(:,:), INTENT(in) ::   k_top, k_bot   ! top & bottom ocean level indices
+      INTEGER                 , INTENT(in) ::   k_mbkuvf       ! flag to recompute mbku, mbkv, mbkf
       !
       INTEGER ::   ji, jj   ! dummy loop indices
@@ -365 +385 @@
          mikv(ji,jj) = MAX(  mikt(ji  ,jj+1) , mikt(ji,jj)  )
          mikf(ji,jj) = MAX(  mikt(ji  ,jj+1) , mikt(ji,jj), mikt(ji+1,jj  ), mikt(ji+1,jj+1)  )
-         !
-         mbku(ji,jj) = MIN(  mbkt(ji+1,jj  ) , mbkt(ji,jj)  )
-         mbkv(ji,jj) = MIN(  mbkt(ji  ,jj+1) , mbkt(ji,jj)  )
-      END_2D
-      ! converte into REAL to use lbc_lnk ; impose a min value of 1 as a zero can be set in lbclnk 
+      END_2D
+
+      IF ( k_mbkuvf==0 ) THEN
+         IF(lwp) WRITE(numout,*) '         mbku, mbkv, mbkf computed from mbkt'
+         DO_2D( 0, 0, 0, 0 )
+            mbku(ji,jj) = MIN(  mbkt(ji+1,jj  ) , mbkt(ji,jj)  )
+            mbkv(ji,jj) = MIN(  mbkt(ji  ,jj+1) , mbkt(ji,jj)  )
+            mbkf(ji,jj) = MIN(  mbkt(ji  ,jj+1) , mbkt(ji,jj), mbkt(ji+1,jj  ), mbkt(ji+1,jj+1)  )
+         END_2D
+      ELSE
+         IF(lwp) WRITE(numout,*) '         mbku, mbkv, mbkf read from file'
+         ! Use mbku, mbkv, mbkf from file
+         ! Ensure these are lower than expected bottom level deduced from mbkt
+         DO_2D( 0, 0, 0, 0 )
+            mbku(ji,jj) = MIN(  mbku(ji,jj), mbkt(ji+1,jj  ) , mbkt(ji,jj)  )
+            mbkv(ji,jj) = MIN(  mbkv(ji,jj), mbkt(ji  ,jj+1) , mbkt(ji,jj)  )
+            mbkf(ji,jj) = MIN(  mbkf(ji,jj), mbkt(ji  ,jj+1) , mbkt(ji,jj), mbkt(ji+1,jj  ), mbkt(ji+1,jj+1)  )
+         END_2D
+      ENDIF
+      ! convert into REAL to use lbc_lnk ; impose a min value of 1 as a zero can be set in lbclnk 
       DO_2D( 0, 0, 0, 0 )
          zk(ji,jj) = REAL( miku(ji,jj), wp )
@@ -399 +434 @@
       CALL lbc_lnk( 'domzgr', zk, 'V', 1.0_wp )
       mbkv(:,:) = MAX( NINT( zk(:,:) ), 1 )
+
+      DO_2D( 0, 0, 0, 0 )
+         zk(ji,jj) = REAL( mbkf(ji,jj), wp )
+      END_2D
+      CALL lbc_lnk( 'domzgr', zk, 'F', 1.0_wp )
+      mbkf(:,:) = MAX( NINT( zk(:,:) ), 1 )
       !
    END SUBROUTINE zgr_top_bot