!!---------------------------------------------------------------------
   !!                     ***  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)         &
#if !defined key_agrif
      , PARAMETER             &
#endif
      ::                      &  
      cp_cfg  =  "gyre"          !: Name of the configuration
   INTEGER                    &
#if !defined key_agrif
      , PARAMETER             &
#endif
      ::                      &
      jp_cfg  =  1     ,      &  !: Resolution of the configuration (degrees)

      ! Data domain size         !!! *  Size of all input files  *
      jpidta  =  30*jp_cfg+2, &  !: 1st lateral dimension ( >= jpiglo )
      jpjdta  =  20*jp_cfg+2, &  !: 2nd lateral dimension ( >= jpjglo )
      jpkdta  =  31         , &  !: Number of levels      ( >= jpk    )

#if defined key_c1d
      ! Zoom domain size         !!! *  C1D zoom  *
      jpiglo  =  3     ,      &  !: 1st dimension of global domain --> i
      jpjglo  =  3     ,      &  !: 2nd dimension of global domain --> j
      ! Domain characteristics
      jperio  =  0               !: Lateral cond. type (between 0 and 6)

   INTEGER                    &
      ::                      &
      ! Starting position of the zoom
      jpizoom =  1     ,      &  !: Left bottom (i,j) indices of the zoom
      jpjzoom =  1               !: in data domain indices
#else
      ! Global domain size       !!! *  Global domain  *
      jpiglo  =  jpidta,      &  !: 1st dimension of global domain --> i
      jpjglo  =  jpjdta,      &  !: 2nd dimension of global domain --> j
      ! Starting position of the zoom
      jpizoom =  1     ,      &  !: Left bottom (i,j) indices of the zoom
      jpjzoom =  1     ,      &  !: in data domain indices
      ! Domain characteristics
      jperio  =  0               !: Lateral cond. type (between 0 and 6)
#endif


   !!  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

   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 higher zacr, the smaller
      !                                           !: the stretching.

      !  If both ppa0 ppa1 and ppsur are specified to pp_to_be_computed, then
      !  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr

   REAL(wp), PARAMETER ::                      &
      &     ppdzmin = pp_not_used           ,  &  !: (meters) vertical thickness of the top layer
      &     pphmax  = pp_not_used                 !: (meters) Maximum depth of the ocean gdepw(jpk)
   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                 !:
   !!---------------------------------------------------------------------