Changeset 5956 for branches/2015/dev_r5151_UKMO_ISF/NEMOGCM/CONFIG/SHARED
- Timestamp:
- 2015-11-30T20:55:41+01:00 (9 years ago)
- Location:
- branches/2015/dev_r5151_UKMO_ISF/NEMOGCM/CONFIG/SHARED
- Files:
-
- 4 edited
- 1 copied
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5151_UKMO_ISF/NEMOGCM/CONFIG/SHARED/domain_def.xml
r5621 r5956 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/2015/dev_r5151_UKMO_ISF/NEMOGCM/CONFIG/SHARED/field_def.xml
r5621 r5956 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/2015/dev_r5151_UKMO_ISF/NEMOGCM/CONFIG/SHARED/namelist_ref
r5953 r5956 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 ! … … 543 529 !!====================================================================== 544 530 !! namlbc lateral momentum boundary condition 545 !! namcla cross land advection546 531 !! namobc open boundaries parameters ("key_obc") 547 532 !! namagrif agrif nested grid ( read by child model only ) ("key_agrif") … … 558 543 / 559 544 !----------------------------------------------------------------------- 560 &namcla ! cross land advection561 !-----------------------------------------------------------------------562 nn_cla = 0 ! advection between 2 ocean pts separates by land563 /564 !-----------------------------------------------------------------------565 &namobc ! open boundaries parameters ("key_obc")566 !-----------------------------------------------------------------------567 ln_obc_clim = .false. ! climatological obc data files (T) or not (F)568 ln_vol_cst = .true. ! impose the total volume conservation (T) or not (F)569 ln_obc_fla = .false. ! Flather open boundary condition570 nn_obcdta = 1 ! = 0 the obc data are equal to the initial state571 ! = 1 the obc data are read in 'obc.dta' files572 cn_obcdta = 'annual' ! set to annual if obc datafile hold 1 year of data573 ! set to monthly if obc datafile hold 1 month of data574 rn_dpein = 1. ! damping time scale for inflow at east open boundary575 rn_dpwin = 1. ! - - - west - -576 rn_dpnin = 1. ! - - - north - -577 rn_dpsin = 1. ! - - - south - -578 rn_dpeob = 3000. ! time relaxation (days) for the east open boundary579 rn_dpwob = 15. ! - - - west - -580 rn_dpnob = 3000. ! - - - north - -581 rn_dpsob = 15. ! - - - south - -582 rn_volemp = 1. ! = 0 the total volume change with the surface flux (E-P-R)583 ! = 1 the total volume remains constant584 /585 !-----------------------------------------------------------------------586 545 &namagrif ! AGRIF zoom ("key_agrif") 587 546 !----------------------------------------------------------------------- … … 590 549 rn_sponge_tra = 2880. ! coefficient for tracer sponge layer [m2/s] 591 550 rn_sponge_dyn = 2880. ! coefficient for dynamics sponge layer [m2/s] 551 ln_chk_bathy = .FALSE. ! 592 552 / 593 553 !----------------------------------------------------------------------- … … 708 668 ! = 2 variable flux (read in geothermal_heating.nc in mW/m2) 709 669 rn_geoflx_cst = 86.4e-3 ! Constant value of geothermal heat flux [W/m2] 710 711 670 / 712 671 !----------------------------------------------------------------------- … … 722 681 !! Tracer (T & S ) namelists 723 682 !!====================================================================== 724 !! nameos equation of state725 !! namtra_adv advection scheme683 !! nameos equation of state 684 !! namtra_adv advection scheme 726 685 !! namtra_adv_mle mixed layer eddy param. (Fox-Kemper param.) 727 !! namtra_ldf lateral diffusion scheme 728 !! namtra_dmp T & S newtonian damping 686 !! namtra_ldf lateral diffusion scheme 687 !! namtra_ldfeiv eddy induced velocity param. 688 !! namtra_dmp T & S newtonian damping 729 689 !!====================================================================== 730 690 ! … … 737 697 ! = 1, S-EOS (simplified eos) 738 698 ln_useCT = .true. ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 739 !!699 ! 740 700 ! ! S-EOS coefficients : 741 !! rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS701 ! rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 742 702 rn_a0 = 1.6550e-1 ! thermal expension coefficient (nn_eos= 1) 743 703 rn_b0 = 7.6554e-1 ! saline expension coefficient (nn_eos= 1) … … 751 711 &namtra_adv ! advection scheme for tracer 752 712 !----------------------------------------------------------------------- 753 ln_traadv_cen2 = .false. ! 2nd order centered scheme 754 ln_traadv_tvd = .true. ! TVD scheme 755 ln_traadv_muscl = .false. ! MUSCL scheme 756 ln_traadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 757 ln_traadv_ubs = .false. ! UBS scheme 758 ln_traadv_qck = .false. ! QUICKEST scheme 759 ln_traadv_msc_ups= .false. ! use upstream scheme within muscl 760 ln_traadv_tvd_zts= .false. ! TVD scheme with sub-timestepping of vertical tracer advection 713 ln_traadv_cen = .false. ! 2nd order centered scheme 714 nn_cen_h = 4 ! =2/4, horizontal 2nd order CEN / 4th order CEN 715 nn_cen_v = 4 ! =2/4, vertical 2nd order CEN / 4th order COMPACT 716 ln_traadv_fct = .false. ! FCT scheme 717 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 718 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 719 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping 720 ! ! (number of sub-timestep = nn_fct_zts) 721 ln_traadv_mus = .false. ! MUSCL scheme 722 ln_mus_ups = .false. ! use upstream scheme near river mouths 723 ln_traadv_ubs = .false. ! UBS scheme 724 nn_ubs_v = 2 ! =2 , vertical 2nd order FCT 725 ln_traadv_qck = .false. ! QUICKEST scheme 761 726 / 762 727 !----------------------------------------------------------------------- 763 728 &namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param) 764 729 !----------------------------------------------------------------------- 765 ln_mle = . true. ! (T) use the Mixed Layer Eddy (MLE) parameterisation730 ln_mle = .false. ! (T) use the Mixed Layer Eddy (MLE) parameterisation 766 731 rn_ce = 0.06 ! magnitude of the MLE (typical value: 0.06 to 0.08) 767 732 nn_mle = 1 ! MLE type: =0 standard Fox-Kemper ; =1 new formulation … … 777 742 !---------------------------------------------------------------------------------- 778 743 ! ! Operator type: 779 ln_traldf_lap = .true. ! laplacian operator 780 ln_traldf_bilap = .false. ! bilaplacian operator 744 ! ! no diffusion: set ln_traldf_lap=..._blp=F 745 ln_traldf_lap = .false. ! laplacian operator 746 ln_traldf_blp = .false. ! bilaplacian operator 781 747 ! ! Direction of action: 782 ln_traldf_level = .false. ! iso-level 783 ln_traldf_hor = .false. ! horizontal (geopotential) (needs "key_ldfslp" when ln_sco=T) 784 ln_traldf_iso = .true. ! iso-neutral (needs "key_ldfslp") 785 ! ! Griffies parameters (all need "key_ldfslp") 786 ln_traldf_grif = .false. ! use griffies triads 787 ln_traldf_gdia = .false. ! output griffies eddy velocities 788 ln_triad_iso = .false. ! pure lateral mixing in ML 789 ln_botmix_grif = .false. ! lateral mixing on bottom 790 ! ! Coefficients 791 ! Eddy-induced (GM) advection always used with Griffies; otherwise needs "key_traldf_eiv" 792 ! Value rn_aeiv_0 is ignored unless = 0 with Held-Larichev spatially varying aeiv 793 ! (key_traldf_c2d & key_traldf_eiv & key_orca_r2, _r1 or _r05) 794 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] 795 rn_aht_0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] 796 rn_ahtb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 797 ! (normally=0; not used with Griffies) 798 rn_slpmax = 0.01 ! slope limit 799 rn_chsmag = 1. ! multiplicative factor in Smagorinsky diffusivity 800 rn_smsh = 1. ! Smagorinsky diffusivity: = 0 - use only sheer 801 rn_aht_m = 2000. ! upper limit or stability criteria for lateral eddy diffusivity (m2/s) 748 ln_traldf_lev = .false. ! iso-level 749 ln_traldf_hor = .false. ! horizontal (geopotential) 750 ln_traldf_iso = .false. ! iso-neutral (standard operator) 751 ln_traldf_triad = .false. ! iso-neutral (triad operator) 752 ! 753 ! ! iso-neutral options: 754 ln_traldf_msc = .false. ! Method of Stabilizing Correction (both operators) 755 rn_slpmax = 0.01 ! slope limit (both operators) 756 ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only) 757 rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only) 758 ln_botmix_triad = .false. ! lateral mixing on bottom (triad only) 759 ! 760 ! ! Coefficients: 761 nn_aht_ijk_t = 0 ! space/time variation of eddy coef 762 ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file 763 ! ! = 0 constant 764 ! ! = 10 F(k) =ldf_c1d 765 ! ! = 20 F(i,j) =ldf_c2d 766 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 767 ! ! = 30 F(i,j,k) =ldf_c2d * ldf_c1d 768 ! ! = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 769 rn_aht_0 = 2000. ! lateral eddy diffusivity (lap. operator) [m2/s] 770 rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] 771 / 772 !---------------------------------------------------------------------------------- 773 &namtra_ldfeiv ! eddy induced velocity param. 774 !---------------------------------------------------------------------------------- 775 ln_ldfeiv =.false. ! use eddy induced velocity parameterization 776 ln_ldfeiv_dia =.false. ! diagnose eiv stream function and velocities 777 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] 778 nn_aei_ijk_t = 21 ! space/time variation of the eiv coeficient 779 ! ! =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file 780 ! ! = 0 constant 781 ! ! = 10 F(k) =ldf_c1d 782 ! ! = 20 F(i,j) =ldf_c2d 783 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 784 ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d 802 785 / 803 786 !----------------------------------------------------------------------- … … 817 800 !! namdyn_vor advection scheme 818 801 !! namdyn_hpg hydrostatic pressure gradient 819 !! namdyn_spg surface pressure gradient (CPP key only)802 !! namdyn_spg surface pressure gradient 820 803 !! namdyn_ldf lateral diffusion scheme 821 804 !!====================================================================== … … 850 833 ln_dynvor_ens = .false. ! energy conserving scheme 851 834 ln_dynvor_mix = .false. ! mixed scheme 852 ln_dynvor_een = .true. ! energy & enstrophy scheme 853 ln_dynvor_een_old = .false. ! energy & enstrophy scheme - original formulation 835 ln_dynvor_een = .false. ! energy & enstrophy scheme 836 nn_een_e3f = 1 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 837 ln_dynvor_msk = .false. ! vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) ! PLEASE DO NOT USE 854 838 / 855 839 !----------------------------------------------------------------------- … … 862 846 ln_hpg_djc = .false. ! s-coordinate (Density Jacobian with Cubic polynomial) 863 847 ln_hpg_prj = .false. ! s-coordinate (Pressure Jacobian scheme) 864 ln_dynhpg_imp = .false. ! time stepping: semi-implicit time scheme (T) 865 ! centered time scheme (F) 866 / 867 !----------------------------------------------------------------------- 868 !namdyn_spg ! surface pressure gradient (CPP key only) 869 !----------------------------------------------------------------------- 870 ! ! explicit free surface ("key_dynspg_exp") 871 ! ! filtered free surface ("key_dynspg_flt") 872 ! ! split-explicit free surface ("key_dynspg_ts") 873 848 / 849 !----------------------------------------------------------------------- 850 &namdyn_spg ! surface pressure gradient 851 !----------------------------------------------------------------------- 852 ln_dynspg_exp = .false. ! explicit free surface 853 ln_dynspg_ts = .false. ! split-explicit free surface 854 ln_bt_fw = .true. ! Forward integration of barotropic Eqs. 855 ln_bt_av = .true. ! Time filtering of barotropic variables 856 nn_bt_flt = 1 ! Time filter choice = 0 None 857 ! ! = 1 Boxcar over nn_baro sub-steps 858 ! ! = 2 Boxcar over 2*nn_baro " " 859 ln_bt_auto = .true. ! Number of sub-step defined from: 860 rn_bt_cmax = 0.8 ! =T : the Maximum Courant Number allowed 861 nn_baro = 30 ! =F : the number of sub-step in rn_rdt seconds 862 / 874 863 !----------------------------------------------------------------------- 875 864 &namdyn_ldf ! lateral diffusion on momentum 876 865 !----------------------------------------------------------------------- 877 866 ! ! Type of the operator : 878 ln_dynldf_lap = .true. ! laplacian operator 879 ln_dynldf_bilap = .false. ! bilaplacian operator 867 ! ! no diffusion: set ln_dynldf_lap=..._blp=F 868 ln_dynldf_lap = .false. ! laplacian operator 869 ln_dynldf_blp = .false. ! bilaplacian operator 880 870 ! ! Direction of action : 881 ln_dynldf_lev el = .false.! iso-level882 ln_dynldf_hor = .true. ! horizontal (geopotential) (require "key_ldfslp" in s-coord.)883 ln_dynldf_iso = .false. ! iso-neutral (require "key_ldfslp")871 ln_dynldf_lev = .false. ! iso-level 872 ln_dynldf_hor = .false. ! horizontal (geopotential) 873 ln_dynldf_iso = .false. ! iso-neutral 884 874 ! ! Coefficient 885 rn_ahm_0_lap = 40000. ! horizontal laplacian eddy viscosity [m2/s] 886 rn_ahmb_0 = 0. ! background eddy viscosity for ldf_iso [m2/s] 887 rn_ahm_0_blp = 0. ! horizontal bilaplacian eddy viscosity [m4/s] 888 rn_cmsmag_1 = 3. ! constant in laplacian Smagorinsky viscosity 889 rn_cmsmag_2 = 3 ! constant in bilaplacian Smagorinsky viscosity 890 rn_cmsh = 1. ! 1 or 0 , if 0 -use only shear for Smagorinsky viscosity 891 rn_ahm_m_blp = -1.e12 ! upper limit for bilap abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp) 892 rn_ahm_m_lap = 40000. ! upper limit for lap ahm < min(dx^2/16rdt, rn_ahm_m_lap) 875 nn_ahm_ijk_t = 0 ! space/time variation of eddy coef 876 ! ! =-30 read in eddy_viscosity_3D.nc file 877 ! ! =-20 read in eddy_viscosity_2D.nc file 878 ! ! = 0 constant 879 ! ! = 10 F(k)=c1d 880 ! ! = 20 F(i,j)=F(grid spacing)=c2d 881 ! ! = 30 F(i,j,k)=c2d*c1d 882 ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) 883 rn_ahm_0 = 40000. ! horizontal laplacian eddy viscosity [m2/s] 884 rn_ahm_b = 0. ! background eddy viscosity for ldf_iso [m2/s] 885 rn_bhm_0 = 1.e+12 ! horizontal bilaplacian eddy viscosity [m4/s] 886 ! 887 ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 893 888 / 894 889 … … 899 894 !! namzdf_ric richardson number dependent vertical mixing ("key_zdfric") 900 895 !! namzdf_tke TKE dependent vertical mixing ("key_zdftke") 901 !! namzdf_kpp KPP dependent vertical mixing ("key_zdfkpp")902 896 !! namzdf_ddm double diffusive mixing parameterization ("key_zdfddm") 903 897 !! namzdf_tmx tidal mixing parameterization ("key_zdftmx") … … 961 955 ! = 1 0.5m at the equator to 30m poleward of 40 degrees 962 956 / 963 !------------------------------------------------------------------------964 &namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionally:965 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb")966 ln_kpprimix = .true. ! shear instability mixing967 rn_difmiw = 1.0e-04 ! constant internal wave viscosity [m2/s]968 rn_difsiw = 0.1e-04 ! constant internal wave diffusivity [m2/s]969 rn_riinfty = 0.8 ! local Richardson Number limit for shear instability970 rn_difri = 0.0050 ! maximum shear mixing at Rig = 0 [m2/s]971 rn_bvsqcon = -0.01e-07 ! Brunt-Vaisala squared for maximum convection [1/s2]972 rn_difcon = 1. ! maximum mixing in interior convection [m2/s]973 nn_avb = 0 ! horizontal averaged (=1) or not (=0) on avt and amv974 nn_ave = 1 ! constant (=0) or profile (=1) background on avt975 /976 957 !----------------------------------------------------------------------- 977 958 &namzdf_gls ! GLS vertical diffusion ("key_zdfgls") … … 1012 993 !! *** Miscellaneous namelists *** 1013 994 !!====================================================================== 1014 !! namsol elliptic solver / island / free surface1015 995 !! nammpp Massively Parallel Processing ("key_mpp_mpi) 1016 996 !! namctl Control prints & Benchmark … … 1021 1001 !!====================================================================== 1022 1002 ! 1023 !-----------------------------------------------------------------------1024 &namsol ! elliptic solver / island / free surface1025 !-----------------------------------------------------------------------1026 nn_solv = 1 ! elliptic solver: =1 preconditioned conjugate gradient (pcg)1027 ! =2 successive-over-relaxation (sor)1028 nn_sol_arp = 0 ! absolute/relative (0/1) precision convergence test1029 rn_eps = 1.e-6 ! absolute precision of the solver1030 nn_nmin = 300 ! minimum of iterations for the SOR solver1031 nn_nmax = 800 ! maximum of iterations for the SOR solver1032 nn_nmod = 10 ! frequency of test for the SOR solver1033 rn_resmax = 1.e-10 ! absolute precision for the SOR solver1034 rn_sor = 1.92 ! optimal coefficient for SOR solver (to be adjusted with the domain)1035 /1036 1003 !----------------------------------------------------------------------- 1037 1004 &nammpp ! Massively Parallel Processing ("key_mpp_mpi) … … 1194 1161 ln_s3d = .false. ! Logical switch for S profile observations 1195 1162 ln_ena = .false. ! Logical switch for ENACT insitu data set 1196 ! ! ln_corLogical switch for Coriolis insitu data set1163 ln_cor = .false. ! Logical switch for Coriolis insitu data set 1197 1164 ln_profb = .false. ! Logical switch for feedback insitu data set 1198 1165 ln_sla = .false. ! Logical switch for SLA observations 1199 1200 1166 ln_sladt = .false. ! Logical switch for AVISO SLA data 1201 1202 1167 ln_slafb = .false. ! Logical switch for feedback SLA data 1203 ! ln_ssh Logical switch for SSH observations 1204 1205 ln_sst = .false. ! Logical switch for SST observations 1206 ln_reysst = .false. ! ln_reysst Logical switch for Reynolds observations 1207 ln_ghrsst = .false. ! ln_ghrsst Logical switch for GHRSST observations 1208 1168 ln_ssh = .false. ! Logical switch for SSH observations 1169 ln_sst = .false. ! Logical switch for SST observations 1170 ln_reysst = .false. ! Logical switch for Reynolds observations 1171 ln_ghrsst = .false. ! Logical switch for GHRSST observations 1209 1172 ln_sstfb = .false. ! Logical switch for feedback SST data 1210 ! ln_sssLogical switch for SSS observations1173 ln_sss = .false. ! Logical switch for SSS observations 1211 1174 ln_seaice = .false. ! Logical switch for Sea Ice observations 1212 ! ln_vel3d Logical switch for velocity observations 1213 ! ln_velavcur Logical switch for velocity daily av. cur. 1214 ! ln_velhrcur Logical switch for velocity high freq. cur. 1215 ! ln_velavadcp Logical switch for velocity daily av. ADCP 1216 ! ln_velhradcp Logical switch for velocity high freq. ADCP 1217 ! ln_velfb Logical switch for feedback velocity data 1218 ! ln_grid_global Global distribtion of observations 1219 ! ln_grid_search_lookup Logical switch for obs grid search w/lookup table 1220 ! grid_search_file Grid search lookup file header 1221 ! enactfiles ENACT input observation file names 1222 ! coriofiles Coriolis input observation file name 1223 ! ! profbfiles: Profile feedback input observation file name 1224 profbfiles = 'profiles_01.nc' 1225 ! ln_profb_enatim Enact feedback input time setting switch 1226 ! slafilesact Active SLA input observation file name 1227 ! slafilespas Passive SLA input observation file name 1228 ! ! slafbfiles: Feedback SLA input observation file name 1229 slafbfiles = 'sla_01.nc' 1230 ! sstfiles GHRSST input observation file name 1231 ! ! sstfbfiles: Feedback SST input observation file name 1232 sstfbfiles = 'sst_01.nc' 1233 ! seaicefiles Sea Ice input observation file names 1234 seaicefiles = 'seaice_01.nc' 1235 ! velavcurfiles Vel. cur. daily av. input file name 1236 ! velhvcurfiles Vel. cur. high freq. input file name 1237 ! velavadcpfiles Vel. ADCP daily av. input file name 1238 ! velhvadcpfiles Vel. ADCP high freq. input file name 1239 ! velfbfiles Vel. feedback input observation file name 1240 ! dobsini Initial date in window YYYYMMDD.HHMMSS 1241 ! dobsend Final date in window YYYYMMDD.HHMMSS 1242 ! n1dint Type of vertical interpolation method 1243 ! n2dint Type of horizontal interpolation method 1244 ! ln_nea Rejection of observations near land switch 1245 nmsshc = 0 ! MSSH correction scheme 1246 ! mdtcorr MDT correction 1247 ! mdtcutoff MDT cutoff for computed correction 1175 ln_vel3d = .false. ! Logical switch for velocity observations 1176 ln_velavcur= .false ! Logical switch for velocity daily av. cur. 1177 ln_velhrcur= .false ! Logical switch for velocity high freq. cur. 1178 ln_velavadcp = .false. ! Logical switch for velocity daily av. ADCP 1179 ln_velhradcp = .false. ! Logical switch for velocity high freq. ADCP 1180 ln_velfb = .false. ! Logical switch for feedback velocity data 1181 ln_grid_global = .false. ! Global distribtion of observations 1182 ln_grid_search_lookup = .false. ! Logical switch for obs grid search w/lookup table 1183 grid_search_file = 'grid_search' ! Grid search lookup file header 1184 ! All of the *files* variables below are arrays. Use namelist_cfg to add more files 1185 enactfiles = 'enact.nc' ! ENACT input observation file names (specify full array in namelist_cfg) 1186 coriofiles = 'corio.nc' ! Coriolis input observation file name 1187 profbfiles = 'profiles_01.nc' ! Profile feedback input observation file name 1188 ln_profb_enatim = .false ! Enact feedback input time setting switch 1189 slafilesact = 'sla_act.nc' ! Active SLA input observation file names 1190 slafilespas = 'sla_pass.nc' ! Passive SLA input observation file names 1191 slafbfiles = 'sla_01.nc' ! slafbfiles: Feedback SLA input observation file names 1192 sstfiles = 'ghrsst.nc' ! GHRSST input observation file names 1193 sstfbfiles = 'sst_01.nc' ! Feedback SST input observation file names 1194 seaicefiles = 'seaice_01.nc' ! Sea Ice input observation file names 1195 velavcurfiles = 'velavcurfile.nc' ! Vel. cur. daily av. input file name 1196 velhrcurfiles = 'velhrcurfile.nc' ! Vel. cur. high freq. input file name 1197 velavadcpfiles = 'velavadcpfile.nc' ! Vel. ADCP daily av. input file name 1198 velhradcpfiles = 'velhradcpfile.nc' ! Vel. ADCP high freq. input file name 1199 velfbfiles = 'velfbfile.nc' ! Vel. feedback input observation file name 1200 dobsini = 20000101.000000 ! Initial date in window YYYYMMDD.HHMMSS 1201 dobsend = 20010101.000000 ! Final date in window YYYYMMDD.HHMMSS 1202 n1dint = 0 ! Type of vertical interpolation method 1203 n2dint = 0 ! Type of horizontal interpolation method 1204 ln_nea = .false. ! Rejection of observations near land switch 1205 nmsshc = 0 ! MSSH correction scheme 1206 mdtcorr = 1.61 ! MDT correction 1207 mdtcutoff = 65.0 ! MDT cutoff for computed correction 1248 1208 ln_altbias = .false. ! Logical switch for alt bias 1249 1209 ln_ignmis = .true. ! Logical switch for ignoring missing files 1250 ! endailyavtypes ENACT daily average types1210 endailyavtypes = 820 ! ENACT daily average types - array (use namelist_cfg to set more values) 1251 1211 ln_grid_global = .true. 1252 1212 ln_grid_search_lookup = .false. -
branches/2015/dev_r5151_UKMO_ISF/NEMOGCM/CONFIG/SHARED/namelist_top_ref
r5621 r5956 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.