Changeset 5989 for branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/CONFIG/SHARED/namelist_ref
- Timestamp:
- 2015-12-03T09:10:32+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/CONFIG/SHARED/namelist_ref
r5260 r5989 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)12 !! 9 - diagnostics (namnc4, namtrd, namspr, namflo, namhsb )13 !! 10 - miscellaneous (nam sol, nammpp, namctl)13 !! 8 - Verical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_ddm, namzdf_tmx) 14 !! 9 - diagnostics (namnc4, namtrd, namspr, namflo, namhsb, namsto) 15 !! 10 - miscellaneous (nammpp, namctl) 14 16 !! 11 - Obs & Assim (namobs, nam_asminc) 15 17 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 37 39 ! = 2 nn_date0 read in restart ; nn_it000 : check consistancy between namelist and restart 38 40 cn_ocerst_in = "restart" ! suffix of ocean restart name (input) 41 cn_ocerst_indir = "." ! directory from which to read input ocean restarts 39 42 cn_ocerst_out = "restart" ! suffix of ocean restart name (output) 43 cn_ocerst_outdir = "." ! directory in which to write output ocean restarts 40 44 nn_istate = 0 ! output the initial state (1) or not (0) 45 ln_rst_list = .false. ! output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F) 41 46 nn_stock = 5475 ! frequency of creation of a restart file (modulo referenced to 1) 47 nn_stocklist = 0,0,0,0,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written 42 48 nn_write = 5475 ! frequency of write in the output file (modulo referenced to nn_it000) 43 49 ln_dimgnnn = .false. ! DIMG file format: 1 file for all processors (F) or by processor (T) 44 50 ln_mskland = .false. ! mask land points in NetCDF outputs (costly: + ~15%) 45 ln_clobber = .false. ! clobber (overwrite) an existing file 51 ln_cfmeta = .false. ! output additional data to netCDF files required for compliance with the CF metadata standard 52 ln_clobber = .true. ! clobber (overwrite) an existing file 46 53 nn_chunksz = 0 ! chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) 47 54 / … … 50 57 !! *** Domain namelists *** 51 58 !!====================================================================== 52 !! namcfg parameters of the configuration 59 !! namcfg parameters of the configuration 53 60 !! namzgr vertical coordinate 54 61 !! namzgr_sco s-coordinate or hybrid z-s-coordinate … … 58 65 ! 59 66 !----------------------------------------------------------------------- 60 &namcfg ! parameters of the configuration 67 &namcfg ! parameters of the configuration 61 68 !----------------------------------------------------------------------- 62 69 cp_cfg = "default" ! name of the configuration … … 72 79 jperio = 0 ! lateral cond. type (between 0 and 6) 73 80 ! = 0 closed ; = 1 cyclic East-West 74 ! = 2 equatorial symmetric ; = 3 North fold T-point pivot 81 ! = 2 equatorial symmetric ; = 3 North fold T-point pivot 75 82 ! = 4 cyclic East-West AND North fold T-point pivot 76 83 ! = 5 North fold F-point pivot 77 84 ! = 6 cyclic East-West AND North fold F-point pivot 78 ln_use_jattr = .false. ! use (T) the file attribute: open_ocean_jstart, if present 85 ln_use_jattr = .false. ! use (T) the file attribute: open_ocean_jstart, if present 79 86 ! in netcdf input files, as the start j-row for reading 80 87 / … … 101 108 !!!!!!! SH94 stretching coefficients (ln_s_sh94 = .true.) 102 109 rn_theta = 6.0 ! surface control parameter (0<=theta<=20) 103 rn_bb = 0.8 ! stretching with SH94 s-sigma 110 rn_bb = 0.8 ! stretching with SH94 s-sigma 104 111 !!!!!!! SF12 stretching coefficient (ln_s_sf12 = .true.) 105 112 rn_alpha = 4.4 ! stretching with SF12 s-sigma … … 110 117 rn_zb_b = -0.2 ! offset for calculating Zb 111 118 !!!!!!!! Other stretching (not SH94 or SF12) [also uses rn_theta above] 112 rn_thetb = 1.0 ! bottom control parameter (0<=thetb<= 1) 119 rn_thetb = 1.0 ! bottom control parameter (0<=thetb<= 1) 113 120 / 114 121 !----------------------------------------------------------------------- … … 118 125 rn_bathy = 0. ! value of the bathymetry. if (=0) bottom flat at jpkm1 119 126 nn_closea = 0 ! remove (=0) or keep (=1) closed seas and lakes (ORCA) 120 nn_msh = 0! create (=1) a mesh file or not (=0)127 nn_msh = 1 ! create (=1) a mesh file or not (=0) 121 128 rn_hmin = -3. ! min depth of the ocean (>0) or min number of ocean level (<0) 122 129 rn_e3zps_min= 20. ! partial step thickness is set larger than the minimum of … … 156 163 / 157 164 !----------------------------------------------------------------------- 158 &namsplit ! time splitting parameters ("key_dynspg_ts")159 !-----------------------------------------------------------------------160 ln_bt_fw = .TRUE. ! Forward integration of barotropic equations161 ln_bt_av = .TRUE. ! Time filtering of barotropic variables162 ln_bt_nn_auto = .TRUE. ! Set nn_baro automatically to be just below163 ! a user defined maximum courant number (rn_bt_cmax)164 nn_baro = 30 ! Number of iterations of barotropic mode165 ! during rn_rdt seconds. Only used if ln_bt_nn_auto=F166 rn_bt_cmax = 0.8 ! Maximum courant number allowed if ln_bt_nn_auto=T167 nn_bt_flt = 1 ! Time filter choice168 ! = 0 None169 ! = 1 Boxcar over nn_baro barotropic steps170 ! = 2 Boxcar over 2*nn_baro " "171 /172 !-----------------------------------------------------------------------173 165 &namcrs ! Grid coarsening for dynamics output and/or 174 166 ! passive tracer coarsened online simulations … … 196 188 &namtsd ! data : Temperature & Salinity 197 189 !----------------------------------------------------------------------- 198 !-----------------------------------------------------------------------199 190 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 200 191 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! … … 215 206 !! namsbc_core CORE bulk formulae formulation 216 207 !! namsbc_mfs MFS bulk formulae formulation 217 !! namsbc_cpl CouPLed formulation ("key_ coupled")208 !! namsbc_cpl CouPLed formulation ("key_oasis3") 218 209 !! namsbc_sas StAndalone Surface module 219 210 !! namtra_qsr penetrative solar radiation … … 235 226 ln_blk_core = .true. ! CORE bulk formulation (T => fill namsbc_core) 236 227 ln_blk_mfs = .false. ! MFS bulk formulation (T => fill namsbc_mfs ) 228 ln_cpl = .false. ! atmosphere coupled formulation ( requires key_oasis3 ) 229 ln_mixcpl = .false. ! forced-coupled mixed formulation ( requires key_oasis3 ) 230 nn_components = 0 ! configuration of the opa-sas OASIS coupling 231 ! =0 no opa-sas OASIS coupling: default single executable configuration 232 ! =1 opa-sas OASIS coupling: multi executable configuration, OPA component 233 ! =2 opa-sas OASIS coupling: multi executable configuration, SAS component 237 234 ln_apr_dyn = .false. ! Patm gradient added in ocean & ice Eqs. (T => fill namsbc_apr ) 238 235 nn_ice = 2 ! =0 no ice boundary condition , … … 245 242 ln_rnf = .true. ! runoffs (T => fill namsbc_rnf) 246 243 nn_isf = 0 ! ice shelf melting/freezing (/=0 => fill namsbc_isf) 247 ! 0 =no isf 1 = presence of ISF 248 ! 2 = bg03 parametrisation 3 = rnf file for isf 244 ! 0 =no isf 1 = presence of ISF 245 ! 2 = bg03 parametrisation 3 = rnf file for isf 249 246 ! 4 = ISF fwf specified 250 247 ! option 1 and 4 need ln_isfcav = .true. (domzgr) … … 277 274 &namsbc_flx ! surface boundary condition : flux formulation 278 275 !----------------------------------------------------------------------- 279 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 276 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 280 277 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 281 278 sn_utau = 'utau' , 24 , 'utau' , .false. , .false., 'yearly' , '' , '' , '' … … 320 317 ln_taudif = .false. ! HF tau contribution: use "mean of stress module - module of the mean stress" data 321 318 rn_zqt = 10. ! Air temperature and humidity reference height (m) 322 rn_zu = 10. ! Wind vector reference height (m) 319 rn_zu = 10. ! Wind vector reference height (m) 323 320 rn_pfac = 1. ! multiplicative factor for precipitation (total & snow) 324 321 rn_efac = 1. ! multiplicative factor for evaporation (0. or 1.) 325 rn_vfac = 0. ! multiplicative factor for ocean/ice velocity 322 rn_vfac = 0. ! multiplicative factor for ocean/ice velocity 326 323 ! in the calculation of the wind stress (0.=absolute winds or 1.=relative winds) 327 324 / … … 342 339 / 343 340 !----------------------------------------------------------------------- 344 &namsbc_cpl ! coupled ocean/atmosphere model ("key_ coupled")341 &namsbc_cpl ! coupled ocean/atmosphere model ("key_oasis3") 345 342 !----------------------------------------------------------------------- 346 343 ! ! description ! multiple ! vector ! vector ! vector ! … … 373 370 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 374 371 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 375 sn_usp = 'sas_grid_U' , 120 , 'vozocrtx' , .true. , .true. , 'yearly' , '' , '' , '' 372 sn_usp = 'sas_grid_U' , 120 , 'vozocrtx' , .true. , .true. , 'yearly' , '' , '' , '' 376 373 sn_vsp = 'sas_grid_V' , 120 , 'vomecrty' , .true. , .true. , 'yearly' , '' , '' , '' 377 374 sn_tem = 'sas_grid_T' , 120 , 'sosstsst' , .true. , .true. , 'yearly' , '' , '' , '' 378 375 sn_sal = 'sas_grid_T' , 120 , 'sosaline' , .true. , .true. , 'yearly' , '' , '' , '' 379 376 sn_ssh = 'sas_grid_T' , 120 , 'sossheig' , .true. , .true. , 'yearly' , '' , '' , '' 380 381 ln_3d_uv = .true. ! specify whether we are supplying a 3D u,v field 377 sn_e3t = 'sas_grid_T' , 120 , 'e3t_m' , .true. , .true. , 'yearly' , '' , '' , '' 378 sn_frq = 'sas_grid_T' , 120 , 'frq_m' , .true. , .true. , 'yearly' , '' , '' , '' 379 380 ln_3d_uve = .true. ! specify whether we are supplying a 3D u,v and e3 field 381 ln_read_frq = .false. ! specify whether we must read frq or not 382 382 cn_dir = './' ! root directory for the location of the bulk files are 383 383 / … … 412 412 413 413 cn_dir = './' ! root directory for the location of the runoff files 414 ln_rnf_emp = .false. ! runoffs included into precipitation field (T) or into a file (F)415 414 ln_rnf_mouth = .true. ! specific treatment at rivers mouths 416 415 rn_hrnf = 15.e0 ! depth over which enhanced vertical mixing is used … … 420 419 ln_rnf_tem = .false. ! read in temperature information for runoff 421 420 ln_rnf_sal = .false. ! read in salinity information for runoff 422 / 423 !----------------------------------------------------------------------- 424 &namsbc_isf ! Top boundary layer (ISF) 421 ln_rnf_depth_ini = .false. ! compute depth at initialisation from runoff file 422 rn_rnf_max = 5.735e-4 ! max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true ) 423 rn_dep_max = 150. ! depth over which runoffs is spread ( ln_rnf_depth_ini = .true ) 424 nn_rnf_depth_file = 0 ! create (=1) a runoff depth file or not (=0) 425 / 426 !----------------------------------------------------------------------- 427 &namsbc_isf ! Top boundary layer (ISF) 425 428 !----------------------------------------------------------------------- 426 429 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! … … 459 462 460 463 cn_dir = './' ! root directory for the location of the bulk files 461 rn_pref = 101000. _wp! reference atmospheric pressure [N/m2]/464 rn_pref = 101000. ! reference atmospheric pressure [N/m2]/ 462 465 ln_ref_apr = .false. ! ref. pressure: global mean Patm (T) or a constant (F) 463 466 ln_apr_obc = .false. ! inverse barometer added to OBC ssh data … … 499 502 ! Initial mass required for an iceberg of each class 500 503 rn_initial_mass = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11 501 ! Proportion of calving mass to apportion to each class 504 ! Proportion of calving mass to apportion to each class 502 505 rn_distribution = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 503 506 ! Ratio between effective and real iceberg mass (non-dim) 504 ! i.e. number of icebergs represented at a point 507 ! i.e. number of icebergs represented at a point 505 508 rn_mass_scaling = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1 506 509 ! thickness of newly calved bergs (m) … … 511 514 rn_bits_erosion_fraction = 0. ! Fraction of erosion melt flux to divert to bergy bits 512 515 rn_sicn_shift = 0. ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1) 513 ln_passive_mode = .false. ! iceberg - ocean decoupling 516 ln_passive_mode = .false. ! iceberg - ocean decoupling 514 517 nn_test_icebergs = 10 ! Create test icebergs of this class (-1 = no) 515 518 ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2) 516 519 rn_test_box = 108.0, 116.0, -66.0, -58.0 517 rn_speed_limit = 0. ! CFL speed limit for a berg 520 rn_speed_limit = 0. ! CFL speed limit for a berg 518 521 519 522 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 520 523 ! ! ! (if <0 months) ! name ! (logical) ! (T/F ) ! 'monthly' ! filename ! pairing ! filename ! 521 524 sn_icb = 'calving' , -1 , 'calvingmask', .true. , .true. , 'yearly' , '' , '' , '' 522 523 cn_dir = './' 525 526 cn_dir = './' 524 527 / 525 528 … … 528 531 !!====================================================================== 529 532 !! namlbc lateral momentum boundary condition 530 !! namcla cross land advection531 533 !! namobc open boundaries parameters ("key_obc") 532 534 !! namagrif agrif nested grid ( read by child model only ) ("key_agrif") … … 543 545 / 544 546 !----------------------------------------------------------------------- 545 &namcla ! cross land advection546 !-----------------------------------------------------------------------547 nn_cla = 0 ! advection between 2 ocean pts separates by land548 /549 !-----------------------------------------------------------------------550 &namobc ! open boundaries parameters ("key_obc")551 !-----------------------------------------------------------------------552 ln_obc_clim = .false. ! climatological obc data files (T) or not (F)553 ln_vol_cst = .true. ! impose the total volume conservation (T) or not (F)554 ln_obc_fla = .false. ! Flather open boundary condition555 nn_obcdta = 1 ! = 0 the obc data are equal to the initial state556 ! = 1 the obc data are read in 'obc.dta' files557 cn_obcdta = 'annual' ! set to annual if obc datafile hold 1 year of data558 ! set to monthly if obc datafile hold 1 month of data559 rn_dpein = 1. ! damping time scale for inflow at east open boundary560 rn_dpwin = 1. ! - - - west - -561 rn_dpnin = 1. ! - - - north - -562 rn_dpsin = 1. ! - - - south - -563 rn_dpeob = 3000. ! time relaxation (days) for the east open boundary564 rn_dpwob = 15. ! - - - west - -565 rn_dpnob = 3000. ! - - - north - -566 rn_dpsob = 15. ! - - - south - -567 rn_volemp = 1. ! = 0 the total volume change with the surface flux (E-P-R)568 ! = 1 the total volume remains constant569 /570 !-----------------------------------------------------------------------571 547 &namagrif ! AGRIF zoom ("key_agrif") 572 548 !----------------------------------------------------------------------- … … 575 551 rn_sponge_tra = 2880. ! coefficient for tracer sponge layer [m2/s] 576 552 rn_sponge_dyn = 2880. ! coefficient for dynamics sponge layer [m2/s] 553 ln_chk_bathy = .FALSE. ! 577 554 / 578 555 !----------------------------------------------------------------------- … … 597 574 ! = 2, use tidal harmonic forcing data from files 598 575 ! = 3, use external data AND tidal harmonic forcing 599 cn_dyn3d = 'none' ! 576 cn_dyn3d = 'none' ! 600 577 nn_dyn3d_dta = 0 ! = 0, bdy data are equal to the initial state 601 578 ! = 1, bdy data are read in 'bdydata .nc' files 602 cn_tra = 'none' ! 579 cn_tra = 'none' ! 603 580 nn_tra_dta = 0 ! = 0, bdy data are equal to the initial state 604 581 ! = 1, bdy data are read in 'bdydata .nc' files 605 cn_ice_lim = 'none' ! 582 cn_ice_lim = 'none' ! 606 583 nn_ice_lim_dta = 0 ! = 0, bdy data are equal to the initial state 607 584 ! = 1, bdy data are read in 'bdydata .nc' files … … 612 589 ln_tra_dmp =.false. ! open boudaries conditions for tracers 613 590 ln_dyn3d_dmp =.false. ! open boundary condition for baroclinic velocities 614 rn_time_dmp = 1. ! Damping time scale in days 591 rn_time_dmp = 1. ! Damping time scale in days 615 592 rn_time_dmp_out = 1. ! Outflow damping time scale 616 593 nn_rimwidth = 10 ! width of the relaxation zone … … 665 642 rn_bfri2_max = 1.e-1 ! max. bottom drag coefficient (non linear case and ln_loglayer=T) 666 643 rn_bfeb2 = 2.5e-3 ! bottom turbulent kinetic energy background (m2/s2) 667 rn_bfrz0 = 3.e-3 ! bottom roughness [m] if ln_loglayer=T 644 rn_bfrz0 = 3.e-3 ! bottom roughness [m] if ln_loglayer=T 668 645 ln_bfr2d = .false. ! horizontal variation of the bottom friction coef (read a 2D mask file ) 669 646 rn_bfrien = 50. ! local multiplying factor of bfr (ln_bfr2d=T) … … 682 659 &nambbc ! bottom temperature boundary condition 683 660 !----------------------------------------------------------------------- 661 ! ! ! (if <0 months) ! 662 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 663 ! ! ! (if <0 months) ! name ! (logical) ! (T/F ) ! 'monthly' ! filename ! pairing ! filename ! 664 sn_qgh ='geothermal_heating.nc', -12. , 'heatflow' , .false. , .true. , 'yearly' , '' , '' , '' 665 ! 666 cn_dir = './' ! root directory for the location of the runoff files 684 667 ln_trabbc = .true. ! Apply a geothermal heating at the ocean bottom 685 668 nn_geoflx = 2 ! geothermal heat flux: = 0 no flux … … 700 683 !! Tracer (T & S ) namelists 701 684 !!====================================================================== 702 !! nameos equation of state703 !! namtra_adv advection scheme685 !! nameos equation of state 686 !! namtra_adv advection scheme 704 687 !! namtra_adv_mle mixed layer eddy param. (Fox-Kemper param.) 705 !! namtra_ldf lateral diffusion scheme 706 !! 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 707 691 !!====================================================================== 708 692 ! … … 711 695 !----------------------------------------------------------------------- 712 696 nn_eos = -1 ! type of equation of state and Brunt-Vaisala frequency 713 ! =-1, TEOS-10 714 ! = 0, EOS-80 697 ! =-1, TEOS-10 698 ! = 0, EOS-80 715 699 ! = 1, S-EOS (simplified eos) 716 700 ln_useCT = .true. ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 717 !!701 ! 718 702 ! ! S-EOS coefficients : 719 !! 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 720 704 rn_a0 = 1.6550e-1 ! thermal expension coefficient (nn_eos= 1) 721 705 rn_b0 = 7.6554e-1 ! saline expension coefficient (nn_eos= 1) … … 729 713 &namtra_adv ! advection scheme for tracer 730 714 !----------------------------------------------------------------------- 731 ln_traadv_cen2 = .false. ! 2nd order centered scheme 732 ln_traadv_tvd = .true. ! TVD scheme 733 ln_traadv_muscl = .false. ! MUSCL scheme 734 ln_traadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 735 ln_traadv_ubs = .false. ! UBS scheme 736 ln_traadv_qck = .false. ! QUICKEST scheme 737 ln_traadv_msc_ups= .false. ! use upstream scheme within muscl 738 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 739 728 / 740 729 !----------------------------------------------------------------------- 741 730 &namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param) 742 731 !----------------------------------------------------------------------- 743 ln_mle = . true. ! (T) use the Mixed Layer Eddy (MLE) parameterisation732 ln_mle = .false. ! (T) use the Mixed Layer Eddy (MLE) parameterisation 744 733 rn_ce = 0.06 ! magnitude of the MLE (typical value: 0.06 to 0.08) 745 734 nn_mle = 1 ! MLE type: =0 standard Fox-Kemper ; =1 new formulation … … 755 744 !---------------------------------------------------------------------------------- 756 745 ! ! Operator type: 757 ln_traldf_lap = .true. ! laplacian operator 758 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 759 749 ! ! Direction of action: 760 ln_traldf_level = .false. ! iso-level 761 ln_traldf_hor = .false. ! horizontal (geopotential) (needs "key_ldfslp" when ln_sco=T) 762 ln_traldf_iso = .true. ! iso-neutral (needs "key_ldfslp") 763 ! ! Griffies parameters (all need "key_ldfslp") 764 ln_traldf_grif = .false. ! use griffies triads 765 ln_traldf_gdia = .false. ! output griffies eddy velocities 766 ln_triad_iso = .false. ! pure lateral mixing in ML 767 ln_botmix_grif = .false. ! lateral mixing on bottom 768 ! ! Coefficients 769 ! Eddy-induced (GM) advection always used with Griffies; otherwise needs "key_traldf_eiv" 770 ! Value rn_aeiv_0 is ignored unless = 0 with Held-Larichev spatially varying aeiv 771 ! (key_traldf_c2d & key_traldf_eiv & key_orca_r2, _r1 or _r05) 772 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] 773 rn_aht_0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] 774 rn_ahtb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 775 ! (normally=0; not used with Griffies) 776 rn_slpmax = 0.01 ! slope limit 777 rn_chsmag = 1. ! multiplicative factor in Smagorinsky diffusivity 778 rn_smsh = 1. ! Smagorinsky diffusivity: = 0 - use only sheer 779 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 780 787 / 781 788 !----------------------------------------------------------------------- … … 795 802 !! namdyn_vor advection scheme 796 803 !! namdyn_hpg hydrostatic pressure gradient 797 !! namdyn_spg surface pressure gradient (CPP key only)804 !! namdyn_spg surface pressure gradient 798 805 !! namdyn_ldf lateral diffusion scheme 799 806 !!====================================================================== … … 803 810 !----------------------------------------------------------------------- 804 811 ln_dynadv_vec = .true. ! vector form (T) or flux form (F) 812 nn_dynkeg = 0 ! scheme for grad(KE): =0 C2 ; =1 Hollingsworth correction 805 813 ln_dynadv_cen2= .false. ! flux form - 2nd order centered scheme 806 814 ln_dynadv_ubs = .false. ! flux form - 3rd order UBS scheme … … 810 818 &nam_vvl ! vertical coordinate options 811 819 !----------------------------------------------------------------------- 812 ln_vvl_zstar = .true. ! zstar vertical coordinate 820 ln_vvl_zstar = .true. ! zstar vertical coordinate 813 821 ln_vvl_ztilde = .false. ! ztilde vertical coordinate: only high frequency variations 814 822 ln_vvl_layer = .false. ! full layer vertical coordinate … … 827 835 ln_dynvor_ens = .false. ! energy conserving scheme 828 836 ln_dynvor_mix = .false. ! mixed scheme 829 ln_dynvor_een = .true. ! energy & enstrophy scheme 830 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 831 840 / 832 841 !----------------------------------------------------------------------- … … 839 848 ln_hpg_djc = .false. ! s-coordinate (Density Jacobian with Cubic polynomial) 840 849 ln_hpg_prj = .false. ! s-coordinate (Pressure Jacobian scheme) 841 ln_dynhpg_imp = .false. ! time stepping: semi-implicit time scheme (T) 842 ! centered time scheme (F) 843 / 844 !----------------------------------------------------------------------- 845 !namdyn_spg ! surface pressure gradient (CPP key only) 846 !----------------------------------------------------------------------- 847 ! ! explicit free surface ("key_dynspg_exp") 848 ! ! filtered free surface ("key_dynspg_flt") 849 ! ! split-explicit free surface ("key_dynspg_ts") 850 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 / 851 865 !----------------------------------------------------------------------- 852 866 &namdyn_ldf ! lateral diffusion on momentum 853 867 !----------------------------------------------------------------------- 854 868 ! ! Type of the operator : 855 ln_dynldf_lap = .true. ! laplacian operator 856 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 857 872 ! ! Direction of action : 858 ln_dynldf_lev el = .false.! iso-level859 ln_dynldf_hor = .true. ! horizontal (geopotential) (require "key_ldfslp" in s-coord.)860 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 861 876 ! ! Coefficient 862 rn_ahm_0_lap = 40000. ! horizontal laplacian eddy viscosity [m2/s] 863 rn_ahmb_0 = 0. ! background eddy viscosity for ldf_iso [m2/s] 864 rn_ahm_0_blp = 0. ! horizontal bilaplacian eddy viscosity [m4/s] 865 rn_cmsmag_1 = 3. ! constant in laplacian Smagorinsky viscosity 866 rn_cmsmag_2 = 3 ! constant in bilaplacian Smagorinsky viscosity 867 rn_cmsh = 1. ! 1 or 0 , if 0 -use only shear for Smagorinsky viscosity 868 rn_ahm_m_blp = -1.e12 ! upper limit for bilap abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp) 869 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) 870 890 / 871 891 … … 876 896 !! namzdf_ric richardson number dependent vertical mixing ("key_zdfric") 877 897 !! namzdf_tke TKE dependent vertical mixing ("key_zdftke") 878 !! namzdf_kpp KPP dependent vertical mixing ("key_zdfkpp")879 898 !! namzdf_ddm double diffusive mixing parameterization ("key_zdfddm") 880 899 !! namzdf_tmx tidal mixing parameterization ("key_zdftmx") … … 932 951 ! = 1 add a tke source below the ML 933 952 ! = 2 add a tke source just at the base of the ML 934 ! = 3 as = 1 applied on HF part of the stress ("key_ coupled")953 ! = 3 as = 1 applied on HF part of the stress ("key_oasis3") 935 954 rn_efr = 0.05 ! fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2) 936 955 nn_htau = 1 ! type of exponential decrease of tke penetration below the ML 937 956 ! = 0 constant 10 m length scale 938 957 ! = 1 0.5m at the equator to 30m poleward of 40 degrees 939 /940 !------------------------------------------------------------------------941 &namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionally:942 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb")943 ln_kpprimix = .true. ! shear instability mixing944 rn_difmiw = 1.0e-04 ! constant internal wave viscosity [m2/s]945 rn_difsiw = 0.1e-04 ! constant internal wave diffusivity [m2/s]946 rn_riinfty = 0.8 ! local Richardson Number limit for shear instability947 rn_difri = 0.0050 ! maximum shear mixing at Rig = 0 [m2/s]948 rn_bvsqcon = -0.01e-07 ! Brunt-Vaisala squared for maximum convection [1/s2]949 rn_difcon = 1. ! maximum mixing in interior convection [m2/s]950 nn_avb = 0 ! horizontal averaged (=1) or not (=0) on avt and amv951 nn_ave = 1 ! constant (=0) or profile (=1) background on avt952 958 / 953 959 !----------------------------------------------------------------------- … … 989 995 !! *** Miscellaneous namelists *** 990 996 !!====================================================================== 991 !! namsol elliptic solver / island / free surface992 997 !! nammpp Massively Parallel Processing ("key_mpp_mpi) 993 998 !! namctl Control prints & Benchmark … … 995 1000 !! namc1d_uvd data: U & V currents ("key_c1d") 996 1001 !! namc1d_dyndmp U & V newtonian damping ("key_c1d") 997 !!====================================================================== 998 ! 999 !----------------------------------------------------------------------- 1000 &namsol ! elliptic solver / island / free surface 1001 !----------------------------------------------------------------------- 1002 nn_solv = 1 ! elliptic solver: =1 preconditioned conjugate gradient (pcg) 1003 ! =2 successive-over-relaxation (sor) 1004 nn_sol_arp = 0 ! absolute/relative (0/1) precision convergence test 1005 rn_eps = 1.e-6 ! absolute precision of the solver 1006 nn_nmin = 300 ! minimum of iterations for the SOR solver 1007 nn_nmax = 800 ! maximum of iterations for the SOR solver 1008 nn_nmod = 10 ! frequency of test for the SOR solver 1009 rn_resmax = 1.e-10 ! absolute precision for the SOR solver 1010 rn_sor = 1.92 ! optimal coefficient for SOR solver (to be adjusted with the domain) 1011 / 1002 !! namsto Stochastic parametrization of EOS 1003 !!====================================================================== 1004 ! 1012 1005 !----------------------------------------------------------------------- 1013 1006 &nammpp ! Massively Parallel Processing ("key_mpp_mpi) … … 1053 1046 ln_dyndmp = .false. ! add a damping term (T) or not (F) 1054 1047 / 1048 !----------------------------------------------------------------------- 1049 &namsto ! Stochastic parametrization of EOS 1050 !----------------------------------------------------------------------- 1051 ln_rststo = .false. ! start from mean parameter (F) or from restart file (T) 1052 ln_rstseed = .true. ! read seed of RNG from restart file 1053 cn_storst_in = "restart_sto" ! suffix of stochastic parameter restart file (input) 1054 cn_storst_out = "restart_sto" ! suffix of stochastic parameter restart file (output) 1055 1056 ln_sto_eos = .false. ! stochastic equation of state 1057 nn_sto_eos = 1 ! number of independent random walks 1058 rn_eos_stdxy = 1.4 ! random walk horz. standard deviation (in grid points) 1059 rn_eos_stdz = 0.7 ! random walk vert. standard deviation (in grid points) 1060 rn_eos_tcor = 1440.0 ! random walk time correlation (in timesteps) 1061 nn_eos_ord = 1 ! order of autoregressive processes 1062 nn_eos_flt = 0 ! passes of Laplacian filter 1063 rn_eos_lim = 2.0 ! limitation factor (default = 3.0) 1064 / 1055 1065 1056 1066 !!====================================================================== … … 1153 1163 ln_s3d = .false. ! Logical switch for S profile observations 1154 1164 ln_ena = .false. ! Logical switch for ENACT insitu data set 1155 ! ! ln_corLogical switch for Coriolis insitu data set1165 ln_cor = .false. ! Logical switch for Coriolis insitu data set 1156 1166 ln_profb = .false. ! Logical switch for feedback insitu data set 1157 1167 ln_sla = .false. ! Logical switch for SLA observations 1158 1159 1168 ln_sladt = .false. ! Logical switch for AVISO SLA data 1160 1161 1169 ln_slafb = .false. ! Logical switch for feedback SLA data 1162 ! ln_ssh Logical switch for SSH observations 1163 1164 ln_sst = .false. ! Logical switch for SST observations 1165 ln_reysst = .false. ! ln_reysst Logical switch for Reynolds observations 1166 ln_ghrsst = .false. ! ln_ghrsst Logical switch for GHRSST observations 1167 1170 ln_ssh = .false. ! Logical switch for SSH observations 1171 ln_sst = .false. ! Logical switch for SST observations 1172 ln_reysst = .false. ! Logical switch for Reynolds observations 1173 ln_ghrsst = .false. ! Logical switch for GHRSST observations 1168 1174 ln_sstfb = .false. ! Logical switch for feedback SST data 1169 ! ln_sssLogical switch for SSS observations1175 ln_sss = .false. ! Logical switch for SSS observations 1170 1176 ln_seaice = .false. ! Logical switch for Sea Ice observations 1171 ! ln_vel3d Logical switch for velocity observations 1172 ! ln_velavcur Logical switch for velocity daily av. cur. 1173 ! ln_velhrcur Logical switch for velocity high freq. cur. 1174 ! ln_velavadcp Logical switch for velocity daily av. ADCP 1175 ! ln_velhradcp Logical switch for velocity high freq. ADCP 1176 ! ln_velfb Logical switch for feedback velocity data 1177 ! ln_grid_global Global distribtion of observations 1178 ! ln_grid_search_lookup Logical switch for obs grid search w/lookup table 1179 ! grid_search_file Grid search lookup file header 1180 ! enactfiles ENACT input observation file names 1181 ! coriofiles Coriolis input observation file name 1182 ! ! profbfiles: Profile feedback input observation file name 1183 profbfiles = 'profiles_01.nc' 1184 ! ln_profb_enatim Enact feedback input time setting switch 1185 ! slafilesact Active SLA input observation file name 1186 ! slafilespas Passive SLA input observation file name 1187 ! ! slafbfiles: Feedback SLA input observation file name 1188 slafbfiles = 'sla_01.nc' 1189 ! sstfiles GHRSST input observation file name 1190 ! ! sstfbfiles: Feedback SST input observation file name 1191 sstfbfiles = 'sst_01.nc' 1192 ! seaicefiles Sea Ice input observation file names 1193 seaicefiles = 'seaice_01.nc' 1194 ! velavcurfiles Vel. cur. daily av. input file name 1195 ! velhvcurfiles Vel. cur. high freq. input file name 1196 ! velavadcpfiles Vel. ADCP daily av. input file name 1197 ! velhvadcpfiles Vel. ADCP high freq. input file name 1198 ! velfbfiles Vel. feedback input observation file name 1199 ! dobsini Initial date in window YYYYMMDD.HHMMSS 1200 ! dobsend Final date in window YYYYMMDD.HHMMSS 1201 ! n1dint Type of vertical interpolation method 1202 ! n2dint Type of horizontal interpolation method 1203 ! ln_nea Rejection of observations near land switch 1204 nmsshc = 0 ! MSSH correction scheme 1205 ! mdtcorr MDT correction 1206 ! mdtcutoff MDT cutoff for computed correction 1177 ln_vel3d = .false. ! Logical switch for velocity observations 1178 ln_velavcur= .false ! Logical switch for velocity daily av. cur. 1179 ln_velhrcur= .false ! Logical switch for velocity high freq. cur. 1180 ln_velavadcp = .false. ! Logical switch for velocity daily av. ADCP 1181 ln_velhradcp = .false. ! Logical switch for velocity high freq. ADCP 1182 ln_velfb = .false. ! Logical switch for feedback velocity data 1183 ln_grid_global = .false. ! Global distribtion of observations 1184 ln_grid_search_lookup = .false. ! Logical switch for obs grid search w/lookup table 1185 grid_search_file = 'grid_search' ! Grid search lookup file header 1186 ! All of the *files* variables below are arrays. Use namelist_cfg to add more files 1187 enactfiles = 'enact.nc' ! ENACT input observation file names (specify full array in namelist_cfg) 1188 coriofiles = 'corio.nc' ! Coriolis input observation file name 1189 profbfiles = 'profiles_01.nc' ! Profile feedback input observation file name 1190 ln_profb_enatim = .false ! Enact feedback input time setting switch 1191 slafilesact = 'sla_act.nc' ! Active SLA input observation file names 1192 slafilespas = 'sla_pass.nc' ! Passive SLA input observation file names 1193 slafbfiles = 'sla_01.nc' ! slafbfiles: Feedback SLA input observation file names 1194 sstfiles = 'ghrsst.nc' ! GHRSST input observation file names 1195 sstfbfiles = 'sst_01.nc' ! Feedback SST input observation file names 1196 seaicefiles = 'seaice_01.nc' ! Sea Ice input observation file names 1197 velavcurfiles = 'velavcurfile.nc' ! Vel. cur. daily av. input file name 1198 velhrcurfiles = 'velhrcurfile.nc' ! Vel. cur. high freq. input file name 1199 velavadcpfiles = 'velavadcpfile.nc' ! Vel. ADCP daily av. input file name 1200 velhradcpfiles = 'velhradcpfile.nc' ! Vel. ADCP high freq. input file name 1201 velfbfiles = 'velfbfile.nc' ! Vel. feedback input observation file name 1202 dobsini = 20000101.000000 ! Initial date in window YYYYMMDD.HHMMSS 1203 dobsend = 20010101.000000 ! Final date in window YYYYMMDD.HHMMSS 1204 n1dint = 0 ! Type of vertical interpolation method 1205 n2dint = 0 ! Type of horizontal interpolation method 1206 ln_nea = .false. ! Rejection of observations near land switch 1207 nmsshc = 0 ! MSSH correction scheme 1208 mdtcorr = 1.61 ! MDT correction 1209 mdtcutoff = 65.0 ! MDT cutoff for computed correction 1207 1210 ln_altbias = .false. ! Logical switch for alt bias 1208 1211 ln_ignmis = .true. ! Logical switch for ignoring missing files 1209 ! endailyavtypes ENACT daily average types1212 endailyavtypes = 820 ! ENACT daily average types - array (use namelist_cfg to set more values) 1210 1213 ln_grid_global = .true. 1211 1214 ln_grid_search_lookup = .false.
Note: See TracChangeset
for help on using the changeset viewer.