source: branches/2013/dev_r3940_CNRS4_IOCRS/NEMOGCM/NEMO/OPA_SRC/par_GYRE.h90 @ 4015

   !!---------------------------------------------------------------------
   !!                     ***  par_GYRE.h90  ***
   !!   Ocean Domain : GYRE configuration at 1/jp_cfg degree resolution 
   !!---------------------------------------------------------------------
   !!----------------------------------------------------------------------
   !! NEMO/OPA 3.3 , NEMO Consortium (2010)
   !! $Id$ 
   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
   !!----------------------------------------------------------------------
   CHARACTER (len=16) :: cp_cfg = "gyre"  !: name of the configuration
   INTEGER            :: jp_cfg = 1       !: resolution of the configuration 
   INTEGER, PARAMETER ::     &
      ! data size              !!! * size of all the input files *
      jpidta  = 30*jp_cfg+2, &  !: 1st horizontal dimension ( >= jpi )
      jpjdta  = 20*jp_cfg+2, &  !: 2nd    "            "    ( >= jpj )
      jpkdta  = 31              !: number of levels         ( >= jpk )

   INTEGER            ::     &
      ! global domain size     !!! * full domain *
      jpiglo  = jpidta,      &  !: 1st dimension of global domain --> i
      jpjglo  = jpjdta,      &  !: 2nd    "                  "    --> j
      ! zoom starting position
      jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom
      jpjzoom =   1   ,      &  !: in data indices

      ! Domain characteristics
      jperio  =     0           !: lateral cond. type (between 0 and 6)

   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config.
   !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr
   REAL(wp), PARAMETER ::   &  !:
      pp_not_used       = 999999._wp  , & !: ???
      pp_to_be_computed =      0._wp      !: ???
   !!
   !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 )
   !!
   INTEGER,PARAMETER   ::    & !:
      jphgr_msh = 5            !: type of horizontal mesh
      !                        ! = 0 curvilinear coordinate on the sphere
      !                        !     read in coordinate.nc file
      !                        ! = 1 geographical mesh on the sphere
      !                        !     with regular grid-spacing
      !                        ! = 2 f-plane with regular grid-spacing
      !                        ! = 3 beta-plane with regular grid-spacing
      !                        ! = 4 Mercator grid with T/U point at the equator  with
      !                        !     isotropic resolution (e1_deg)
      !                        ! =5  beta-plane with regular grid-spacing and rotated domain (GYRE configuration)

      !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid.
      !   The mercator grid starts only approximately at gphi0 because
      !   of the constraint that the equator be a T point.
   REAL(wp) &
#if !defined key_agrif
      , PARAMETER  &
#endif
      ::     &  !:
      ppglam0  =    0.0_wp,   &  !: longitude of first raw and column T-point (jphgr_msh = 1)
      ppgphi0  =   29.0_wp,   &  !: latitude  of first raw and column T-point (jphgr_msh = 1)
      !                          ! latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3)
      ppe1_deg = pp_not_used ,   &  !: zonal      grid-spacing (degrees)
      ppe2_deg = pp_not_used ,   &  !: meridional grid-spacing (degrees)
      !
      ppe1_m   = pp_not_used,    &  !: zonal      grid-spacing (meters )
      ppe2_m   = pp_not_used        !: meridional grid-spacing (meters )
   !!
   !!  Coefficients associated with the vertical coordinate system
   !!

   REAL(wp), PARAMETER  ::       &  !:
      &     ppsur = -2033.194295283385_wp   ,  &  !: 
      &     ppa0  =  155.8325369664153_wp   ,  &  !:
      &     ppa1  =  146.3615918601890_wp   ,  &  !:
      !
      &     ppkth =  17.28520372419791_wp   ,  &  !: (non dimensional): gives the approximate
      !                                           !    layer number above which  stretching will
      !                                           !    be maximum. Usually of order jpk/2.
      &     ppacr =  5.000000000000000_wp         !: (non dimensional): stretching factor
      !                                           !    for the grid. The highest zacr, the smallest
      !                                           !    the stretching.

   !!
   !!  If all ppa0 ppa1 and ppsur are specified to 0, then
   !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr
   !!
   REAL(wp), PARAMETER ::        &  !:
      &     ppdzmin = pp_not_used   ,  &  !: (meters): depth of the top (first) model layer
      !                             !            depth of second "w" level
      &     pphmax  = pp_not_used         !: (meters): maximum depth of the ocean
      !                             !            depth of the last "w" level
   LOGICAL,  PARAMETER ::        &
      &     ldbletanh = .FALSE.                   !: Use/do not use double tanf function for vertical coordinates
   REAL(wp), PARAMETER ::        &
      &     ppa2    = pp_not_used           ,  &  !: Double tanh function parameters
      &     ppkth2  = pp_not_used           ,  &  !:
      &     ppacr2  = pp_not_used                 !:

   !!---------------------------------------------------------------------