New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
par_GYRE.h90 in branches/2013/dev_r3987_UKMO6_C1D/NEMOGCM/NEMO/OPA_SRC – NEMO

source: branches/2013/dev_r3987_UKMO6_C1D/NEMOGCM/NEMO/OPA_SRC/par_GYRE.h90 @ 4144

Last change on this file since 4144 was 4144, checked in by rfurner, 11 years ago

Commit for 2013 changes; see #1085

  • Property svn:keywords set to Id
File size: 5.9 KB
RevLine 
[93]1   !!---------------------------------------------------------------------
2   !!                     ***  par_GYRE.h90  ***
3   !!   Ocean Domain : GYRE configuration at 1/jp_cfg degree resolution
4   !!---------------------------------------------------------------------
[247]5   !!----------------------------------------------------------------------
[2528]6   !! NEMO/OPA 3.3 , NEMO Consortium (2010)
[1152]7   !! $Id$
[2528]8   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
[247]9   !!----------------------------------------------------------------------
[4144]10   CHARACTER (len=16)         &
[392]11#if !defined key_agrif
[4144]12      , PARAMETER             &
[389]13#endif
[4144]14      ::                      & 
15      cp_cfg  =  "gyre"          !: Name of the configuration
16   INTEGER                    &
[392]17#if !defined key_agrif
[4144]18      , PARAMETER             &
[389]19#endif
[4144]20      ::                      &
21      jp_cfg  =  1     ,      &  !: Resolution of the configuration (degrees)
[93]22
[4144]23      ! Data domain size         !!! *  Size of all input files  *
24      jpidta  =  30*jp_cfg+2, &  !: 1st lateral dimension ( >= jpiglo )
25      jpjdta  =  20*jp_cfg+2, &  !: 2nd lateral dimension ( >= jpjglo )
26      jpkdta  =  31         , &  !: Number of levels      ( >= jpk    )
[93]27
[4144]28#if defined key_c1d
29      ! Zoom domain size         !!! *  C1D zoom  *
30      jpiglo  =  3     ,      &  !: 1st dimension of global domain --> i
31      jpjglo  =  3     ,      &  !: 2nd dimension of global domain --> j
32      ! Domain characteristics
33      jperio  =  0               !: Lateral cond. type (between 0 and 6)
[93]34
[4144]35   INTEGER                    &
36      ::                      &
37      ! Starting position of the zoom
38      jpizoom =  1     ,      &  !: Left bottom (i,j) indices of the zoom
39      jpjzoom =  1               !: in data domain indices
40#else
41      ! Global domain size       !!! *  Global domain  *
42      jpiglo  =  jpidta,      &  !: 1st dimension of global domain --> i
43      jpjglo  =  jpjdta,      &  !: 2nd dimension of global domain --> j
44      ! Starting position of the zoom
45      jpizoom =  1     ,      &  !: Left bottom (i,j) indices of the zoom
46      jpjzoom =  1     ,      &  !: in data domain indices
[93]47      ! Domain characteristics
[4144]48      jperio  =  0               !: Lateral cond. type (between 0 and 6)
49#endif
[93]50
[4144]51
[93]52   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config.
53   !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr
54
[4144]55   REAL(wp), PARAMETER ::              &
56      pp_not_used       = 999999._wp , &  !:
57      pp_to_be_computed = 0._wp           !:
58
59
60   !! Coefficients associated with the horizontal coordinate system
61
62   INTEGER, PARAMETER  ::     &
63      jphgr_msh = 5              !: type of horizontal mesh
64      !                          !: = 0 curvilinear coordinate on the sphere
65      !                          !:     read in coordinate.nc file
66      !                          !: = 1 geographical mesh on the sphere
67      !                          !:     with regular grid-spacing
68      !                          !: = 2 f-plane with regular grid-spacing
69      !                          !: = 3 beta-plane with regular grid-spacing
70      !                          !: = 4 Mercator grid with T/U point at the equator  with
71      !                          !:     isotropic resolution (e1_deg)
72      !                          !: = 5 beta-plane with regular grid-spacing and rotated domain (GYRE configuration)
73
[93]74      !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid.
75      !   The mercator grid starts only approximately at gphi0 because
76      !   of the constraint that the equator be a T point.
[4144]77
78   REAL(wp)                   &
[443]79#if !defined key_agrif
[4144]80      , PARAMETER             &
[443]81#endif
[4144]82      ::                      &
83      ppglam0  = 0.0_wp     , &  !: longitude of first raw and column T-point   (jphgr_msh = 1)
84      ppgphi0  = 29.0_wp    , &  !: latitude  of first raw and column T-point   (jphgr_msh = 1)
85      !                          !: latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3)
86      ppe1_deg = pp_not_used, &  !: zonal      grid-spacing (degrees)
87      ppe2_deg = pp_not_used, &  !: meridional grid-spacing (degrees)
[93]88      !
[4144]89      ppe1_m   = pp_not_used, &  !: zonal      grid-spacing (meters )
90      ppe2_m   = pp_not_used     !: meridional grid-spacing (meters )
91
92
[93]93   !!  Coefficients associated with the vertical coordinate system
94
[4144]95   REAL(wp), PARAMETER  ::                     &
96      &     ppsur =   -2033.194295283385_wp ,  &  !:
97      &     ppa0  =   155.8325369664153_wp  ,  &  !:
98      &     ppa1  =   146.3615918601890_wp  ,  &  !:
99      &     ppkth =   17.28520372419791_wp  ,  &  !: (non dimensional): gives the approximate
100      !                                           !: layer number above which  stretching will
101      !                                           !: be maximum. Usually of order jpk/2.
102      &     ppacr =   5.000000000000000_wp        !: (non dimensional): stretching factor
103      !                                           !: for the grid. The higher zacr, the smaller
104      !                                           !: the stretching.
[93]105
[4144]106      !  If both ppa0 ppa1 and ppsur are specified to pp_to_be_computed, then
107      !  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr
108
109   REAL(wp), PARAMETER ::                      &
110      &     ppdzmin = pp_not_used           ,  &  !: (meters) vertical thickness of the top layer
111      &     pphmax  = pp_not_used                 !: (meters) Maximum depth of the ocean gdepw(jpk)
112   LOGICAL,  PARAMETER ::                      &
[2528]113      &     ldbletanh = .FALSE.                   !: Use/do not use double tanf function for vertical coordinates
[4144]114   REAL(wp), PARAMETER ::                      &
[2528]115      &     ppa2    = pp_not_used           ,  &  !: Double tanh function parameters
116      &     ppkth2  = pp_not_used           ,  &  !:
117      &     ppacr2  = pp_not_used                 !:
[93]118   !!---------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.