Changeset 907 for trunk/CONFIG/ORCA2_LIM/EXP00/namelist
- Timestamp:
- 2008-04-25T09:35:26+02:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/CONFIG/ORCA2_LIM/EXP00/namelist
r889 r907 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 … … 40 55 nstock = 5475 41 56 nwrite = 5475 42 ln_dimgnnn = .FALSE. 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 ! ln_dm2dc Daily mean to Diurnal Cycle short wave (qsr) 203 ! ln_ssr Sea Surface Restoring on T and/or S (fill namsbc_ssr) 204 ! nn_fwb FreshWater Budget: =0 no control , 205 ! =1 annual global mean of e-p-r set to zero , 206 ! =2 global mean of e-p-r set to zero at each nn_fsbc time step 207 &namsbc 208 nn_fsbc = 5 209 ln_ana = .false. 210 ln_flx = .false. 211 ln_blk_clio = .true. 212 ln_blk_core = .false. 213 ln_cpl = .false. 214 nn_ice = 2 215 ln_dm2dc = .false. 216 ln_rnf = .false. 217 ln_ssr = .false. 218 nn_fwb = 0 219 / 220 !----------------------------------------------------------------------- 221 ! namsbc_ana analytical surface boundary condition 222 !----------------------------------------------------------------------- 223 ! nn_tau000 gently increase the stress over the first ntau_rst time-steps 224 ! rn_utau0 uniform value for the i-stress 225 ! rn_vtau0 uniform value for the j-stress 226 ! rn_q0 uniform value for the total heat flux 227 ! rn_qsr0 uniform value for the solar radiation 228 ! rn_emp0 uniform value for the freswater budget (E-P) 229 &namsbc_ana 230 nn_tau000 = 0 231 rn_utau0 = 0.5 232 rn_vtau0 = 0.e0 233 rn_q0 = 0.e0 234 rn_qsr0 = 0.e0 235 rn_emp0 = 0.e0 236 / 237 !----------------------------------------------------------------------- 238 ! namsbc_flx surface boundary condition : flux formulation (#ifdef "key_sbc_flux") 239 !----------------------------------------------------------------------- 240 ! cn_dir directory in which the model is executed 241 ! 242 ! THE ORDER OF THE FILES MATTER: 243 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 244 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 245 &namsbc_flx 246 cn_dir = './' 247 sn_utau = 'utau.nc' , 24. , 'utau' , .false. , 0 , 0 248 sn_vtau = 'vtau.nc' , 24. , 'vtau' , .false. , 0 , 0 249 sn_qtot = 'qtot.nc' , 24. , 'qtot' , .false. , 0 , 0 250 sn_qsr = 'qsr.nc' , 24. , 'qsr' , .false. , 0 , 0 251 sn_emp = 'emp.nc' , 24. , 'emp' , .false. , 0 , 0 252 / 253 !----------------------------------------------------------------------- 254 ! namsbc_clio CLIO bulk formulea 255 !----------------------------------------------------------------------- 256 ! cn_dir directory in which the model is executed 257 ! 258 ! THE ORDER OF THE FILES MATTER: 259 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 260 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 261 &namsbc_clio 262 cn_dir = './' 263 sn_utau = 'taux_1m' , -12. , 'sozotaux', .false. , 1 , 0 264 sn_vtau = 'tauy_1m' , -12. , 'sometauy', .false. , 1 , 0 265 sn_wndm = 'flx' , -12. , 'socliowi', .false. , 1 , 0 266 sn_tair = 'flx' , -12. , 'socliot1', .false. , 1 , 0 267 sn_humi = 'flx' , -12. , 'socliohu', .false. , 1 , 0 268 sn_ccov = 'flx' , -12. , 'socliocl', .false. , 1 , 0 269 sn_prec = 'flx' , -12. , 'socliopl', .false. , 1 , 0 270 / 271 !----------------------------------------------------------------------- 272 ! namsbc_core CORE bulk formulea 273 !----------------------------------------------------------------------- 274 ! cn_dir directory in which the model is executed 275 ! ln_2m logical flag to use air temp. and hum referenced at 2m instead 10m 276 ! alpha_precip multiplication factor for precipitation (total & snow) 277 ! 278 ! THE ORDER OF THE FILES MATTER: 279 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 280 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 281 &namsbc_core 282 cn_dir = './' 283 ln_2m = .false. 284 alpha_precip = 1. 285 sn_wndi = 'u10' , 24. , 'U_10_MOD', .false. , 1 , 0 286 sn_wndj = 'v10' , 24. , 'V_10_MOD', .false. , 1 , 0 287 sn_qsr = 'rad' , 24. , 'SWDN_MOD', .false. , 1 , 0 288 sn_qlw = 'rad' , 24. , 'LWDN_MOD', .false. , 1 , 0 289 sn_tair = 't10' , 24. , 'T_10_MOD', .false. , 1 , 0 290 sn_humi = 'q10' , 24. , 'Q_10_MOD', .false. , 1 , 0 291 sn_prec = 'precip' , -12. , 'RAIN' , .false. , 1 , 0 292 sn_snow = 'precip' , -12. , 'SNOW' , .false. , 1 , 0 293 / 294 !----------------------------------------------------------------------- 295 ! namsbc_cpl coupled ocean/atmosphere model ("key_coupled") 296 !----------------------------------------------------------------------- 297 &namsbc_cpl 476 298 / 477 299 !----------------------------------------------------------------------- … … 489 311 / 490 312 !----------------------------------------------------------------------- 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 313 ! namsbc_rnf runoffs namelist surface boundary condition 314 !----------------------------------------------------------------------- 315 !nn_runoff =0 no, 1 runoff, 2 runoff+river mouth ups adv 316 !rn_hrnf runoffs, depth over which enhanced vertical mixing is used 317 !rn_avt_rnf runoffs, value of the additional vertical mixing coef. [m2/s] 318 !cn_dir Root directory for location of ssr files 319 !sn_rnf informations about the runoff file to be read 320 !sn_cnf informations about the runoff mouth file to be read 321 ! THE ORDER OF THE FILES MATTER: 322 ! ! file name ! frequency (hours) ! variable ! time interpolation ! 323 ! ! ! (if <0 months) ! name ! (logical) ! 324 &namsbc_rnf 325 nn_runoff = 2 326 rn_hrnf = 0.e0 327 rn_avt_rnf = 1.e-3 328 cn_dir = './' 329 sn_rnf = 'runoff_1m_nomask.nc' , -12. , 'sorunoff', .true. , 1 , 0 330 sn_cnf = 'runoff_1m_nomask.nc' , 0. , 'socoefr' , .false. , 1 , 0 331 / 332 !----------------------------------------------------------------------- 333 ! namsbc_ssr surface boundary condition : sea surface restoring 334 !----------------------------------------------------------------------- 335 ! nn_sstr =0/1 add a retroaction term in the surface heat flux 336 ! nn_sssr =O/1/2 add a damping term in the surface freshwater flux 337 ! (=1, salt flux, concentration/dillution emps only) 338 ! (=2, volume flux, both emp and emps are updated) 339 ! dqdt magnitude of the retroaction on temperature [W/m2/K] 340 ! deds magnitude of the damping on salinity [mm/day] 341 ! THE ORDER OF THE FILES MATTER: 342 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 343 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 344 &namsbc_ssr 345 nn_sstr = 0 346 nn_sssr = 0 347 dqdt = -40. 348 deds = 27.7 349 sn_sst = 'sst_data.nc' , 24. , 'sst' , .false. , 0 , 0 350 sn_sss = 'sss_data.nc' , -12. , 'sss' , .true. , 0 , 0 351 / 352 !----------------------------------------------------------------------- 353 ! namalb albedo parameters 354 !----------------------------------------------------------------------- 355 ! cgren correction of the snow or ice albedo to take into account 356 ! albice albedo of melting ice in the arctic and antarctic 357 ! alphd coefficients for linear interpolation used to compute albedo 358 ! between two extremes values (Pyane, 1972) 359 ! alphc " " 360 ! alphdi " " 361 &namalb 362 cgren = 0.06 363 albice = 0.5 364 alphd = 0.80 365 alphc = 0.65 366 alphdi = 0.72 367 / 368 369 !!====================================================================== 370 !! *** Lateral boundary condition *** 371 !!====================================================================== 372 !! namlbc lateral momentum boundary condition 373 !! namcla cross land advection 374 !! namobc open boundaries parameters ("key_obc") 375 !!====================================================================== 376 377 !----------------------------------------------------------------------- 378 ! namlbc lateral momentum boundary condition 379 !----------------------------------------------------------------------- 380 ! shlat shlat = 0 , free slip ; 0 < shlat < 2 , partial slip 381 ! shlat = 2 , no slip ; 2 < shlat , strong slip 382 &namlbc 383 shlat = 2. 384 / 385 !----------------------------------------------------------------------- 386 ! namcla cross land advection 387 !----------------------------------------------------------------------- 388 ! n_cla advection between 2 ocean pts separates by land 389 &namcla 390 n_cla = 1 609 391 / 610 392 !----------------------------------------------------------------------- … … 638 420 ln_vol_cst = .false. 639 421 / 640 !----------------------------------------------------------------------- 641 ! namflo float parameters (#ifdef key_float) 422 423 !!====================================================================== 424 !! *** Bottom boundary condition *** 425 !!====================================================================== 426 !! nambfr bottom friction 427 !! nambbc bottom temperature boundary condition ("key_trabbc") 428 !! nambbl bottom boundary layer scheme ("key_trabbl_dif","key_trabbl_adv") 429 !!====================================================================== 430 !----------------------------------------------------------------------- 431 ! nambfr bottom friction 432 !----------------------------------------------------------------------- 433 ! nbotfr type of bottom friction : = 0 : no slip ; = 1 : linear friction 434 ! = 3 : free slip ; = 2 : non linear friction 435 ! bfri1 bottom drag coefficient (linear case) 436 ! bfri2 bottom drag coefficient (non linear case) 437 ! bfeb2 bottom turbulent kinetic energy (m2/s2) 438 &nambfr 439 nbotfr = 1 440 bfri1 = 4.e-4 441 bfri2 = 1.e-3 442 bfeb2 = 2.5e-3 443 / 444 !----------------------------------------------------------------------- 445 ! nambbc bottom temperature boundary condition 446 !----------------------------------------------------------------------- 447 ! ngeo_flux = 0 no geothermal heat flux 448 ! = 1 constant geothermal heat flux 449 ! = 2 variable geothermal heat flux (read in geothermal_heating.nc in mW/m2) 450 ! ngeo_flux_const Constant value of geothermal heat flux (W/m2) 451 &nambbc 452 ngeo_flux = 2 453 ngeo_flux_const = 86.4e-3 454 / 455 !----------------------------------------------------------------------- 456 ! nambbl bottom boundary layer scheme 457 !----------------------------------------------------------------------- 458 ! ! diffusive bbl ("key_trabbl") 459 ! ! advective bbl ("key_trabbl_adv") 460 ! atrbbl mixing coefficient of the bottom boundary layer scheme (m2/s) 461 &nambbl 462 atrbbl = 10000. 463 / 464 465 !!====================================================================== 466 !! Tracer (T & S ) namelists 467 !!====================================================================== 468 !! nameos equation of state 469 !! nam_traadv advection scheme 470 !! nam_traldf lateral diffusion scheme 471 !! namtdp tracer newtonian damping ("key_tradmp") 472 !!====================================================================== 473 474 !----------------------------------------------------------------------- 475 ! nameos ocean physical parameters 476 !----------------------------------------------------------------------- 477 ! neos type of equation of state and Brunt-Vaisala frequency 478 ! ! = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 479 ! ! = 1, linear: rho(T) = rau0 * ( 1.028 - ralpha * T ) 480 ! ! = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 481 ! ralpha thermal expension coefficient (neos= 1 or 2) 482 ! rbeta saline expension coefficient (neos= 2) 483 &nameos 484 neos = 0 485 ralpha = 2.e-4 486 rbeta = 0.001 487 / 488 !----------------------------------------------------------------------- 489 ! nam_traadv advection scheme for tracer 490 !----------------------------------------------------------------------- 491 ! ln_traadv_cen2 2nd order centered scheme (default T) 492 ! ln_traadv_tvd TVD scheme (default F) 493 ! ln_traadv_muscl MUSCL scheme (default F) 494 ! ln_traadv_muscl2 MUSCL scheme + cen2 at boundaries (default F) 495 ! ln_traadv_ubs UBS scheme (default F) 496 ! ln_traadv_qck QUICKEST scheme (default F) 497 &nam_traadv 498 ln_traadv_cen2 = .true. 499 ln_traadv_tvd = .false. 500 ln_traadv_muscl = .false. 501 ln_traadv_muscl2 = .false. 502 ln_traadv_ubs = .false. 503 ln_traadv_qck = .false. 504 / 505 !----------------------------------------------------------------------- 506 ! nam_traldf lateral diffusion scheme for tracer 507 !----------------------------------------------------------------------- 508 ! Type of the operator : 509 ! ln_traldf_lap laplacian operator (default T) 510 ! ln_traldf_bilap bilaplacian operator (default F) 511 ! Direction of action : 512 ! ln_traldf_level iso-level (default F) 513 ! ln_traldf_hor horizontal (geopotential) (default F) (require "key_ldfslp" when ln_sco=T) 514 ! ln_traldf_iso iso-neutral (default T) (require "key_ldfslp") 515 ! Coefficient 516 ! aht0 horizontal eddy diffusivity for tracers (m2/s) 517 ! ahtb0 background eddy diffusivity for isopycnal diffusion (m2/s) 518 ! aeiv0 eddy induced velocity coefficient (m2/s) (require "key_traldf_eiv") 519 &nam_traldf 520 ln_traldf_lap = .true. 521 ln_traldf_bilap = .false. 522 ln_traldf_level = .false. 523 ln_traldf_hor = .false. 524 ln_traldf_iso = .true. 525 aht0 = 2000. 526 ahtb0 = 0. 527 aeiv0 = 2000. 528 / 529 !----------------------------------------------------------------------- 530 ! namtdp tracer newtonian damping ('key_tradmp') 531 !----------------------------------------------------------------------- 532 ! ndmp type of damping in temperature and salinity 533 ! ! ='latitude', damping poleward of 'ndmp' degrees and function 534 ! ! of the distance-to-coast. Red and Med Seas as ndmp=-1 535 ! ! =-1 damping only in Med and Red Seas 536 ! ndmpf =1 create a damping.coeff NetCDF file (the 3D damping array) 537 ! nmldmp type of damping in the mixed layer 538 ! ! =0 damping throughout the water column 539 ! ! =1 no damping in the mixed layer defined by avt >5cm2/s 540 ! ! =2 no damping in the mixed layer defined rho<rho(surf)+.01 541 ! sdmp surface time scale for internal damping (days) 542 ! bdmp bottom time scale for internal damping (days) 543 ! hdmp depth of transition between sdmp and bdmp (meters) 544 &namtdp 545 ndmp = -1 546 ndmpf = 1 547 nmldmp = 1 548 sdmp = 50. 549 bdmp = 360. 550 hdmp = 800. 551 / 552 553 !!====================================================================== 554 !! *** Dynamics namelists *** 555 !!====================================================================== 556 !! nam_dynadv formulation of the momentum advection 557 !! nam_dynvor advection scheme 558 !! nam_dynhpg hydrostatic pressure gradient 559 !! namflg hydrostatic pressure gradient time stepping 560 !! nam_dynspg surface pressure gradient (CPP key only) 561 !! nam_dynldf lateral diffusion scheme 562 !!====================================================================== 563 564 !----------------------------------------------------------------------- 565 ! nam_dynadv formulation of the momentum advection 566 !----------------------------------------------------------------------- 567 ! ln_dynadv_vec vector form (T) or flux form (F) (default T) 568 ! ln_dynadv_cen2 flux form - 2nd order centered scheme (default T) 569 ! ln_dynadv_ubs flux form - 3rd order UBS scheme (default F) 570 &nam_dynadv 571 ln_dynadv_vec = .true. 572 ln_dynadv_cen2 = .false. 573 ln_dynadv_ubs = .false. 574 / 575 !----------------------------------------------------------------------- 576 ! nam_dynvor option of physics/algorithm (not control by CPP keys) 577 !----------------------------------------------------------------------- 578 ! ln_dynvor_ens vorticity trends: enstrophy conserving scheme (default T) 579 ! ln_dynvor_ene " " : energy conserving scheme (default F) 580 ! ln_dynvor_mix " " : mixed scheme (default F) 581 ! ln_dynvor_een " " : energy & enstrophy scheme (default F) 582 &nam_dynvor 583 ln_dynvor_ene = .false. 584 ln_dynvor_ens = .true. 585 ln_dynvor_mix = .false. 586 ln_dynvor_een = .false. 587 / 588 !----------------------------------------------------------------------- 589 ! nam_dynhpg Hydrostatic pressure gradient option 590 !----------------------------------------------------------------------- 591 ! type of pressure gradient scheme (choose one only!) 592 ! ln_hpg_zco z-coordinate - full steps (default T) 593 ! ln_hpg_zps z-coordinate - partial steps (interpolation) 594 ! ln_hpg_sco s-coordinate (standard jacobian formulation) 595 ! ln_hpg_hel s-coordinate (helsinki modification) 596 ! ln_hpg_wdj s-coordinate (weighted density jacobian) 597 ! ln_hpg_djc s-coordinate (Density Jacobian with Cubic polynomial) 598 ! ln_hpg_rot s-coordinate (ROTated axes scheme) 599 ! parameters 600 ! gamm weighting coefficient (wdj scheme) 601 &nam_dynhpg 602 ln_hpg_zco = .false. 603 ln_hpg_zps = .true. 604 ln_hpg_sco = .false. 605 ln_hpg_hel = .false. 606 ln_hpg_wdj = .false. 607 ln_hpg_djc = .false. 608 ln_hpg_rot = .false. 609 gamm = 0.e0 610 / 611 !----------------------------------------------------------------------- 612 ! namflg algorithm flags (algorithm not control by CPP keys) 613 !----------------------------------------------------------------------- 614 ! ln_dynhpg_imp hydrostatic pressure gradient: semi-implicit time scheme (T) 615 ! centered time scheme (F) 616 ! nn_dynhpg_rst add dynhpg implicit variables in restart ot not (1/0) 617 &namflg 618 ln_dynhpg_imp = .false. 619 nn_dynhpg_rst = 0 620 / 621 !----------------------------------------------------------------------- 622 ! nam_dynspg surface pressure gradient (CPP key only) 623 !----------------------------------------------------------------------- 624 ! ! explicit free surface ("key_dynspg_exp") 625 ! ! filtered free surface ("key_dynspg_flt") 626 ! ! split-explicit free surface ("key_dynspg_ts") 627 ! ! rigid-lid ("key_dynspg_rl") 628 629 !----------------------------------------------------------------------- 630 ! nam_dynldf lateral diffusion on momentum 631 !----------------------------------------------------------------------- 632 ! Type of the operator : 633 ! ln_dynldf_lap laplacian operator (default T) 634 ! ln_dynldf_bilap bilaplacian operator (default F) 635 ! Direction of action : 636 ! ln_dynldf_level iso-level (default F) 637 ! ln_dynldf_hor horizontal (geopotential) (default F) (require "key_ldfslp" if ln_sco=T) 638 ! ln_dynldf_iso iso-neutral (default T) (require "key_ldfslp") 639 ! Coefficient 640 ! ahm0 horizontal eddy viscosity for the dynamics (m2/s) 641 ! ahmb0 background eddy viscosity for isopycnal diffusion (m2/s) 642 &nam_dynldf 643 ln_dynldf_lap = .true. 644 ln_dynldf_bilap = .false. 645 ln_dynldf_level = .false. 646 ln_dynldf_hor = .true. 647 ln_dynldf_iso = .false. 648 ahm0 = 40000. 649 ahmb0 = 0. 650 / 651 652 !!====================================================================== 653 !! Tracers & Dynamics vertical physics namelists 654 !!====================================================================== 655 !! namzdf vertical physics 656 !! namnpc non penetrative convection 657 !! namric richardson number dependent vertical mixing ("key_zdfric" ) 658 !! namtke TKE dependent vertical mixing ("key_zdftke" ) 659 !! namkpp KPP dependent vertical mixing ("key_zdfkpp" ) 660 !! namddm double diffusive mixing parameterization ("key_zdfddm" ) 661 !!====================================================================== 662 663 !----------------------------------------------------------------------- 664 ! namzdf vertical physics 665 !----------------------------------------------------------------------- 666 ! ln_zdfevd enhanced vertical diffusion (default T) 667 ! ln_zdfnpc Non-Penetrative Convection (default T) 668 ! avm0 vertical eddy viscosity for the dynamic (m2/s) (also background Kz if not "key_zdfcst") 669 ! avt0 vertical eddy diffusivity for tracers (m2/s) (also background Kz if not "key_zdfcst") 670 ! avevd vertical coefficient for enhanced diffusion scheme (m2/s) 671 ! n_evdm = 0 apply enhanced mixing on tracer only 672 ! = 1 apply enhanced mixing on both tracer and momentum 673 ! ln_zdfexp vertical physics: (=T) time splitting (T) (Default=F) 674 ! (=F) euler backward (F) 675 ! n_zdfexp number of sub-timestep for time splitting scheme 676 &namzdf 677 ln_zdfevd = .true. 678 ln_zdfnpc = .false. 679 avm0 = 1.2e-4 680 avt0 = 1.2e-5 681 avevd = 100. 682 n_evdm = 1 683 ln_zdfexp = .false. 684 n_zdfexp = 3 685 / 686 !----------------------------------------------------------------------- 687 ! namnpc non penetrative convection 688 !----------------------------------------------------------------------- 689 ! nnpc1 non penetrative convective scheme frequency 690 ! nnpc2 non penetrative convective scheme print frequency 691 &namnpc 692 nnpc1 = 1 693 nnpc2 = 365 694 / 695 !----------------------------------------------------------------------- 696 ! namric richardson number dependent vertical diffusion ("key_zdfric" ) 697 !----------------------------------------------------------------------- 698 ! avmri maximum value of the vertical viscosity 699 ! alp coefficient of the parameterization 700 ! nric coefficient of the parameterization 701 &namric 702 avmri = 100.e-4 703 alp = 5. 704 nric = 2 705 / 706 !----------------------------------------------------------------------- 707 ! namtke turbulent eddy kinetic dependent vertical diffusion ("key_zdftke") 708 !----------------------------------------------------------------------- 709 ! ln_rstke flag to restart with tke from a run without tke (default F) 710 ! ediff coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) ) 711 ! ediss coef. of the Kolmogoroff dissipation 712 ! ebb coef. of the surface input of tke 713 ! efave coef. to applied to the tke diffusion ( avtke=efave*avm ) 714 ! emin minimum value of tke (m^2/s^2) 715 ! emin0 surface minimum value of tke (m^2/s^2) 716 ! nitke number of restart iterative loops 717 ! ri_c critic richardson number = 2/9 = 0.22222222 (hard coded) 718 ! nmxl length used = 0 bounded by the distance to surface and bottom 719 ! ! = 1 bounded by the local vertical scale factor 720 ! ! = 2 first vertical derivative of mixing length bounded by 1 721 ! npdl Prandtl number = 0 no vertical prandtl number (avt=avm) 722 ! ! = 1 prandtl number function of richarson number (avt=pdl*avm) 723 ! ! = 2 same as = 1 but a shapiro filter is applied on pdl 724 ! nave = horizontal averaged (=1) or not (=0) of avt (default =1) 725 ! navb = 0 cst background avt0, avm0 / =1 profile used on avtb 726 &namtke 727 ln_rstke = .false. 728 ediff = 0.1 729 ediss = 0.7 730 ebb = 3.75 731 efave = 1. 732 emin = 1.e-6 733 emin0 = 1.e-4 734 nitke = 50 735 nmxl = 2 736 npdl = 1 737 navb = 0 738 / 739 !----------------------------------------------------------------------- 740 ! namkpp K-Profile Parameterization dependent vertical mixing ("key_zdfkpp" ) 741 !----------------------------------------------------------------------- 742 ! ln_kpprimix shear instability mixing (default T) 743 ! difmiw constant internal wave viscosity (m2/s) 744 ! difsiw constant internal wave diffusivity (m2/s) 745 ! Riinfty local Richardson Number limit for shear instability 746 ! difri maximum shear mixing at Rig = 0 (m2/s) 747 ! bvsqcon Brunt-Vaisala squared (1/s**2) for maximum convection 748 ! difcon maximum mixing in interior convection (m2/s) 749 ! nave = 0/1 flag for horizontal average on avt, avmu, avmv 750 ! navb = 0/1 flag for constant or profile background avt 751 &namkpp 752 ln_kpprimix = .true. 753 difmiw = 1.e-04 754 difsiw = 0.1e-04 755 Riinfty = 0.8 756 difri = 0.0050 757 bvsqcon = -0.01e-07 758 difcon = 1. 759 navb = 0 760 nave = 1 761 / 762 !----------------------------------------------------------------------- 763 ! namddm double diffusive mixing parameterization ("key_zdfddm") 764 !----------------------------------------------------------------------- 765 ! avts maximum avs for dd mixing 766 ! hsbfr heat/salt buoyancy flux ratio 767 &namddm 768 avts = 1.e-4 769 hsbfr = 1.6 770 / 771 772 !!====================================================================== 773 !! *** Miscelaneous namelists *** 774 !!====================================================================== 775 !! namsol elliptic solver / island / free surface 776 !!====================================================================== 777 778 !----------------------------------------------------------------------- 779 ! namsol elliptic solver / island / free surface 780 !----------------------------------------------------------------------- 781 ! nsolv elliptic solver (=1 preconditioned conjugate gradient: pcg) 782 ! (=2 successive-over-relaxation: sor) 783 ! (=3 FETI currently it does not work! ("key_feti") 784 ! nsol_arp absolute/relative (0/1) precision convergence test 785 ! nmin minimum of iterations for the SOR solver 786 ! nmax maximum of iterations for the SOR solver 787 ! nmod frequency of test for the SOR solver 788 ! eps absolute precision of the solver 789 ! resmax absolute precision for the SOR solver 790 ! sor optimal coefficient for SOR solver 791 ! epsisl absolute precision on stream function solver 792 ! nmisl maximum pcg iterations for island 793 ! rnu strength of the additional force used in free surface b.c. 794 &namsol 795 nsolv = 1 796 nsol_arp = 0 797 nmin = 300 798 nmax = 800 799 nmod = 10 800 eps = 1.E-6 801 resmax = 1.E-10 802 sor = 1.92 803 epsisl = 1.e-10 804 nmisl = 4000 805 rnu = 1. 806 / 807 !!====================================================================== 808 !! *** Diagnostics namelists *** 809 !!====================================================================== 810 !! namtrd dynamics and/or tracer trends ("key_trddyn","key_trdtra","key_trdmld") 811 !! namgap level mean model-data gap ("key_diagap") 812 !! namspr surface pressure diagnosed in rigid-lid ("key_diaspr") 813 !! namflo float parameters ("key_float") 814 !! namptr Poleward Transport Diagnostics 815 !!====================================================================== 816 817 !----------------------------------------------------------------------- 818 ! namtrd diagnostics on dynamics and/or tracer trends ("key_trddyn" and/or "key_trdtra") 819 ! or mixed-layer trends ('key_trdmld') 820 ! or barotropic vorticity ("key_trdvor") 821 !----------------------------------------------------------------------- 822 ! ntrd time step frequency dynamics and tracers trends 823 ! nctls control surface type in mixed-layer trends (0,1 or n<jpk) 824 ! ln_trdmld_restart restart for ML diagnostics 825 ! ucf unit conversion factor (=1 -> /seconds | =86400. -> /day) 826 ! ln_trdmld_instant flag to diagnose trends of instantantaneous or mean ML T/S 827 &namtrd 828 ntrd = 365 829 nctls = 0 830 ln_trdmld_restart = .false. 831 ucf = 1. 832 ln_trdmld_instant = .false. 833 / 834 !----------------------------------------------------------------------- 835 ! namgap level mean model-data gap ('key_diagap') 836 !----------------------------------------------------------------------- 837 ! ngap time-step frequency of model-data gap computation 838 ! nprg time-step frequency of gap print in model output 839 &namgap 840 ngap = 15 841 nprg = 10 842 / 843 !----------------------------------------------------------------------- 844 ! namspr surface pressure diagnostic 845 !----------------------------------------------------------------------- 846 ! nmaxp maximum of iterations for the solver 847 ! epsp absolute precision of the solver 848 ! niterp number of iteration done by the solver 849 &namspr 850 nmaxp = 1000 851 epsp = 1.e-3 852 niterp = 400 853 / 854 !----------------------------------------------------------------------- 855 ! namflo float parameters ("key_float") 642 856 !----------------------------------------------------------------------- 643 857 ! ln_rstflo boolean term for float restart (true or false) … … 654 868 ln_flork4 = .false. 655 869 / 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 / 870 !----------------------------------------------------------------------- 871 ! namptr Poleward Transport Diagnostic 872 !----------------------------------------------------------------------- 873 ! ln_diaptr logical flag for Poleward transport computation 874 ! ln_subbas logical flag for Atlantic/Pacific/Indian basins computation 875 ! need input basins mask file named "subbasins.nc" 876 ! nf_ptr Frequency of computation 877 &namptr 878 ln_diaptr = .false. 879 ln_subbas = .false. 880 nf_ptr = 15 881 /
Note: See TracChangeset
for help on using the changeset viewer.