-
Notifications
You must be signed in to change notification settings - Fork 1
/
par_ORCA_R2.h90
168 lines (157 loc) · 8.24 KB
/
par_ORCA_R2.h90
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
!!---------------------------------------------------------------------
!! *** par_ORCA_R2.h90 ***
!! Ocean Domain : 2 degrees resolution global ocean
!! (0RCA_R2 configuration)
!!---------------------------------------------------------------------
!!----------------------------------------------------------------------
!! NEMO/OPA 3.3 , NEMO Consortium (2010)
!! $Id: par_ORCA_R2.h90 2715 2011-03-30 15:58:35Z rblod $
!! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
!!----------------------------------------------------------------------
CHARACTER (len=16) &
#if !defined key_agrif
, PARAMETER &
#endif
:: &
cp_cfg = "orca" !: name of the configuration
INTEGER &
#if !defined key_agrif
, PARAMETER &
#endif
:: &
jp_cfg = 2, & !: resolution of the configuration (degrees)
! data size !!! * size of all input files *
jpidta = 182, & !: 1st lateral dimension ( >= jpiglo )
jpjdta = 149, & !: 2nd " " ( >= jpjglo )
jpkdta = 75 !: number of levels ( >= jpk )
#if defined key_antarctic
! zoom domain size !!! * antarctic zoom *
INTEGER &
#if !defined key_agrif
, PARAMETER &
#endif
:: &
jpiglo = jpidta, & !: 1st dimension of global domain --> i
jpjglo = 50, & !: 2nd " " --> j
! zoom starting position
jpizoom = 1 , & !: left bottom (i,j) indices of the zoom
jpjzoom = 1 , & !: in data domain indices
! Domain characteristics
jperio = 1 !: lateral cond. type (between 0 and 6)
#elif defined key_arctic
! zoom domain size !!! * arctic zoom *
INTEGER &
#if !defined key_agrif
, PARAMETER &
#endif
:: &
jpiglo = 142 , & !: 1st dimension of global domain --> i
jpjglo = jpjdta-97+1, & !: 2nd " " --> j
! zoom starting position
jpizoom = 21 , & !: left bottom (i,j) indices of the zoom
jpjzoom = 97 , & !: in data domain indices
! Domain characteristics
jperio = 3 !: lateral cond. type (between 0 and 6)
#elif defined key_c1d
! global domain size !!! * global domain *
INTEGER &
#if !defined key_agrif
, PARAMETER &
#endif
:: &
jpiglo = 3 , & !: 1st dimension of global domain --> i
jpjglo = 3 , & !: 2nd " " --> j
! starting position of the zoom
! jpizoom = 61 , & !: left bottom (i,j) indices of the zoom
! jpjzoom = 133 , & !: in data domain indices (160W,75N)
! jpizoom = 61 , & !: left bottom (i,j) indices of the zoom
! jpjzoom = 110 , & !: in data domain indices (160W,50N)
! jpizoom = 61 , & !: left bottom (i,j) indices of the zoom
! jpjzoom = 97 , & !: in data domain indices (160W,30N)
! jpizoom = 61 , & !: left bottom (i,j) indices of the zoom
! jpjzoom = 86 , & !: in data domain indices (160W,10N)
! jpizoom = 61 , & !: left bottom (i,j) indices of the zoom
! jpjzoom = 49 , & !: in data domain indices (160W,30S)
! jpizoom = 61 , & !: left bottom (i,j) indices of the zoom
! jpjzoom = 27 , & !: in data domain indices (160W,60S)
! jpizoom = 61 , & !: left bottom (i,j) indices of the zoom
! jpjzoom = 7 , & !: in data domain indices (160W,75S)
! jpizoom = 110 , & !: left bottom (i,j) indices of the zoom
! jpjzoom = 97 , & !: in data domain indices (64W,31.5N) BATS site
jpizoom = 175 , & !: left bottom (i,j) indices of the zoom
jpjzoom = 36 , & !: in data domain indices KERFIX site
! jpizoom = 179 , & !: left bottom (i,j) indices of the zoom
! jpjzoom = 38 , & !: in data domain indices PLUME site
! Domain characteristics
jperio = 0 !: lateral cond. type (between 0 and 6)
#else
! global domain size !!! * global domain *
INTEGER &
#if !defined key_agrif
, PARAMETER &
#endif
:: &
jpiglo = jpidta, & !: 1st dimension of global domain --> i
jpjglo = jpjdta, & !: 2nd " " --> 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 = 4 !: 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 (jphgr_msh /= 0 )
INTEGER,PARAMETER :: & !
jphgr_msh = 0 !: 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)
! 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) ,PARAMETER :: & !
ppglam0 = pp_not_used, & !: longitude of first raw and column T-point (jphgr_msh = 1)
ppgphi0 = pp_not_used, & !: 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 )
!!
!! Vertical grid parameter for domzgr
!! ==================================
!!
REAL(wp), PARAMETER :: &
& ppsur = -4762.96143546300_wp , & !: ORCA r4, r2 and r05 coefficients
& ppa0 = 255.58049070440_wp , & !: (default coefficients)
& ppa1 = 245.58132232490_wp , & !:
& ppkth = 21.43336197938_wp , & !: (non dimensional): gives the approximate
! !: layer number above which stretching will
! !: be maximum. Usually of order jpk/2.
& ppacr = 3.00000000000_wp !: (non dimensional): stretching factor
! !: for the grid. The highest zacr, the smallest
! !: the stretching.
!!
!! If both 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) 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 !:
!!---------------------------------------------------------------------