Ignore:
Timestamp:
2016-09-01T12:17:31+02:00 (4 years ago)
Author:
gm
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/NEMO/OPA_SRC/USR/usrdef_zgr.F90

    r6717 r6904  
    8282      CALL zgr_z  ( pdept_1d, pdepw_1d, pe3t_1d , pe3w_1d )  ! Reference z-coordinate system 
    8383      ! 
    84       CALL zgr_top_bot( k_top , k_bot )                      ! top and bottom ocean level indices 
     84      CALL zgr_msk_top_bot( k_top , k_bot )                  ! masked top and bottom ocean t-level indices 
    8585      ! 
    8686      !                                                      ! z-coordinate (3D arrays) from the 1D z-coord. 
     
    9898      !! 
    9999      !! ** Purpose :   set the depth of model levels and the resulting  
    100       !!      vertical scale factors. 
     100      !!              vertical scale factors. 
    101101      !! 
    102102      !! ** Method  :   z-coordinate system (use in all type of coordinate) 
    103       !!      The depth of model levels is defined from an analytical 
    104       !!      function the derivative of which gives the scale factors. 
    105       !!      both depth and scale factors only depend on k (1d arrays). 
    106       !!              w-level: pdepw_1d  = pdep(k) 
    107       !!                       pe3w_1d(k) = dk(pdep)(k)     = e3(k) 
    108       !!              t-level: pdept_1d  = pdep(k+0.5) 
    109       !!                       pe3t_1d(k) = dk(pdep)(k+0.5) = e3(k+0.5) 
     103      !!              The depth of model levels is defined from an analytical 
     104      !!              function the derivative of which gives the scale factors. 
     105      !!              both depth and scale factors only depend on k (1d arrays). 
     106      !!                 w-level: pdepw_1d  = pdep(k) 
     107      !!                          pe3w_1d(k) = dk(pdep)(k)     = e3(k) 
     108      !!                 t-level: pdept_1d  = pdep(k+0.5) 
     109      !!                          pe3t_1d(k) = dk(pdep)(k+0.5) = e3(k+0.5) 
    110110      !! 
    111111      !!      Here the Madec & Imbard (1996) function is used 
     
    193193 
    194194 
    195    SUBROUTINE zgr_top_bot( k_top , k_bot ) 
    196       !!---------------------------------------------------------------------- 
    197       !!                    ***  ROUTINE zgr_top_bot  *** 
    198       !! 
    199       !! ** Purpose :   set the top and bottom ocean levels 
     195   SUBROUTINE zgr_msk_top_bot( k_top , k_bot ) 
     196      !!---------------------------------------------------------------------- 
     197      !!                    ***  ROUTINE zgr_msk_top_bot  *** 
     198      !! 
     199      !! ** Purpose :   set the masked top and bottom ocean t-levels 
    200200      !! 
    201201      !! ** Method  :   GYRE case = closed flat box ocean without ocean cavities 
     
    208208      INTEGER , DIMENSION(:,:), INTENT(out) ::   k_top , k_bot   ! first & last ocean level 
    209209      ! 
    210       INTEGER  ::   ji, jj               ! dummy loop indices 
    211       INTEGER  ::   ii0, ii1, ij0, ij1   ! local indices 
    212210      REAL(wp), DIMENSION(jpi,jpj) ::   z2d   ! 2D local workspace 
    213211      !!---------------------------------------------------------------------- 
     
    220218      z2d(:,:) = REAL( jpkm1 , wp )          ! flat bottom 
    221219      ! 
    222       CALL lbc_lnk( z2d, 'T', 1. )           ! applied the lateral boundary condition (here jperio=0 ==>> closed) 
     220      CALL lbc_lnk( z2d, 'T', 1. )           ! set surrounding land to zero (here jperio=0 ==>> closed) 
    223221      ! 
    224222      k_bot(:,:) = INT( z2d(:,:) )           ! =jpkm1 over the ocean point, =0 elsewhere 
    225223      ! 
    226       k_top(:,:) = MIN( 1 , k_bot(:,:) )     ! =1     over the ocean point, =0 elsewhere 
    227       ! 
    228    END SUBROUTINE zgr_top_bot 
     224      k_top(:,:) = MIN( 1 , k_bot(:,:) )     ! = 1    over the ocean point, =0 elsewhere 
     225      ! 
     226   END SUBROUTINE zgr_msk_top_bot 
    229227    
    230228 
Note: See TracChangeset for help on using the changeset viewer.