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_ORCA_R2.h90 in branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC – NEMO

source: branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/par_ORCA_R2.h90 @ 2379

Last change on this file since 2379 was 2379, checked in by acc, 13 years ago

nemo_v3_3beta. First stage of adding ORCA1 support. See ticket #758. Added optional double tanh function for vertical grid

  • Property svn:keywords set to Id
File size: 8.0 KB
Line 
1   !!---------------------------------------------------------------------
2   !!                     ***  par_ORCA_R2.h90  *** 
3   !!   Ocean Domain : 2 degrees resolution global ocean
4   !!                  (0RCA_R2 configuration)
5   !!---------------------------------------------------------------------
6   !!----------------------------------------------------------------------
7   !! NEMO/OPA 3.3 , NEMO Consortium (2010)
8   !! $Id$
9   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
10   !!----------------------------------------------------------------------
11   CHARACTER (len=16)      &
12#if !defined key_agrif
13      , PARAMETER  &
14#endif
15      ::    & 
16      cp_cfg = "orca"           !: name of the configuration
17   INTEGER     &
18#if !defined key_agrif
19      , PARAMETER  &
20#endif
21      :: &
22      jp_cfg = 2,            &  !: resolution of the configuration (degrees)
23
24      ! data size              !!! * size of all input files *
25      jpidta  = 182,         &  !: 1st lateral dimension ( >= jpiglo )
26      jpjdta  = 149,         &  !: 2nd    "       "      ( >= jpjglo )
27      jpkdta  = 31              !: number of levels      ( >= jpk    )
28
29#if defined key_antarctic
30      ! zoom domain size       !!! *  antarctic zoom  *
31   INTEGER     &
32#if !defined key_agrif
33      , PARAMETER  &
34#endif
35      :: &
36      jpiglo  = jpidta,      &  !: 1st dimension of global domain --> i
37      jpjglo  = 50,          &  !: 2nd    "                  "    --> j
38      jpk     = jpkdta,      &  !: number of vertical levels
39      ! zoom starting position
40      jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom
41      jpjzoom =   1   ,      &  !: in data domain indices
42      ! Domain characteristics
43      jperio  =   1             !: lateral cond. type (between 0 and 6)
44
45#elif defined key_arctic
46      ! zoom domain size       !!! *  arctic zoom  *
47   INTEGER    &
48#if !defined key_agrif
49      , PARAMETER  &
50#endif
51      :: &
52      jpiglo  = 142   ,      &  !: 1st dimension of global domain --> i
53      jpjglo  = jpjdta-97+1, &  !: 2nd    "                  "    --> j
54      jpk     = jpkdta,      &  !: number of vertical levels
55      ! zoom starting position
56      jpizoom =  21   ,      &  !: left bottom (i,j) indices of the zoom
57      jpjzoom =  97   ,      &  !: in data domain indices
58      ! Domain characteristics
59      jperio  =   3             !: lateral cond. type (between 0 and 6)
60
61#elif defined key_c1d
62      ! global domain size     !!! *  global domain  *
63   INTEGER    &
64#if !defined key_agrif
65      , PARAMETER  &
66#endif
67      :: &
68      jpiglo  = 3     ,      &  !: 1st dimension of global domain --> i
69      jpjglo  = 3     ,      &  !: 2nd    "                  "    --> j
70      jpk     = jpkdta,      &  !: number of vertical levels
71      ! starting position of the zoom
72      ! jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom
73      ! jpjzoom =   133  ,    &  !: in data domain indices (160W,75N)
74      ! jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom
75      ! jpjzoom =   110  ,    &  !: in data domain indices (160W,50N)
76      ! jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom
77      ! jpjzoom =   97   ,    &  !: in data domain indices (160W,30N)
78      ! jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom
79      ! jpjzoom =   86   ,    &  !: in data domain indices (160W,10N)
80      ! jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom
81      ! jpjzoom =   49   ,    &  !: in data domain indices (160W,30S)
82      ! jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom
83      ! jpjzoom =   27   ,    &  !: in data domain indices (160W,60S)
84      jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom
85      jpjzoom =    7   ,    &  !: in data domain indices (160W,75S)
86      ! Domain characteristics
87      jperio  =   0            !: lateral cond. type (between 0 and 6)
88#else
89      ! global domain size     !!! *  global domain  *
90   INTEGER    &
91#if !defined key_agrif
92      , PARAMETER  &
93#endif
94      :: &
95      jpiglo  = jpidta,      &  !: 1st dimension of global domain --> i
96      jpjglo  = jpjdta,      &  !: 2nd    "                  "    --> j
97      jpk     = jpkdta,      &  !: number of vertical levels
98      ! starting position of the zoom
99      jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom
100      jpjzoom =   1   ,      &  !: in data domain indices
101      ! Domain characteristics
102      jperio  =   4             !: lateral cond. type (between 0 and 6)
103
104#endif
105
106   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config.
107   !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr
108   REAL(wp), PARAMETER ::   &
109      pp_not_used       = 999999_wp , &  !:
110      pp_to_be_computed = 0._wp          !:
111
112   !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 )
113
114   INTEGER,PARAMETER   ::    & !
115      jphgr_msh = 0            !: type of horizontal mesh
116      !                        !  = 0 curvilinear coordinate on the sphere
117      !                        !      read in coordinate.nc file
118      !                        !  = 1 geographical mesh on the sphere
119      !                        !      with regular grid-spacing
120      !                        !  = 2 f-plane with regular grid-spacing
121      !                        !  = 3 beta-plane with regular grid-spacing
122      !                        !  = 4 Mercator grid with T/U point at the equator  with
123      !                        !      isotropic resolution (e1_deg)
124
125      !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid.
126      !   The mercator grid starts only approximately at gphi0 because
127      !   of the constraint that the equator be a T point.
128   REAL(wp) ,PARAMETER ::       &  !
129      ppglam0  = pp_not_used,   &  !: longitude of first raw and column T-point (jphgr_msh = 1)
130      ppgphi0  = pp_not_used,   &  !: latitude  of first raw and column T-point (jphgr_msh = 1)
131      !                            !  latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3)
132      ppe1_deg = pp_not_used,   &  !: zonal      grid-spacing (degrees)
133      ppe2_deg = pp_not_used,   &  !: meridional grid-spacing (degrees)
134      !
135      ppe1_m   = pp_not_used,   &  !: zonal      grid-spacing (meters )
136      ppe2_m   = pp_not_used       !: meridional grid-spacing (meters )
137
138   !!
139   !! Vertical grid parameter for domzgr
140   !! ==================================
141   !!
142   REAL(wp), PARAMETER  ::       &
143      &     ppsur = -4762.96143546300_wp    ,  &  !: ORCA r4, r2 and r05 coefficients
144      &     ppa0  =   255.58049070440_wp    ,  &  !: (default coefficients)
145      &     ppa1  =   245.58132232490_wp    ,  &  !:
146      &     ppkth =    21.43336197938_wp    ,  &  !: (non dimensional): gives the approximate
147      !                                           !: layer number above which  stretching will
148      !                                           !: be maximum. Usually of order jpk/2.
149      &     ppacr =     3.00000000000_wp          !: (non dimensional): stretching factor
150      !                                           !: for the grid. The highest zacr, the smallest
151      !                                           !: the stretching.
152
153   !!
154   !!  If both ppa0 ppa1 and ppsur are specified to 0, then
155   !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr
156   !!
157   REAL(wp), PARAMETER ::        &
158      &     ppdzmin = pp_not_used           ,  &  !: (meters) vertical thickness of the top layer
159      &     pphmax  = pp_not_used                 !: (meters) Maximum depth of the ocean gdepw(jpk)
160   LOGICAL,  PARAMETER ::        &
161      &     ldbletanh = .FALSE.                   !: Use/do not use double tanf function for vertical coordinates
162   REAL(wp), PARAMETER ::        &
163      &     ppa2    = pp_not_used           ,  &  !: Double tanh function parameters
164      &     ppkth2  = pp_not_used           ,  &  !:
165      &     ppacr2  = pp_not_used                 !:
166   !!---------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.