- Timestamp:
- 2015-12-14T10:27:28+01:00 (8 years ago)
- Location:
- branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/SHARED
- Files:
-
- 4 edited
- 1 copied
-
1_namelist_ref (copied) (copied from trunk/NEMOGCM/CONFIG/SHARED/1_namelist_ref)
-
domain_def.xml (modified) (2 diffs)
-
field_def.xml (modified) (11 diffs)
-
namelist_ref (modified) (20 diffs)
-
namelist_top_ref (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/SHARED/domain_def.xml
r5600 r6043 163 163 164 164 <domain id="grid_V" long_name="grid V"/> 165 165 166 <domain_group id="grid_W"> 166 167 <domain id="grid_W" long_name="grid W"/> … … 168 169 <domain id="EqW" zoom_ibegin="1" zoom_jbegin="0000" zoom_ni="0000" zoom_nj="1" /> 169 170 </domain_group> 171 172 <domain id="grid_F" long_name="grid F"/> 170 173 171 174 <domain_group id="scalarpoint"> -
branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/SHARED/field_def.xml
r5600 r6043 25 25 26 26 <field id="toce" long_name="temperature" standard_name="sea_water_potential_temperature" unit="degC" grid_ref="grid_T_3D"/> 27 <field id="toce_e3t" long_name="temperature * e3t" unit="degC*m"grid_ref="grid_T_3D" > toce * e3t </field >27 <field id="toce_e3t" long_name="temperature (thickness weighted)" unit="degC" grid_ref="grid_T_3D" > toce * e3t </field > 28 28 <field id="soce" long_name="salinity" standard_name="sea_water_practical_salinity" unit="1e-3" grid_ref="grid_T_3D"/> 29 <field id="soce_e3t" long_name="salinity * e3t" unit="1e-3*m" grid_ref="grid_T_3D" > soce * e3t </field > 29 <field id="soce_e3t" long_name="salinity (thickness weighted)" unit="1e-3" grid_ref="grid_T_3D" > soce * e3t </field > 30 31 <!-- t-eddy viscosity coefficients (ldfdyn) --> 32 <field id="ahmt_2d" long_name=" surface t-eddy viscosity coefficient" unit="m2/s or m4/s" /> 33 <field id="ahmt_3d" long_name=" 3D t-eddy viscosity coefficient" unit="m2/s or m4/s" grid_ref="grid_T_3D"/> 30 34 31 35 <field id="sst" long_name="sea surface temperature" standard_name="sea_surface_temperature" unit="degC" /> … … 367 371 <field id="utau" long_name="Wind Stress along i-axis" standard_name="surface_downward_x_stress" unit="N/m2" /> 368 372 <field id="uoce" long_name="ocean current along i-axis" standard_name="sea_water_x_velocity" unit="m/s" grid_ref="grid_U_3D" /> 369 <field id="uoce_e3u" long_name="ocean current along i-axis * e3u" unit="m2/s"grid_ref="grid_U_3D" > uoce * e3u </field>373 <field id="uoce_e3u" long_name="ocean current along i-axis (thickness weighted)" unit="m/s" grid_ref="grid_U_3D" > uoce * e3u </field> 370 374 <field id="ssu" long_name="ocean surface current along i-axis" unit="m/s" /> 371 375 <field id="sbu" long_name="ocean bottom current along i-axis" unit="m/s" /> … … 375 379 <field id="uoces" long_name="ocean transport along i-axis times salinity (CRS)" unit="1e-3*m/s" grid_ref="grid_U_3D" /> 376 380 381 <!-- u-eddy coefficients (ldftra) --> 382 <field id="ahtu_2d" long_name=" surface u-eddy diffusivity coefficient" unit="m2/s or m4/s" /> 383 <field id="ahtu_3d" long_name=" 3D u-EIV coefficient" unit="m2/s or m4/s" grid_ref="grid_U_3D"/> 384 <field id="aeiu_2d" long_name=" surface u-EIV coefficient" unit="m2/s" /> 385 <field id="aeiu_3d" long_name=" 3D u-EIV coefficient" unit="m2/s" grid_ref="grid_U_3D"/> 386 377 387 <!-- variables available with MLE --> 378 388 <field id="psiu_mle" long_name="MLE streamfunction along i-axis" unit="m3/s" grid_ref="grid_U_3D" /> 379 389 380 <!-- uoce_eiv: available with key_traldf_eiv and key_diaeiv-->390 <!-- uoce_eiv: available EIV --> 381 391 <field id="uoce_eiv" long_name="EIV ocean current along i-axis" standard_name="bolus_sea_water_x_velocity" unit="m/s" grid_ref="grid_U_3D" /> 382 392 … … 402 412 <field id="vtau" long_name="Wind Stress along j-axis" standard_name="surface_downward_y_stress" unit="N/m2" /> 403 413 <field id="voce" long_name="ocean current along j-axis" standard_name="sea_water_y_velocity" unit="m/s" grid_ref="grid_V_3D" /> 404 <field id="voce_e3v" long_name="ocean current along j-axis * e3v" unit="m2/s"grid_ref="grid_V_3D" > voce * e3v </field>414 <field id="voce_e3v" long_name="ocean current along j-axis (thickness weighted)" unit="m/s" grid_ref="grid_V_3D" > voce * e3v </field> 405 415 <field id="ssv" long_name="ocean surface current along j-axis" unit="m/s" /> 406 416 <field id="sbv" long_name="ocean bottom current along j-axis" unit="m/s" /> … … 410 420 <field id="voces" long_name="ocean transport along j-axis times salinity (CRS)" unit="1e-3*m/s" grid_ref="grid_V_3D" /> 411 421 422 <!-- v-eddy coefficients (ldftra, ldfdyn) --> 423 <field id="ahtv_2d" long_name=" surface v-eddy diffusivity coefficient" unit="m2/s or (m4/s)^1/2" /> 424 <field id="ahtv_3d" long_name=" 3D v-eddy diffusivity coefficient" unit="m2/s or (m4/s)^1/2" grid_ref="grid_V_3D"/> 425 <field id="aeiv_2d" long_name=" surface v-EIV coefficient" unit="m2/s" /> 426 <field id="aeiv_3d" long_name=" 3D v-EIV coefficient" unit="m2/s" grid_ref="grid_V_3D" /> 427 412 428 <!-- variables available with MLE --> 413 429 <field id="psiv_mle" long_name="MLE streamfunction along j-axis" unit="m3/s" grid_ref="grid_V_3D" /> 414 430 415 <!-- voce_eiv: available with key_traldf_eiv and key_diaeiv-->431 <!-- voce_eiv: available with EIV --> 416 432 <field id="voce_eiv" long_name="EIV ocean current along j-axis" standard_name="bolus_sea_water_y_velocity" unit="m/s" grid_ref="grid_V_3D" /> 417 433 … … 438 454 <field id="wocetr_eff" long_name="effective ocean vertical transport" unit="m3/s" /> 439 455 440 <!-- woce_eiv: available with key_traldf_eiv and key_diaeiv-->456 <!-- woce_eiv: available with EIV --> 441 457 <field id="woce_eiv" long_name="EIV ocean vertical velocity" standard_name="bolus_upward_sea_water_velocity" unit="m/s" /> 442 458 … … 449 465 450 466 <!-- avt_evd and avm_evd: available with ln_zdfevd --> 451 <field id="avt_evd" long_name="convective enhancement tovertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_convection" unit="m2/s" />452 <field id="avm_evd" long_name="convective enhancement tovertical viscosity" standard_name="ocean_vertical_momentum_diffusivity_due_to_convection" unit="m2/s" />467 <field id="avt_evd" long_name="convective enhancement of vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_convection" unit="m2/s" /> 468 <field id="avm_evd" long_name="convective enhancement of vertical viscosity" standard_name="ocean_vertical_momentum_diffusivity_due_to_convection" unit="m2/s" /> 453 469 454 470 <!-- avt_tide: available with key_zdftmx --> … … 459 475 <field id="w_masstr2" long_name="square of vertical mass transport" standard_name="square_of_upward_ocean_mass_transport" unit="kg2/s2" /> 460 476 461 <!-- aht2d and aht2d_eiv : available with key_traldf_eiv and key_traldf_c2d-->477 <!-- aht2d and aht2d_eiv --> 462 478 <field id="aht2d" long_name="lateral eddy diffusivity" standard_name="ocean_tracer_xy_laplacian_diffusivity" unit="m2/s" grid_ref="grid_W_2D" /> 463 479 <field id="aht2d_eiv" long_name="EIV lateral eddy diffusivity" standard_name="ocean_tracer_bolus_laplacian_diffusivity" unit="m2/s" grid_ref="grid_W_2D" /> 464 480 </field_group> 465 481 482 <!-- F grid --> 483 <!-- f-eddy viscosity coefficients (ldfdyn) --> 484 <field id="ahmf_2d" long_name=" surface f-eddy viscosity coefficient" unit="m2/s or m4/s" /> 485 <field id="ahmf_3d" long_name=" 3D f-eddy viscosity coefficient" unit="m2/s or m4/s" grid_ref="grid_T_3D"/> 486 466 487 <!-- scalar variables available with key_diaar5 --> 467 488 … … 651 672 <field id="ketrd_hpg" long_name="ke-trend: hydrostatic pressure gradient" unit="W/s^3" /> 652 673 <field id="ketrd_spg" long_name="ke-trend: surface pressure gradient" unit="W/s^3" /> 653 <field id="ketrd_spgexp" long_name="ke-trend: surface pressure gradient (explicit)" unit="W/s^3" />654 <field id="ketrd_spgflt" long_name="ke-trend: surface pressure gradient (filter)" unit="W/s^3" />655 <field id="ssh_flt" long_name="filtered contribution to ssh (dynspg_flt)" unit="m" grid_ref="grid_T_2D" />656 674 <field id="w0" long_name="surface vertical velocity" unit="m/s" grid_ref="grid_T_2D" /> 657 675 <field id="pw0_exp" long_name="surface pressure flux due to ssh" unit="W/s^2" grid_ref="grid_T_2D" /> … … 695 713 <field id="utrd_hpg" long_name="i-trend: hydrostatic pressure gradient" unit="m/s^2" /> 696 714 <field id="utrd_spg" long_name="i-trend: surface pressure gradient" unit="m/s^2" /> 697 <field id="utrd_spgexp" long_name="i-trend: surface pressure gradient (explicit)" unit="m/s^2" />698 <field id="utrd_spgflt" long_name="i-trend: surface pressure gradient (filtered)" unit="m/s^2" />699 715 <field id="utrd_keg" long_name="i-trend: KE gradient or hor. adv." unit="m/s^2" /> 700 716 <field id="utrd_rvo" long_name="i-trend: relative vorticity or metric term" unit="m/s^2" /> … … 715 731 <field id="vtrd_hpg" long_name="j-trend: hydrostatic pressure gradient" unit="m/s^2" /> 716 732 <field id="vtrd_spg" long_name="j-trend: surface pressure gradient" unit="m/s^2" /> 717 <field id="vtrd_spgexp" long_name="j-trend: surface pressure gradient (explicit)" unit="m/s^2" />718 <field id="vtrd_spgflt" long_name="j-trend: surface pressure gradient (filtered)" unit="m/s^2" />719 733 <field id="vtrd_keg" long_name="j-trend: KE gradient or hor. adv." unit="m/s^2" /> 720 734 <field id="vtrd_rvo" long_name="j-trend: relative vorticity or metric term" unit="m/s^2" /> -
branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/SHARED/namelist_ref
r5600 r6043 1 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 !! namelist_ref 1 3 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 4 !! NEMO/OPA : 1 - run manager (namrun) … … 5 7 !! namsbc_cpl, namtra_qsr, namsbc_rnf, 6 8 !! namsbc_apr, namsbc_ssr, namsbc_alb) 7 !! 4 - lateral boundary (namlbc, nam cla, namobc, namagrif, nambdy, nambdy_tide)9 !! 4 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) 8 10 !! 5 - bottom boundary (nambfr, nambbc, nambbl) 9 !! 6 - Tracer (nameos, namtra_adv, namtra_ldf, namtra_ dmp)11 !! 6 - Tracer (nameos, namtra_adv, namtra_ldf, namtra_ldfeiv, namtra_dmp) 10 12 !! 7 - dynamics (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) 11 !! 8 - Verical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_ kpp, namzdf_ddm, namzdf_tmx)13 !! 8 - Verical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_ddm, namzdf_tmx) 12 14 !! 9 - diagnostics (namnc4, namtrd, namspr, namflo, namhsb, namsto) 13 !! 10 - miscellaneous (nam sol, nammpp, namctl)15 !! 10 - miscellaneous (nammpp, namctl) 14 16 !! 11 - Obs & Assim (namobs, nam_asminc) 15 17 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 48 50 ln_mskland = .false. ! mask land points in NetCDF outputs (costly: + ~15%) 49 51 ln_cfmeta = .false. ! output additional data to netCDF files required for compliance with the CF metadata standard 50 ln_clobber = . false. ! clobber (overwrite) an existing file52 ln_clobber = .true. ! clobber (overwrite) an existing file 51 53 nn_chunksz = 0 ! chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) 52 54 / … … 161 163 / 162 164 !----------------------------------------------------------------------- 163 &namsplit ! time splitting parameters ("key_dynspg_ts")164 !-----------------------------------------------------------------------165 ln_bt_fw = .TRUE. ! Forward integration of barotropic equations166 ln_bt_av = .TRUE. ! Time filtering of barotropic variables167 ln_bt_nn_auto = .TRUE. ! Set nn_baro automatically to be just below168 ! a user defined maximum courant number (rn_bt_cmax)169 nn_baro = 30 ! Number of iterations of barotropic mode170 ! during rn_rdt seconds. Only used if ln_bt_nn_auto=F171 rn_bt_cmax = 0.8 ! Maximum courant number allowed if ln_bt_nn_auto=T172 nn_bt_flt = 1 ! Time filter choice173 ! = 0 None174 ! = 1 Boxcar over nn_baro barotropic steps175 ! = 2 Boxcar over 2*nn_baro " "176 /177 !-----------------------------------------------------------------------178 165 &namcrs ! Grid coarsening for dynamics output and/or 179 166 ! passive tracer coarsened online simulations … … 200 187 !----------------------------------------------------------------------- 201 188 &namtsd ! data : Temperature & Salinity 202 !-----------------------------------------------------------------------203 189 !----------------------------------------------------------------------- 204 190 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! … … 545 531 !!====================================================================== 546 532 !! namlbc lateral momentum boundary condition 547 !! namcla cross land advection548 533 !! namobc open boundaries parameters ("key_obc") 549 534 !! namagrif agrif nested grid ( read by child model only ) ("key_agrif") … … 560 545 / 561 546 !----------------------------------------------------------------------- 562 &namcla ! cross land advection563 !-----------------------------------------------------------------------564 nn_cla = 0 ! advection between 2 ocean pts separates by land565 /566 !-----------------------------------------------------------------------567 &namobc ! open boundaries parameters ("key_obc")568 !-----------------------------------------------------------------------569 ln_obc_clim = .false. ! climatological obc data files (T) or not (F)570 ln_vol_cst = .true. ! impose the total volume conservation (T) or not (F)571 ln_obc_fla = .false. ! Flather open boundary condition572 nn_obcdta = 1 ! = 0 the obc data are equal to the initial state573 ! = 1 the obc data are read in 'obc.dta' files574 cn_obcdta = 'annual' ! set to annual if obc datafile hold 1 year of data575 ! set to monthly if obc datafile hold 1 month of data576 rn_dpein = 1. ! damping time scale for inflow at east open boundary577 rn_dpwin = 1. ! - - - west - -578 rn_dpnin = 1. ! - - - north - -579 rn_dpsin = 1. ! - - - south - -580 rn_dpeob = 3000. ! time relaxation (days) for the east open boundary581 rn_dpwob = 15. ! - - - west - -582 rn_dpnob = 3000. ! - - - north - -583 rn_dpsob = 15. ! - - - south - -584 rn_volemp = 1. ! = 0 the total volume change with the surface flux (E-P-R)585 ! = 1 the total volume remains constant586 /587 !-----------------------------------------------------------------------588 547 &namagrif ! AGRIF zoom ("key_agrif") 589 548 !----------------------------------------------------------------------- … … 592 551 rn_sponge_tra = 2880. ! coefficient for tracer sponge layer [m2/s] 593 552 rn_sponge_dyn = 2880. ! coefficient for dynamics sponge layer [m2/s] 553 ln_chk_bathy = .FALSE. ! 594 554 / 595 555 !----------------------------------------------------------------------- … … 710 670 ! = 2 variable flux (read in geothermal_heating.nc in mW/m2) 711 671 rn_geoflx_cst = 86.4e-3 ! Constant value of geothermal heat flux [W/m2] 712 713 672 / 714 673 !----------------------------------------------------------------------- … … 724 683 !! Tracer (T & S ) namelists 725 684 !!====================================================================== 726 !! nameos equation of state727 !! namtra_adv advection scheme685 !! nameos equation of state 686 !! namtra_adv advection scheme 728 687 !! namtra_adv_mle mixed layer eddy param. (Fox-Kemper param.) 729 !! namtra_ldf lateral diffusion scheme 730 !! namtra_dmp T & S newtonian damping 688 !! namtra_ldf lateral diffusion scheme 689 !! namtra_ldfeiv eddy induced velocity param. 690 !! namtra_dmp T & S newtonian damping 731 691 !!====================================================================== 732 692 ! … … 739 699 ! = 1, S-EOS (simplified eos) 740 700 ln_useCT = .true. ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 741 !!701 ! 742 702 ! ! S-EOS coefficients : 743 !! rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS703 ! rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 744 704 rn_a0 = 1.6550e-1 ! thermal expension coefficient (nn_eos= 1) 745 705 rn_b0 = 7.6554e-1 ! saline expension coefficient (nn_eos= 1) … … 753 713 &namtra_adv ! advection scheme for tracer 754 714 !----------------------------------------------------------------------- 755 ln_traadv_cen2 = .false. ! 2nd order centered scheme 756 ln_traadv_tvd = .true. ! TVD scheme 757 ln_traadv_muscl = .false. ! MUSCL scheme 758 ln_traadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 759 ln_traadv_ubs = .false. ! UBS scheme 760 ln_traadv_qck = .false. ! QUICKEST scheme 761 ln_traadv_msc_ups= .false. ! use upstream scheme within muscl 762 ln_traadv_tvd_zts= .false. ! TVD scheme with sub-timestepping of vertical tracer advection 715 ln_traadv_cen = .false. ! 2nd order centered scheme 716 nn_cen_h = 4 ! =2/4, horizontal 2nd order CEN / 4th order CEN 717 nn_cen_v = 4 ! =2/4, vertical 2nd order CEN / 4th order COMPACT 718 ln_traadv_fct = .false. ! FCT scheme 719 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 720 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 721 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping 722 ! ! (number of sub-timestep = nn_fct_zts) 723 ln_traadv_mus = .false. ! MUSCL scheme 724 ln_mus_ups = .false. ! use upstream scheme near river mouths 725 ln_traadv_ubs = .false. ! UBS scheme 726 nn_ubs_v = 2 ! =2 , vertical 2nd order FCT 727 ln_traadv_qck = .false. ! QUICKEST scheme 763 728 / 764 729 !----------------------------------------------------------------------- 765 730 &namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param) 766 731 !----------------------------------------------------------------------- 767 ln_mle = . true. ! (T) use the Mixed Layer Eddy (MLE) parameterisation732 ln_mle = .false. ! (T) use the Mixed Layer Eddy (MLE) parameterisation 768 733 rn_ce = 0.06 ! magnitude of the MLE (typical value: 0.06 to 0.08) 769 734 nn_mle = 1 ! MLE type: =0 standard Fox-Kemper ; =1 new formulation … … 779 744 !---------------------------------------------------------------------------------- 780 745 ! ! Operator type: 781 ln_traldf_lap = .true. ! laplacian operator 782 ln_traldf_bilap = .false. ! bilaplacian operator 746 ! ! no diffusion: set ln_traldf_lap=..._blp=F 747 ln_traldf_lap = .false. ! laplacian operator 748 ln_traldf_blp = .false. ! bilaplacian operator 783 749 ! ! Direction of action: 784 ln_traldf_level = .false. ! iso-level 785 ln_traldf_hor = .false. ! horizontal (geopotential) (needs "key_ldfslp" when ln_sco=T) 786 ln_traldf_iso = .true. ! iso-neutral (needs "key_ldfslp") 787 ! ! Griffies parameters (all need "key_ldfslp") 788 ln_traldf_grif = .false. ! use griffies triads 789 ln_traldf_gdia = .false. ! output griffies eddy velocities 790 ln_triad_iso = .false. ! pure lateral mixing in ML 791 ln_botmix_grif = .false. ! lateral mixing on bottom 792 ! ! Coefficients 793 ! Eddy-induced (GM) advection always used with Griffies; otherwise needs "key_traldf_eiv" 794 ! Value rn_aeiv_0 is ignored unless = 0 with Held-Larichev spatially varying aeiv 795 ! (key_traldf_c2d & key_traldf_eiv & key_orca_r2, _r1 or _r05) 796 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] 797 rn_aht_0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] 798 rn_ahtb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 799 ! (normally=0; not used with Griffies) 800 rn_slpmax = 0.01 ! slope limit 801 rn_chsmag = 1. ! multiplicative factor in Smagorinsky diffusivity 802 rn_smsh = 1. ! Smagorinsky diffusivity: = 0 - use only sheer 803 rn_aht_m = 2000. ! upper limit or stability criteria for lateral eddy diffusivity (m2/s) 750 ln_traldf_lev = .false. ! iso-level 751 ln_traldf_hor = .false. ! horizontal (geopotential) 752 ln_traldf_iso = .false. ! iso-neutral (standard operator) 753 ln_traldf_triad = .false. ! iso-neutral (triad operator) 754 ! 755 ! ! iso-neutral options: 756 ln_traldf_msc = .false. ! Method of Stabilizing Correction (both operators) 757 rn_slpmax = 0.01 ! slope limit (both operators) 758 ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only) 759 rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only) 760 ln_botmix_triad = .false. ! lateral mixing on bottom (triad only) 761 ! 762 ! ! Coefficients: 763 nn_aht_ijk_t = 0 ! space/time variation of eddy coef 764 ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file 765 ! ! = 0 constant 766 ! ! = 10 F(k) =ldf_c1d 767 ! ! = 20 F(i,j) =ldf_c2d 768 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 769 ! ! = 30 F(i,j,k) =ldf_c2d * ldf_c1d 770 ! ! = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 771 rn_aht_0 = 2000. ! lateral eddy diffusivity (lap. operator) [m2/s] 772 rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] 773 / 774 !---------------------------------------------------------------------------------- 775 &namtra_ldfeiv ! eddy induced velocity param. 776 !---------------------------------------------------------------------------------- 777 ln_ldfeiv =.false. ! use eddy induced velocity parameterization 778 ln_ldfeiv_dia =.false. ! diagnose eiv stream function and velocities 779 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] 780 nn_aei_ijk_t = 21 ! space/time variation of the eiv coeficient 781 ! ! =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file 782 ! ! = 0 constant 783 ! ! = 10 F(k) =ldf_c1d 784 ! ! = 20 F(i,j) =ldf_c2d 785 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 786 ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d 804 787 / 805 788 !----------------------------------------------------------------------- … … 819 802 !! namdyn_vor advection scheme 820 803 !! namdyn_hpg hydrostatic pressure gradient 821 !! namdyn_spg surface pressure gradient (CPP key only)804 !! namdyn_spg surface pressure gradient 822 805 !! namdyn_ldf lateral diffusion scheme 823 806 !!====================================================================== … … 852 835 ln_dynvor_ens = .false. ! energy conserving scheme 853 836 ln_dynvor_mix = .false. ! mixed scheme 854 ln_dynvor_een = .true. ! energy & enstrophy scheme 855 ln_dynvor_een_old = .false. ! energy & enstrophy scheme - original formulation 837 ln_dynvor_een = .false. ! energy & enstrophy scheme 838 nn_een_e3f = 1 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 839 ln_dynvor_msk = .false. ! vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) ! PLEASE DO NOT USE 856 840 / 857 841 !----------------------------------------------------------------------- … … 864 848 ln_hpg_djc = .false. ! s-coordinate (Density Jacobian with Cubic polynomial) 865 849 ln_hpg_prj = .false. ! s-coordinate (Pressure Jacobian scheme) 866 ln_dynhpg_imp = .false. ! time stepping: semi-implicit time scheme (T) 867 ! centered time scheme (F) 868 / 869 !----------------------------------------------------------------------- 870 !namdyn_spg ! surface pressure gradient (CPP key only) 871 !----------------------------------------------------------------------- 872 ! ! explicit free surface ("key_dynspg_exp") 873 ! ! filtered free surface ("key_dynspg_flt") 874 ! ! split-explicit free surface ("key_dynspg_ts") 875 850 / 851 !----------------------------------------------------------------------- 852 &namdyn_spg ! surface pressure gradient 853 !----------------------------------------------------------------------- 854 ln_dynspg_exp = .false. ! explicit free surface 855 ln_dynspg_ts = .false. ! split-explicit free surface 856 ln_bt_fw = .true. ! Forward integration of barotropic Eqs. 857 ln_bt_av = .true. ! Time filtering of barotropic variables 858 nn_bt_flt = 1 ! Time filter choice = 0 None 859 ! ! = 1 Boxcar over nn_baro sub-steps 860 ! ! = 2 Boxcar over 2*nn_baro " " 861 ln_bt_auto = .true. ! Number of sub-step defined from: 862 rn_bt_cmax = 0.8 ! =T : the Maximum Courant Number allowed 863 nn_baro = 30 ! =F : the number of sub-step in rn_rdt seconds 864 / 876 865 !----------------------------------------------------------------------- 877 866 &namdyn_ldf ! lateral diffusion on momentum 878 867 !----------------------------------------------------------------------- 879 868 ! ! Type of the operator : 880 ln_dynldf_lap = .true. ! laplacian operator 881 ln_dynldf_bilap = .false. ! bilaplacian operator 869 ! ! no diffusion: set ln_dynldf_lap=..._blp=F 870 ln_dynldf_lap = .false. ! laplacian operator 871 ln_dynldf_blp = .false. ! bilaplacian operator 882 872 ! ! Direction of action : 883 ln_dynldf_lev el = .false.! iso-level884 ln_dynldf_hor = .true. ! horizontal (geopotential) (require "key_ldfslp" in s-coord.)885 ln_dynldf_iso = .false. ! iso-neutral (require "key_ldfslp")873 ln_dynldf_lev = .false. ! iso-level 874 ln_dynldf_hor = .false. ! horizontal (geopotential) 875 ln_dynldf_iso = .false. ! iso-neutral 886 876 ! ! Coefficient 887 rn_ahm_0_lap = 40000. ! horizontal laplacian eddy viscosity [m2/s] 888 rn_ahmb_0 = 0. ! background eddy viscosity for ldf_iso [m2/s] 889 rn_ahm_0_blp = 0. ! horizontal bilaplacian eddy viscosity [m4/s] 890 rn_cmsmag_1 = 3. ! constant in laplacian Smagorinsky viscosity 891 rn_cmsmag_2 = 3 ! constant in bilaplacian Smagorinsky viscosity 892 rn_cmsh = 1. ! 1 or 0 , if 0 -use only shear for Smagorinsky viscosity 893 rn_ahm_m_blp = -1.e12 ! upper limit for bilap abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp) 894 rn_ahm_m_lap = 40000. ! upper limit for lap ahm < min(dx^2/16rdt, rn_ahm_m_lap) 877 nn_ahm_ijk_t = 0 ! space/time variation of eddy coef 878 ! ! =-30 read in eddy_viscosity_3D.nc file 879 ! ! =-20 read in eddy_viscosity_2D.nc file 880 ! ! = 0 constant 881 ! ! = 10 F(k)=c1d 882 ! ! = 20 F(i,j)=F(grid spacing)=c2d 883 ! ! = 30 F(i,j,k)=c2d*c1d 884 ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) 885 rn_ahm_0 = 40000. ! horizontal laplacian eddy viscosity [m2/s] 886 rn_ahm_b = 0. ! background eddy viscosity for ldf_iso [m2/s] 887 rn_bhm_0 = 1.e+12 ! horizontal bilaplacian eddy viscosity [m4/s] 888 ! 889 ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 895 890 / 896 891 … … 901 896 !! namzdf_ric richardson number dependent vertical mixing ("key_zdfric") 902 897 !! namzdf_tke TKE dependent vertical mixing ("key_zdftke") 903 !! namzdf_kpp KPP dependent vertical mixing ("key_zdfkpp")904 898 !! namzdf_ddm double diffusive mixing parameterization ("key_zdfddm") 905 899 !! namzdf_tmx tidal mixing parameterization ("key_zdftmx") … … 963 957 ! = 1 0.5m at the equator to 30m poleward of 40 degrees 964 958 / 965 !------------------------------------------------------------------------966 &namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionally:967 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb")968 ln_kpprimix = .true. ! shear instability mixing969 rn_difmiw = 1.0e-04 ! constant internal wave viscosity [m2/s]970 rn_difsiw = 0.1e-04 ! constant internal wave diffusivity [m2/s]971 rn_riinfty = 0.8 ! local Richardson Number limit for shear instability972 rn_difri = 0.0050 ! maximum shear mixing at Rig = 0 [m2/s]973 rn_bvsqcon = -0.01e-07 ! Brunt-Vaisala squared for maximum convection [1/s2]974 rn_difcon = 1. ! maximum mixing in interior convection [m2/s]975 nn_avb = 0 ! horizontal averaged (=1) or not (=0) on avt and amv976 nn_ave = 1 ! constant (=0) or profile (=1) background on avt977 /978 959 !----------------------------------------------------------------------- 979 960 &namzdf_gls ! GLS vertical diffusion ("key_zdfgls") … … 1014 995 !! *** Miscellaneous namelists *** 1015 996 !!====================================================================== 1016 !! namsol elliptic solver / island / free surface1017 997 !! nammpp Massively Parallel Processing ("key_mpp_mpi) 1018 998 !! namctl Control prints & Benchmark … … 1023 1003 !!====================================================================== 1024 1004 ! 1025 !-----------------------------------------------------------------------1026 &namsol ! elliptic solver / island / free surface1027 !-----------------------------------------------------------------------1028 nn_solv = 1 ! elliptic solver: =1 preconditioned conjugate gradient (pcg)1029 ! =2 successive-over-relaxation (sor)1030 nn_sol_arp = 0 ! absolute/relative (0/1) precision convergence test1031 rn_eps = 1.e-6 ! absolute precision of the solver1032 nn_nmin = 300 ! minimum of iterations for the SOR solver1033 nn_nmax = 800 ! maximum of iterations for the SOR solver1034 nn_nmod = 10 ! frequency of test for the SOR solver1035 rn_resmax = 1.e-10 ! absolute precision for the SOR solver1036 rn_sor = 1.92 ! optimal coefficient for SOR solver (to be adjusted with the domain)1037 /1038 1005 !----------------------------------------------------------------------- 1039 1006 &nammpp ! Massively Parallel Processing ("key_mpp_mpi) -
branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/SHARED/namelist_top_ref
r5600 r6043 41 41 &namtrc_adv ! advection scheme for passive tracer 42 42 !----------------------------------------------------------------------- 43 ln_trcadv_cen2 = .false. ! 2nd order centered scheme 44 ln_trcadv_tvd = .true. ! TVD scheme 45 ln_trcadv_muscl = .false. ! MUSCL scheme 46 ln_trcadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 47 ln_trcadv_ubs = .false. ! UBS scheme 48 ln_trcadv_qck = .false. ! QUICKEST scheme 49 ln_trcadv_msc_ups = .false. ! use upstream scheme within muscl 43 ln_trcadv_cen = .false. ! 2nd order centered scheme 44 nn_cen_h = 4 ! =2/4, horizontal 2nd order CEN / 4th order CEN 45 nn_cen_v = 4 ! =2/4, vertical 2nd order CEN / 4th order COMPACT 46 ln_trcadv_fct = .false. ! FCT scheme 47 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 48 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 49 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping 50 ! ! (number of sub-timestep = nn_fct_zts) 51 ln_trcadv_mus = .false. ! MUSCL scheme 52 ln_mus_ups = .false. ! use upstream scheme near river mouths 53 ln_trcadv_ubs = .false. ! UBS scheme 54 nn_ubs_v = 2 ! =2 , vertical 2nd order FCT 55 ln_trcadv_qck = .false. ! QUICKEST scheme 50 56 / 51 57 !----------------------------------------------------------------------- 52 58 &namtrc_ldf ! lateral diffusion scheme for passive tracer 53 59 !----------------------------------------------------------------------- 54 ! ! Type of the operator : 55 ln_trcldf_lap = .true. ! laplacian operator 56 ln_trcldf_bilap = .false. ! bilaplacian operator 57 ! Direction of action : 58 ln_trcldf_level = .false. ! iso-level 59 ln_trcldf_hor = .false. ! horizontal (geopotential) (require "key_ldfslp" when ln_sco=T) 60 ln_trcldf_iso = .true. ! iso-neutral (require "key_ldfslp") 61 ! ! Coefficient 62 rn_ahtrc_0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] 63 rn_ahtrb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 60 ! ! Type of the operator: 61 ln_trcldf_lap = .true. ! laplacian operator 62 ln_trcldf_blp = .false. ! bilaplacian operator 63 ! ! Direction of action: 64 ln_trcldf_lev = .false. ! iso-level 65 ln_trcldf_hor = .false. ! horizontal (geopotential) 66 ln_trcldf_iso = .true. ! iso-neutral (standard operator) 67 ln_trcldf_triad = .false. ! iso-neutral (triad operator) 68 ! ! Coefficient 69 rn_ahtrc_0 = 2000. ! lateral eddy diffusivity (lap. operator) [m2/s] 70 rn_bhtrc_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] 64 71 / 65 72 !-----------------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.
