Changeset 990 for branches/dev_003_CPL/CONFIG/ORCA2_LIM/EXP00
- Timestamp:
- 2008-05-23T16:38:21+02:00 (16 years ago)
- Location:
- branches/dev_003_CPL/CONFIG/ORCA2_LIM/EXP00
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_003_CPL/CONFIG/ORCA2_LIM/EXP00/1_namelist
r889 r990 1 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 ! OPA namelist : model option and parameter input 3 ! ------------- 4 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 5 ! 1 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 !! NEMO/OPA : 1 - run manager (namrun, nam_ctl, nam_mpp, nam_mpp_dyndist, nam_ctl) 3 !! namelists 2 - Domain (nam_zgr, nam_zgr_sco, namdom) 4 !! 3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core 5 !! namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, namalb) 6 !! 4 - lateral boundary (namlbc, namcla, namobc) 7 !! 5 - bottom boundary (nambfr, nambbc, nambbl) 8 !! 6 - Tracer (nameos, nam_traadv, nam_traldf, namtdp) 9 !! 7 - dynamics (nam_dynadv, nam_dynvor, nam_dynhpg, namflg, nam_dynspg, nam_dynldf) 10 !! 8 - Verical physics (namzdf, namnpc, namric, namtke, namkpp, namddm) 11 !! 9 - diagnostics (namtrd, namgap, namspr, namflo, namptr) 12 !! 9 - miscellaneous (namsol) 13 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 14 ! CAUTION: some scripts does not support CAPITALs for logical definition 15 ! ******* use .true. or .false. and NOT .TRUE. or .FALSE. 16 17 !!====================================================================== 18 !! *** Run management namelists *** 19 !!====================================================================== 20 !! namrun parameters of the run 21 !! nam_mpp Massively Parallel Processing 22 !! nam_mpp_dyndist 23 !! nam_ctl Control prints & Benchmark 24 !!====================================================================== 25 6 26 !----------------------------------------------------------------------- 7 27 ! namrun parameters of the run … … 10 30 ! cexper experience name for vairmer format 11 31 ! ln_rstart boolean term for restart (true or false) 12 ! nrstdt control of the restart timestep: 13 ! = 0 restart, do not control nit000 in the restart file. 14 ! = 1 restart, control nit000 in the restart file. Do not 15 ! use the date in the restart file (use ndate0 in namelist) 16 ! = 2 restart, control nit000 in the restart file, use the date 17 ! in the restart file. ndate0 in the namelist is ignored. 32 ! nrstdt restart control = 0 restart, do not control nit000 in the restart file. 33 ! ! = 1 restart, control nit000 in the restart file. Do not 34 ! ! use the date in the restart file (use ndate0 in namelist) 35 ! ! = 2 restart, control nit000 in the restart file, use the date 36 ! ! in the restart file. ndate0 in the namelist is ignored. 18 37 ! nit000 number of the first time step 19 38 ! nitend number of the last time step … … 24 43 ! nwrite frequency of OUTPUT file 25 44 ! ln_dimgnnn (F/T) 1 DIMG file - (for all proc/per proc) 26 !27 ! CAUTION: for usual run scripts, logical value of28 ! ******* ln_rstart must be .true. or .false.29 ! and NOT .TRUE. or .FALSE.30 45 &namrun 31 46 no = 0 … … 38 53 nleapy = 0 39 54 ninist = 0 40 nstock = 10950 41 nwrite = 10950 42 ln_dimgnnn = .FALSE. 55 nstock = 10950 56 nwrite = 10950 57 ln_dimgnnn = .false. 58 / 59 !----------------------------------------------------------------------- 60 ! nam_mpp Massively Parallel Processing 61 !----------------------------------------------------------------------- 62 ! c_mpi_send mpi send/recieve type 63 ! = 'S' : standard blocking send 64 ! = 'B' : buffer blocking send 65 ! = 'I' : immediate non-blocking send 66 &nam_mpp 67 c_mpi_send = 'S' 68 / 69 !----------------------------------------------------------------------- 70 ! nam_mpp_dyndist Massively Parallel Distribution ("key_agrif" && "key_mpp_dyndist") 71 !----------------------------------------------------------------------- 72 ! jpni number of processors following i 73 ! jpnj number of processors following j 74 ! jpnij number of local domains 75 &nam_mpp_dyndist 76 jpni = 1 77 jpnj = 1 78 jpnij = 1 43 79 / 44 80 !----------------------------------------------------------------------- … … 69 105 nbit_cmp = 0 70 106 / 71 !----------------------------------------------------------------------- 72 ! nam_mpp Massively Parallel Processing 73 !----------------------------------------------------------------------- 74 ! c_mpi_send mpi send/recieve type 75 ! = 'S' : standard blocking send 76 ! = 'B' : buffer blocking send 77 ! = 'I' : immediate non-blocking send 78 &nam_mpp 79 c_mpi_send = 'S' 80 / 81 !----------------------------------------------------------------------- 82 ! nam_mpp_dyndist Massively Parallel Distribution 83 ! #ifdef ( key_agrif && key_mpp_dyndist ) 84 !----------------------------------------------------------------------- 85 ! jpni number of processors following i 86 ! jpnj number of processors following j 87 ! jpnij number of local domains 88 &nam_mpp_dyndist 89 jpni = 1 90 jpnj = 1 91 jpnij = 1 92 / 107 108 !!====================================================================== 109 !! *** Domain namelists *** 110 !!====================================================================== 111 !! nam_zgr vertical coordinate 112 !! nam_zgr_sco s-coordinate or hybrid z-s-coordinate 113 !! namdom space and time domain (bathymetry, mesh, timestep) 114 !!====================================================================== 115 93 116 !----------------------------------------------------------------------- 94 117 ! nam_zgr vertical coordinate 95 118 !----------------------------------------------------------------------- 96 ! ln_zco z-coordinate - full steps (T/F) 119 ! ln_zco z-coordinate - full steps (T/F) ("key_zco" may also be defined) 97 120 ! ln_zps z-coordinate - partial steps (T/F) 98 121 ! ln_sco s- or hybrid z-s-coordinate (T/F) … … 118 141 / 119 142 !----------------------------------------------------------------------- 120 ! nam_traadv advection scheme for tracer (option not control by CPP keys)121 !-----------------------------------------------------------------------122 ! ln_traadv_cen2 2nd order centered scheme (default T)123 ! ln_traadv_tvd TVD scheme (default F)124 ! ln_traadv_muscl MUSCL scheme (default F)125 ! ln_traadv_muscl2 MUSCL2 scheme (default F)126 ! ln_traadv_ubs UBS scheme (default F)127 ! ln_traadv_qck QUICKEST scheme (default F)128 &nam_traadv129 ln_traadv_cen2 = .false.130 ln_traadv_tvd = .true.131 ln_traadv_muscl = .false.132 ln_traadv_muscl2 = .false.133 ln_traadv_ubs = .false.134 ln_traadv_qck = .false.135 /136 !-----------------------------------------------------------------------137 ! nam_traldf lateral diffusion scheme for tracer (option not control by CPP keys)138 !-----------------------------------------------------------------------139 ! Type of the operator :140 ! ln_traldf_lap laplacian operator (default T)141 ! ln_traldf_bilap bilaplacian operator (default F)142 ! Direction of action :143 ! ln_traldf_level iso-level (default F)144 ! ln_traldf_hor horizontal (geopotential) (default F)^**145 ! ln_traldf_iso iso-neutral (default T)^*146 ! Coefficient147 ! aht0 horizontal eddy diffusivity for tracers (m2/s)148 ! ahtb0 background eddy diffusivity for isopycnal diffusion (m2/s)149 ! aeiv0 eddy induced velocity coefficient (m2/s)150 ! ^* require key_ldfslp to compute the direction of the lateral diffusion151 ! ^** require key_ldfslp in s-coordinate152 &nam_traldf153 ln_traldf_lap = .true.154 ln_traldf_bilap = .false.155 ln_traldf_level = .false.156 ln_traldf_hor = .false.157 ln_traldf_iso = .true.158 aht0 = 1000.159 ahtb0 = 0.160 aeiv0 = 0.161 /162 !-----------------------------------------------------------------------163 ! nam_dynldf lateral diffusion on momentum164 !-----------------------------------------------------------------------165 ! Type of the operator :166 ! ln_dynldf_lap laplacian operator (default T)167 ! ln_dynldf_bilap bilaplacian operator (default F)168 ! Direction of action :169 ! ln_dynldf_level iso-level (default F)170 ! ln_dynldf_hor horizontal (geopotential) (default F)^**171 ! ln_dynldf_iso iso-neutral (default T)^*172 ! Coefficient173 ! ahm0 horizontal eddy viscosity for the dynamics (m2/s)174 ! ahmb0 background eddy viscosity for isopycnal diffusion (m2/s)175 &nam_dynldf176 ln_dynldf_lap = .false.177 ln_dynldf_bilap = .true.178 ln_dynldf_level = .false.179 ln_dynldf_hor = .true.180 ln_dynldf_iso = .false.181 ahm0 = -8.5e+11182 ahmb0 = 0.183 /184 !-----------------------------------------------------------------------185 ! namflg algorithm flags (algorithm not control by CPP keys)186 !-----------------------------------------------------------------------187 ! ln_dynhpg_imp hydrostatic pressure gradient: semi-implicit time scheme (T)188 ! centered time scheme (F)189 ! nn_dynhpg_rst add dynhpg implicit variables in restart ot not (1/0)190 &namflg191 ln_dynhpg_imp = .false.192 nn_dynhpg_rst = 0193 /194 !-----------------------------------------------------------------------195 ! nam_dynhpg Hydrostatic pressure gradient option196 !-----------------------------------------------------------------------197 ! type of pressure gradient scheme (choose one only!)198 ! ln_hpg_zco z-coordinate - full steps (default T)199 ! ln_hpg_zps z-coordinate - partial steps (interpolation)200 ! ln_hpg_sco s-coordinate (standard jacobian formulation)201 ! ln_hpg_hel s-coordinate (helsinki modification)202 ! ln_hpg_wdj s-coordinate (weighted density jacobian)203 ! ln_hpg_djc s-coordinate (Density Jacobian with Cubic polynomial)204 ! ln_hpg_rot s-coordinate (ROTated axes scheme)205 ! parameters206 ! gamm weighting coefficient (wdj scheme)207 &nam_dynhpg208 ln_hpg_zco = .false.209 ln_hpg_zps = .true.210 ln_hpg_sco = .false.211 ln_hpg_hel = .false.212 ln_hpg_wdj = .false.213 ln_hpg_djc = .false.214 ln_hpg_rot = .false.215 gamm = 0.e0216 /217 !-----------------------------------------------------------------------218 ! nam_dynadv option of physics/algorithm (not control by CPP keys)219 !-----------------------------------------------------------------------220 ! ln_dynadv_vec vector form flag221 ! ln_dynadv_cen2 flux form - 2nd order centered scheme (default T)222 ! ln_dynadv_ubs flux form - 3rd order UBS scheme (default F)223 &nam_dynadv224 ln_dynadv_vec = .TRUE.225 ln_dynadv_cen2 = .FALSE.226 ln_dynadv_ubs = .FALSE.227 /228 !-----------------------------------------------------------------------229 ! nam_dynvor option of physics/algorithm (not control by CPP keys)230 !-----------------------------------------------------------------------231 ! ln_dynvor_ens vorticity trends: enstrophy conserving scheme (default T)232 ! ln_dynvor_ene " " : energy conserving scheme (default F)233 ! ln_dynvor_mix " " : mixed scheme (default F)234 ! ln_dynvor_een " " : energy & enstrophy scheme (default F)235 &nam_dynvor236 ln_dynvor_ene = .FALSE.237 ln_dynvor_ens = .FALSE.238 ln_dynvor_mix = .FALSE.239 ln_dynvor_een = .TRUE.240 /241 !-----------------------------------------------------------------------242 ! namalb albedo parameters243 !-----------------------------------------------------------------------244 ! cgren correction of the snow or ice albedo to take into account245 ! albice albedo of melting ice in the arctic and antarctic246 ! alphd coefficients for linear interpolation used to compute albedo247 ! between two extremes values (Pyane, 1972)248 ! alphc " "249 ! alphdi " "250 &namalb251 cgren = 0.06252 albice = 0.5253 alphd = 0.80254 alphc = 0.65255 alphdi = 0.72256 /257 !-----------------------------------------------------------------------258 143 ! namdom space and time domain (bathymetry, mesh, timestep) 259 144 !----------------------------------------------------------------------- 260 ! ntopo = 0/1 ,compute/read the bathymetry file 261 ! (mbathy, nb of T-ocean levels) 145 ! ntopo = 0/1 ,compute/read the bathymetry file (mbathy, nb of T-ocean levels) 262 146 ! e3zps_min the thickness of the partial step is set larger than the 263 ! e3zps_rat the minimum of e3zps_min and e3zps_rat * e3t 264 ! (N.B. 0<e3zps_rat<1) 265 ! (coordinates, scale factors) 147 ! e3zps_rat the minimum of e3zps_min and e3zps_rat * e3t (N.B. 0<e3zps_rat<1) 266 148 ! nmsh =1 create a mesh file (coordinates, scale factors, masks) 267 ! nacc the acceleration of convergence method 268 ! = 0, no acceleration, rdt = rdttra 269 ! = 1, acceleration used, rdt < rdttra(k) 149 ! nacc acceleration of convergence method = 0, not used, rdt = rdttra 150 ! = 1, used, rdt < rdttra(k) 270 151 ! atfp asselin time filter parameter 271 152 ! rdt time step for the dynamics (and tracer if nacc=0) … … 274 155 ! rdth depth variation of tracer time step 275 156 ! rdtbt barotropic time step (for the time splitting algorithm) 276 ! nclosea 0 no closed sea277 ! 1 closed sea (Black Sea, Caspian Sea, Great US Lakes...)157 ! nclosea =0 remove the closed sea from the global domain (orca configuration) 158 ! =1 closed sea (Black Sea, Caspian Sea, Great US Lakes...) 278 159 &namdom 279 160 ntopo = 1 … … 290 171 nclosea = 0 291 172 / 292 !----------------------------------------------------------------------- 293 ! namfwb freshwater budget correction 294 !----------------------------------------------------------------------- 295 ! ln_fwb logical flag for freshwater budget correction (0 annual mean) 296 &namfwb 297 ln_fwb = .false. 298 / 299 !----------------------------------------------------------------------- 300 ! namptr Poleward Transport Diagnostic 301 !----------------------------------------------------------------------- 302 ! ln_diaptr logical flag for Poleward transport computation 303 ! ln_subbas logical flag for Atlantic/Pacific/Indian basins computation 304 ! need input basins mask file named "subbasins.nc" 305 ! nf_ptr Frequency of computation 306 &namptr 307 ln_diaptr = .false. 308 ln_subbas = .false. 309 nf_ptr = 15 310 / 311 !----------------------------------------------------------------------- 312 ! namcro cross land advection 313 !----------------------------------------------------------------------- 314 ! n_cla advection between 2 ocean pts separates by land 315 &namcla 316 n_cla = 0 317 / 318 !----------------------------------------------------------------------- 319 ! namzdf vertical physics 320 !----------------------------------------------------------------------- 321 ! ln_zdfevd enhanced vertical diffusion (default T) 322 ! ln_zdfnpc Non-Penetrative Convection (default T) 323 ! avm0 vertical eddy viscosity for the dynamic (m2/s) 324 ! avt0 vertical eddy diffusivity for tracers (m2/s) 325 ! avevd vertical coefficient for enhanced diffusion scheme (m2/s) 326 ! nevdm = 0 apply enhanced mixing on tracer only 327 ! = 1 apply enhanced mixing on both tracer and momentum 328 ! ln_zdfexp vertical physics: (=T) time splitting (T) (Default=F) 329 ! (=F) euler backward (F) 330 ! n_zdfexp number of sub-timestep for time splitting scheme 331 &namzdf 332 ln_zdfevd = .true. 333 ln_zdfnpc = .false. 334 avm0 = 1.2e-4 335 avt0 = 1.2e-5 336 avevd = 100. 337 n_evdm = 1 338 ln_zdfexp = .false. 339 n_zdfexp = 3 340 / 341 !----------------------------------------------------------------------- 342 ! namnpc vnon penetrative convection 343 !----------------------------------------------------------------------- 344 ! nnpc1 non penetrative convective scheme frequency 345 ! nnpc2 non penetrative convective scheme print frequency 346 &namnpc 347 nnpc1 = 1 348 nnpc2 = 365 349 / 350 !----------------------------------------------------------------------- 351 ! nambbl bottom boundary layer scheme 352 !----------------------------------------------------------------------- 353 ! atrbbl lateral tracer coeff. for bottom boundary layer scheme(m2/s) 354 &nambbl 355 atrbbl = 10000. 356 / 357 !----------------------------------------------------------------------- 358 ! namric richardson number dependent vertical diffusion 359 ! ( #ifdef "key_zdfrichardson" ) 360 !----------------------------------------------------------------------- 361 ! avmri maximum value of the vertical viscosity 362 ! alp coefficient of the parameterization 363 ! nric coefficient of the parameterization 364 &namric 365 avmri = 100.e-4 366 alp = 5. 367 nric = 2 368 / 369 !----------------------------------------------------------------------- 370 ! namtke turbulent eddy kinetic dependent vertical diffusion 371 ! ( #ifdef "key_zdftke" ) 372 !----------------------------------------------------------------------- 373 ! ln_rstke flag to restart with tke from a run without tke (default F) 374 ! ediff coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) ) 375 ! ediss coef. of the Kolmogoroff dissipation 376 ! ebb coef. of the surface input of tke 377 ! efave coef. to applied to the tke diffusion ( avtke=efave*avm ) 378 ! emin minimum value of tke (m^2/s^2) 379 ! emin0 surface minimum value of tke (m^2/s^2) 380 ! nitke number of restart iterative loops 381 ! ri_c critic richardson number 382 ! nmxl flag on mixing length used 383 ! = 0 bounded by the distance to surface and bottom 384 ! = 1 bounded by the local vertical scale factor 385 ! = 2 first vertical derivative of mixing length bounded by 1 386 ! npdl flag on prandtl number 387 ! = 0 no vertical prandtl number (avt=avm) 388 ! = 1 prandtl number function of richarson number (avt=pdl*avm) 389 ! = 2 same as = 1 but a shapiro filter is applied on pdl 390 ! nave = horizontal averaged (=1) or not (=0) of avt (default =1) 391 ! navb = 0 cst background avt0, avm0 / =1 profile used on avtb 392 &namtke 393 ln_rstke = .false. 394 ediff = 0.1 395 ediss = 0.7 396 ebb = 60. 397 efave = 1. 398 emin = 1.e-6 399 emin0 = 1.e-4 400 nitke = 50 401 nmxl = 2 402 npdl = 1 403 navb = 0 404 / 405 !----------------------------------------------------------------------- 406 ! namkpp K-Profile Parameterization dependent vertical diffusion 407 ! ( #ifdef "key_zdfkpp" ) 408 !----------------------------------------------------------------------- 409 ! ln_kpprimix shear instability mixing (default T) 410 ! difmiw constant internal wave viscosity (m2/s) 411 ! difsiw constant internal wave diffusivity (m2/s) 412 ! Riinfty local Richardson Number limit for shear instability 413 ! difri maximum shear mixing at Rig = 0 (m2/s) 414 ! bvsqcon Brunt-Vaisala squared (1/s**2) for maximum convection 415 ! difcon maximum mixing in interior convection (m2/s) 416 ! nave = 0/1 flag for horizontal average on avt, avmu, avmv 417 ! navb = 0/1 flag for constant or profile background avt 418 &namkpp 419 ln_kpprimix = .true. 420 difmiw = 1.e-04 421 difsiw = 0.1e-04 422 Riinfty = 0.8 423 difri = 0.0050 424 bvsqcon = -0.01e-07 425 difcon = 1. 426 navb = 0 427 nave = 1 428 / 429 !----------------------------------------------------------------------- 430 ! namddm double diffusive mixing parameterization 431 !----------------------------------------------------------------------- 432 ! avts maximum avs for dd mixing 433 ! hsbfr heat/salt buoyancy flux ratio 434 &namddm 435 avts = 1.e-4 436 hsbfr = 1.6 437 / 438 !----------------------------------------------------------------------- 439 ! namlbc lateral momentum boundary condition 440 !----------------------------------------------------------------------- 441 ! shlat lateral boundary condition on velocity 442 ! shlat = 0 , free slip 443 ! 0 < shlat < 2 , partial slip 444 ! shlat = 2 , no slip 445 ! 2 < shlat , strong slip 446 &namlbc 447 shlat = 2. 448 / 449 !----------------------------------------------------------------------- 450 ! nambfr bottom friction 451 !----------------------------------------------------------------------- 452 ! nbotfr type of bottom friction 453 ! nbotfr = 0 , no slip 454 ! nbotfr = 1 , linear friction 455 ! nbotfr = 2 , nonlinear friction 456 ! nbotfr = 3 , free slip 457 ! bfri1 bottom drag coefficient (linear case) 458 ! bfri2 bottom drag coefficient (non linear case) 459 ! bfeb2 bottom turbulent kinetic energy (m^2/s^2) 460 &nambfr 461 nbotfr = 1 462 bfri1 = 4.e-4 463 bfri2 = 1.e-3 464 bfeb2 = 2.5e-3 465 / 466 !----------------------------------------------------------------------- 467 ! nambbc bottom temperature boundary condition 468 !----------------------------------------------------------------------- 469 ! ngeo_flux = 0 no geothermal heat flux 470 ! = 1 constant geothermal heat flux 471 ! = 2 variable geothermal heat flux (read in geothermal_heating.nc) 472 ! ( C A U T I O N : flux in mW/m2 in the NetCDF file ) 473 ! ngeo_flux_const Constant value of geothermal heat flux (W/m2) 474 &nambbc 475 ngeo_flux = 2 476 ngeo_flux_const = 86.4e-3 173 174 !!====================================================================== 175 !! *** Surface Boundary Condition namelists *** 176 !!====================================================================== 177 !! namsbc surface boundary condition 178 !! namsbc_ana analytical formulation 179 !! namsbc_flx flux formulation 180 !! namsbc_clio CLIO bulk formulea formulation 181 !! namsbc_core CORE bulk formulea formulation 182 !! namsbc_cpl CouPLed formulation ("key_coupled") 183 !! namqsr penetrative solar radiation 184 !! namsbc_rnf river runoffs 185 !! namsbc_ssr sea surface restoring term (for T and/or S) 186 !! namalb albedo parameters 187 !!====================================================================== 188 189 !----------------------------------------------------------------------- 190 ! namsbc Surface Boundary Condition (surface module) 191 !----------------------------------------------------------------------- 192 ! nn_fsbc frequency of surface boundary condition computation 193 ! (= the frequency of sea-ice model call) 194 ! ln_ana analytical formulation (fill namsbc_ana ) 195 ! ln_flx flux formulation (fill namsbc_flx ) 196 ! ln_blk_clio CLIO bulk formulation (fill namsbc_core) 197 ! ln_blk_core CORE bulk formulation (fill namsbc_clio) 198 ! ln_cpl Coupled formulation (fill namsbc_cpl ) 199 ! nn_ice =0 no ice boundary condition , 200 ! =1 observed ice-cover , 201 ! =2 ice-model used ("key_lim3" or "key_lim2) 202 ! nn_ico_cpl ice-ocean coupling : =0 LIM-3 old case 203 ! =1 stresses computed using now ocean velocity 204 ! =2 combination of 0 and 1 cases 205 ! ln_dm2dc Daily mean to Diurnal Cycle short wave (qsr) 206 ! ln_ssr Sea Surface Restoring on T and/or S (fill namsbc_ssr) 207 ! nn_fwb FreshWater Budget: =0 no control , 208 ! =1 annual global mean of e-p-r set to zero , 209 ! =2 global mean of e-p-r set to zero at each nn_fsbc time step 210 &namsbc 211 nn_fsbc = 5 212 ln_ana = .false. 213 ln_flx = .false. 214 ln_blk_clio = .true. 215 ln_blk_core = .false. 216 ln_cpl = .false. 217 nn_ice = 2 218 nn_ico_cpl = 0 219 ln_dm2dc = .false. 220 ln_rnf = .false. 221 ln_ssr = .false. 222 nn_fwb = 0 223 / 224 !----------------------------------------------------------------------- 225 ! namsbc_ana analytical surface boundary condition 226 !----------------------------------------------------------------------- 227 ! nn_tau000 gently increase the stress over the first ntau_rst time-steps 228 ! rn_utau0 uniform value for the i-stress 229 ! rn_vtau0 uniform value for the j-stress 230 ! rn_q0 uniform value for the total heat flux 231 ! rn_qsr0 uniform value for the solar radiation 232 ! rn_emp0 uniform value for the freswater budget (E-P) 233 &namsbc_ana 234 nn_tau000 = 0 235 rn_utau0 = 0.5 236 rn_vtau0 = 0.e0 237 rn_q0 = 0.e0 238 rn_qsr0 = 0.e0 239 rn_emp0 = 0.e0 240 / 241 !----------------------------------------------------------------------- 242 ! namsbc_flx surface boundary condition : flux formulation (#ifdef "key_sbc_flux") 243 !----------------------------------------------------------------------- 244 ! cn_dir directory in which the model is executed 245 ! 246 ! THE ORDER OF THE FILES MATTER: 247 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 248 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 249 &namsbc_flx 250 cn_dir = './' 251 sn_utau = 'utau.nc' , 24. , 'utau' , .false. , 0 , 0 252 sn_vtau = 'vtau.nc' , 24. , 'vtau' , .false. , 0 , 0 253 sn_qtot = 'qtot.nc' , 24. , 'qtot' , .false. , 0 , 0 254 sn_qsr = 'qsr.nc' , 24. , 'qsr' , .false. , 0 , 0 255 sn_emp = 'emp.nc' , 24. , 'emp' , .false. , 0 , 0 256 / 257 !----------------------------------------------------------------------- 258 ! namsbc_clio CLIO bulk formulea 259 !----------------------------------------------------------------------- 260 ! cn_dir directory in which the model is executed 261 ! 262 ! THE ORDER OF THE FILES MATTER: 263 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 264 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 265 &namsbc_clio 266 cn_dir = './' 267 sn_utau = 'taux_1m' , -12. , 'sozotaux', .false. , 1 , 0 268 sn_vtau = 'tauy_1m' , -12. , 'sometauy', .false. , 1 , 0 269 sn_wndm = 'flx' , -12. , 'socliowi', .false. , 1 , 0 270 sn_tair = 'flx' , -12. , 'socliot1', .false. , 1 , 0 271 sn_humi = 'flx' , -12. , 'socliohu', .false. , 1 , 0 272 sn_ccov = 'flx' , -12. , 'socliocl', .false. , 1 , 0 273 sn_prec = 'flx' , -12. , 'socliopl', .false. , 1 , 0 274 / 275 !----------------------------------------------------------------------- 276 ! namsbc_core CORE bulk formulea 277 !----------------------------------------------------------------------- 278 ! cn_dir directory in which the model is executed 279 ! ln_2m logical flag to use air temp. and hum referenced at 2m instead 10m 280 ! alpha_precip multiplication factor for precipitation (total & snow) 281 ! 282 ! THE ORDER OF THE FILES MATTER: 283 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 284 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 285 &namsbc_core 286 cn_dir = './' 287 ln_2m = .false. 288 alpha_precip = 1. 289 sn_wndi = 'u10' , 24. , 'U_10_MOD', .false. , 1 , 0 290 sn_wndj = 'v10' , 24. , 'V_10_MOD', .false. , 1 , 0 291 sn_qsr = 'rad' , 24. , 'SWDN_MOD', .false. , 1 , 0 292 sn_qlw = 'rad' , 24. , 'LWDN_MOD', .false. , 1 , 0 293 sn_tair = 't10' , 24. , 'T_10_MOD', .false. , 1 , 0 294 sn_humi = 'q10' , 24. , 'Q_10_MOD', .false. , 1 , 0 295 sn_prec = 'precip' , -12. , 'RAIN' , .false. , 1 , 0 296 sn_snow = 'precip' , -12. , 'SNOW' , .false. , 1 , 0 297 / 298 !----------------------------------------------------------------------- 299 ! namsbc_cpl coupled ocean/atmosphere model ("key_coupled") 300 !----------------------------------------------------------------------- 301 &namsbc_cpl 477 302 / 478 303 !----------------------------------------------------------------------- … … 490 315 / 491 316 !----------------------------------------------------------------------- 492 ! namtdp tracer newtonian damping ('key_tradmp') 493 !----------------------------------------------------------------------- 494 ! ndmp type of damping in temperature and salinity 495 ! (='latitude', damping poleward of 'ndmp' degrees and function 496 ! of the distance-to-coast. Red and Med Seas as ndmp=-1) 497 ! (=-1 damping only in Med and Red Seas) 498 ! ndmpf =1 create a damping.coeff NetCDF file (the 3D damping array) 499 ! nmldmp type of damping in the mixed layer 500 ! (=0 damping throughout the water column) 501 ! (=1 no damping in the mixed layer defined by avt >5cm2/s ) 502 ! (=2 no damping in the mixed layer defined rho<rho(surf)+.01 ) 503 ! sdmp surface time scale for internal damping (days) 504 ! bdmp bottom time scale for internal damping (days) 505 ! hdmp depth of transition between sdmp and bdmp (meters) 506 &namtdp 507 ndmp = 90 508 ndmpf = 1 509 nmldmp = 1 510 sdmp = 50. 511 bdmp = 360. 512 hdmp = 800. 513 / 514 !----------------------------------------------------------------------- 515 ! nameos ocean physical parameters 516 !----------------------------------------------------------------------- 517 ! neos type of equation of state and Brunt-Vaisala frequency 518 ! = 0, UNESCO (formulation of Jackett and McDougall (1994) 519 ! and of McDougall (1987) ) 520 ! = 1, linear: rho(T) = rau0 * ( 1.028 - ralpha * T ) 521 ! = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 522 ! with rau0=1020 set in parcst routine 523 ! ralpha thermal expension coefficient (linear equation of state) 524 ! rbeta saline expension coefficient (linear equation of state) 525 &nameos 526 neos = 0 527 ralpha = 2.e-4 528 rbeta = 0.001 529 / 530 !----------------------------------------------------------------------- 531 ! namsol elliptic solver / island / free surface 532 !----------------------------------------------------------------------- 533 ! nsolv elliptic solver (=1 preconditioned conjugate gradient: pcg) 534 ! (=2 successive-over-relaxation: sor) 535 ! (=3 FETI: fet, all require "key_feti" defined) 536 ! nsol_arp absolute/relative (0/1) precision convergence test 537 ! nmin minimum of iterations for the SOR solver 538 ! nmax maximum of iterations for the SOR solver 539 ! nmod frequency of test for the SOR solver 540 ! eps absolute precision of the solver 541 ! resmax absolute precision for the SOR solver 542 ! sor optimal coefficient for SOR solver 543 ! epsisl absolute precision on stream function solver 544 ! nmisl maximum pcg iterations for island 545 ! rnu strength of the additional force used in free surface b.c. 546 &namsol 547 nsolv = 1 548 nsol_arp = 0 549 nmin = 300 550 nmax = 800 551 nmod = 10 552 eps = 1.E-6 553 resmax = 1.E-10 554 sor = 1.92 555 epsisl = 1.e-10 556 nmisl = 4000 557 rnu = 1. 558 / 559 !======================================================================= 560 ! Diagnostics namelists 561 ! namtrd dynamics and/or tracer trends 562 ! namgap level mean model-data gap 563 ! namznl zonal mean heat & freshwater fluxes computation 564 ! namspr surface pressure in rigid-lid 565 !======================================================================= 566 !----------------------------------------------------------------------- 567 ! namtrd diagnostics on dynamics and/or tracer trends 568 ! ('key_trdyn' and/or 'key_trdtra') 569 ! or mixed-layer trends ('key_trdmld') 570 !----------------------------------------------------------------------- 571 ! ntrd time step frequency dynamics and tracers trends 572 ! nctls control surface type in mixed-layer trends (0,1 or n<jpk) 573 ! ln_trdmld_restart restart for ML diagnostics 574 ! ucf unit conversion factor (=1 -> /seconds | =86400. -> /day) 575 ! ln_trdmld_instant flag to diagnose trends of instantantaneous or mean ML T/S 576 &namtrd 577 ntrd = 365 578 nctls = 0 579 ln_trdmld_restart = .false. 580 ucf = 1. 581 ln_trdmld_instant = .false. 582 / 583 !----------------------------------------------------------------------- 584 ! namgap level mean model-data gap ('key_diagap') 585 !----------------------------------------------------------------------- 586 ! ngap time-step frequency of model-data gap computation 587 ! nprg time-step frequency of gap print in model output 588 &namgap 589 ngap = 15 590 nprg = 10 591 / 592 !----------------------------------------------------------------------- 593 ! namznl zonal mean heat & freshwater fluxes computation 594 ! (#ifdef "key_diaznl") 595 !----------------------------------------------------------------------- 596 ! nfznl time-step frequency of zonal mean fluxes computation 597 &namznl 598 nfznl = 15 599 / 600 !----------------------------------------------------------------------- 601 ! namspr surface pressure diagnostic 602 !----------------------------------------------------------------------- 603 ! nmaxp maximum of iterations for the solver 604 ! epsp absolute precision of the solver 605 ! niterp number of iteration done by the solver 606 &namspr 607 nmaxp = 1000 608 epsp = 1.e-3 609 niterp = 400 317 ! namsbc_rnf runoffs namelist surface boundary condition 318 !----------------------------------------------------------------------- 319 !cn_dir Root directory for location of ssr files 320 !nn_runoff =0 no, 1 runoff, 2 runoff+river mouth ups adv 321 !rn_hrnf runoffs, depth over which enhanced vertical mixing is used 322 !rn_avt_rnf runoffs, value of the additional vertical mixing coef. [m2/s] 323 !sn_rnf informations about the runoff file to be read 324 !sn_cnf informations about the runoff mouth file to be read 325 ! THE ORDER OF THE FILES MATTER: 326 ! ! file name ! frequency (hours) ! variable ! time interpolation ! 327 ! ! ! (if <0 months) ! name ! (logical) ! 328 &namsbc_rnf 329 cn_dir = './' 330 nn_runoff = 2 331 rn_hrnf = 0.e0 332 rn_avt_rnf = 1.e-3 333 sn_rnf = 'runoff_1m_nomask.nc' , -12. , 'sorunoff', .true. , 1 , 0 334 sn_cnf = 'runoff_1m_nomask.nc' , 0. , 'socoefr' , .false. , 1 , 0 335 / 336 !----------------------------------------------------------------------- 337 ! namsbc_ssr surface boundary condition : sea surface restoring 338 !----------------------------------------------------------------------- 339 ! cn_dir directory in which the model is executed 340 ! nn_sstr =0/1 add a retroaction term in the surface heat flux 341 ! nn_sssr =O/1/2 add a damping term in the surface freshwater flux 342 ! (=1, salt flux, concentration/dillution emps only) 343 ! (=2, volume flux, both emp and emps are updated) 344 ! dqdt magnitude of the retroaction on temperature [W/m2/K] 345 ! deds magnitude of the damping on salinity [mm/day] 346 ! THE ORDER OF THE FILES MATTER: 347 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 348 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 349 &namsbc_ssr 350 cn_dir = './' 351 nn_sstr = 0 352 nn_sssr = 0 353 dqdt = -40. 354 deds = 27.7 355 sn_sst = 'sst_data.nc' , 24. , 'sst' , .false. , 0 , 0 356 sn_sss = 'sss_data.nc' , -12. , 'sss' , .true. , 0 , 0 357 / 358 !----------------------------------------------------------------------- 359 ! namalb albedo parameters 360 !----------------------------------------------------------------------- 361 ! cgren correction of the snow or ice albedo to take into account 362 ! albice albedo of melting ice in the arctic and antarctic 363 ! alphd coefficients for linear interpolation used to compute albedo 364 ! between two extremes values (Pyane, 1972) 365 ! alphc " " 366 ! alphdi " " 367 &namalb 368 cgren = 0.06 369 albice = 0.5 370 alphd = 0.80 371 alphc = 0.65 372 alphdi = 0.72 373 / 374 375 !!====================================================================== 376 !! *** Lateral boundary condition *** 377 !!====================================================================== 378 !! namlbc lateral momentum boundary condition 379 !! namcla cross land advection 380 !! namobc open boundaries parameters ("key_obc") 381 !!====================================================================== 382 383 !----------------------------------------------------------------------- 384 ! namlbc lateral momentum boundary condition 385 !----------------------------------------------------------------------- 386 ! shlat shlat = 0 , free slip ; 0 < shlat < 2 , partial slip 387 ! shlat = 2 , no slip ; 2 < shlat , strong slip 388 &namlbc 389 shlat = 2. 390 / 391 !----------------------------------------------------------------------- 392 ! namcla cross land advection 393 !----------------------------------------------------------------------- 394 ! n_cla advection between 2 ocean pts separates by land 395 &namcla 396 n_cla = 0 610 397 / 611 398 !----------------------------------------------------------------------- … … 639 426 ln_vol_cst = .false. 640 427 / 641 !----------------------------------------------------------------------- 642 ! namflo float parameters (#ifdef key_float) 428 429 !!====================================================================== 430 !! *** Bottom boundary condition *** 431 !!====================================================================== 432 !! nambfr bottom friction 433 !! nambbc bottom temperature boundary condition ("key_trabbc") 434 !! nambbl bottom boundary layer scheme ("key_trabbl_dif","key_trabbl_adv") 435 !!====================================================================== 436 !----------------------------------------------------------------------- 437 ! nambfr bottom friction 438 !----------------------------------------------------------------------- 439 ! nbotfr type of bottom friction : = 0 : no slip ; = 1 : linear friction 440 ! = 3 : free slip ; = 2 : non linear friction 441 ! bfri1 bottom drag coefficient (linear case) 442 ! bfri2 bottom drag coefficient (non linear case) 443 ! bfeb2 bottom turbulent kinetic energy (m2/s2) 444 &nambfr 445 nbotfr = 1 446 bfri1 = 4.e-4 447 bfri2 = 1.e-3 448 bfeb2 = 2.5e-3 449 / 450 !----------------------------------------------------------------------- 451 ! nambbc bottom temperature boundary condition 452 !----------------------------------------------------------------------- 453 ! ngeo_flux = 0 no geothermal heat flux 454 ! = 1 constant geothermal heat flux 455 ! = 2 variable geothermal heat flux (read in geothermal_heating.nc in mW/m2) 456 ! ngeo_flux_const Constant value of geothermal heat flux (W/m2) 457 &nambbc 458 ngeo_flux = 2 459 ngeo_flux_const = 86.4e-3 460 / 461 !----------------------------------------------------------------------- 462 ! nambbl bottom boundary layer scheme 463 !----------------------------------------------------------------------- 464 ! ! diffusive bbl ("key_trabbl") 465 ! ! advective bbl ("key_trabbl_adv") 466 ! atrbbl mixing coefficient of the bottom boundary layer scheme (m2/s) 467 &nambbl 468 atrbbl = 10000. 469 / 470 471 !!====================================================================== 472 !! Tracer (T & S ) namelists 473 !!====================================================================== 474 !! nameos equation of state 475 !! nam_traadv advection scheme 476 !! nam_traldf lateral diffusion scheme 477 !! namtdp tracer newtonian damping ("key_tradmp") 478 !!====================================================================== 479 480 !----------------------------------------------------------------------- 481 ! nameos ocean physical parameters 482 !----------------------------------------------------------------------- 483 ! neos type of equation of state and Brunt-Vaisala frequency 484 ! ! = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 485 ! ! = 1, linear: rho(T) = rau0 * ( 1.028 - ralpha * T ) 486 ! ! = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 487 ! ralpha thermal expension coefficient (neos= 1 or 2) 488 ! rbeta saline expension coefficient (neos= 2) 489 &nameos 490 neos = 0 491 ralpha = 2.e-4 492 rbeta = 0.001 493 / 494 !----------------------------------------------------------------------- 495 ! nam_traadv advection scheme for tracer 496 !----------------------------------------------------------------------- 497 ! ln_traadv_cen2 2nd order centered scheme (default T) 498 ! ln_traadv_tvd TVD scheme (default F) 499 ! ln_traadv_muscl MUSCL scheme (default F) 500 ! ln_traadv_muscl2 MUSCL scheme + cen2 at boundaries (default F) 501 ! ln_traadv_ubs UBS scheme (default F) 502 ! ln_traadv_qck QUICKEST scheme (default F) 503 &nam_traadv 504 ln_traadv_cen2 = .false. 505 ln_traadv_tvd = .true. 506 ln_traadv_muscl = .false. 507 ln_traadv_muscl2 = .false. 508 ln_traadv_ubs = .false. 509 ln_traadv_qck = .false. 510 / 511 !----------------------------------------------------------------------- 512 ! nam_traldf lateral diffusion scheme for tracer 513 !----------------------------------------------------------------------- 514 ! Type of the operator : 515 ! ln_traldf_lap laplacian operator (default T) 516 ! ln_traldf_bilap bilaplacian operator (default F) 517 ! Direction of action : 518 ! ln_traldf_level iso-level (default F) 519 ! ln_traldf_hor horizontal (geopotential) (default F) (require "key_ldfslp" when ln_sco=T) 520 ! ln_traldf_iso iso-neutral (default T) (require "key_ldfslp") 521 ! Coefficient 522 ! aht0 horizontal eddy diffusivity for tracers (m2/s) 523 ! ahtb0 background eddy diffusivity for isopycnal diffusion (m2/s) 524 ! aeiv0 eddy induced velocity coefficient (m2/s) (require "key_traldf_eiv") 525 &nam_traldf 526 ln_traldf_lap = .true. 527 ln_traldf_bilap = .false. 528 ln_traldf_level = .false. 529 ln_traldf_hor = .false. 530 ln_traldf_iso = .true. 531 aht0 = 1000. 532 ahtb0 = 0. 533 aeiv0 = 0. 534 / 535 !----------------------------------------------------------------------- 536 ! namtdp tracer newtonian damping ('key_tradmp') 537 !----------------------------------------------------------------------- 538 ! ndmp type of damping in temperature and salinity 539 ! ! ='latitude', damping poleward of 'ndmp' degrees and function 540 ! ! of the distance-to-coast. Red and Med Seas as ndmp=-1 541 ! ! =-1 damping only in Med and Red Seas 542 ! ndmpf =1 create a damping.coeff NetCDF file (the 3D damping array) 543 ! nmldmp type of damping in the mixed layer 544 ! ! =0 damping throughout the water column 545 ! ! =1 no damping in the mixed layer defined by avt >5cm2/s 546 ! ! =2 no damping in the mixed layer defined rho<rho(surf)+.01 547 ! sdmp surface time scale for internal damping (days) 548 ! bdmp bottom time scale for internal damping (days) 549 ! hdmp depth of transition between sdmp and bdmp (meters) 550 &namtdp 551 ndmp = 90 552 ndmpf = 1 553 nmldmp = 1 554 sdmp = 50. 555 bdmp = 360. 556 hdmp = 800. 557 / 558 559 !!====================================================================== 560 !! *** Dynamics namelists *** 561 !!====================================================================== 562 !! nam_dynadv formulation of the momentum advection 563 !! nam_dynvor advection scheme 564 !! nam_dynhpg hydrostatic pressure gradient 565 !! namflg hydrostatic pressure gradient time stepping 566 !! nam_dynspg surface pressure gradient (CPP key only) 567 !! nam_dynldf lateral diffusion scheme 568 !!====================================================================== 569 570 !----------------------------------------------------------------------- 571 ! nam_dynadv formulation of the momentum advection 572 !----------------------------------------------------------------------- 573 ! ln_dynadv_vec vector form (T) or flux form (F) (default T) 574 ! ln_dynadv_cen2 flux form - 2nd order centered scheme (default T) 575 ! ln_dynadv_ubs flux form - 3rd order UBS scheme (default F) 576 &nam_dynadv 577 ln_dynadv_vec = .true. 578 ln_dynadv_cen2 = .false. 579 ln_dynadv_ubs = .false. 580 / 581 !----------------------------------------------------------------------- 582 ! nam_dynvor option of physics/algorithm (not control by CPP keys) 583 !----------------------------------------------------------------------- 584 ! ln_dynvor_ens vorticity trends: enstrophy conserving scheme (default T) 585 ! ln_dynvor_ene " " : energy conserving scheme (default F) 586 ! ln_dynvor_mix " " : mixed scheme (default F) 587 ! ln_dynvor_een " " : energy & enstrophy scheme (default F) 588 &nam_dynvor 589 ln_dynvor_ene = .false. 590 ln_dynvor_ens = .false. 591 ln_dynvor_mix = .false. 592 ln_dynvor_een = .true. 593 / 594 !----------------------------------------------------------------------- 595 ! nam_dynhpg Hydrostatic pressure gradient option 596 !----------------------------------------------------------------------- 597 ! type of pressure gradient scheme (choose one only!) 598 ! ln_hpg_zco z-coordinate - full steps (default T) 599 ! ln_hpg_zps z-coordinate - partial steps (interpolation) 600 ! ln_hpg_sco s-coordinate (standard jacobian formulation) 601 ! ln_hpg_hel s-coordinate (helsinki modification) 602 ! ln_hpg_wdj s-coordinate (weighted density jacobian) 603 ! ln_hpg_djc s-coordinate (Density Jacobian with Cubic polynomial) 604 ! ln_hpg_rot s-coordinate (ROTated axes scheme) 605 ! parameters 606 ! gamm weighting coefficient (wdj scheme) 607 &nam_dynhpg 608 ln_hpg_zco = .false. 609 ln_hpg_zps = .true. 610 ln_hpg_sco = .false. 611 ln_hpg_hel = .false. 612 ln_hpg_wdj = .false. 613 ln_hpg_djc = .false. 614 ln_hpg_rot = .false. 615 gamm = 0.e0 616 / 617 !----------------------------------------------------------------------- 618 ! namflg algorithm flags (algorithm not control by CPP keys) 619 !----------------------------------------------------------------------- 620 ! ln_dynhpg_imp hydrostatic pressure gradient: semi-implicit time scheme (T) 621 ! centered time scheme (F) 622 ! nn_dynhpg_rst add dynhpg implicit variables in restart ot not (1/0) 623 &namflg 624 ln_dynhpg_imp = .false. 625 nn_dynhpg_rst = 0 626 / 627 !----------------------------------------------------------------------- 628 ! nam_dynspg surface pressure gradient (CPP key only) 629 !----------------------------------------------------------------------- 630 ! ! explicit free surface ("key_dynspg_exp") 631 ! ! filtered free surface ("key_dynspg_flt") 632 ! ! split-explicit free surface ("key_dynspg_ts") 633 ! ! rigid-lid ("key_dynspg_rl") 634 635 !----------------------------------------------------------------------- 636 ! nam_dynldf lateral diffusion on momentum 637 !----------------------------------------------------------------------- 638 ! Type of the operator : 639 ! ln_dynldf_lap laplacian operator (default T) 640 ! ln_dynldf_bilap bilaplacian operator (default F) 641 ! Direction of action : 642 ! ln_dynldf_level iso-level (default F) 643 ! ln_dynldf_hor horizontal (geopotential) (default F) (require "key_ldfslp" if ln_sco=T) 644 ! ln_dynldf_iso iso-neutral (default T) (require "key_ldfslp") 645 ! Coefficient 646 ! ahm0 horizontal eddy viscosity for the dynamics (m2/s) 647 ! ahmb0 background eddy viscosity for isopycnal diffusion (m2/s) 648 &nam_dynldf 649 ln_dynldf_lap = .false. 650 ln_dynldf_bilap = .true. 651 ln_dynldf_level = .false. 652 ln_dynldf_hor = .true. 653 ln_dynldf_iso = .false. 654 ahm0 = -8.5e+11 655 ahmb0 = 0. 656 / 657 658 !!====================================================================== 659 !! Tracers & Dynamics vertical physics namelists 660 !!====================================================================== 661 !! namzdf vertical physics 662 !! namnpc non penetrative convection 663 !! namric richardson number dependent vertical mixing ("key_zdfric" ) 664 !! namtke TKE dependent vertical mixing ("key_zdftke" ) 665 !! namkpp KPP dependent vertical mixing ("key_zdfkpp" ) 666 !! namddm double diffusive mixing parameterization ("key_zdfddm" ) 667 !!====================================================================== 668 669 !----------------------------------------------------------------------- 670 ! namzdf vertical physics 671 !----------------------------------------------------------------------- 672 ! ln_zdfevd enhanced vertical diffusion (default T) 673 ! ln_zdfnpc Non-Penetrative Convection (default T) 674 ! avm0 vertical eddy viscosity for the dynamic (m2/s) (also background Kz if not "key_zdfcst") 675 ! avt0 vertical eddy diffusivity for tracers (m2/s) (also background Kz if not "key_zdfcst") 676 ! avevd vertical coefficient for enhanced diffusion scheme (m2/s) 677 ! n_evdm = 0 apply enhanced mixing on tracer only 678 ! = 1 apply enhanced mixing on both tracer and momentum 679 ! ln_zdfexp vertical physics: (=T) time splitting (T) (Default=F) 680 ! (=F) euler backward (F) 681 ! n_zdfexp number of sub-timestep for time splitting scheme 682 &namzdf 683 ln_zdfevd = .true. 684 ln_zdfnpc = .false. 685 avm0 = 1.2e-4 686 avt0 = 1.2e-5 687 avevd = 100. 688 n_evdm = 1 689 ln_zdfexp = .false. 690 n_zdfexp = 3 691 / 692 !----------------------------------------------------------------------- 693 ! namnpc non penetrative convection 694 !----------------------------------------------------------------------- 695 ! nnpc1 non penetrative convective scheme frequency 696 ! nnpc2 non penetrative convective scheme print frequency 697 &namnpc 698 nnpc1 = 1 699 nnpc2 = 365 700 / 701 !----------------------------------------------------------------------- 702 ! namric richardson number dependent vertical diffusion ("key_zdfric" ) 703 !----------------------------------------------------------------------- 704 ! avmri maximum value of the vertical viscosity 705 ! alp coefficient of the parameterization 706 ! nric coefficient of the parameterization 707 &namric 708 avmri = 100.e-4 709 alp = 5. 710 nric = 2 711 / 712 !----------------------------------------------------------------------- 713 ! namtke turbulent eddy kinetic dependent vertical diffusion ("key_zdftke") 714 !----------------------------------------------------------------------- 715 ! ln_rstke flag to restart with tke from a run without tke (default F) 716 ! ediff coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) ) 717 ! ediss coef. of the Kolmogoroff dissipation 718 ! ebb coef. of the surface input of tke 719 ! efave coef. to applied to the tke diffusion ( avtke=efave*avm ) 720 ! emin minimum value of tke (m^2/s^2) 721 ! emin0 surface minimum value of tke (m^2/s^2) 722 ! nitke number of restart iterative loops 723 ! ri_c critic richardson number = 2/9 = 0.22222222 (hard coded) 724 ! nmxl length used = 0 bounded by the distance to surface and bottom 725 ! ! = 1 bounded by the local vertical scale factor 726 ! ! = 2 first vertical derivative of mixing length bounded by 1 727 ! npdl Prandtl number = 0 no vertical prandtl number (avt=avm) 728 ! ! = 1 prandtl number function of richarson number (avt=pdl*avm) 729 ! ! = 2 same as = 1 but a shapiro filter is applied on pdl 730 ! nave = horizontal averaged (=1) or not (=0) of avt (default =1) 731 ! navb = 0 cst background avt0, avm0 / =1 profile used on avtb 732 &namtke 733 ln_rstke = .false. 734 ediff = 0.1 735 ediss = 0.7 736 ebb = 60. 737 efave = 1. 738 emin = 1.e-6 739 emin0 = 1.e-4 740 nitke = 50 741 nmxl = 2 742 npdl = 1 743 navb = 0 744 / 745 !----------------------------------------------------------------------- 746 ! namkpp K-Profile Parameterization dependent vertical mixing ("key_zdfkpp" ) 747 !----------------------------------------------------------------------- 748 ! ln_kpprimix shear instability mixing (default T) 749 ! difmiw constant internal wave viscosity (m2/s) 750 ! difsiw constant internal wave diffusivity (m2/s) 751 ! Riinfty local Richardson Number limit for shear instability 752 ! difri maximum shear mixing at Rig = 0 (m2/s) 753 ! bvsqcon Brunt-Vaisala squared (1/s**2) for maximum convection 754 ! difcon maximum mixing in interior convection (m2/s) 755 ! nave = 0/1 flag for horizontal average on avt, avmu, avmv 756 ! navb = 0/1 flag for constant or profile background avt 757 &namkpp 758 ln_kpprimix = .true. 759 difmiw = 1.e-04 760 difsiw = 0.1e-04 761 Riinfty = 0.8 762 difri = 0.0050 763 bvsqcon = -0.01e-07 764 difcon = 1. 765 navb = 0 766 nave = 1 767 / 768 !----------------------------------------------------------------------- 769 ! namddm double diffusive mixing parameterization ("key_zdfddm") 770 !----------------------------------------------------------------------- 771 ! avts maximum avs for dd mixing 772 ! hsbfr heat/salt buoyancy flux ratio 773 &namddm 774 avts = 1.e-4 775 hsbfr = 1.6 776 / 777 778 !!====================================================================== 779 !! *** Miscelaneous namelists *** 780 !!====================================================================== 781 !! namsol elliptic solver / island / free surface 782 !!====================================================================== 783 784 !----------------------------------------------------------------------- 785 ! namsol elliptic solver / island / free surface 786 !----------------------------------------------------------------------- 787 ! nsolv elliptic solver (=1 preconditioned conjugate gradient: pcg) 788 ! (=2 successive-over-relaxation: sor) 789 ! (=3 FETI currently it does not work! ("key_feti") 790 ! nsol_arp absolute/relative (0/1) precision convergence test 791 ! nmin minimum of iterations for the SOR solver 792 ! nmax maximum of iterations for the SOR solver 793 ! nmod frequency of test for the SOR solver 794 ! eps absolute precision of the solver 795 ! resmax absolute precision for the SOR solver 796 ! sor optimal coefficient for SOR solver 797 ! epsisl absolute precision on stream function solver 798 ! nmisl maximum pcg iterations for island 799 ! rnu strength of the additional force used in free surface b.c. 800 &namsol 801 nsolv = 1 802 nsol_arp = 0 803 nmin = 300 804 nmax = 800 805 nmod = 10 806 eps = 1.E-6 807 resmax = 1.E-10 808 sor = 1.92 809 epsisl = 1.e-10 810 nmisl = 4000 811 rnu = 1. 812 / 813 !!====================================================================== 814 !! *** Diagnostics namelists *** 815 !!====================================================================== 816 !! namtrd dynamics and/or tracer trends ("key_trddyn","key_trdtra","key_trdmld") 817 !! namgap level mean model-data gap ("key_diagap") 818 !! namspr surface pressure diagnosed in rigid-lid ("key_diaspr") 819 !! namflo float parameters ("key_float") 820 !! namptr Poleward Transport Diagnostics 821 !!====================================================================== 822 823 !----------------------------------------------------------------------- 824 ! namtrd diagnostics on dynamics and/or tracer trends ("key_trddyn" and/or "key_trdtra") 825 ! or mixed-layer trends ('key_trdmld') 826 ! or barotropic vorticity ("key_trdvor") 827 !----------------------------------------------------------------------- 828 ! ntrd time step frequency dynamics and tracers trends 829 ! nctls control surface type in mixed-layer trends (0,1 or n<jpk) 830 ! ln_trdmld_restart restart for ML diagnostics 831 ! ucf unit conversion factor (=1 -> /seconds | =86400. -> /day) 832 ! ln_trdmld_instant flag to diagnose trends of instantantaneous or mean ML T/S 833 &namtrd 834 ntrd = 365 835 nctls = 0 836 ln_trdmld_restart = .false. 837 ucf = 1. 838 ln_trdmld_instant = .false. 839 / 840 !----------------------------------------------------------------------- 841 ! namgap level mean model-data gap ('key_diagap') 842 !----------------------------------------------------------------------- 843 ! ngap time-step frequency of model-data gap computation 844 ! nprg time-step frequency of gap print in model output 845 &namgap 846 ngap = 15 847 nprg = 10 848 / 849 !----------------------------------------------------------------------- 850 ! namspr surface pressure diagnostic 851 !----------------------------------------------------------------------- 852 ! nmaxp maximum of iterations for the solver 853 ! epsp absolute precision of the solver 854 ! niterp number of iteration done by the solver 855 &namspr 856 nmaxp = 1000 857 epsp = 1.e-3 858 niterp = 400 859 / 860 !----------------------------------------------------------------------- 861 ! namflo float parameters ("key_float") 643 862 !----------------------------------------------------------------------- 644 863 ! ln_rstflo boolean term for float restart (true or false) … … 655 874 ln_flork4 = .false. 656 875 / 657 !!====================================================================== 658 !! Surface Boundary Condition namelists 659 !! 660 !! namsbc surface boundary condition 661 !! namsbc_ana analytical formulation of sbc 662 !! namsbc_flx flux formulation of sbc 663 !! namsbc_core CORE bulk formulea formulation of sbc 664 !! namsbc_clio CLIO bulk formulea formulation of sbc 665 !! namsbc_cpl CouPLed formulation of sbc 666 !! namsbc_ssr sea surface restoring term (for T and/or S) 667 !!====================================================================== 668 !----------------------------------------------------------------------- 669 ! namsbc Surface Boundary Condition (surface module) 670 !----------------------------------------------------------------------- 671 ! nn_fsbc frequency of surface boundary condition computation 672 ! (also correspond to the frequency of sea-ice model call) 673 ! 674 ! ln_ana analytical formulation (fill namsbc_ana ) 675 ! ln_flx flux formulation (fill namsbc_flx ) (overwritten by key_sbc_flux ) 676 ! ln_blk_clio CLIO bulk formulation (fill namsbc_core) (overwritten by key_sbc_bulk_clio) 677 ! ln_blk_core CORE bulk formulation (fill namsbc_clio) (overwritten by key_sbc_bulk_core) 678 ! ln_cpl Coupled formulation (fill namsbc_cpl ) (overwritten by key_sbc_coupled ) 679 ! 680 ! nn_ice =0 no ice boundary condition , 681 ! =1 observed ice-cover , 682 ! =2 ice-model ice-cover (overwritten by key_ice_lim) 683 ! 684 ! ln_dm2dc Daily mean to Diurnal Cycle short wave (qsr) 685 ! 686 ! ln_ssr Sea Surface restoring on SST and/or SSS (fill namsbc_ssr) 687 ! (overwritten by key_dta_sst and/or key_dtasss) 688 ! 689 ! nn_fwb FreshWater Budget: =0 no control , 690 ! =1 annual global mean of emp set to zero , 691 ! =2 global mean of emp set to zero at each nn_fsbc time step 692 &namsbc 693 nn_fsbc = 5 694 ln_ana = .false. 695 ln_flx = .false. 696 ln_blk_clio = .false. 697 ln_blk_core = .true. 698 ln_cpl = .false. 699 nn_ice = 2 700 ln_dm2dc = .false. 701 ln_rnf = .false. 702 ln_ssr = .false. 703 nn_fwb = 0 704 / 705 !----------------------------------------------------------------------- 706 ! namsbc_rnf runoffs namelist surface boundary condition 707 !----------------------------------------------------------------------- 708 !nn_runoff =0 no, 1 runoff, 2 runoff+river mouth ups adv 709 !rn_hrnf runoffs, depth over which enhanced vertical mixing is used 710 !rn_avt_rnf runoffs, value of the additional vertical mixing coef. [m2/s] 711 !cn_dir Root directory for location of ssr files 712 !sn_rnf informations about the runoff file to be read 713 !sn_cnf informations about the runoff mouth file to be read 714 ! THE ORDER OF THE FILES MATTER: 715 ! ! file name ! frequency (hours) ! variable ! time interpolation ! 716 ! ! ! (if <0 months) ! name ! (logical) ! 717 &namsbc_rnf 718 nn_runoff = 2 719 rn_hrnf = 0.e0 720 rn_avt_rnf = 1.e-3 721 cn_dir = './' 722 sn_rnf = 'runoff_1m_nomask.nc' , -12. , 'sorunoff', .TRUE. , 1 , 0 723 sn_cnf = 'runoff_1m_nomask.nc' , 0. , 'socoefr' , .FALSE. , 1 , 0 724 / 725 !----------------------------------------------------------------------- 726 ! namsbc_ana analytical surface boundary condition 727 !----------------------------------------------------------------------- 728 ! nn_tau000 gently increase the stress over the first ntau_rst time-steps 729 ! rn_utau0 uniform value used as default surface heat flux 730 ! rn_vtau0 uniform value used as default solar radiation flux 731 ! rn_q0 uniform value used as default surface heat flux 732 ! rn_qsr0 uniform value used as default solar radiation flux 733 ! rn_emp0 uniform value used as default surface freswater budget (E-P) 734 &namsbc_ana 735 nn_tau000 = 0 736 rn_utau0 = 0.5 737 rn_vtau0 = 0.e0 738 rn_q0 = 0.e0 739 rn_qsr0 = 0.e0 740 rn_emp0 = 0.e0 741 / 742 !----------------------------------------------------------------------- 743 ! namsbc_flx surface boundary condition : flux formulation (#ifdef "key_sbc_flux") 744 !----------------------------------------------------------------------- 745 ! cn_dir directory in which the model is executed 746 ! 747 ! THE ORDER OF THE FILES MATTER: 748 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 749 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 750 &namsbc_flx 751 cn_dir = './' 752 sn_utau = 'utau.nc' , 24. , 'utau' , .FALSE. , 0 , 0 753 sn_vtau = 'vtau.nc' , 24. , 'vtau' , .FALSE. , 0 , 0 754 sn_qtot = 'qtot.nc' , 24. , 'qtot' , .FALSE. , 0 , 0 755 sn_qsr = 'qsr.nc' , 24. , 'qsr' , .FALSE. , 0 , 0 756 sn_emp = 'emp.nc' , 24. , 'emp' , .FALSE. , 0 , 0 757 758 / 759 !----------------------------------------------------------------------- 760 ! namsbc_clio CLIO bulk formulea 761 !----------------------------------------------------------------------- 762 ! cn_dir directory in which the model is executed 763 ! 764 ! THE ORDER OF THE FILES MATTER: 765 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 766 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 767 &namsbc_clio 768 cn_dir = './' 769 sn_utau = 'taux_1m' , -12. , 'sozotaux', .FALSE. , 1 , 0 770 sn_vtau = 'tauy_1m' , -12. , 'sometauy', .FALSE. , 1 , 0 771 sn_wndm = 'flx' , -12. , 'socliowi', .FALSE. , 1 , 0 772 sn_tair = 'flx' , -12. , 'socliot1', .FALSE. , 1 , 0 773 sn_humi = 'flx' , -12. , 'socliohu', .FALSE. , 1 , 0 774 sn_ccov = 'flx' , -12. , 'socliocl', .FALSE. , 1 , 0 775 sn_prec = 'flx' , -12. , 'socliopl', .FALSE. , 1 , 0 776 / 777 !----------------------------------------------------------------------- 778 ! namsbc_core CORE bulk formulea 779 !----------------------------------------------------------------------- 780 ! cn_dir directory in which the model is executed 781 ! ln_2m logical flag to use air temp. and hum referenced at 2m instead 10m 782 ! alpha_precip multiplication factor for precipitation (total & snow) 783 ! 784 ! THE ORDER OF THE FILES MATTER: 785 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 786 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 787 &namsbc_core 788 cn_dir = './' 789 ln_2m = .FALSE. 790 alpha_precip = 1. 791 sn_wndi = 'u10' , 24. , 'U_10_MOD', .FALSE. , 1 , 0 792 sn_wndj = 'v10' , 24. , 'V_10_MOD', .FALSE. , 1 , 0 793 sn_qsr = 'rad' , 24. , 'SWDN_MOD', .FALSE. , 1 , 0 794 sn_qlw = 'rad' , 24. , 'LWDN_MOD', .FALSE. , 1 , 0 795 sn_tair = 't10' , 24. , 'T_10_MOD', .FALSE. , 1 , 0 796 sn_humi = 'q10' , 24. , 'Q_10_MOD', .FALSE. , 1 , 0 797 sn_prec = 'precip' , -12. , 'RAIN' , .FALSE. , 1 , 0 798 sn_snow = 'precip' , -12. , 'SNOW' , .FALSE. , 1 , 0 799 / 800 !----------------------------------------------------------------------- 801 ! namsbc_cpl coupled ocean/atmosphere model (#ifdef "key_coupled") 802 !----------------------------------------------------------------------- 803 &namsbc_cpl 804 / 805 !----------------------------------------------------------------------- 806 ! namsbc_ssr surface boundary condition : sea surface restoring 807 !----------------------------------------------------------------------- 808 ! nn_sstr =0/1 add a retroaction term in the surface heat flux 809 ! nn_sssr =O/1/2 add a damping term in the surface freshwater flux 810 ! (=1, salt flux, concentration/dillution emps only) 811 ! (=2, volume flux, both emp and emps are updated) 812 ! dqdt magnitude of the retroaction on temperature [W/m2/K] 813 ! deds magnitude of the damping on salinity [mm/day] 814 ! THE ORDER OF THE FILES MATTER: 815 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 816 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 817 &namsbc_ssr 818 nn_sstr = 0 819 nn_sssr = 0 820 dqdt = -40. 821 deds = 27.7 822 sn_sst = 'sst_data.nc' , 24. , 'sst' , .FALSE. , 0 , 0 823 sn_sss = 'sss_data.nc' , -12. , 'sss' , .TRUE. , 0 , 0 824 / 825 !----------------------------------------------------------------------- 826 ! namagrif agrif parameters (#ifdef key_agrif) 827 !----------------------------------------------------------------------- 828 ! nbclineupdate baroclinic update frequency 829 ! ln_spc_dyn use 0 as special value for dynamics 830 ! visc_tra viscosity coeeficient for tracers sponge layer 831 ! visc_dyn viscosity coeeficient for dynamics sponge layer 832 &namagrif 833 nbclineupdate = 3 834 ln_spc_dyn = .TRUE. 835 visc_tra = 2880. 836 visc_dyn = 2880. 837 / 876 !----------------------------------------------------------------------- 877 ! namptr Poleward Transport Diagnostic 878 !----------------------------------------------------------------------- 879 ! ln_diaptr logical flag for Poleward transport computation 880 ! ln_subbas logical flag for Atlantic/Pacific/Indian basins computation 881 ! need input basins mask file named "subbasins.nc" 882 ! nf_ptr Frequency of computation 883 &namptr 884 ln_diaptr = .false. 885 ln_subbas = .false. 886 nf_ptr = 15 887 / -
branches/dev_003_CPL/CONFIG/ORCA2_LIM/EXP00/namelist
r889 r990 1 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 ! OPA namelist : model option and parameter input 3 ! ------------- 4 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 5 ! 1 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 !! NEMO/OPA : 1 - run manager (namrun, nam_ctl, nam_mpp, nam_mpp_dyndist, nam_ctl) 3 !! namelists 2 - Domain (nam_zgr, nam_zgr_sco, namdom) 4 !! 3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core 5 !! namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, namalb) 6 !! 4 - lateral boundary (namlbc, namcla, namobc) 7 !! 5 - bottom boundary (nambfr, nambbc, nambbl) 8 !! 6 - Tracer (nameos, nam_traadv, nam_traldf, namtdp) 9 !! 7 - dynamics (nam_dynadv, nam_dynvor, nam_dynhpg, namflg, nam_dynspg, nam_dynldf) 10 !! 8 - Verical physics (namzdf, namnpc, namric, namtke, namkpp, namddm) 11 !! 9 - diagnostics (namtrd, namgap, namspr, namflo, namptr) 12 !! 9 - miscellaneous (namsol) 13 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 14 ! CAUTION: some scripts does not support CAPITALs for logical definition 15 ! ******* use .true. or .false. and NOT .TRUE. or .FALSE. 16 17 !!====================================================================== 18 !! *** Run management namelists *** 19 !!====================================================================== 20 !! namrun parameters of the run 21 !! nam_mpp Massively Parallel Processing 22 !! nam_mpp_dyndist 23 !! nam_ctl Control prints & Benchmark 24 !!====================================================================== 25 6 26 !----------------------------------------------------------------------- 7 27 ! namrun parameters of the run … … 10 30 ! cexper experience name for vairmer format 11 31 ! ln_rstart boolean term for restart (true or false) 12 ! nrstdt control of the restart timestep: 13 ! = 0 restart, do not control nit000 in the restart file. 14 ! = 1 restart, control nit000 in the restart file. Do not 15 ! use the date in the restart file (use ndate0 in namelist) 16 ! = 2 restart, control nit000 in the restart file, use the date 17 ! in the restart file. ndate0 in the namelist is ignored. 32 ! nrstdt restart control = 0 restart, do not control nit000 in the restart file. 33 ! ! = 1 restart, control nit000 in the restart file. Do not 34 ! ! use the date in the restart file (use ndate0 in namelist) 35 ! ! = 2 restart, control nit000 in the restart file, use the date 36 ! ! in the restart file. ndate0 in the namelist is ignored. 18 37 ! nit000 number of the first time step 19 38 ! nitend number of the last time step … … 24 43 ! nwrite frequency of OUTPUT file 25 44 ! ln_dimgnnn (F/T) 1 DIMG file - (for all proc/per proc) 26 !27 ! CAUTION: for usual run scripts, logical value of28 ! ******* ln_rstart must be .true. or .false.29 ! and NOT .TRUE. or .FALSE.30 45 &namrun 31 46 no = 0 … … 38 53 nleapy = 0 39 54 ninist = 0 40 nstock = 5475 41 nwrite = 5475 42 ln_dimgnnn = .FALSE. 55 nstock = 5475 56 nwrite = 5475 57 ln_dimgnnn = .false. 58 / 59 !----------------------------------------------------------------------- 60 ! nam_mpp Massively Parallel Processing 61 !----------------------------------------------------------------------- 62 ! c_mpi_send mpi send/recieve type 63 ! = 'S' : standard blocking send 64 ! = 'B' : buffer blocking send 65 ! = 'I' : immediate non-blocking send 66 &nam_mpp 67 c_mpi_send = 'S' 68 / 69 !----------------------------------------------------------------------- 70 ! nam_mpp_dyndist Massively Parallel Distribution ("key_agrif" && "key_mpp_dyndist") 71 !----------------------------------------------------------------------- 72 ! jpni number of processors following i 73 ! jpnj number of processors following j 74 ! jpnij number of local domains 75 &nam_mpp_dyndist 76 jpni = 1 77 jpnj = 1 78 jpnij = 1 43 79 / 44 80 !----------------------------------------------------------------------- … … 69 105 nbit_cmp = 0 70 106 / 71 !----------------------------------------------------------------------- 72 ! nam_mpp Massively Parallel Processing 73 !----------------------------------------------------------------------- 74 ! c_mpi_send mpi send/recieve type 75 ! = 'S' : standard blocking send 76 ! = 'B' : buffer blocking send 77 ! = 'I' : immediate non-blocking send 78 &nam_mpp 79 c_mpi_send = 'S' 80 / 81 !----------------------------------------------------------------------- 82 ! nam_mpp_dyndist Massively Parallel Distribution 83 ! #ifdef ( key_agrif && key_mpp_dyndist ) 84 !----------------------------------------------------------------------- 85 ! jpni number of processors following i 86 ! jpnj number of processors following j 87 ! jpnij number of local domains 88 &nam_mpp_dyndist 89 jpni = 1 90 jpnj = 1 91 jpnij = 1 92 / 107 108 !!====================================================================== 109 !! *** Domain namelists *** 110 !!====================================================================== 111 !! nam_zgr vertical coordinate 112 !! nam_zgr_sco s-coordinate or hybrid z-s-coordinate 113 !! namdom space and time domain (bathymetry, mesh, timestep) 114 !!====================================================================== 115 93 116 !----------------------------------------------------------------------- 94 117 ! nam_zgr vertical coordinate 95 118 !----------------------------------------------------------------------- 96 ! ln_zco z-coordinate - full steps (T/F) 119 ! ln_zco z-coordinate - full steps (T/F) ("key_zco" may also be defined) 97 120 ! ln_zps z-coordinate - partial steps (T/F) 98 121 ! ln_sco s- or hybrid z-s-coordinate (T/F) … … 118 141 / 119 142 !----------------------------------------------------------------------- 120 ! nam_traadv advection scheme for tracer (option not control by CPP keys)121 !-----------------------------------------------------------------------122 ! ln_traadv_cen2 2nd order centered scheme (default T)123 ! ln_traadv_tvd TVD scheme (default F)124 ! ln_traadv_muscl MUSCL scheme (default F)125 ! ln_traadv_muscl2 MUSCL2 scheme (default F)126 ! ln_traadv_ubs UBS scheme (default F)127 ! ln_traadv_qck QUICKEST scheme (default F)128 &nam_traadv129 ln_traadv_cen2 = .true.130 ln_traadv_tvd = .false.131 ln_traadv_muscl = .false.132 ln_traadv_muscl2 = .false.133 ln_traadv_ubs = .false.134 ln_traadv_qck = .false.135 /136 !-----------------------------------------------------------------------137 ! nam_traldf lateral diffusion scheme for tracer (option not control by CPP keys)138 !-----------------------------------------------------------------------139 ! Type of the operator :140 ! ln_traldf_lap laplacian operator (default T)141 ! ln_traldf_bilap bilaplacian operator (default F)142 ! Direction of action :143 ! ln_traldf_level iso-level (default F)144 ! ln_traldf_hor horizontal (geopotential) (default F)^**145 ! ln_traldf_iso iso-neutral (default T)^*146 ! Coefficient147 ! aht0 horizontal eddy diffusivity for tracers (m2/s)148 ! ahtb0 background eddy diffusivity for isopycnal diffusion (m2/s)149 ! aeiv0 eddy induced velocity coefficient (m2/s)150 ! ^* require key_ldfslp to compute the direction of the lateral diffusion151 ! ^** require key_ldfslp in s-coordinate152 &nam_traldf153 ln_traldf_lap = .true.154 ln_traldf_bilap = .false.155 ln_traldf_level = .false.156 ln_traldf_hor = .false.157 ln_traldf_iso = .true.158 aht0 = 2000.159 ahtb0 = 0.160 aeiv0 = 2000.161 /162 !-----------------------------------------------------------------------163 ! nam_dynldf lateral diffusion on momentum164 !-----------------------------------------------------------------------165 ! Type of the operator :166 ! ln_dynldf_lap laplacian operator (default T)167 ! ln_dynldf_bilap bilaplacian operator (default F)168 ! Direction of action :169 ! ln_dynldf_level iso-level (default F)170 ! ln_dynldf_hor horizontal (geopotential) (default F)^**171 ! ln_dynldf_iso iso-neutral (default T)^*172 ! Coefficient173 ! ahm0 horizontal eddy viscosity for the dynamics (m2/s)174 ! ahmb0 background eddy viscosity for isopycnal diffusion (m2/s)175 &nam_dynldf176 ln_dynldf_lap = .true.177 ln_dynldf_bilap = .false.178 ln_dynldf_level = .false.179 ln_dynldf_hor = .true.180 ln_dynldf_iso = .false.181 ahm0 = 40000.182 ahmb0 = 0.183 /184 !-----------------------------------------------------------------------185 ! namflg algorithm flags (algorithm not control by CPP keys)186 !-----------------------------------------------------------------------187 ! ln_dynhpg_imp hydrostatic pressure gradient: semi-implicit time scheme (T)188 ! centered time scheme (F)189 ! nn_dynhpg_rst add dynhpg implicit variables in restart ot not (1/0)190 &namflg191 ln_dynhpg_imp = .false.192 nn_dynhpg_rst = 0193 /194 !-----------------------------------------------------------------------195 ! nam_dynhpg Hydrostatic pressure gradient option196 !-----------------------------------------------------------------------197 ! type of pressure gradient scheme (choose one only!)198 ! ln_hpg_zco z-coordinate - full steps (default T)199 ! ln_hpg_zps z-coordinate - partial steps (interpolation)200 ! ln_hpg_sco s-coordinate (standard jacobian formulation)201 ! ln_hpg_hel s-coordinate (helsinki modification)202 ! ln_hpg_wdj s-coordinate (weighted density jacobian)203 ! ln_hpg_djc s-coordinate (Density Jacobian with Cubic polynomial)204 ! ln_hpg_rot s-coordinate (ROTated axes scheme)205 ! parameters206 ! gamm weighting coefficient (wdj scheme)207 &nam_dynhpg208 ln_hpg_zco = .false.209 ln_hpg_zps = .true.210 ln_hpg_sco = .false.211 ln_hpg_hel = .false.212 ln_hpg_wdj = .false.213 ln_hpg_djc = .false.214 ln_hpg_rot = .false.215 gamm = 0.e0216 /217 !-----------------------------------------------------------------------218 ! nam_dynadv option of physics/algorithm (not control by CPP keys)219 !-----------------------------------------------------------------------220 ! ln_dynadv_vec vector form flag221 ! ln_dynadv_cen2 flux form - 2nd order centered scheme (default T)222 ! ln_dynadv_ubs flux form - 3rd order UBS scheme (default F)223 &nam_dynadv224 ln_dynadv_vec = .TRUE.225 ln_dynadv_cen2 = .FALSE.226 ln_dynadv_ubs = .FALSE.227 /228 !-----------------------------------------------------------------------229 ! nam_dynvor option of physics/algorithm (not control by CPP keys)230 !-----------------------------------------------------------------------231 ! ln_dynvor_ens vorticity trends: enstrophy conserving scheme (default T)232 ! ln_dynvor_ene " " : energy conserving scheme (default F)233 ! ln_dynvor_mix " " : mixed scheme (default F)234 ! ln_dynvor_een " " : energy & enstrophy scheme (default F)235 &nam_dynvor236 ln_dynvor_ene = .FALSE.237 ln_dynvor_ens = .TRUE.238 ln_dynvor_mix = .FALSE.239 ln_dynvor_een = .FALSE.240 /241 !-----------------------------------------------------------------------242 ! namalb albedo parameters243 !-----------------------------------------------------------------------244 ! cgren correction of the snow or ice albedo to take into account245 ! albice albedo of melting ice in the arctic and antarctic246 ! alphd coefficients for linear interpolation used to compute albedo247 ! between two extremes values (Pyane, 1972)248 ! alphc " "249 ! alphdi " "250 &namalb251 cgren = 0.06252 albice = 0.5253 alphd = 0.80254 alphc = 0.65255 alphdi = 0.72256 /257 !-----------------------------------------------------------------------258 143 ! namdom space and time domain (bathymetry, mesh, timestep) 259 144 !----------------------------------------------------------------------- 260 ! ntopo = 0/1 ,compute/read the bathymetry file 261 ! (mbathy, nb of T-ocean levels) 145 ! ntopo = 0/1 ,compute/read the bathymetry file (mbathy, nb of T-ocean levels) 262 146 ! e3zps_min the thickness of the partial step is set larger than the 263 ! e3zps_rat the minimum of e3zps_min and e3zps_rat * e3t 264 ! (N.B. 0<e3zps_rat<1) 147 ! e3zps_rat the minimum of e3zps_min and e3zps_rat * e3t (N.B. 0<e3zps_rat<1) 265 148 ! nmsh =1 create a mesh file (coordinates, scale factors, masks) 266 ! nacc the acceleration of convergence method 267 ! = 0, no acceleration, rdt = rdttra 268 ! = 1, acceleration used, rdt < rdttra(k) 149 ! nacc acceleration of convergence method = 0, not used, rdt = rdttra 150 ! = 1, used, rdt < rdttra(k) 269 151 ! atfp asselin time filter parameter 270 152 ! rdt time step for the dynamics (and tracer if nacc=0) … … 273 155 ! rdth depth variation of tracer time step 274 156 ! rdtbt barotropic time step (for the time splitting algorithm) 275 ! nclosea 0 no closed sea276 ! 1 closed sea (Black Sea, Caspian Sea, Great US Lakes...)157 ! nclosea =0 remove the closed sea from the global domain (orca configuration) 158 ! =1 closed sea (Black Sea, Caspian Sea, Great US Lakes...) 277 159 &namdom 278 160 ntopo = 1 … … 289 171 nclosea = 0 290 172 / 291 !----------------------------------------------------------------------- 292 ! namfwb freshwater budget correction 293 !----------------------------------------------------------------------- 294 ! ln_fwb logical flag for freshwater budget correction (0 annual mean) 295 &namfwb 296 ln_fwb = .true. 297 / 298 !----------------------------------------------------------------------- 299 ! namptr Poleward Transport Diagnostic 300 !----------------------------------------------------------------------- 301 ! ln_diaptr logical flag for Poleward transport computation 302 ! ln_subbas logical flag for Atlantic/Pacific/Indian basins computation 303 ! need input basins mask file named "subbasins.nc" 304 ! nf_ptr Frequency of computation 305 &namptr 306 ln_diaptr = .false. 307 ln_subbas = .false. 308 nf_ptr = 15 309 / 310 !----------------------------------------------------------------------- 311 ! namcro cross land advection 312 !----------------------------------------------------------------------- 313 ! n_cla advection between 2 ocean pts separates by land 314 &namcla 315 n_cla = 1 316 / 317 !----------------------------------------------------------------------- 318 ! namzdf vertical physics 319 !----------------------------------------------------------------------- 320 ! ln_zdfevd enhanced vertical diffusion (default T) 321 ! ln_zdfnpc Non-Penetrative Convection (default T) 322 ! avm0 vertical eddy viscosity for the dynamic (m2/s) 323 ! avt0 vertical eddy diffusivity for tracers (m2/s) 324 ! avevd vertical coefficient for enhanced diffusion scheme (m2/s) 325 ! nevdm = 0 apply enhanced mixing on tracer only 326 ! = 1 apply enhanced mixing on both tracer and momentum 327 ! ln_zdfexp vertical physics: (=T) time splitting (T) (Default=F) 328 ! (=F) euler backward (F) 329 ! n_zdfexp number of sub-timestep for time splitting scheme 330 &namzdf 331 ln_zdfevd = .true. 332 ln_zdfnpc = .false. 333 avm0 = 1.2e-4 334 avt0 = 1.2e-5 335 avevd = 100. 336 n_evdm = 1 337 ln_zdfexp = .false. 338 n_zdfexp = 3 339 / 340 !----------------------------------------------------------------------- 341 ! namnpc vnon penetrative convection 342 !----------------------------------------------------------------------- 343 ! nnpc1 non penetrative convective scheme frequency 344 ! nnpc2 non penetrative convective scheme print frequency 345 &namnpc 346 nnpc1 = 1 347 nnpc2 = 365 348 / 349 !----------------------------------------------------------------------- 350 ! nambbl bottom boundary layer scheme 351 !----------------------------------------------------------------------- 352 ! atrbbl lateral tracer coeff. for bottom boundary layer scheme(m2/s) 353 &nambbl 354 atrbbl = 10000. 355 / 356 !----------------------------------------------------------------------- 357 ! namric richardson number dependent vertical diffusion 358 ! ( #ifdef "key_zdfrichardson" ) 359 !----------------------------------------------------------------------- 360 ! avmri maximum value of the vertical viscosity 361 ! alp coefficient of the parameterization 362 ! nric coefficient of the parameterization 363 &namric 364 avmri = 100.e-4 365 alp = 5. 366 nric = 2 367 / 368 !----------------------------------------------------------------------- 369 ! namtke turbulent eddy kinetic dependent vertical diffusion 370 ! ( #ifdef "key_zdftke" ) 371 !----------------------------------------------------------------------- 372 ! ln_rstke flag to restart with tke from a run without tke (default F) 373 ! ediff coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) ) 374 ! ediss coef. of the Kolmogoroff dissipation 375 ! ebb coef. of the surface input of tke 376 ! efave coef. to applied to the tke diffusion ( avtke=efave*avm ) 377 ! emin minimum value of tke (m^2/s^2) 378 ! emin0 surface minimum value of tke (m^2/s^2) 379 ! nitke number of restart iterative loops 380 ! ri_c critic richardson number 381 ! nmxl flag on mixing length used 382 ! = 0 bounded by the distance to surface and bottom 383 ! = 1 bounded by the local vertical scale factor 384 ! = 2 first vertical derivative of mixing length bounded by 1 385 ! npdl flag on prandtl number 386 ! = 0 no vertical prandtl number (avt=avm) 387 ! = 1 prandtl number function of richarson number (avt=pdl*avm) 388 ! = 2 same as = 1 but a shapiro filter is applied on pdl 389 ! nave = horizontal averaged (=1) or not (=0) of avt (default =1) 390 ! navb = 0 cst background avt0, avm0 / =1 profile used on avtb 391 &namtke 392 ln_rstke = .false. 393 ediff = 0.1 394 ediss = 0.7 395 ebb = 3.75 396 efave = 1. 397 emin = 1.e-6 398 emin0 = 1.e-4 399 nitke = 50 400 nmxl = 2 401 npdl = 1 402 navb = 0 403 / 404 !----------------------------------------------------------------------- 405 ! namkpp K-Profile Parameterization dependent vertical diffusion 406 ! ( #ifdef "key_zdfkpp" ) 407 !----------------------------------------------------------------------- 408 ! ln_kpprimix shear instability mixing (default T) 409 ! difmiw constant internal wave viscosity (m2/s) 410 ! difsiw constant internal wave diffusivity (m2/s) 411 ! Riinfty local Richardson Number limit for shear instability 412 ! difri maximum shear mixing at Rig = 0 (m2/s) 413 ! bvsqcon Brunt-Vaisala squared (1/s**2) for maximum convection 414 ! difcon maximum mixing in interior convection (m2/s) 415 ! nave = 0/1 flag for horizontal average on avt, avmu, avmv 416 ! navb = 0/1 flag for constant or profile background avt 417 &namkpp 418 ln_kpprimix = .true. 419 difmiw = 1.e-04 420 difsiw = 0.1e-04 421 Riinfty = 0.8 422 difri = 0.0050 423 bvsqcon = -0.01e-07 424 difcon = 1. 425 navb = 0 426 nave = 1 427 / 428 !----------------------------------------------------------------------- 429 ! namddm double diffusive mixing parameterization 430 !----------------------------------------------------------------------- 431 ! avts maximum avs for dd mixing 432 ! hsbfr heat/salt buoyancy flux ratio 433 &namddm 434 avts = 1.e-4 435 hsbfr = 1.6 436 / 437 !----------------------------------------------------------------------- 438 ! namlbc lateral momentum boundary condition 439 !----------------------------------------------------------------------- 440 ! shlat lateral boundary condition on velocity 441 ! shlat = 0 , free slip 442 ! 0 < shlat < 2 , partial slip 443 ! shlat = 2 , no slip 444 ! 2 < shlat , strong slip 445 &namlbc 446 shlat = 2. 447 / 448 !----------------------------------------------------------------------- 449 ! nambfr bottom friction 450 !----------------------------------------------------------------------- 451 ! nbotfr type of bottom friction 452 ! nbotfr = 0 , no slip 453 ! nbotfr = 1 , linear friction 454 ! nbotfr = 2 , nonlinear friction 455 ! nbotfr = 3 , free slip 456 ! bfri1 bottom drag coefficient (linear case) 457 ! bfri2 bottom drag coefficient (non linear case) 458 ! bfeb2 bottom turbulent kinetic energy (m^2/s^2) 459 &nambfr 460 nbotfr = 1 461 bfri1 = 4.e-4 462 bfri2 = 1.e-3 463 bfeb2 = 2.5e-3 464 / 465 !----------------------------------------------------------------------- 466 ! nambbc bottom temperature boundary condition 467 !----------------------------------------------------------------------- 468 ! ngeo_flux = 0 no geothermal heat flux 469 ! = 1 constant geothermal heat flux 470 ! = 2 variable geothermal heat flux (read in geothermal_heating.nc) 471 ! ( C A U T I O N : flux in mW/m2 in the NetCDF file ) 472 ! ngeo_flux_const Constant value of geothermal heat flux (W/m2) 473 &nambbc 474 ngeo_flux = 2 475 ngeo_flux_const = 86.4e-3 173 174 !!====================================================================== 175 !! *** Surface Boundary Condition namelists *** 176 !!====================================================================== 177 !! namsbc surface boundary condition 178 !! namsbc_ana analytical formulation 179 !! namsbc_flx flux formulation 180 !! namsbc_clio CLIO bulk formulea formulation 181 !! namsbc_core CORE bulk formulea formulation 182 !! namsbc_cpl CouPLed formulation ("key_coupled") 183 !! namqsr penetrative solar radiation 184 !! namsbc_rnf river runoffs 185 !! namsbc_ssr sea surface restoring term (for T and/or S) 186 !! namalb albedo parameters 187 !!====================================================================== 188 189 !----------------------------------------------------------------------- 190 ! namsbc Surface Boundary Condition (surface module) 191 !----------------------------------------------------------------------- 192 ! nn_fsbc frequency of surface boundary condition computation 193 ! (= the frequency of sea-ice model call) 194 ! ln_ana analytical formulation (fill namsbc_ana ) 195 ! ln_flx flux formulation (fill namsbc_flx ) 196 ! ln_blk_clio CLIO bulk formulation (fill namsbc_core) 197 ! ln_blk_core CORE bulk formulation (fill namsbc_clio) 198 ! ln_cpl Coupled formulation (fill namsbc_cpl ) 199 ! nn_ice =0 no ice boundary condition , 200 ! =1 observed ice-cover , 201 ! =2 ice-model used ("key_lim3" or "key_lim2) 202 ! nn_ico_cpl ice-ocean coupling : =0 LIM-3 old case 203 ! =1 stresses computed using now ocean velocity 204 ! =2 combination of 0 and 1 cases 205 ! ln_dm2dc Daily mean to Diurnal Cycle short wave (qsr) 206 ! ln_ssr Sea Surface Restoring on T and/or S (fill namsbc_ssr) 207 ! nn_fwb FreshWater Budget: =0 no control , 208 ! =1 annual global mean of e-p-r set to zero , 209 ! =2 global mean of e-p-r set to zero at each nn_fsbc time step 210 &namsbc 211 nn_fsbc = 5 212 ln_ana = .false. 213 ln_flx = .false. 214 ln_blk_clio = .true. 215 ln_blk_core = .false. 216 ln_cpl = .false. 217 nn_ice = 2 218 nn_ico_cpl = 0 219 ln_dm2dc = .false. 220 ln_rnf = .false. 221 ln_ssr = .false. 222 nn_fwb = 0 223 / 224 !----------------------------------------------------------------------- 225 ! namsbc_ana analytical surface boundary condition 226 !----------------------------------------------------------------------- 227 ! nn_tau000 gently increase the stress over the first ntau_rst time-steps 228 ! rn_utau0 uniform value for the i-stress 229 ! rn_vtau0 uniform value for the j-stress 230 ! rn_q0 uniform value for the total heat flux 231 ! rn_qsr0 uniform value for the solar radiation 232 ! rn_emp0 uniform value for the freswater budget (E-P) 233 &namsbc_ana 234 nn_tau000 = 0 235 rn_utau0 = 0.5 236 rn_vtau0 = 0.e0 237 rn_q0 = 0.e0 238 rn_qsr0 = 0.e0 239 rn_emp0 = 0.e0 240 / 241 !----------------------------------------------------------------------- 242 ! namsbc_flx surface boundary condition : flux formulation (#ifdef "key_sbc_flux") 243 !----------------------------------------------------------------------- 244 ! cn_dir directory in which the model is executed 245 ! 246 ! THE ORDER OF THE FILES MATTER: 247 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 248 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 249 &namsbc_flx 250 cn_dir = './' 251 sn_utau = 'utau.nc' , 24. , 'utau' , .false. , 0 , 0 252 sn_vtau = 'vtau.nc' , 24. , 'vtau' , .false. , 0 , 0 253 sn_qtot = 'qtot.nc' , 24. , 'qtot' , .false. , 0 , 0 254 sn_qsr = 'qsr.nc' , 24. , 'qsr' , .false. , 0 , 0 255 sn_emp = 'emp.nc' , 24. , 'emp' , .false. , 0 , 0 256 / 257 !----------------------------------------------------------------------- 258 ! namsbc_clio CLIO bulk formulea 259 !----------------------------------------------------------------------- 260 ! cn_dir directory in which the model is executed 261 ! 262 ! THE ORDER OF THE FILES MATTER: 263 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 264 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 265 &namsbc_clio 266 cn_dir = './' 267 sn_utau = 'taux_1m' , -12. , 'sozotaux', .false. , 1 , 0 268 sn_vtau = 'tauy_1m' , -12. , 'sometauy', .false. , 1 , 0 269 sn_wndm = 'flx' , -12. , 'socliowi', .false. , 1 , 0 270 sn_tair = 'flx' , -12. , 'socliot1', .false. , 1 , 0 271 sn_humi = 'flx' , -12. , 'socliohu', .false. , 1 , 0 272 sn_ccov = 'flx' , -12. , 'socliocl', .false. , 1 , 0 273 sn_prec = 'flx' , -12. , 'socliopl', .false. , 1 , 0 274 / 275 !----------------------------------------------------------------------- 276 ! namsbc_core CORE bulk formulea 277 !----------------------------------------------------------------------- 278 ! cn_dir directory in which the model is executed 279 ! ln_2m logical flag to use air temp. and hum referenced at 2m instead 10m 280 ! alpha_precip multiplication factor for precipitation (total & snow) 281 ! 282 ! THE ORDER OF THE FILES MATTER: 283 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 284 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 285 &namsbc_core 286 cn_dir = './' 287 ln_2m = .false. 288 alpha_precip = 1. 289 sn_wndi = 'u10' , 24. , 'U_10_MOD', .false. , 1 , 0 290 sn_wndj = 'v10' , 24. , 'V_10_MOD', .false. , 1 , 0 291 sn_qsr = 'rad' , 24. , 'SWDN_MOD', .false. , 1 , 0 292 sn_qlw = 'rad' , 24. , 'LWDN_MOD', .false. , 1 , 0 293 sn_tair = 't10' , 24. , 'T_10_MOD', .false. , 1 , 0 294 sn_humi = 'q10' , 24. , 'Q_10_MOD', .false. , 1 , 0 295 sn_prec = 'precip' , -12. , 'RAIN' , .false. , 1 , 0 296 sn_snow = 'precip' , -12. , 'SNOW' , .false. , 1 , 0 297 / 298 !----------------------------------------------------------------------- 299 ! namsbc_cpl coupled ocean/atmosphere model ("key_coupled") 300 !----------------------------------------------------------------------- 301 &namsbc_cpl 476 302 / 477 303 !----------------------------------------------------------------------- … … 489 315 / 490 316 !----------------------------------------------------------------------- 491 ! namtdp tracer newtonian damping ('key_tradmp') 492 !----------------------------------------------------------------------- 493 ! ndmp type of damping in temperature and salinity 494 ! (='latitude', damping poleward of 'ndmp' degrees and function 495 ! of the distance-to-coast. Red and Med Seas as ndmp=-1) 496 ! (=-1 damping only in Med and Red Seas) 497 ! ndmpf =1 create a damping.coeff NetCDF file (the 3D damping array) 498 ! nmldmp type of damping in the mixed layer 499 ! (=0 damping throughout the water column) 500 ! (=1 no damping in the mixed layer defined by avt >5cm2/s ) 501 ! (=2 no damping in the mixed layer defined rho<rho(surf)+.01 ) 502 ! sdmp surface time scale for internal damping (days) 503 ! bdmp bottom time scale for internal damping (days) 504 ! hdmp depth of transition between sdmp and bdmp (meters) 505 &namtdp 506 ndmp = -1 507 ndmpf = 1 508 nmldmp = 1 509 sdmp = 50. 510 bdmp = 360. 511 hdmp = 800. 512 / 513 !----------------------------------------------------------------------- 514 ! nameos ocean physical parameters 515 !----------------------------------------------------------------------- 516 ! neos type of equation of state and Brunt-Vaisala frequency 517 ! = 0, UNESCO (formulation of Jackett and McDougall (1994) 518 ! and of McDougall (1987) ) 519 ! = 1, linear: rho(T) = rau0 * ( 1.028 - ralpha * T ) 520 ! = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 521 ! with rau0=1020 set in parcst routine 522 ! ralpha thermal expension coefficient (linear equation of state) 523 ! rbeta saline expension coefficient (linear equation of state) 524 &nameos 525 neos = 0 526 ralpha = 2.e-4 527 rbeta = 0.001 528 / 529 !----------------------------------------------------------------------- 530 ! namsol elliptic solver / island / free surface 531 !----------------------------------------------------------------------- 532 ! nsolv elliptic solver (=1 preconditioned conjugate gradient: pcg) 533 ! (=2 successive-over-relaxation: sor) 534 ! (=3 FETI: fet, all require "key_feti" defined) 535 ! nsol_arp absolute/relative (0/1) precision convergence test 536 ! nmin minimum of iterations for the SOR solver 537 ! nmax maximum of iterations for the SOR solver 538 ! nmod frequency of test for the SOR solver 539 ! eps absolute precision of the solver 540 ! resmax absolute precision for the SOR solver 541 ! sor optimal coefficient for SOR solver 542 ! epsisl absolute precision on stream function solver 543 ! nmisl maximum pcg iterations for island 544 ! rnu strength of the additional force used in free surface b.c. 545 &namsol 546 nsolv = 1 547 nsol_arp = 0 548 nmin = 300 549 nmax = 800 550 nmod = 10 551 eps = 1.E-6 552 resmax = 1.E-10 553 sor = 1.92 554 epsisl = 1.e-10 555 nmisl = 4000 556 rnu = 1. 557 / 558 !======================================================================= 559 ! Diagnostics namelists 560 ! namtrd dynamics and/or tracer trends 561 ! namgap level mean model-data gap 562 ! namznl zonal mean heat & freshwater fluxes computation 563 ! namspr surface pressure in rigid-lid 564 !======================================================================= 565 !----------------------------------------------------------------------- 566 ! namtrd diagnostics on dynamics and/or tracer trends 567 ! ('key_trdyn' and/or 'key_trdtra') 568 ! or mixed-layer trends ('key_trdmld') 569 !----------------------------------------------------------------------- 570 ! ntrd time step frequency dynamics and tracers trends 571 ! nctls control surface type in mixed-layer trends (0,1 or n<jpk) 572 ! ln_trdmld_restart restart for ML diagnostics 573 ! ucf unit conversion factor (=1 -> /seconds | =86400. -> /day) 574 ! ln_trdmld_instant flag to diagnose trends of instantantaneous or mean ML T/S 575 &namtrd 576 ntrd = 365 577 nctls = 0 578 ln_trdmld_restart = .false. 579 ucf = 1. 580 ln_trdmld_instant = .false. 581 / 582 !----------------------------------------------------------------------- 583 ! namgap level mean model-data gap ('key_diagap') 584 !----------------------------------------------------------------------- 585 ! ngap time-step frequency of model-data gap computation 586 ! nprg time-step frequency of gap print in model output 587 &namgap 588 ngap = 15 589 nprg = 10 590 / 591 !----------------------------------------------------------------------- 592 ! namznl zonal mean heat & freshwater fluxes computation 593 ! (#ifdef "key_diaznl") 594 !----------------------------------------------------------------------- 595 ! nfznl time-step frequency of zonal mean fluxes computation 596 &namznl 597 nfznl = 15 598 / 599 !----------------------------------------------------------------------- 600 ! namspr surface pressure diagnostic 601 !----------------------------------------------------------------------- 602 ! nmaxp maximum of iterations for the solver 603 ! epsp absolute precision of the solver 604 ! niterp number of iteration done by the solver 605 &namspr 606 nmaxp = 1000 607 epsp = 1.e-3 608 niterp = 400 317 ! namsbc_rnf runoffs namelist surface boundary condition 318 !----------------------------------------------------------------------- 319 !cn_dir Root directory for location of ssr files 320 !nn_runoff =0 no, 1 runoff, 2 runoff+river mouth ups adv 321 !rn_hrnf runoffs, depth over which enhanced vertical mixing is used 322 !rn_avt_rnf runoffs, value of the additional vertical mixing coef. [m2/s] 323 !sn_rnf informations about the runoff file to be read 324 !sn_cnf informations about the runoff mouth file to be read 325 ! THE ORDER OF THE FILES MATTER: 326 ! ! file name ! frequency (hours) ! variable ! time interpolation ! 327 ! ! ! (if <0 months) ! name ! (logical) ! 328 &namsbc_rnf 329 cn_dir = './' 330 nn_runoff = 2 331 rn_hrnf = 0.e0 332 rn_avt_rnf = 1.e-3 333 sn_rnf = 'runoff_1m_nomask.nc' , -12. , 'sorunoff', .true. , 1 , 0 334 sn_cnf = 'runoff_1m_nomask.nc' , 0. , 'socoefr' , .false. , 1 , 0 335 / 336 !----------------------------------------------------------------------- 337 ! namsbc_ssr surface boundary condition : sea surface restoring 338 !----------------------------------------------------------------------- 339 ! cn_dir directory in which the model is executed 340 ! nn_sstr =0/1 add a retroaction term in the surface heat flux 341 ! nn_sssr =O/1/2 add a damping term in the surface freshwater flux 342 ! (=1, salt flux, concentration/dillution emps only) 343 ! (=2, volume flux, both emp and emps are updated) 344 ! dqdt magnitude of the retroaction on temperature [W/m2/K] 345 ! deds magnitude of the damping on salinity [mm/day] 346 ! THE ORDER OF THE FILES MATTER: 347 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 348 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 349 &namsbc_ssr 350 cn_dir = './' 351 nn_sstr = 0 352 nn_sssr = 0 353 dqdt = -40. 354 deds = 27.7 355 sn_sst = 'sst_data.nc' , 24. , 'sst' , .false. , 0 , 0 356 sn_sss = 'sss_data.nc' , -12. , 'sss' , .true. , 0 , 0 357 / 358 !----------------------------------------------------------------------- 359 ! namalb albedo parameters 360 !----------------------------------------------------------------------- 361 ! cgren correction of the snow or ice albedo to take into account 362 ! albice albedo of melting ice in the arctic and antarctic 363 ! alphd coefficients for linear interpolation used to compute albedo 364 ! between two extremes values (Pyane, 1972) 365 ! alphc " " 366 ! alphdi " " 367 &namalb 368 cgren = 0.06 369 albice = 0.5 370 alphd = 0.80 371 alphc = 0.65 372 alphdi = 0.72 373 / 374 375 !!====================================================================== 376 !! *** Lateral boundary condition *** 377 !!====================================================================== 378 !! namlbc lateral momentum boundary condition 379 !! namcla cross land advection 380 !! namobc open boundaries parameters ("key_obc") 381 !!====================================================================== 382 383 !----------------------------------------------------------------------- 384 ! namlbc lateral momentum boundary condition 385 !----------------------------------------------------------------------- 386 ! shlat shlat = 0 , free slip ; 0 < shlat < 2 , partial slip 387 ! shlat = 2 , no slip ; 2 < shlat , strong slip 388 &namlbc 389 shlat = 2. 390 / 391 !----------------------------------------------------------------------- 392 ! namcla cross land advection 393 !----------------------------------------------------------------------- 394 ! n_cla advection between 2 ocean pts separates by land 395 &namcla 396 n_cla = 1 609 397 / 610 398 !----------------------------------------------------------------------- … … 638 426 ln_vol_cst = .false. 639 427 / 640 !----------------------------------------------------------------------- 641 ! namflo float parameters (#ifdef key_float) 428 429 !!====================================================================== 430 !! *** Bottom boundary condition *** 431 !!====================================================================== 432 !! nambfr bottom friction 433 !! nambbc bottom temperature boundary condition ("key_trabbc") 434 !! nambbl bottom boundary layer scheme ("key_trabbl_dif","key_trabbl_adv") 435 !!====================================================================== 436 !----------------------------------------------------------------------- 437 ! nambfr bottom friction 438 !----------------------------------------------------------------------- 439 ! nbotfr type of bottom friction : = 0 : no slip ; = 1 : linear friction 440 ! = 3 : free slip ; = 2 : non linear friction 441 ! bfri1 bottom drag coefficient (linear case) 442 ! bfri2 bottom drag coefficient (non linear case) 443 ! bfeb2 bottom turbulent kinetic energy (m2/s2) 444 &nambfr 445 nbotfr = 1 446 bfri1 = 4.e-4 447 bfri2 = 1.e-3 448 bfeb2 = 2.5e-3 449 / 450 !----------------------------------------------------------------------- 451 ! nambbc bottom temperature boundary condition 452 !----------------------------------------------------------------------- 453 ! ngeo_flux = 0 no geothermal heat flux 454 ! = 1 constant geothermal heat flux 455 ! = 2 variable geothermal heat flux (read in geothermal_heating.nc in mW/m2) 456 ! ngeo_flux_const Constant value of geothermal heat flux (W/m2) 457 &nambbc 458 ngeo_flux = 2 459 ngeo_flux_const = 86.4e-3 460 / 461 !----------------------------------------------------------------------- 462 ! nambbl bottom boundary layer scheme 463 !----------------------------------------------------------------------- 464 ! ! diffusive bbl ("key_trabbl") 465 ! ! advective bbl ("key_trabbl_adv") 466 ! atrbbl mixing coefficient of the bottom boundary layer scheme (m2/s) 467 &nambbl 468 atrbbl = 10000. 469 / 470 471 !!====================================================================== 472 !! Tracer (T & S ) namelists 473 !!====================================================================== 474 !! nameos equation of state 475 !! nam_traadv advection scheme 476 !! nam_traldf lateral diffusion scheme 477 !! namtdp tracer newtonian damping ("key_tradmp") 478 !!====================================================================== 479 480 !----------------------------------------------------------------------- 481 ! nameos ocean physical parameters 482 !----------------------------------------------------------------------- 483 ! neos type of equation of state and Brunt-Vaisala frequency 484 ! ! = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 485 ! ! = 1, linear: rho(T) = rau0 * ( 1.028 - ralpha * T ) 486 ! ! = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 487 ! ralpha thermal expension coefficient (neos= 1 or 2) 488 ! rbeta saline expension coefficient (neos= 2) 489 &nameos 490 neos = 0 491 ralpha = 2.e-4 492 rbeta = 0.001 493 / 494 !----------------------------------------------------------------------- 495 ! nam_traadv advection scheme for tracer 496 !----------------------------------------------------------------------- 497 ! ln_traadv_cen2 2nd order centered scheme (default T) 498 ! ln_traadv_tvd TVD scheme (default F) 499 ! ln_traadv_muscl MUSCL scheme (default F) 500 ! ln_traadv_muscl2 MUSCL scheme + cen2 at boundaries (default F) 501 ! ln_traadv_ubs UBS scheme (default F) 502 ! ln_traadv_qck QUICKEST scheme (default F) 503 &nam_traadv 504 ln_traadv_cen2 = .true. 505 ln_traadv_tvd = .false. 506 ln_traadv_muscl = .false. 507 ln_traadv_muscl2 = .false. 508 ln_traadv_ubs = .false. 509 ln_traadv_qck = .false. 510 / 511 !----------------------------------------------------------------------- 512 ! nam_traldf lateral diffusion scheme for tracer 513 !----------------------------------------------------------------------- 514 ! Type of the operator : 515 ! ln_traldf_lap laplacian operator (default T) 516 ! ln_traldf_bilap bilaplacian operator (default F) 517 ! Direction of action : 518 ! ln_traldf_level iso-level (default F) 519 ! ln_traldf_hor horizontal (geopotential) (default F) (require "key_ldfslp" when ln_sco=T) 520 ! ln_traldf_iso iso-neutral (default T) (require "key_ldfslp") 521 ! Coefficient 522 ! aht0 horizontal eddy diffusivity for tracers (m2/s) 523 ! ahtb0 background eddy diffusivity for isopycnal diffusion (m2/s) 524 ! aeiv0 eddy induced velocity coefficient (m2/s) (require "key_traldf_eiv") 525 &nam_traldf 526 ln_traldf_lap = .true. 527 ln_traldf_bilap = .false. 528 ln_traldf_level = .false. 529 ln_traldf_hor = .false. 530 ln_traldf_iso = .true. 531 aht0 = 2000. 532 ahtb0 = 0. 533 aeiv0 = 2000. 534 / 535 !----------------------------------------------------------------------- 536 ! namtdp tracer newtonian damping ('key_tradmp') 537 !----------------------------------------------------------------------- 538 ! ndmp type of damping in temperature and salinity 539 ! ! ='latitude', damping poleward of 'ndmp' degrees and function 540 ! ! of the distance-to-coast. Red and Med Seas as ndmp=-1 541 ! ! =-1 damping only in Med and Red Seas 542 ! ndmpf =1 create a damping.coeff NetCDF file (the 3D damping array) 543 ! nmldmp type of damping in the mixed layer 544 ! ! =0 damping throughout the water column 545 ! ! =1 no damping in the mixed layer defined by avt >5cm2/s 546 ! ! =2 no damping in the mixed layer defined rho<rho(surf)+.01 547 ! sdmp surface time scale for internal damping (days) 548 ! bdmp bottom time scale for internal damping (days) 549 ! hdmp depth of transition between sdmp and bdmp (meters) 550 &namtdp 551 ndmp = -1 552 ndmpf = 1 553 nmldmp = 1 554 sdmp = 50. 555 bdmp = 360. 556 hdmp = 800. 557 / 558 559 !!====================================================================== 560 !! *** Dynamics namelists *** 561 !!====================================================================== 562 !! nam_dynadv formulation of the momentum advection 563 !! nam_dynvor advection scheme 564 !! nam_dynhpg hydrostatic pressure gradient 565 !! namflg hydrostatic pressure gradient time stepping 566 !! nam_dynspg surface pressure gradient (CPP key only) 567 !! nam_dynldf lateral diffusion scheme 568 !!====================================================================== 569 570 !----------------------------------------------------------------------- 571 ! nam_dynadv formulation of the momentum advection 572 !----------------------------------------------------------------------- 573 ! ln_dynadv_vec vector form (T) or flux form (F) (default T) 574 ! ln_dynadv_cen2 flux form - 2nd order centered scheme (default T) 575 ! ln_dynadv_ubs flux form - 3rd order UBS scheme (default F) 576 &nam_dynadv 577 ln_dynadv_vec = .true. 578 ln_dynadv_cen2 = .false. 579 ln_dynadv_ubs = .false. 580 / 581 !----------------------------------------------------------------------- 582 ! nam_dynvor option of physics/algorithm (not control by CPP keys) 583 !----------------------------------------------------------------------- 584 ! ln_dynvor_ens vorticity trends: enstrophy conserving scheme (default T) 585 ! ln_dynvor_ene " " : energy conserving scheme (default F) 586 ! ln_dynvor_mix " " : mixed scheme (default F) 587 ! ln_dynvor_een " " : energy & enstrophy scheme (default F) 588 &nam_dynvor 589 ln_dynvor_ene = .false. 590 ln_dynvor_ens = .true. 591 ln_dynvor_mix = .false. 592 ln_dynvor_een = .false. 593 / 594 !----------------------------------------------------------------------- 595 ! nam_dynhpg Hydrostatic pressure gradient option 596 !----------------------------------------------------------------------- 597 ! type of pressure gradient scheme (choose one only!) 598 ! ln_hpg_zco z-coordinate - full steps (default T) 599 ! ln_hpg_zps z-coordinate - partial steps (interpolation) 600 ! ln_hpg_sco s-coordinate (standard jacobian formulation) 601 ! ln_hpg_hel s-coordinate (helsinki modification) 602 ! ln_hpg_wdj s-coordinate (weighted density jacobian) 603 ! ln_hpg_djc s-coordinate (Density Jacobian with Cubic polynomial) 604 ! ln_hpg_rot s-coordinate (ROTated axes scheme) 605 ! parameters 606 ! gamm weighting coefficient (wdj scheme) 607 &nam_dynhpg 608 ln_hpg_zco = .false. 609 ln_hpg_zps = .true. 610 ln_hpg_sco = .false. 611 ln_hpg_hel = .false. 612 ln_hpg_wdj = .false. 613 ln_hpg_djc = .false. 614 ln_hpg_rot = .false. 615 gamm = 0.e0 616 / 617 !----------------------------------------------------------------------- 618 ! namflg algorithm flags (algorithm not control by CPP keys) 619 !----------------------------------------------------------------------- 620 ! ln_dynhpg_imp hydrostatic pressure gradient: semi-implicit time scheme (T) 621 ! centered time scheme (F) 622 ! nn_dynhpg_rst add dynhpg implicit variables in restart ot not (1/0) 623 &namflg 624 ln_dynhpg_imp = .false. 625 nn_dynhpg_rst = 0 626 / 627 !----------------------------------------------------------------------- 628 ! nam_dynspg surface pressure gradient (CPP key only) 629 !----------------------------------------------------------------------- 630 ! ! explicit free surface ("key_dynspg_exp") 631 ! ! filtered free surface ("key_dynspg_flt") 632 ! ! split-explicit free surface ("key_dynspg_ts") 633 ! ! rigid-lid ("key_dynspg_rl") 634 635 !----------------------------------------------------------------------- 636 ! nam_dynldf lateral diffusion on momentum 637 !----------------------------------------------------------------------- 638 ! Type of the operator : 639 ! ln_dynldf_lap laplacian operator (default T) 640 ! ln_dynldf_bilap bilaplacian operator (default F) 641 ! Direction of action : 642 ! ln_dynldf_level iso-level (default F) 643 ! ln_dynldf_hor horizontal (geopotential) (default F) (require "key_ldfslp" if ln_sco=T) 644 ! ln_dynldf_iso iso-neutral (default T) (require "key_ldfslp") 645 ! Coefficient 646 ! ahm0 horizontal eddy viscosity for the dynamics (m2/s) 647 ! ahmb0 background eddy viscosity for isopycnal diffusion (m2/s) 648 &nam_dynldf 649 ln_dynldf_lap = .true. 650 ln_dynldf_bilap = .false. 651 ln_dynldf_level = .false. 652 ln_dynldf_hor = .true. 653 ln_dynldf_iso = .false. 654 ahm0 = 40000. 655 ahmb0 = 0. 656 / 657 658 !!====================================================================== 659 !! Tracers & Dynamics vertical physics namelists 660 !!====================================================================== 661 !! namzdf vertical physics 662 !! namnpc non penetrative convection 663 !! namric richardson number dependent vertical mixing ("key_zdfric" ) 664 !! namtke TKE dependent vertical mixing ("key_zdftke" ) 665 !! namkpp KPP dependent vertical mixing ("key_zdfkpp" ) 666 !! namddm double diffusive mixing parameterization ("key_zdfddm" ) 667 !!====================================================================== 668 669 !----------------------------------------------------------------------- 670 ! namzdf vertical physics 671 !----------------------------------------------------------------------- 672 ! ln_zdfevd enhanced vertical diffusion (default T) 673 ! ln_zdfnpc Non-Penetrative Convection (default T) 674 ! avm0 vertical eddy viscosity for the dynamic (m2/s) (also background Kz if not "key_zdfcst") 675 ! avt0 vertical eddy diffusivity for tracers (m2/s) (also background Kz if not "key_zdfcst") 676 ! avevd vertical coefficient for enhanced diffusion scheme (m2/s) 677 ! n_evdm = 0 apply enhanced mixing on tracer only 678 ! = 1 apply enhanced mixing on both tracer and momentum 679 ! ln_zdfexp vertical physics: (=T) time splitting (T) (Default=F) 680 ! (=F) euler backward (F) 681 ! n_zdfexp number of sub-timestep for time splitting scheme 682 &namzdf 683 ln_zdfevd = .true. 684 ln_zdfnpc = .false. 685 avm0 = 1.2e-4 686 avt0 = 1.2e-5 687 avevd = 100. 688 n_evdm = 1 689 ln_zdfexp = .false. 690 n_zdfexp = 3 691 / 692 !----------------------------------------------------------------------- 693 ! namnpc non penetrative convection 694 !----------------------------------------------------------------------- 695 ! nnpc1 non penetrative convective scheme frequency 696 ! nnpc2 non penetrative convective scheme print frequency 697 &namnpc 698 nnpc1 = 1 699 nnpc2 = 365 700 / 701 !----------------------------------------------------------------------- 702 ! namric richardson number dependent vertical diffusion ("key_zdfric" ) 703 !----------------------------------------------------------------------- 704 ! avmri maximum value of the vertical viscosity 705 ! alp coefficient of the parameterization 706 ! nric coefficient of the parameterization 707 &namric 708 avmri = 100.e-4 709 alp = 5. 710 nric = 2 711 / 712 !----------------------------------------------------------------------- 713 ! namtke turbulent eddy kinetic dependent vertical diffusion ("key_zdftke") 714 !----------------------------------------------------------------------- 715 ! ln_rstke flag to restart with tke from a run without tke (default F) 716 ! ediff coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) ) 717 ! ediss coef. of the Kolmogoroff dissipation 718 ! ebb coef. of the surface input of tke 719 ! efave coef. to applied to the tke diffusion ( avtke=efave*avm ) 720 ! emin minimum value of tke (m^2/s^2) 721 ! emin0 surface minimum value of tke (m^2/s^2) 722 ! nitke number of restart iterative loops 723 ! ri_c critic richardson number = 2/9 = 0.22222222 (hard coded) 724 ! nmxl length used = 0 bounded by the distance to surface and bottom 725 ! ! = 1 bounded by the local vertical scale factor 726 ! ! = 2 first vertical derivative of mixing length bounded by 1 727 ! npdl Prandtl number = 0 no vertical prandtl number (avt=avm) 728 ! ! = 1 prandtl number function of richarson number (avt=pdl*avm) 729 ! ! = 2 same as = 1 but a shapiro filter is applied on pdl 730 ! nave = horizontal averaged (=1) or not (=0) of avt (default =1) 731 ! navb = 0 cst background avt0, avm0 / =1 profile used on avtb 732 &namtke 733 ln_rstke = .false. 734 ediff = 0.1 735 ediss = 0.7 736 ebb = 3.75 737 efave = 1. 738 emin = 1.e-6 739 emin0 = 1.e-4 740 nitke = 50 741 nmxl = 2 742 npdl = 1 743 navb = 0 744 / 745 !----------------------------------------------------------------------- 746 ! namkpp K-Profile Parameterization dependent vertical mixing ("key_zdfkpp" ) 747 !----------------------------------------------------------------------- 748 ! ln_kpprimix shear instability mixing (default T) 749 ! difmiw constant internal wave viscosity (m2/s) 750 ! difsiw constant internal wave diffusivity (m2/s) 751 ! Riinfty local Richardson Number limit for shear instability 752 ! difri maximum shear mixing at Rig = 0 (m2/s) 753 ! bvsqcon Brunt-Vaisala squared (1/s**2) for maximum convection 754 ! difcon maximum mixing in interior convection (m2/s) 755 ! nave = 0/1 flag for horizontal average on avt, avmu, avmv 756 ! navb = 0/1 flag for constant or profile background avt 757 &namkpp 758 ln_kpprimix = .true. 759 difmiw = 1.e-04 760 difsiw = 0.1e-04 761 Riinfty = 0.8 762 difri = 0.0050 763 bvsqcon = -0.01e-07 764 difcon = 1. 765 navb = 0 766 nave = 1 767 / 768 !----------------------------------------------------------------------- 769 ! namddm double diffusive mixing parameterization ("key_zdfddm") 770 !----------------------------------------------------------------------- 771 ! avts maximum avs for dd mixing 772 ! hsbfr heat/salt buoyancy flux ratio 773 &namddm 774 avts = 1.e-4 775 hsbfr = 1.6 776 / 777 778 !!====================================================================== 779 !! *** Miscelaneous namelists *** 780 !!====================================================================== 781 !! namsol elliptic solver / island / free surface 782 !!====================================================================== 783 784 !----------------------------------------------------------------------- 785 ! namsol elliptic solver / island / free surface 786 !----------------------------------------------------------------------- 787 ! nsolv elliptic solver (=1 preconditioned conjugate gradient: pcg) 788 ! (=2 successive-over-relaxation: sor) 789 ! (=3 FETI currently it does not work! ("key_feti") 790 ! nsol_arp absolute/relative (0/1) precision convergence test 791 ! nmin minimum of iterations for the SOR solver 792 ! nmax maximum of iterations for the SOR solver 793 ! nmod frequency of test for the SOR solver 794 ! eps absolute precision of the solver 795 ! resmax absolute precision for the SOR solver 796 ! sor optimal coefficient for SOR solver 797 ! epsisl absolute precision on stream function solver 798 ! nmisl maximum pcg iterations for island 799 ! rnu strength of the additional force used in free surface b.c. 800 &namsol 801 nsolv = 1 802 nsol_arp = 0 803 nmin = 300 804 nmax = 800 805 nmod = 10 806 eps = 1.E-6 807 resmax = 1.E-10 808 sor = 1.92 809 epsisl = 1.e-10 810 nmisl = 4000 811 rnu = 1. 812 / 813 !!====================================================================== 814 !! *** Diagnostics namelists *** 815 !!====================================================================== 816 !! namtrd dynamics and/or tracer trends ("key_trddyn","key_trdtra","key_trdmld") 817 !! namgap level mean model-data gap ("key_diagap") 818 !! namspr surface pressure diagnosed in rigid-lid ("key_diaspr") 819 !! namflo float parameters ("key_float") 820 !! namptr Poleward Transport Diagnostics 821 !!====================================================================== 822 823 !----------------------------------------------------------------------- 824 ! namtrd diagnostics on dynamics and/or tracer trends ("key_trddyn" and/or "key_trdtra") 825 ! or mixed-layer trends ('key_trdmld') 826 ! or barotropic vorticity ("key_trdvor") 827 !----------------------------------------------------------------------- 828 ! ntrd time step frequency dynamics and tracers trends 829 ! nctls control surface type in mixed-layer trends (0,1 or n<jpk) 830 ! ln_trdmld_restart restart for ML diagnostics 831 ! ucf unit conversion factor (=1 -> /seconds | =86400. -> /day) 832 ! ln_trdmld_instant flag to diagnose trends of instantantaneous or mean ML T/S 833 &namtrd 834 ntrd = 365 835 nctls = 0 836 ln_trdmld_restart = .false. 837 ucf = 1. 838 ln_trdmld_instant = .false. 839 / 840 !----------------------------------------------------------------------- 841 ! namgap level mean model-data gap ('key_diagap') 842 !----------------------------------------------------------------------- 843 ! ngap time-step frequency of model-data gap computation 844 ! nprg time-step frequency of gap print in model output 845 &namgap 846 ngap = 15 847 nprg = 10 848 / 849 !----------------------------------------------------------------------- 850 ! namspr surface pressure diagnostic 851 !----------------------------------------------------------------------- 852 ! nmaxp maximum of iterations for the solver 853 ! epsp absolute precision of the solver 854 ! niterp number of iteration done by the solver 855 &namspr 856 nmaxp = 1000 857 epsp = 1.e-3 858 niterp = 400 859 / 860 !----------------------------------------------------------------------- 861 ! namflo float parameters ("key_float") 642 862 !----------------------------------------------------------------------- 643 863 ! ln_rstflo boolean term for float restart (true or false) … … 654 874 ln_flork4 = .false. 655 875 / 656 !!====================================================================== 657 !! Surface Boundary Condition namelists 658 !! 659 !! namsbc surface boundary condition 660 !! namsbc_ana analytical formulation of sbc 661 !! namsbc_flx flux formulation of sbc 662 !! namsbc_core CORE bulk formulea formulation of sbc 663 !! namsbc_clio CLIO bulk formulea formulation of sbc 664 !! namsbc_cpl CouPLed formulation of sbc 665 !! namsbc_ssr sea surface restoring term (for T and/or S) 666 !!====================================================================== 667 !----------------------------------------------------------------------- 668 ! namsbc Surface Boundary Condition (surface module) 669 !----------------------------------------------------------------------- 670 ! nn_fsbc frequency of surface boundary condition computation 671 ! (also correspond to the frequency of sea-ice model call) 672 ! 673 ! ln_ana analytical formulation (fill namsbc_ana ) 674 ! ln_flx flux formulation (fill namsbc_flx ) (overwritten by key_sbc_flux ) 675 ! ln_blk_clio CLIO bulk formulation (fill namsbc_core) (overwritten by key_sbc_bulk_clio) 676 ! ln_blk_core CORE bulk formulation (fill namsbc_clio) (overwritten by key_sbc_bulk_core) 677 ! ln_cpl Coupled formulation (fill namsbc_cpl ) (overwritten by key_sbc_coupled ) 678 ! 679 ! nn_ice =0 no ice boundary condition , 680 ! =1 observed ice-cover , 681 ! =2 ice-model ice-cover (overwritten by key_ice_lim) 682 ! 683 ! ln_dm2dc Daily mean to Diurnal Cycle short wave (qsr) 684 ! 685 ! ln_ssr Sea Surface restoring on SST and/or SSS (fill namsbc_ssr) 686 ! (overwritten by key_dta_sst and/or key_dtasss) 687 ! 688 ! nn_fwb FreshWater Budget: =0 no control , 689 ! =1 annual global mean of emp set to zero , 690 ! =2 global mean of emp set to zero at each nn_fsbc time step 691 &namsbc 692 nn_fsbc = 5 693 ln_ana = .false. 694 ln_flx = .false. 695 ln_blk_clio = .false. 696 ln_blk_core = .true. 697 ln_cpl = .false. 698 nn_ice = 2 699 ln_dm2dc = .false. 700 ln_rnf = .false. 701 ln_ssr = .false. 702 nn_fwb = 0 703 / 704 !----------------------------------------------------------------------- 705 ! namsbc_rnf runoffs namelist surface boundary condition 706 !----------------------------------------------------------------------- 707 !nn_runoff =0 no, 1 runoff, 2 runoff+river mouth ups adv 708 !rn_hrnf runoffs, depth over which enhanced vertical mixing is used 709 !rn_avt_rnf runoffs, value of the additional vertical mixing coef. [m2/s] 710 !cn_dir Root directory for location of ssr files 711 !sn_rnf informations about the runoff file to be read 712 !sn_cnf informations about the runoff mouth file to be read 713 ! THE ORDER OF THE FILES MATTER: 714 ! ! file name ! frequency (hours) ! variable ! time interpolation ! 715 ! ! ! (if <0 months) ! name ! (logical) ! 716 &namsbc_rnf 717 nn_runoff = 2 718 rn_hrnf = 0.e0 719 rn_avt_rnf = 1.e-3 720 cn_dir = './' 721 sn_rnf = 'runoff_1m_nomask.nc' , -12. , 'sorunoff', .TRUE. , 1 , 0 722 sn_cnf = 'runoff_1m_nomask.nc' , 0. , 'socoefr' , .FALSE. , 1 , 0 723 / 724 !----------------------------------------------------------------------- 725 ! namsbc_ana analytical surface boundary condition 726 !----------------------------------------------------------------------- 727 ! nn_tau000 gently increase the stress over the first ntau_rst time-steps 728 ! rn_utau0 uniform value used as default surface heat flux 729 ! rn_vtau0 uniform value used as default solar radiation flux 730 ! rn_q0 uniform value used as default surface heat flux 731 ! rn_qsr0 uniform value used as default solar radiation flux 732 ! rn_emp0 uniform value used as default surface freswater budget (E-P) 733 &namsbc_ana 734 nn_tau000 = 0 735 rn_utau0 = 0.5 736 rn_vtau0 = 0.e0 737 rn_q0 = 0.e0 738 rn_qsr0 = 0.e0 739 rn_emp0 = 0.e0 740 / 741 !----------------------------------------------------------------------- 742 ! namsbc_flx surface boundary condition : flux formulation (#ifdef "key_sbc_flux") 743 !----------------------------------------------------------------------- 744 ! cn_dir directory in which the model is executed 745 ! 746 ! THE ORDER OF THE FILES MATTER: 747 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 748 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 749 &namsbc_flx 750 cn_dir = './' 751 sn_utau = 'utau.nc' , 24. , 'utau' , .FALSE. , 0 , 0 752 sn_vtau = 'vtau.nc' , 24. , 'vtau' , .FALSE. , 0 , 0 753 sn_qtot = 'qtot.nc' , 24. , 'qtot' , .FALSE. , 0 , 0 754 sn_qsr = 'qsr.nc' , 24. , 'qsr' , .FALSE. , 0 , 0 755 sn_emp = 'emp.nc' , 24. , 'emp' , .FALSE. , 0 , 0 756 757 / 758 !----------------------------------------------------------------------- 759 ! namsbc_clio CLIO bulk formulea 760 !----------------------------------------------------------------------- 761 ! cn_dir directory in which the model is executed 762 ! 763 ! THE ORDER OF THE FILES MATTER: 764 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 765 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 766 &namsbc_clio 767 cn_dir = './' 768 sn_utau = 'taux_1m' , -12. , 'sozotaux', .FALSE. , 1 , 0 769 sn_vtau = 'tauy_1m' , -12. , 'sometauy', .FALSE. , 1 , 0 770 sn_wndm = 'flx' , -12. , 'socliowi', .FALSE. , 1 , 0 771 sn_tair = 'flx' , -12. , 'socliot1', .FALSE. , 1 , 0 772 sn_humi = 'flx' , -12. , 'socliohu', .FALSE. , 1 , 0 773 sn_ccov = 'flx' , -12. , 'socliocl', .FALSE. , 1 , 0 774 sn_prec = 'flx' , -12. , 'socliopl', .FALSE. , 1 , 0 775 / 776 !----------------------------------------------------------------------- 777 ! namsbc_core CORE bulk formulea 778 !----------------------------------------------------------------------- 779 ! cn_dir directory in which the model is executed 780 ! ln_2m logical flag to use air temp. and hum referenced at 2m instead 10m 781 ! alpha_precip multiplication factor for precipitation (total & snow) 782 ! 783 ! THE ORDER OF THE FILES MATTER: 784 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 785 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 786 &namsbc_core 787 cn_dir = './' 788 ln_2m = .FALSE. 789 alpha_precip = 1. 790 sn_wndi = 'u10' , 24. , 'U_10_MOD', .FALSE. , 1 , 0 791 sn_wndj = 'v10' , 24. , 'V_10_MOD', .FALSE. , 1 , 0 792 sn_qsr = 'rad' , 24. , 'SWDN_MOD', .FALSE. , 1 , 0 793 sn_qlw = 'rad' , 24. , 'LWDN_MOD', .FALSE. , 1 , 0 794 sn_tair = 't10' , 24. , 'T_10_MOD', .FALSE. , 1 , 0 795 sn_humi = 'q10' , 24. , 'Q_10_MOD', .FALSE. , 1 , 0 796 sn_prec = 'precip' , -12. , 'RAIN' , .FALSE. , 1 , 0 797 sn_snow = 'precip' , -12. , 'SNOW' , .FALSE. , 1 , 0 798 / 799 !----------------------------------------------------------------------- 800 ! namsbc_cpl coupled ocean/atmosphere model (#ifdef "key_coupled") 801 !----------------------------------------------------------------------- 802 &namsbc_cpl 803 / 804 !----------------------------------------------------------------------- 805 ! namsbc_ssr surface boundary condition : sea surface restoring 806 !----------------------------------------------------------------------- 807 ! nn_sstr =0/1 add a retroaction term in the surface heat flux 808 ! nn_sssr =O/1/2 add a damping term in the surface freshwater flux 809 ! (=1, salt flux, concentration/dillution emps only) 810 ! (=2, volume flux, both emp and emps are updated) 811 ! dqdt magnitude of the retroaction on temperature [W/m2/K] 812 ! deds magnitude of the damping on salinity [mm/day] 813 ! THE ORDER OF THE FILES MATTER: 814 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 815 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 816 &namsbc_ssr 817 nn_sstr = 0 818 nn_sssr = 0 819 dqdt = -40. 820 deds = 27.7 821 sn_sst = 'sst_data.nc' , 24. , 'sst' , .FALSE. , 0 , 0 822 sn_sss = 'sss_data.nc' , -12. , 'sss' , .TRUE. , 0 , 0 823 / 876 !----------------------------------------------------------------------- 877 ! namptr Poleward Transport Diagnostic 878 !----------------------------------------------------------------------- 879 ! ln_diaptr logical flag for Poleward transport computation 880 ! ln_subbas logical flag for Atlantic/Pacific/Indian basins computation 881 ! need input basins mask file named "subbasins.nc" 882 ! nf_ptr Frequency of computation 883 &namptr 884 ln_diaptr = .false. 885 ln_subbas = .false. 886 nf_ptr = 15 887 /
Note: See TracChangeset
for help on using the changeset viewer.