Newer
Older
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
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!! NEMO/OCE Configuration namelist : used to overwrite defaults values defined in SHARED/namelist_ref
!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
!-----------------------------------------------------------------------
&namusr_def ! User defined : OVERFLOW configuration
!-----------------------------------------------------------------------
! ! type of vertical coordinate
ln_zco = .false. ! z-coordinate
ln_zps = .true. ! z-partial-step coordinate
ln_sco = .false. ! s-coordinate
rn_dx = 5000. ! horizontal resolution [meters]
rn_dz = 60. ! vertical resolution [meters]
rn_f0 = 1.e-4 ! coriolis [s-1]
/
!-----------------------------------------------------------------------
&namrun ! parameters of the run
!-----------------------------------------------------------------------
nn_no = 0 ! job number (no more used...)
cn_exp = "DOME" ! experience name
nn_it000 = 1 ! first time step
nn_itend = 11520 ! here 16h of simulation (=5760 time-step) abort after 5802 for zps: pb of physics conditions
nn_stock = 0 ! frequency of creation of a restart file (modulo referenced to 1)
/
!-----------------------------------------------------------------------
&namcfg ! parameters of the configuration
!-----------------------------------------------------------------------
ln_read_cfg = .true. ! (=T) read the domain configuration file
! ! (=F) user defined configuration (F => create/check namusr_def)
cn_domcfg = "DOME_domcfg" ! domain configuration filename
!
/
!-----------------------------------------------------------------------
&namtile ! parameters of the tiling
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namdom ! space and time domain (bathymetry, mesh, timestep)
!-----------------------------------------------------------------------
ln_linssh = .false.
rn_Dt = 300. ! time step for the dynamics (and tracer if nn_acc=0)
rn_atfp = 0.1 ! asselin time filter parameter
ln_meshmask = .false. ! =T create a mesh file
/
!-----------------------------------------------------------------------
&namtsd ! data : Temperature & Salinity
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namsbc ! Surface Boundary Condition (surface module)
!-----------------------------------------------------------------------
nn_fsbc = 1 ! frequency of surface boundary condition computation
! (also = the frequency of sea-ice & iceberg model call)
ln_usr = .true. ! user defined formulation (T => check usrdef_sbc)
ln_blk = .false. ! Bulk formulation (T => fill namsbc_blk )
nn_ice = 0 ! =0 no ice boundary condition
ln_traqsr = .false. ! Light penetration in the ocean (T => fill namtra_qsr )
ln_rnf = .false. ! runoffs (T => fill namsbc_rnf)
ln_ssr = .false. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr)
nn_fwb = 0 ! FreshWater Budget: =0 unchecked
/
!-----------------------------------------------------------------------
&namlbc ! lateral momentum boundary condition (default: NO selection)
!-----------------------------------------------------------------------
rn_shlat = 0. ! free slip
/
!!======================================================================
!! *** Top/Bottom boundary condition *** !!
!!======================================================================
!! namdrg top/bottom drag coefficient (default: NO selection)
!! namdrg_top top friction (ln_isfcav=T)
!! namdrg_bot bottom friction
!! nambbc bottom temperature boundary condition (default: OFF)
!! nambbl bottom boundary layer scheme (default: OFF)
!!======================================================================
!
!-----------------------------------------------------------------------
&namdrg ! top/bottom drag coefficient (default: NO selection)
!-----------------------------------------------------------------------
ln_drg_OFF = .false. ! free-slip : Cd = 0 (F => fill namdrg_bot
ln_lin = .false. ! linear drag: Cd = Cd0 Uc0 & namdrg_top)
ln_non_lin = .true. ! non-linear drag: Cd = Cd0 |U|
ln_loglayer = .false. ! logarithmic drag: Cd = vkarmn/log(z/z0) |U|
!
ln_drgimp = .true. ! implicit top/bottom friction flag
/
!-----------------------------------------------------------------------
&namdrg_bot ! BOTTOM friction (ln_OFF =F)
!-----------------------------------------------------------------------
rn_Cd0 = 2.e-3 ! drag coefficient [-]
rn_ke0 = 0. ! background kinetic energy [m2/s2] (non-linear cases)
/
!-----------------------------------------------------------------------
&nameos ! ocean physical parameters
!-----------------------------------------------------------------------
ln_seos = .true. ! = Use simplified equation of state (S-EOS)
! ! rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS
rn_a0 = 0.2 ! thermal expension coefficient (for simplified equation of state)
rn_b0 = 0. ! saline expension coefficient (for simplified equation of state)
rn_lambda1 = 0. ! cabbeling coeff in T^2 (=0 for linear eos)
rn_lambda2 = 0. ! cabbeling coeff in S^2 (=0 for linear eos)
rn_mu1 = 0. ! thermobaric coeff. in T (=0 for linear eos)
rn_mu2 = 0. ! thermobaric coeff. in S (=0 for linear eos)
rn_nu = 0. ! cabbeling coeff in T*S (=0 for linear eos)
/
!-----------------------------------------------------------------------
&namtra_adv ! advection scheme for tracer
!----------------------------------------------------------------------
ln_traadv_fct = .true. ! FCT scheme
nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order
nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order
/
!-----------------------------------------------------------------------
&namtra_ldf ! lateral diffusion scheme for tracers (default: NO selection)
!-----------------------------------------------------------------------
! ! Operator type:
ln_traldf_OFF = .true. ! No explicit diffusion
ln_traldf_lap = .false. ! laplacian operator
ln_traldf_blp = .false. ! bilaplacian operator
!
! ! Direction of action:
ln_traldf_lev = .false. ! iso-level
ln_traldf_hor = .true. ! horizontal (geopotential)
ln_traldf_iso = .false. ! iso-neutral (standard operator)
ln_traldf_triad = .false. ! iso-neutral (triad operator)
!
! ! iso-neutral options:
ln_traldf_msc = .false. ! Method of Stabilizing Correction (both operators)
rn_slpmax = 0.01 ! slope limit (both operators)
ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only)
rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only)
ln_botmix_triad = .false. ! lateral mixing on bottom (triad only)
!
! ! Coefficients:
nn_aht_ijk_t = 20 ! space/time variation of eddy coefficient:
! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file
! ! = 0 constant
! ! = 10 F(k) =ldf_c1d
! ! = 20 F(i,j) =ldf_c2d
! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation
! ! = 30 F(i,j,k) =ldf_c2d * ldf_c1d
! ! = 31 F(i,j,k,t)=F(local velocity and grid-spacing)
! ! time invariant coefficients: aht0 = 1/2 Ud*Ld (lap case)
! ! or = 1/12 Ud*Ld^3 (blp case)
rn_Ud = 0.02 ! lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30)
rn_Ld = 200.e+3 ! lateral diffusive length [m] (nn_aht_ijk_t= 0, 10)
/
!!======================================================================
!! *** Dynamics namelists *** !!
!!======================================================================
!
!-----------------------------------------------------------------------
&namdyn_adv ! formulation of the momentum advection (default: NO selection)
!-----------------------------------------------------------------------
ln_dynadv_OFF = .false. ! linear dynamics (no momentum advection)
ln_dynadv_vec = .false. ! vector form (T) or flux form (F)
nn_dynkeg = 0 ! scheme for grad(KE): =0 C2 ; =1 Hollingsworth correction
ln_dynadv_cen2 = .false. ! flux form - 2nd order centered scheme
ln_dynadv_ubs = .true. ! flux form - 3rd order UBS scheme
/
!-----------------------------------------------------------------------
&nam_vvl ! vertical coordinate options (default: zstar)
!-----------------------------------------------------------------------
ln_vvl_zstar = .true. ! zstar vertical coordinate
/
!-----------------------------------------------------------------------
&namdyn_vor ! option of physics/algorithm
!-----------------------------------------------------------------------
ln_dynvor_een = .true. ! energy & enstrophy scheme
!
nn_e3f_typ = 1 ! type of e3f (EEN, ENE, ENS, MIX only) =0 e3f = mi(mj(e3t))/4
! ! =1 e3f = mi(mj(e3t))/mi(mj( tmask))
/
!-----------------------------------------------------------------------
&namdyn_hpg ! Hydrostatic pressure gradient option
!-----------------------------------------------------------------------
ln_hpg_zps = .false. ! z-coordinate - partial steps (interpolation)
ln_hpg_djc = .true. ! s-coordinate (standard jacobian formulation)
/
!-----------------------------------------------------------------------
&namdyn_spg ! Surface pressure gradient
!-----------------------------------------------------------------------
ln_dynspg_ts = .true. ! split-explicit free surface
ln_bt_fw = .true. ! Forward integration of barotropic Eqs.
ln_bt_av = .true. ! Time filtering of barotropic variables
nn_bt_flt = 1 ! Time filter choice = 0 None
! ! = 1 Boxcar over nn_e sub-steps
! ! = 2 Boxcar over 2*nn_e " "
ln_bt_auto = .false. ! Number of sub-step defined from:
nn_e = 20 ! =F : the number of sub-step in rn_Dt seconds
/
!-----------------------------------------------------------------------
&namdyn_ldf ! lateral diffusion on momentum (default: NO selection)
!-----------------------------------------------------------------------
! ! Type of the operator :
ln_dynldf_OFF = .true. ! No operator (i.e. no explicit diffusion)
/
!!======================================================================
!! vertical physics namelists !!
!!======================================================================
!-----------------------------------------------------------------------
&namzdf ! vertical physics (default: NO selection)
!-----------------------------------------------------------------------
ln_zad_Aimp = .false. ! Courant number dependent scheme (Shchepetkin 2015)
! ! type of vertical closure (required)
ln_zdfcst = .false. ! constant mixing
ln_zdfric = .false. ! local Richardson dependent formulation (T => fill namzdf_ric)
ln_zdftke = .false. ! Turbulent Kinetic Energy closure (T => fill namzdf_tke)
ln_zdfgls = .true. ! Generic Length Scale closure (T => fill namzdf_gls)
ln_zdfosm = .false. ! OSMOSIS BL closure (T => fill namzdf_osm)
!
! ! convection
ln_zdfevd = .false. ! enhanced vertical diffusion
nn_evdm = 0 ! apply on tracer (=0) or on tracer and momentum (=1)
rn_evd = 100. ! mixing coefficient [m2/s]
ln_zdfnpc = .false. ! Non-Penetrative Convective algorithm
nn_npc = 1 ! frequency of application of npc
nn_npcp = 365 ! npc control print frequency
!
ln_zdfddm = .false. ! double diffusive mixing
rn_avts = 1.e-4 ! maximum avs (vertical mixing on salinity)
rn_hsbfr = 1.6 ! heat/salt buoyancy flux ratio
!
! ! gravity wave-driven vertical mixing
ln_zdfiwm = .false. ! internal wave-induced mixing (T => fill namzdf_iwm)
ln_zdfswm = .false. ! surface wave-induced mixing (T => ln_wave=ln_sdw=T )
!
! ! coefficients
rn_avm0 = 1.e-5 ! vertical eddy viscosity [m2/s] (background Kz if ln_zdfcst=F)
rn_avt0 = 1.e-6 ! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F)
nn_avb = 0 ! profile for background avt & avm (=1) or not (=0)
nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0)
/
!-----------------------------------------------------------------------
&nambdy ! unstructured open boundaries (default: OFF)
!-----------------------------------------------------------------------
ln_bdy = .true. ! Use unstructured open boundaries
nb_bdy = 2 ! number of open boundary sets
ln_coords_file = .false.,.false. ! =T : read bdy coordinates from file
ln_mask_file = .false. ! =T : read mask from file
cn_dyn2d = 'frs','frs' !
nn_dyn2d_dta = 0,0 ! = 0, bdy data are equal to the initial state
! ! = 1, bdy data are read in 'bdydata .nc' files
! ! = 2, use tidal harmonic forcing data from files
! ! = 3, use external data AND tidal harmonic forcing
cn_dyn3d = 'frs','frs' !
nn_dyn3d_dta = 0,0 ! = 0, bdy data are equal to the initial state
! ! = 1, bdy data are read in 'bdydata .nc' files
cn_tra = 'frs','frs' !
nn_tra_dta = 0,0 ! = 0, bdy data are equal to the initial state
! ! = 1, bdy data are read in 'bdydata .nc' files
cn_ice = 'none','none' !
nn_ice_dta = 0,0 ! = 0, bdy data are equal to the initial state
! ! = 1, bdy data are read in 'bdydata .nc' files
!
ln_tra_dmp =.false.,.false. ! open boudaries conditions for tracers
ln_dyn3d_dmp =.false.,.false. ! open boundary condition for baroclinic velocities
rn_time_dmp = 1.,1. ! Damping time scale in days
rn_time_dmp_out = 5.,5. ! Outflow damping time scale
nn_rimwidth = 1,10 ! width of the relaxation zone
ln_vol = .true.
/
!-----------------------------------------------------------------------
&nambdy_index ! bdy segment definition
!-----------------------------------------------------------------------
ctypebdy='N'
nbdyind = -1
nbdybeg = -1
nbdyend = -1
/
!-----------------------------------------------------------------------
&nambdy_index ! bdy segment definition
!-----------------------------------------------------------------------
ctypebdy='W'
nbdyind = 2
nbdybeg = 1
nbdyend = 164
/