1 | !----------------------------------------------------------------------- |
---|
2 | &namrun ! parameters of the run |
---|
3 | !----------------------------------------------------------------------- |
---|
4 | nn_no = 0 ! job number (no more used...) |
---|
5 | cn_exp = "eORCA1L3_OSM_N" ! experience name |
---|
6 | nn_it000 = 1 !8993 !33 !8992 ! first time step |
---|
7 | nn_itend = 35040 !72 !35040 !35040 !6 4 !11680 !32 !2464 !35040 !64 !9000 !116800 ! last time step (std 5475) |
---|
8 | nn_date0 = 19800101 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) |
---|
9 | nn_leapy = 0 ! Leap year calendar (1) or not (0) |
---|
10 | ln_rstart = .false. !.false. ! start from rest (F) or from a restart file (T) |
---|
11 | ! |
---|
12 | nn_rstctl = 2 !0 ! 2 ! restart control ==> activated only if ln_rstart=T |
---|
13 | cn_ocerst_in = "restart" ! "eORCA1L3_OSM_no_w_00000032_restart" ! suffix of ocean restart name (input) |
---|
14 | cn_ocerst_indir = "RESTARTS" ! directory from which to read input ocean restarts |
---|
15 | cn_ocerst_outdir = "RESTARTS" ! directory from which to read input ocean restarts |
---|
16 | nn_istate = 0 !1 ! output the initial state (1) or not (0) |
---|
17 | ln_rst_list = .false. ! output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F) |
---|
18 | nn_stock = 11680 !8892 ! 32 !8892 !11680 ! frequency of creation of a restart file (modulo referenced to 1) |
---|
19 | nn_stocklist = 8892,11680,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written |
---|
20 | nn_write = 11680 ! frequency of write in the output file (modulo referenced to nn_it000) |
---|
21 | ln_mskland = .false. ! mask land points in NetCDF outputs (costly: + ~15%) |
---|
22 | ln_cfmeta = .false. ! output additional data to netCDF files required for compliance with the CF metadata standard |
---|
23 | nn_chunksz = 0 ! chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) |
---|
24 | / |
---|
25 | !----------------------------------------------------------------------- |
---|
26 | &namcfg ! parameters of the configuration (default: user defined GYRE) |
---|
27 | !----------------------------------------------------------------------- |
---|
28 | ln_read_cfg = .true. ! (=T) read the domain configuration file |
---|
29 | ! ! (=F) user defined configuration ==>>> see usrdef(_...) modules |
---|
30 | cn_domcfg = "eORCA_R1_zps_domcfg" ! domain configuration filename |
---|
31 | / |
---|
32 | !----------------------------------------------------------------------- |
---|
33 | &namdom ! time and space domain |
---|
34 | !----------------------------------------------------------------------- |
---|
35 | rn_rdt = 2700. ! time step for the dynamics |
---|
36 | ln_meshmask = .false. !.true. ! create meshfile |
---|
37 | / |
---|
38 | !----------------------------------------------------------------------- |
---|
39 | &namtsd ! data : Temperature & Salinity |
---|
40 | !----------------------------------------------------------------------- |
---|
41 | ! ! =T read T-S fields for: |
---|
42 | ln_tsd_init = .true. ! ocean initialisation |
---|
43 | ln_tsd_dmp = .false. ! T-S restoring (see namtra_dmp) |
---|
44 | |
---|
45 | cn_dir = '/hpcdata/scratch/agn/OSMOSIS/FORCING/SHACONEMO_eORCA_LIM3/' ! root directory for the T-S data location |
---|
46 | !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! |
---|
47 | ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! |
---|
48 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! |
---|
49 | sn_tem = 'conservative_temperature_WOA13_decav_Reg1L75_clim', -1 ,'votemper' , .true. , .true. , 'yearly' , 'weights_3D_WOA13d1_2_eorca1_bilinear.nc' , '' , '' |
---|
50 | sn_sal = 'absolute_salinity_WOA13_decav_Reg1L75_clim' , -1 ,'vosaline' , .true. , .true. , 'yearly' , 'weights_3D_WOA13d1_2_eorca1_bilinear.nc' , '' , '' |
---|
51 | / |
---|
52 | !----------------------------------------------------------------------- |
---|
53 | &namwad ! Wetting and drying default is no WAD |
---|
54 | !----------------------------------------------------------------------- |
---|
55 | / |
---|
56 | !----------------------------------------------------------------------- |
---|
57 | &namcrs ! coarsened grid (for outputs and/or TOP) (ln_crs =T) |
---|
58 | !----------------------------------------------------------------------- |
---|
59 | / |
---|
60 | !----------------------------------------------------------------------- |
---|
61 | &namc1d ! 1D configuration options ("key_c1d") |
---|
62 | !----------------------------------------------------------------------- |
---|
63 | / |
---|
64 | !----------------------------------------------------------------------- |
---|
65 | &namc1d_dyndmp ! U & V newtonian damping ("key_c1d") |
---|
66 | !----------------------------------------------------------------------- |
---|
67 | / |
---|
68 | !----------------------------------------------------------------------- |
---|
69 | &namc1d_uvd ! data: U & V currents ("key_c1d") |
---|
70 | !----------------------------------------------------------------------- |
---|
71 | / |
---|
72 | !----------------------------------------------------------------------- |
---|
73 | &namsbc ! Surface Boundary Condition (surface module) |
---|
74 | !----------------------------------------------------------------------- |
---|
75 | nn_fsbc = 1 ! frequency of surface boundary condition computation |
---|
76 | ! (also = the frequency of sea-ice & iceberg model call) |
---|
77 | ! Type of air-sea fluxes |
---|
78 | ln_blk = .true. ! Bulk formulation (T => fill namsbc_blk ) |
---|
79 | nn_ice = 2 ! =0 no ice boundary condition |
---|
80 | ! =1 use observed ice-cover ( => fill namsbc_iif ) |
---|
81 | ! =2 or 3 automatically for LIM3 or CICE ("key_lim3" or "key_cice") |
---|
82 | ! except in AGRIF zoom where it has to be specified |
---|
83 | ln_ice_embd = .false. ! =T embedded sea-ice (pressure + mass and salt exchanges) |
---|
84 | ! =F levitating ice (no pressure, mass and salt exchanges) |
---|
85 | ! Misc. options of sbc : |
---|
86 | ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr) |
---|
87 | ln_dm2dc = .true. ! daily mean to diurnal cycle on short wave |
---|
88 | ln_rnf = .true. ! runoffs (T => fill namsbc_rnf) |
---|
89 | ln_ssr = .true. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) |
---|
90 | nn_fwb = 2 ! FreshWater Budget: =0 unchecked |
---|
91 | ! =1 global mean of e-p-r set to zero at each time step |
---|
92 | ! =2 annual global mean of e-p-r set to zero |
---|
93 | ln_isf = .true. ! ice shelf (T => fill namsbc_isf) |
---|
94 | ln_wave = .true. ! Activate coupling with wave (T => fill namsbc_wave) |
---|
95 | ln_cdgw = .false. ! Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave) |
---|
96 | ln_sdw = .true. !.false. ! Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave) |
---|
97 | nn_sdrift = 1 ! Parameterization for the calculation of 3D-Stokes drift from the surface Stokes drift |
---|
98 | ! = 0 Breivik 2015 parameterization: v_z=v_0*[exp(2*k*z)/(1-8*k*z)] |
---|
99 | ! = 1 Phillips: v_z=v_o*[exp(2*k*z)-beta*sqrt(-2*k*pi*z)*erfc(sqrt(-2*k*z))] |
---|
100 | ! = 2 Phillips as (1) but using the wave frequency from a wave model |
---|
101 | ln_tauwoc = .false. ! Activate ocean stress modified by external wave induced stress (T => ln_wave=.true. & fill namsbc_wave) |
---|
102 | ln_stcor = .true. ! Activate Stokes Coriolis term (T => ln_wave=.true. & ln_sdw=.true. & fill namsbc_wave) |
---|
103 | nn_lsm = 0 ! =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) , |
---|
104 | ! =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field) |
---|
105 | / |
---|
106 | !----------------------------------------------------------------------- |
---|
107 | &namsbc_flx ! surface boundary condition : flux formulation |
---|
108 | !----------------------------------------------------------------------- |
---|
109 | / |
---|
110 | !----------------------------------------------------------------------- |
---|
111 | &namsbc_blk ! namsbc_blk generic Bulk formula (ln_blk =T) |
---|
112 | !----------------------------------------------------------------------- |
---|
113 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! |
---|
114 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! |
---|
115 | sn_wndi = 'u_10_core_oonly' , 6 , 'U_10_MOD', .true. , .false. , 'yearly' , 'weights_coreII_2_eorca1_bicubic.nc' , 'U1' , '' |
---|
116 | sn_wndj = 'v_10_core_oonly' , 6 , 'V_10_MOD', .true. , .false. , 'yearly' , 'weights_coreII_2_eorca1_bicubic.nc' , 'V1' , '' |
---|
117 | sn_qsr = 'rad_core_oonly' , 24 , 'SWDN_MOD', .false. , .false. , 'yearly' , 'weights_coreII_2_eorca1_bilinear.nc' , '' , '' |
---|
118 | sn_qlw = 'rad_core_oonly' , 24 , 'LWDN_MOD', .true. , .false. , 'yearly' , 'weights_coreII_2_eorca1_bilinear.nc' , '' , '' |
---|
119 | sn_tair = 't_10_core_oonly' , 6 , 'T_10_MOD', .false. , .false. , 'yearly' , 'weights_coreII_2_eorca1_bilinear.nc' , '' , '' |
---|
120 | sn_humi = 'q_10_core_oonly' , 6 , 'Q_10_MOD', .false. , .false. , 'yearly' , 'weights_coreII_2_eorca1_bilinear.nc' , '' , '' |
---|
121 | sn_prec = 'tprecip_core_oonly' , -1 , 'TPRECIP' , .true. , .false. , 'yearly' , 'weights_coreII_2_eorca1_bilinear.nc' , '' , '' |
---|
122 | sn_snow = 'tprecip_core_oonly' , -1 , 'SNOW' , .true. , .false. , 'yearly' , 'weights_coreII_2_eorca1_bilinear.nc' , '' , '' |
---|
123 | sn_slp = 'slp_core_oonly' , 6 , 'SLP' , .true. , .false. , 'yearly' , 'weights_coreII_2_eorca1_bilinear.nc' , '' , '' |
---|
124 | ln_NCAR = .true. ! "NCAR" algorithm (Large and Yeager 2008) |
---|
125 | cn_dir = '/hpcdata/scratch/agn/OSMOSIS/FORCING/IAF_CORE/' ! root directory for the location of the flux files |
---|
126 | / |
---|
127 | !----------------------------------------------------------------------- |
---|
128 | &namsbc_cpl ! coupled ocean/atmosphere model ("key_oasis3") |
---|
129 | !----------------------------------------------------------------------- |
---|
130 | / |
---|
131 | !----------------------------------------------------------------------- |
---|
132 | &namsbc_sas ! Stand-Alone Surface boundary condition |
---|
133 | !----------------------------------------------------------------------- |
---|
134 | / |
---|
135 | !----------------------------------------------------------------------- |
---|
136 | &namsbc_iif ! Ice-IF : use observed ice cover (nn_ice = 1) |
---|
137 | !----------------------------------------------------------------------- |
---|
138 | / |
---|
139 | !----------------------------------------------------------------------- |
---|
140 | &namtra_qsr ! penetrative solar radiation (ln_traqsr =T) |
---|
141 | !----------------------------------------------------------------------- |
---|
142 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! |
---|
143 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! |
---|
144 | sn_chl ='merged_ESACCI_BIOMER4V1R1_CHL_REG05', -1 , 'CHLA' , .true. , .true. , 'yearly' , 'weights_reg05_2_eorca1_bilinear.nc' , '' , '' |
---|
145 | cn_dir = '/hpcdata/scratch/agn/OSMOSIS/FORCING/SHACONEMO_eORCA_LIM3/' ! root directory for the location of the flux files |
---|
146 | ln_qsr_rgb = .false. ! RGB (Red-Green-Blue) light penetration |
---|
147 | ln_qsr_2bd = .true. ! 2 bands light penetration |
---|
148 | ln_qsr_bio = .false. ! bio-model light penetration |
---|
149 | nn_chldta = 1 ! RGB : Chl data (=1) or cst value (=0) |
---|
150 | rn_abs = 0.58 ! RGB & 2 bands: fraction of light (rn_si1) |
---|
151 | rn_si0 = 0.35 ! RGB & 2 bands: shortess depth of extinction |
---|
152 | rn_si1 = 23.0 ! 2 bands: longest depth of extinction |
---|
153 | / |
---|
154 | !----------------------------------------------------------------------- |
---|
155 | &namsbc_rnf ! runoffs namelist surface boundary condition (ln_rnf =T) |
---|
156 | !----------------------------------------------------------------------- |
---|
157 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask! |
---|
158 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! |
---|
159 | sn_rnf = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc', -1 , 'sorunoff', .true. , .true. , 'yearly' , '' , '' , '' |
---|
160 | sn_cnf = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc', 0 , 'socoeff' , .false. , .true. , 'yearly' , '' , '' , '' |
---|
161 | |
---|
162 | cn_dir = '/hpcdata/scratch/agn/OSMOSIS/FORCING/SHACONEMO_eORCA_LIM3/' ! root directory for the location of the flux files |
---|
163 | ln_rnf_mouth = .false. ! specific treatment at rivers mouths |
---|
164 | ln_rnf_depth = .false. ! read in depth information for runoff |
---|
165 | ln_rnf_tem = .false. ! read in temperature information for runoff |
---|
166 | ln_rnf_sal = .false. ! read in salinity information for runoff |
---|
167 | ln_rnf_depth_ini = .false. !.true. ! compute depth at initialisation from runoff file |
---|
168 | rn_rnf_max = 0.05 ! max value of the runoff climatology over global domain ( if ln_rnf_depth_ini = .true ) |
---|
169 | rn_dep_max = 150. ! depth over which runoffs is spread ( ln_rnf_depth_ini = .true ) |
---|
170 | nn_rnf_depth_file = 0 ! create (=1) a runoff depth file or not (=0) |
---|
171 | / |
---|
172 | !----------------------------------------------------------------------- |
---|
173 | &namsbc_isf ! Top boundary layer (ISF) (nn_isf >0) |
---|
174 | !----------------------------------------------------------------------- |
---|
175 | ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask! |
---|
176 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! |
---|
177 | sn_rnfisf = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc' , -12 ,'sornfisf' , .false. , .true. , 'yearly' , '' , '' , '' |
---|
178 | sn_depmax_isf = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc' , -12 ,'sodepmax_isf' , .false. , .true. , 'yearly' , '' , '' , '' |
---|
179 | sn_depmin_isf = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc' , -12 ,'sodepmin_isf' , .false. , .true. , 'yearly' , '' , '' , '' |
---|
180 | nn_isf = 3 ! ice shelf melting/freezing |
---|
181 | ! 1 = presence of ISF 2 = bg03 parametrisation |
---|
182 | ! 3 = rnf file for isf 4 = ISF fwf specified |
---|
183 | ! option 1 and 4 need ln_isfcav = .true. (domzgr) |
---|
184 | / |
---|
185 | !----------------------------------------------------------------------- |
---|
186 | &namsbc_iscpl ! land ice / ocean coupling option |
---|
187 | !----------------------------------------------------------------------- |
---|
188 | / |
---|
189 | !----------------------------------------------------------------------- |
---|
190 | &namsbc_apr ! Atmospheric pressure used as ocean forcing (ln_apr_dyn =T) |
---|
191 | !----------------------------------------------------------------------- |
---|
192 | / |
---|
193 | !----------------------------------------------------------------------- |
---|
194 | &namsbc_ssr ! surface boundary condition : sea surface restoring (ln_ssr =T) |
---|
195 | !----------------------------------------------------------------------- |
---|
196 | ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! |
---|
197 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! |
---|
198 | sn_sss = 'sss_absolute_salinity_WOA13_decav_Reg1L75_clim', -1. , 'sosaline', .true. , .true. , 'yearly' , 'weights_WOA13d1_2_eorca1_bilinear.nc' , '' , '' |
---|
199 | cn_dir = '/hpcdata/scratch/agn/OSMOSIS/FORCING/SHACONEMO_eORCA_LIM3/' ! root directory for the location of the flux files |
---|
200 | rn_deds = -166.67 ! magnitude of the damping on salinity [mm/day] |
---|
201 | ln_sssr_bnd = .true. ! flag to bound erp term (associated with nn_sssr=2) |
---|
202 | rn_sssr_bnd = 4.e0 ! ABS(Max/Min) value of the damping erp term [mm/day] |
---|
203 | / |
---|
204 | !----------------------------------------------------------------------- |
---|
205 | &namsbc_wave ! External fields from wave model (ln_wave=T) |
---|
206 | !----------------------------------------------------------------------- |
---|
207 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! |
---|
208 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! |
---|
209 | ! sn_cdg = 'sdw_wave' , 1 , 'drag_coeff' , .true. , .false. , 'yearly' , '' , '' , '' |
---|
210 | sn_usd = 'ust' , 6 , 'ust' , .true. , .false. , 'monthly' , 'weights_grid05_eORCA1_bicubic.nc' , 'U1' , '' |
---|
211 | sn_vsd = 'vst' , 6 , 'vst' , .true. , .false. , 'monthly' , 'weights_grid05_eORCA1_bicubic.nc' , 'V1' , '' |
---|
212 | sn_hsw = 'swh' , 6 , 'swh' , .true. , .false. , 'monthly' , 'weights_grid05_eORCA1_bilin.nc' , '' , '' |
---|
213 | sn_wmp = 'mwp' , 6 , 'mwp' , .true. , .false. , 'monthly' , 'weights_grid05_eORCA1_bilin.nc' , '' , '' |
---|
214 | ! sn_wnum = 'sdw_wave' , 1 , 'wave_num' , .true. , .false. , 'yearly' , '' , '' , '' |
---|
215 | ! sn_tauoc = 'sdw_wave' , 1 , 'wave_stress', .true. , .false. , 'yearly' , '' , '' , '' |
---|
216 | ! |
---|
217 | cn_dir = '/hpcdata/scratch/agn/OSMOSIS/FORCING/ERA_INTERIM/WAVES/' ! root directory for the location of the flux files |
---|
218 | / |
---|
219 | !----------------------------------------------------------------------- |
---|
220 | &namberg ! iceberg parameters (default: No iceberg) |
---|
221 | !----------------------------------------------------------------------- |
---|
222 | / |
---|
223 | !----------------------------------------------------------------------- |
---|
224 | &namlbc ! lateral momentum boundary condition |
---|
225 | !----------------------------------------------------------------------- |
---|
226 | ! ! free slip ! partial slip ! no slip ! strong slip |
---|
227 | rn_shlat = 0. ! shlat = 0 ! 0 < shlat < 2 ! shlat = 2 ! 2 < shlat |
---|
228 | ln_vorlat = .false. ! consistency of vorticity boundary condition with analytical Eqs. |
---|
229 | / |
---|
230 | !----------------------------------------------------------------------- |
---|
231 | &namagrif ! AGRIF zoom ("key_agrif") |
---|
232 | !----------------------------------------------------------------------- |
---|
233 | / |
---|
234 | !----------------------------------------------------------------------- |
---|
235 | &nam_tide ! tide parameters |
---|
236 | !----------------------------------------------------------------------- |
---|
237 | / |
---|
238 | !----------------------------------------------------------------------- |
---|
239 | &nambdy ! unstructured open boundaries |
---|
240 | !----------------------------------------------------------------------- |
---|
241 | / |
---|
242 | !----------------------------------------------------------------------- |
---|
243 | &nambdy_dta ! open boundaries - external data |
---|
244 | !----------------------------------------------------------------------- |
---|
245 | / |
---|
246 | !----------------------------------------------------------------------- |
---|
247 | &nambdy_tide ! tidal forcing at open boundaries |
---|
248 | !----------------------------------------------------------------------- |
---|
249 | / |
---|
250 | !----------------------------------------------------------------------- |
---|
251 | &namdrg ! top/bottom drag coefficient (default: NO selection) |
---|
252 | !----------------------------------------------------------------------- |
---|
253 | ln_lin = .true. ! linear drag: Cd = Cd0 Uc0 |
---|
254 | / |
---|
255 | !----------------------------------------------------------------------- |
---|
256 | &namdrg_top ! TOP friction (ln_isfcav=T) |
---|
257 | !----------------------------------------------------------------------- |
---|
258 | / |
---|
259 | !----------------------------------------------------------------------- |
---|
260 | &namdrg_bot ! BOTTOM friction |
---|
261 | !----------------------------------------------------------------------- |
---|
262 | / |
---|
263 | !----------------------------------------------------------------------- |
---|
264 | &nambbc ! bottom temperature boundary condition (default: NO) |
---|
265 | !----------------------------------------------------------------------- |
---|
266 | ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! |
---|
267 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F ) ! 'monthly' ! filename ! pairing ! filename ! |
---|
268 | sn_qgh ='Goutorbe_ghflux.nc', -12. , 'gh_flux' , .false. , .true. , 'yearly' , 'weights_Goutorbe1_2_eorca1_bilinear.nc' , '' , '' |
---|
269 | ! |
---|
270 | ln_trabbc = .true. ! Apply a geothermal heating at the ocean bottom |
---|
271 | nn_geoflx = 2 ! geothermal heat flux: = 0 no flux |
---|
272 | ! ! = 1 constant flux |
---|
273 | ! ! = 2 variable flux (read in geothermal_heating.nc in mW/m2) |
---|
274 | rn_geoflx_cst = 86.4e-3 ! Constant value of geothermal heat flux [W/m2] |
---|
275 | cn_dir = '/hpcdata/scratch/agn/OSMOSIS/FORCING/SHACONEMO_eORCA_LIM3/' ! root directory for the location of the flux files |
---|
276 | / |
---|
277 | !----------------------------------------------------------------------- |
---|
278 | &nambbl ! bottom boundary layer scheme (default: NO) |
---|
279 | !----------------------------------------------------------------------- |
---|
280 | / |
---|
281 | !----------------------------------------------------------------------- |
---|
282 | &nameos ! ocean Equation Of Seawater (default: NO) |
---|
283 | !----------------------------------------------------------------------- |
---|
284 | ln_teos10 = .true. ! = Use TEOS-10 equation of state |
---|
285 | / |
---|
286 | !----------------------------------------------------------------------- |
---|
287 | &namtra_adv ! advection scheme for tracer (default: NO selection) |
---|
288 | !----------------------------------------------------------------------- |
---|
289 | ln_traadv_fct = .true. ! FCT scheme |
---|
290 | nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order |
---|
291 | nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order |
---|
292 | / |
---|
293 | !----------------------------------------------------------------------- |
---|
294 | &namtra_mle ! mixed layer eddy parametrisation (Fox-Kemper param) (default: NO) |
---|
295 | !----------------------------------------------------------------------- |
---|
296 | ln_mle = .true. ! (T) use the Mixed Layer Eddy (MLE) parameterisation |
---|
297 | rn_ce = 0.04 ! magnitude of the MLE (typical value: 0.06 to 0.08) |
---|
298 | |
---|
299 | / |
---|
300 | !----------------------------------------------------------------------- |
---|
301 | &namtra_ldf ! lateral diffusion scheme for tracers |
---|
302 | !---------------------------------------------------------------------------------- |
---|
303 | ! ! Operator type: |
---|
304 | ln_traldf_lap = .true. ! laplacian operator |
---|
305 | ln_traldf_blp = .false. ! bilaplacian operator |
---|
306 | ! |
---|
307 | ! ! Direction of action: |
---|
308 | ln_traldf_iso = .true. ! iso-neutral (standard operator) |
---|
309 | ln_traldf_triad = .false. ! iso-neutral (triad operator) |
---|
310 | ! |
---|
311 | ! ! iso-neutral options: |
---|
312 | ln_traldf_msc = .true. ! Method of Stabilizing Correction (both operators) |
---|
313 | rn_slpmax = 0.01 ! slope limit (both operators) |
---|
314 | ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only) |
---|
315 | rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only) |
---|
316 | ln_botmix_triad = .false. ! lateral mixing on bottom (triad only) |
---|
317 | ! |
---|
318 | ! ! Coefficients: |
---|
319 | nn_aht_ijk_t = 20 ! space/time variation of eddy coefficient: |
---|
320 | ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file |
---|
321 | ! ! = 0 constant |
---|
322 | ! ! = 10 F(k) =ldf_c1d |
---|
323 | ! ! = 20 F(i,j) =ldf_c2d |
---|
324 | ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation |
---|
325 | ! ! = 30 F(i,j,k) =ldf_c2d * ldf_c1d |
---|
326 | ! ! = 31 F(i,j,k,t)=F(local velocity and grid-spacing) |
---|
327 | ! ! time invariant coefficients: aht0 = 1/2 Ud*Ld (lap case) |
---|
328 | ! ! or = 1/12 Ud*Ld^3 (blp case) |
---|
329 | / |
---|
330 | !----------------------------------------------------------------------- |
---|
331 | &namtra_ldfeiv ! eddy induced velocity param. (default: NO) |
---|
332 | !----------------------------------------------------------------------- |
---|
333 | ln_ldfeiv = .true. ! use eddy induced velocity parameterization |
---|
334 | rn_aeiv_0 = 1000. ! eddy induced velocity coefficient [m2/s] |
---|
335 | nn_aei_ijk_t = 21 ! space/time variation of the eiv coeficient |
---|
336 | ! ! =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file |
---|
337 | ! ! = 0 constant |
---|
338 | ! ! = 10 F(k) =ldf_c1d |
---|
339 | ! ! = 20 F(i,j) =ldf_c2d |
---|
340 | ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation |
---|
341 | ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d |
---|
342 | ln_ldfeiv_dia =.true. ! diagnose eiv stream function and velocities |
---|
343 | / |
---|
344 | !----------------------------------------------------------------------- |
---|
345 | &namtra_dmp ! tracer: T & S newtonian damping (default: NO) |
---|
346 | !----------------------------------------------------------------------- |
---|
347 | ln_tradmp = .false. |
---|
348 | / |
---|
349 | !----------------------------------------------------------------------- |
---|
350 | &nam_vvl ! vertical coordinate options (default: z-star) |
---|
351 | !----------------------------------------------------------------------- |
---|
352 | ln_vvl_zstar = .true. ! z-star vertical coordinate |
---|
353 | / |
---|
354 | !----------------------------------------------------------------------- |
---|
355 | &namdyn_adv ! formulation of the momentum advection (default: NO selection) |
---|
356 | !----------------------------------------------------------------------- |
---|
357 | ln_dynadv_vec = .true. ! vector form - 2nd centered scheme |
---|
358 | nn_dynkeg = 0 ! grad(KE) scheme: =0 C2 ; =1 Hollingsworth correction |
---|
359 | ln_dynadv_cen2 = .false. ! flux form - 2nd order centered scheme |
---|
360 | / |
---|
361 | !----------------------------------------------------------------------- |
---|
362 | &namdyn_vor ! Vorticity / Coriolis scheme (default: NO) |
---|
363 | !----------------------------------------------------------------------- |
---|
364 | ln_dynvor_een = .true. ! energy & enstrophy scheme |
---|
365 | nn_een_e3f = 1 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) |
---|
366 | / |
---|
367 | !----------------------------------------------------------------------- |
---|
368 | &namdyn_hpg ! Hydrostatic pressure gradient option (default: NO selection) |
---|
369 | !----------------------------------------------------------------------- |
---|
370 | ln_hpg_sco = .true. ! s-coordinate (standard jacobian formulation) |
---|
371 | / |
---|
372 | !----------------------------------------------------------------------- |
---|
373 | &namdyn_spg ! surface pressure gradient (default: NO) |
---|
374 | !----------------------------------------------------------------------- |
---|
375 | ln_dynspg_ts = .true. ! split-explicit free surface |
---|
376 | ln_bt_fw = .false. ! Forward integration of barotropic Eqs. |
---|
377 | ln_bt_av = .true. ! Time filtering of barotropic variables |
---|
378 | nn_bt_flt = 1 ! Time filter choice = 0 None |
---|
379 | ! ! = 1 Boxcar over nn_baro sub-steps |
---|
380 | ! ! = 2 Boxcar over 2*nn_baro " " |
---|
381 | ln_bt_auto = .true. ! Number of sub-step defined from: |
---|
382 | rn_bt_cmax = 0.8 ! =T : the Maximum Courant Number allowed |
---|
383 | nn_baro = 30 ! =F : the number of sub-step in rn_rdt seconds |
---|
384 | rn_bt_alpha = 0. ! Temporal diffusion parameter (if ln_bt_av=F) |
---|
385 | / |
---|
386 | !----------------------------------------------------------------------- |
---|
387 | &namdyn_ldf ! lateral diffusion on momentum (default: NO selection) |
---|
388 | !----------------------------------------------------------------------- |
---|
389 | ! ! Type of the operator : |
---|
390 | ! ! no diffusion: set ln_dynldf_lap=..._blp=F |
---|
391 | ln_dynldf_lap = .true. ! laplacian operator |
---|
392 | ! ! Direction of action : |
---|
393 | ln_dynldf_lev = .true. ! iso-level |
---|
394 | ! ! Coefficient |
---|
395 | nn_ahm_ijk_t = 30 ! space/time variation of eddy coef |
---|
396 | ! ! =-30 read in eddy_viscosity_3D.nc file |
---|
397 | ! ! =-20 read in eddy_viscosity_2D.nc file |
---|
398 | ! ! = 0 constant |
---|
399 | ! ! = 10 F(k)=c1d |
---|
400 | ! ! = 20 F(i,j)=F(grid spacing)=c2d |
---|
401 | ! ! = 30 F(i,j,k)=c2d*c1d |
---|
402 | ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) |
---|
403 | ! ! = 32 F(i,j,k)=F(local gridscale and deformation rate) |
---|
404 | ! ! time invariant coefficients : ahm = 1/2 Uv*Lv (lap case) |
---|
405 | ! ! or = 1/12 Uv*Lv^3 (blp case) |
---|
406 | / |
---|
407 | !----------------------------------------------------------------------- |
---|
408 | &namzdf ! vertical physics (default: NO selection) |
---|
409 | !----------------------------------------------------------------------- |
---|
410 | ln_zdfosm = .true. ! OSMOSIS BL closure (T => fill namzdf_osm) |
---|
411 | ln_zdftke = .false. |
---|
412 | ! |
---|
413 | ! ! convection |
---|
414 | ln_zdfevd = .false. ! enhanced vertical diffusion |
---|
415 | rn_avm0 = 1.2e-4 ! vertical eddy viscosity [m2/s] (background Kz if ln_zdfcst=F) |
---|
416 | rn_avt0 = 1.2e-5 ! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F) |
---|
417 | / |
---|
418 | !----------------------------------------------------------------------- |
---|
419 | &namzdf_ric ! richardson number dependent vertical diffusion (ln_zdfric =T) |
---|
420 | !----------------------------------------------------------------------- |
---|
421 | / |
---|
422 | !----------------------------------------------------------------------- |
---|
423 | &namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion (ln_zdftke =T) |
---|
424 | !----------------------------------------------------------------------- |
---|
425 | ! |
---|
426 | nn_etau = 0 ! penetration of tke below the mixed layer (ML) due to NIWs |
---|
427 | ! = 0 none ; = 1 add a tke source below the ML |
---|
428 | ! = 2 add a tke source just at the base of the ML |
---|
429 | ! = 3 as = 1 applied on HF part of the stress (ln_cpl=T) |
---|
430 | / |
---|
431 | !----------------------------------------------------------------------- |
---|
432 | &namzdf_gls ! GLS vertical diffusion (ln_zdfgls =T) |
---|
433 | !----------------------------------------------------------------------- |
---|
434 | / |
---|
435 | !----------------------------------------------------------------------- |
---|
436 | &namzdf_osm ! OSM vertical diffusion (ln_zdfosm =T) |
---|
437 | !----------------------------------------------------------------------- |
---|
438 | ln_use_osm_la = .false. ! Use namelist rn_osm_la |
---|
439 | rn_osm_la = 0.3 ! Turbulent Langmuir number |
---|
440 | rn_osm_dstokes = 5. ! Depth scale of Stokes drift (m) |
---|
441 | nn_ave = 0 ! choice of horizontal averaging on avt, avmu, avmv |
---|
442 | ln_dia_osm = .true. ! output OSMOSIS-OBL variables |
---|
443 | rn_osm_hbl0 = 10. ! initial hbl value |
---|
444 | ln_kpprimix = .true. ! Use KPP-style Ri# mixing below BL |
---|
445 | rn_riinfty = 0.7 ! Highest local Ri_g permitting shear instability |
---|
446 | rn_difri = 0.005 ! max Ri# diffusivity at Ri_g = 0 (m^2/s) |
---|
447 | ln_convmix = .true. ! Use convective instability mixing below BL |
---|
448 | rn_difconv = 1. !0.01 !1. ! diffusivity when unstable below BL (m2/s) |
---|
449 | nn_osm_wave = 2 ! Method used to calculate Stokes drift |
---|
450 | ! = 2: Use ECMWF wave fields |
---|
451 | ! = 1: Pierson Moskowitz wave spectrum |
---|
452 | ! = 0: Constant La# = 0.3 |
---|
453 | !========================only if key_FK_osm is set ================================= |
---|
454 | ln_osm_mle = .true. ! Use integrated FK-OSM model |
---|
455 | / |
---|
456 | !----------------------------------------------------------------------- |
---|
457 | &namosm_mle ! mixed layer eddy parametrisation (Fox-Kemper) (default: OFF) |
---|
458 | !----------------------------------------------------------------------- |
---|
459 | rn_osm_mle_ce = 0.04 ! magnitude of the MLE (typical value: 0.06 to 0.08) |
---|
460 | nn_osm_mle = 0 ! MLE type: =0 standard Fox-Kemper ; =1 new formulation |
---|
461 | rn_osm_mle_lf = 5.e+3 ! typical scale of mixed layer front (meters) (case rn_osm_mle=0) |
---|
462 | rn_osm_mle_time = 43200 !172800. ! time scale for mixing momentum across the mixed layer (seconds) (case rn_osm_mle=0) |
---|
463 | rn_osm_mle_lat = 20. ! reference latitude (degrees) of MLE coef. (case rn_mle=1) |
---|
464 | rn_osm_mle_rho_c = 0.03 ! 0.01 ! delta rho criterion used to calculate MLD for FK |
---|
465 | rn_osm_mle_thresh = 0.0001 ! delta b criterion used for FK MLD criterion |
---|
466 | rn_osm_mle_tau = 172800. ! time scale for FK-OSM (seconds) (case rn_osm_mle=0) |
---|
467 | !========================only if key_FK_osm is set ================================= |
---|
468 | / |
---|
469 | !----------------------------------------------------------------------- |
---|
470 | &namzdf_iwm ! internal wave-driven mixing parameterization (ln_zdfiwm =T) |
---|
471 | !----------------------------------------------------------------------- |
---|
472 | / |
---|
473 | !----------------------------------------------------------------------- |
---|
474 | &nammpp ! Massively Parallel Processing ("key_mpp_mpi") |
---|
475 | !----------------------------------------------------------------------- |
---|
476 | ln_nnogather = .true. ! activate code to avoid mpi_allgather use at the northfold |
---|
477 | jpni = 16 !12 !16 !12 !16 !12 !16 ! 12 jpni number of processors following i (set automatically if < 1) |
---|
478 | jpnj = 15 !6 !15 !6 !15 !6 ! 15 ! 8 jpnj number of processors following j (set automatically if < 1) |
---|
479 | / |
---|
480 | !----------------------------------------------------------------------- |
---|
481 | &namctl ! Control prints |
---|
482 | !----------------------------------------------------------------------- |
---|
483 | ln_ctl = .false. ! Toggle all report printing on/off (T/F); Ignored if sn_cfctl%l_config is T |
---|
484 | / |
---|
485 | !----------------------------------------------------------------------- |
---|
486 | &namsto ! Stochastic parametrization of EOS (default: NO) |
---|
487 | !----------------------------------------------------------------------- |
---|
488 | / |
---|
489 | !----------------------------------------------------------------------- |
---|
490 | &namtrd ! trend diagnostics (default F) |
---|
491 | !----------------------------------------------------------------------- |
---|
492 | / |
---|
493 | !----------------------------------------------------------------------- |
---|
494 | &namptr ! Poleward Transport Diagnostic (default F) |
---|
495 | !----------------------------------------------------------------------- |
---|
496 | / |
---|
497 | !----------------------------------------------------------------------- |
---|
498 | &namhsb ! Heat and salt budgets (default F) |
---|
499 | !----------------------------------------------------------------------- |
---|
500 | / |
---|
501 | !----------------------------------------------------------------------- |
---|
502 | &namhth ! Extra ML depth and thermocline depth diagnostics (default F) |
---|
503 | !----------------------------------------------------------------------- |
---|
504 | ln_diahth = .true. ! produce extra diagnostics of BL and thermocline depth |
---|
505 | / |
---|
506 | !----------------------------------------------------------------------- |
---|
507 | &namdiu ! Cool skin and warm layer models (default F) |
---|
508 | !----------------------------------------------------------------------- |
---|
509 | / |
---|
510 | !----------------------------------------------------------------------- |
---|
511 | &namflo ! float parameters ("key_float") |
---|
512 | !----------------------------------------------------------------------- |
---|
513 | / |
---|
514 | !----------------------------------------------------------------------- |
---|
515 | &nam_diaharm ! Harmonic analysis of tidal constituents ("key_diaharm") |
---|
516 | !----------------------------------------------------------------------- |
---|
517 | / |
---|
518 | !----------------------------------------------------------------------- |
---|
519 | &namdct ! transports through some sections ("key_diadct") |
---|
520 | !----------------------------------------------------------------------- |
---|
521 | / |
---|
522 | !----------------------------------------------------------------------- |
---|
523 | &nam_diatmb ! Top Middle Bottom Output (default F) |
---|
524 | !----------------------------------------------------------------------- |
---|
525 | / |
---|
526 | !----------------------------------------------------------------------- |
---|
527 | &nam_dia25h ! 25h Mean Output (default F) |
---|
528 | !----------------------------------------------------------------------- |
---|
529 | ln_dia25h = .false. ! Choose 25h mean output or not |
---|
530 | / |
---|