Changeset 6904 for branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/NEMO/OPA_SRC/USR/usrdef_zgr.F90

Timestamp:

2016-09-01T12:17:31+02:00 (8 years ago)

Author:

gm

Message:

#1692 - branch SIMPLIF_2_usrdef: OVERFLOW configuration (zco & sco) + small bug corrections

File:

• branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/NEMO/OPA_SRC/USR/usrdef_zgr.F90 (modified) (5 diffs)

branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/NEMO/OPA_SRC/USR/usrdef_zgr.F90

@@ -82 +82 @@
       CALL zgr_z  ( pdept_1d, pdepw_1d, pe3t_1d , pe3w_1d )  ! Reference z-coordinate system
       !
-      CALL zgr_top_bot( k_top , k_bot )                      ! top and bottom ocean level indices
+      CALL zgr_msk_top_bot( k_top , k_bot )                  ! masked top and bottom ocean t-level indices
       !
       !                                                      ! z-coordinate (3D arrays) from the 1D z-coord.
@@ -98 +98 @@
       !!
       !! ** Purpose :   set the depth of model levels and the resulting 
-      !!      vertical scale factors.
+      !!              vertical scale factors.
       !!
       !! ** Method  :   z-coordinate system (use in all type of coordinate)
-      !!      The depth of model levels is defined from an analytical
-      !!      function the derivative of which gives the scale factors.
-      !!      both depth and scale factors only depend on k (1d arrays).
-      !!              w-level: pdepw_1d  = pdep(k)
-      !!                       pe3w_1d(k) = dk(pdep)(k)     = e3(k)
-      !!              t-level: pdept_1d  = pdep(k+0.5)
-      !!                       pe3t_1d(k) = dk(pdep)(k+0.5) = e3(k+0.5)
+      !!              The depth of model levels is defined from an analytical
+      !!              function the derivative of which gives the scale factors.
+      !!              both depth and scale factors only depend on k (1d arrays).
+      !!                 w-level: pdepw_1d  = pdep(k)
+      !!                          pe3w_1d(k) = dk(pdep)(k)     = e3(k)
+      !!                 t-level: pdept_1d  = pdep(k+0.5)
+      !!                          pe3t_1d(k) = dk(pdep)(k+0.5) = e3(k+0.5)
       !!
       !!      Here the Madec & Imbard (1996) function is used
@@ -193 +193 @@
 
 
-   SUBROUTINE zgr_top_bot( k_top , k_bot )
-      !!----------------------------------------------------------------------
-      !!                    ***  ROUTINE zgr_top_bot  ***
-      !!
-      !! ** Purpose :   set the top and bottom ocean levels
+   SUBROUTINE zgr_msk_top_bot( k_top , k_bot )
+      !!----------------------------------------------------------------------
+      !!                    ***  ROUTINE zgr_msk_top_bot  ***
+      !!
+      !! ** Purpose :   set the masked top and bottom ocean t-levels
       !!
       !! ** Method  :   GYRE case = closed flat box ocean without ocean cavities
@@ -208 +208 @@
       INTEGER , DIMENSION(:,:), INTENT(out) ::   k_top , k_bot   ! first & last ocean level
       !
-      INTEGER  ::   ji, jj               ! dummy loop indices
-      INTEGER  ::   ii0, ii1, ij0, ij1   ! local indices
       REAL(wp), DIMENSION(jpi,jpj) ::   z2d   ! 2D local workspace
       !!----------------------------------------------------------------------
@@ -220 +218 @@
       z2d(:,:) = REAL( jpkm1 , wp )          ! flat bottom
       !
-      CALL lbc_lnk( z2d, 'T', 1. )           ! applied the lateral boundary condition (here jperio=0 ==>> closed)
+      CALL lbc_lnk( z2d, 'T', 1. )           ! set surrounding land to zero (here jperio=0 ==>> closed)
       !
       k_bot(:,:) = INT( z2d(:,:) )           ! =jpkm1 over the ocean point, =0 elsewhere
       !
-      k_top(:,:) = MIN( 1 , k_bot(:,:) )     ! =1     over the ocean point, =0 elsewhere
-      !
-   END SUBROUTINE zgr_top_bot
+      k_top(:,:) = MIN( 1 , k_bot(:,:) )     ! = 1    over the ocean point, =0 elsewhere
+      !
+   END SUBROUTINE zgr_msk_top_bot