- Timestamp:
- 2019-12-10T15:44:23+01:00 (5 years ago)
- Location:
- NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/cfgs/SHARED
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/cfgs/SHARED/field_def_nemo-oce.xml
r12109 r12154 41 41 <field id="ahmt_3d" long_name=" 3D t-eddy viscosity coefficient" unit="m2/s or m4/s" grid_ref="grid_T_3D"/> 42 42 43 <field id="sst" long_name="sea surface temperature" standard_name="sea_surface_temperature" unit="degC" /> 43 <field id="sst" long_name="Bulk sea surface temperature" standard_name="bulk_sea_surface_temperature" unit="degC" /> 44 <field id="t_skin" long_name="Skin temperature aka SSST" standard_name="skin_temperature" unit="degC" /> 44 45 <field id="sst2" long_name="square of sea surface temperature" standard_name="square_of_sea_surface_temperature" unit="degC2" > sst * sst </field > 45 46 <field id="sstmax" long_name="max of sea surface temperature" field_ref="sst" operation="maximum" /> … … 300 301 301 302 <!-- *_oce variables available with ln_blk_clio or ln_blk_core --> 303 <field id="rho_air" long_name="Air density at 10m above sea surface" standard_name="rho_air_10m" unit="kg/m3" /> 304 <field id="dt_skin" long_name="SSST-SST temperature difference" standard_name="SSST-SST" unit="K" /> 302 305 <field id="qlw_oce" long_name="Longwave Downward Heat Flux over open ocean" standard_name="surface_net_downward_longwave_flux" unit="W/m2" /> 303 306 <field id="qsb_oce" long_name="Sensible Downward Heat Flux over open ocean" standard_name="surface_downward_sensible_heat_flux" unit="W/m2" /> 304 307 <field id="qla_oce" long_name="Latent Downward Heat Flux over open ocean" standard_name="surface_downward_latent_heat_flux" unit="W/m2" /> 308 <field id="evap_oce" long_name="Evaporation over open ocean" standard_name="evaporation" unit="kg/m2/s" /> 305 309 <field id="qt_oce" long_name="total flux at ocean surface" standard_name="surface_downward_heat_flux_in_sea_water" unit="W/m2" /> 306 310 <field id="qsr_oce" long_name="solar heat flux at ocean surface" standard_name="net_downward_shortwave_flux_at_sea_water_surface" unit="W/m2" /> … … 362 366 </field_group> 363 367 364 <!-- scalar variables --> 365 <field_group id="SBC_0D" grid_ref="grid_1point" > 368 369 </field_group> <!-- SBC --> 370 371 <!-- ABL --> 372 <field_group id="ABL" > <!-- time step automaticaly defined based on nn_fsbc --> 373 374 <!-- variables available with ABL on atmospheric T grid--> 375 <field_group id="grid_ABL3D" grid_ref="grid_TA_3D" > 376 <field id="u_abl" long_name="ABL i-horizontal velocity" standard_name="abl_x_velocity" unit="m/s" /> 377 <field id="v_abl" long_name="ABL j-horizontal velocity" standard_name="abl_y_velocity" unit="m/s" /> 378 <field id="t_abl" long_name="ABL potential temperature" standard_name="abl_theta" unit="K" /> 379 <field id="q_abl" long_name="ABL specific humidity" standard_name="abl_qspe" unit="kg/kg" /> 380 <!-- debug (to be removed) --> 381 <field id="u_dta" long_name="DTA i-horizontal velocity" standard_name="dta_x_velocity" unit="m/s" /> 382 <field id="v_dta" long_name="DTA j-horizontal velocity" standard_name="dta_y_velocity" unit="m/s" /> 383 <field id="t_dta" long_name="DTA potential temperature" standard_name="dta_theta" unit="K" /> 384 <field id="q_dta" long_name="DTA specific humidity" standard_name="dta_qspe" unit="kg/kg" /> 385 <field id="coeft" long_name="ABL nudging coefficient" standard_name="coeft" unit="" /> 386 <field id="tke_abl" long_name="ABL turbulent kinetic energy" standard_name="abl_tke" unit="m2/s2" /> 387 <field id="avm_abl" long_name="ABL turbulent viscosity" standard_name="abl_avm" unit="m2/s" /> 388 <field id="avt_abl" long_name="ABL turbulent diffusivity" standard_name="abl_avt" unit="m2/s" /> 389 <field id="mxl_abl" long_name="ABL mixing length" standard_name="abl_mxl" unit="m" /> 366 390 </field_group> 367 391 368 </field_group> <!-- SBC --> 369 392 <field_group id="grid_ABL2D" grid_ref="grid_TA_2D" > 393 <field id="pblh" long_name="ABL height" standard_name="abl_height" unit="m" /> 394 <field id="uz1_abl" long_name="ABL i-horizontal velocity" standard_name="abl_x_velocity" unit="m/s" /> 395 <field id="vz1_abl" long_name="ABL j-horizontal velocity" standard_name="abl_y_velocity" unit="m/s" /> 396 <field id="uvz1_abl" long_name="ABL wind speed module" standard_name="abl_wind_speed" unit="m/s" > sqrt( uz1_abl^2 + vz1_abl^2 ) </field> 397 <field id="tz1_abl" long_name="ABL potential temperature" standard_name="abl_theta" unit="K" /> 398 <field id="qz1_abl" long_name="ABL specific humidity" standard_name="abl_qspe" unit="kg/kg" /> 399 <field id="uz1_dta" long_name="DTA i-horizontal velocity" standard_name="dta_x_velocity" unit="m/s" /> 400 <field id="vz1_dta" long_name="DTA j-horizontal velocity" standard_name="dta_y_velocity" unit="m/s" /> 401 <field id="uvz1_dta" long_name="DTA wind speed module" standard_name="dta_wind_speed" unit="m/s" > sqrt( uz1_dta^2 + vz1_dta^2 ) </field> 402 <field id="tz1_dta" long_name="DTA potential temperature" standard_name="dta_theta" unit="K" /> 403 <field id="qz1_dta" long_name="DTA specific humidity" standard_name="dta_qspe" unit="kg/kg" /> 404 <!-- debug (to be removed) --> 405 <field id="uz1_geo" long_name="GEO i-horizontal velocity" standard_name="geo_x_velocity" unit="m/s" /> 406 <field id="vz1_geo" long_name="GEO j-horizontal velocity" standard_name="geo_y_velocity" unit="m/s" /> 407 <field id="uvz1_geo" long_name="GEO wind speed module" standard_name="geo_wind_speed" unit="m/s" > sqrt( uz1_geo^2 + vz1_geo^2 ) </field> 408 </field_group> 409 410 </field_group> <!-- ABL --> 411 412 370 413 <!-- U grid --> 371 414 … … 391 434 <field id="uocet" long_name="ocean transport along i-axis times temperature (CRS)" unit="degC*m/s" grid_ref="grid_U_3D" /> 392 435 <field id="uoces" long_name="ocean transport along i-axis times salinity (CRS)" unit="1e-3*m/s" grid_ref="grid_U_3D" /> 436 <field id="ssuww" long_name="ocean surface wind work along i-axis" standard_name="surface_x_wind_work" unit="N/m*s" > utau * ssu </field> 393 437 394 438 <!-- u-eddy diffusivity coefficients (available if ln_traldf_OFF=F) --> … … 448 492 <field id="vocet" long_name="ocean transport along j-axis times temperature (CRS)" unit="degC*m/s" grid_ref="grid_V_3D" /> 449 493 <field id="voces" long_name="ocean transport along j-axis times salinity (CRS)" unit="1e-3*m/s" grid_ref="grid_V_3D" /> 494 <field id="ssvww" long_name="ocean surface wind work along j-axis" standard_name="surface_y_wind_work" unit="N/m*s" > vtau * ssv </field> 450 495 451 496 <!-- v-eddy diffusivity coefficients (available if ln_traldf_OFF=F) --> … … 589 634 590 635 591 <!-- variables available with key_float-->636 <!-- variables available with ln_floats --> 592 637 593 638 <field_group id="floatvar" grid_ref="grid_T_nfloat" operation="instant" > -
NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/cfgs/SHARED/namelist_ref
r12113 r12154 5 5 !! namelists 2 - Surface boundary (namsbc, namsbc_flx, namsbc_blk, namsbc_cpl, 6 6 !! namsbc_sas, namtra_qsr, namsbc_rnf, 7 !! namsbc_isf, namsbc_iscpl, namsbc_apr, 7 !! namsbc_isf, namsbc_iscpl, namsbc_apr, 8 8 !! namsbc_ssr, namsbc_wave, namberg) 9 9 !! 3 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) … … 65 65 ln_clobber = .true. ! clobber (overwrite) an existing file 66 66 nn_chunksz = 0 ! chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) 67 ln_xios_read = . FALSE. ! use XIOS to read restart file (only for a single file restart)67 ln_xios_read = .false. ! use XIOS to read restart file (only for a single file restart) 68 68 nn_wxios = 0 ! use XIOS to write restart file 0 - no, 1 - single file output, 2 - multiple file output 69 69 / … … 88 88 cn_domcfg = "domain_cfg" ! domain configuration filename 89 89 ! 90 ln_closea = .false. ! T => keep closed seas (defined by closea_mask field) in the 90 ln_closea = .false. ! T => keep closed seas (defined by closea_mask field) in the 91 91 ! ! domain and apply special treatment of freshwater fluxes. 92 ! ! F => suppress closed seas (defined by closea_mask field) 92 ! ! F => suppress closed seas (defined by closea_mask field) 93 93 ! ! from the bathymetry at runtime. 94 94 ! ! If closea_mask field doesn't exist in the domain_cfg file … … 106 106 ln_tsd_init = .false. ! ocean initialisation 107 107 ln_tsd_dmp = .false. ! T-S restoring (see namtra_dmp) 108 108 109 109 cn_dir = './' ! root directory for the T-S data location 110 110 !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! … … 195 195 nn_fsbc = 2 ! frequency of SBC module call 196 196 ! ! (control sea-ice & iceberg model call) 197 ! Type of air-sea fluxes 197 ! Type of air-sea fluxes 198 198 ln_usr = .false. ! user defined formulation (T => check usrdef_sbc) 199 199 ln_flx = .false. ! flux formulation (T => fill namsbc_flx ) 200 200 ln_blk = .false. ! Bulk formulation (T => fill namsbc_blk ) 201 ln_abl = .false. ! ABL formulation (T => fill namsbc_abl ) 201 202 ! ! Type of coupling (Ocean/Ice/Atmosphere) : 202 203 ln_cpl = .false. ! atmosphere coupled formulation ( requires key_oasis3 ) … … 205 206 ! ! =0 no opa-sas OASIS coupling: default single executable config. 206 207 ! ! =1 opa-sas OASIS coupling: multi executable config., OPA component 207 ! ! =2 opa-sas OASIS coupling: multi executable config., SAS component 208 ! ! =2 opa-sas OASIS coupling: multi executable config., SAS component 208 209 ! Sea-ice : 209 nn_ice = 0 ! =0 no ice boundary condition 210 nn_ice = 0 ! =0 no ice boundary condition 210 211 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 211 212 ! ! =2 or 3 automatically for SI3 or CICE ("key_si3" or "key_cice") … … 213 214 ln_ice_embd = .false. ! =T embedded sea-ice (pressure + mass and salt exchanges) 214 215 ! ! =F levitating ice (no pressure, mass and salt exchanges) 215 ! Misc. options of sbc : 216 ! Misc. options of sbc : 216 217 ln_traqsr = .false. ! Light penetration in the ocean (T => fill namtra_qsr) 217 218 ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave … … 225 226 ln_wave = .false. ! Activate coupling with wave (T => fill namsbc_wave) 226 227 ln_cdgw = .false. ! Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave) 227 ln_sdw = .false. ! Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave) 228 ln_sdw = .false. ! Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave) 228 229 nn_sdrift = 0 ! Parameterization for the calculation of 3D-Stokes drift from the surface Stokes drift 229 230 ! ! = 0 Breivik 2015 parameterization: v_z=v_0*[exp(2*k*z)/(1-8*k*z)] … … 250 251 / 251 252 !----------------------------------------------------------------------- 252 &namsbc_blk ! namsbc_blk generic Bulk formula 253 &namsbc_blk ! namsbc_blk generic Bulk formula (ln_blk =T) 253 254 !----------------------------------------------------------------------- 254 255 ! ! bulk algorithm : 255 ln_NCAR = .false.! "NCAR" algorithm (Large and Yeager 2008)256 ln_NCAR = .true. ! "NCAR" algorithm (Large and Yeager 2008) 256 257 ln_COARE_3p0 = .false. ! "COARE 3.0" algorithm (Fairall et al. 2003) 257 ln_COARE_3p 5 = .false. ! "COARE 3.5" algorithm (Edson et al. 2013)258 ln_ECMWF = .false. ! "ECMWF" algorithm (IFS cycle 31)258 ln_COARE_3p6 = .false. ! "COARE 3.6" algorithm (Edson et al. 2013) 259 ln_ECMWF = .false. ! "ECMWF" algorithm (IFS cycle 45r1) 259 260 ! 260 rn_zqt = 10. ! Air temperature & humidity reference height (m) 261 rn_zu = 10. ! Wind vector reference height (m) 262 ln_Cd_L12 = .false. ! air-ice drags = F(ice concentration) (Lupkes et al. 2012) 263 ln_Cd_L15 = .false. ! air-ice drags = F(ice concentration) (Lupkes et al. 2015) 264 ln_taudif = .false. ! HF tau contribution: use "mean of stress module - module of the mean stress" data 265 rn_pfac = 1. ! multiplicative factor for precipitation (total & snow) 266 rn_efac = 1. ! multiplicative factor for evaporation (0. or 1.) 267 rn_vfac = 0. ! multiplicative factor for ocean & ice velocity used to 268 ! ! calculate the wind stress (0.=absolute or 1.=relative winds) 269 261 rn_zqt = 10. ! Air temperature & humidity reference height (m) 262 rn_zu = 10. ! Wind vector reference height (m) 263 ln_Cd_L12 = .false. ! air-ice drags = F(ice conc.) (Lupkes et al. 2012) 264 ln_Cd_L15 = .false. ! air-ice drags = F(ice conc.) (Lupkes et al. 2015) 265 ! ! - module of the mean stress" data 266 rn_pfac = 1. ! multipl. factor for precipitation (total & snow) 267 rn_efac = 1. ! multipl. factor for evaporation (0. or 1.) 268 rn_vfac = 0. ! multipl. factor for ocean & ice velocity 269 ! ! used to calculate the wind stress 270 ! ! (0. => absolute or 1. => relative winds) 271 ln_skin_cs = .false. ! use the cool-skin parameterization 272 ln_skin_wl = .false. ! use the warm-layer parameterization 273 ! ! ==> only available in ECMWF and COARE algorithms 274 ln_humi_sph = .true. ! humidity "sn_humi" is specific humidity [kg/kg] 275 ln_humi_dpt = .false. ! humidity "sn_humi" is dew-point temperature [K] 276 ln_humi_rlh = .false. ! humidity "sn_humi" is relative humidity [%] 277 ! 270 278 cn_dir = './' ! root directory for the bulk data location 271 279 !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________! … … 278 286 sn_tair = 't_10.15JUNE2009_fill' , 6. , 'T_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 279 287 sn_humi = 'q_10.15JUNE2009_fill' , 6. , 'Q_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 288 sn_hpgi = 'NONE' , 24. , 'uhpg' , .false. , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'UG' , '' 289 sn_hpgj = 'NONE' , 24. , 'vhpg' , .false. , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'VG' , '' 280 290 sn_prec = 'ncar_precip.15JUNE2009_fill', -1. , 'PRC_MOD1', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 281 291 sn_snow = 'ncar_precip.15JUNE2009_fill', -1. , 'SNOW' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 282 292 sn_slp = 'slp.15JUNE2009_fill' , 6. , 'SLP' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 283 sn_tdif = 'taudif_core' , 24 , 'taudif' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 293 / 294 !----------------------------------------------------------------------- 295 &namsbc_abl ! Atmospheric Boundary Layer formulation (ln_abl = T) 296 !----------------------------------------------------------------------- 297 cn_dir = './' ! root directory for the location of the ABL grid file 298 cn_dom = 'dom_cfg_abl.nc' 299 300 cn_ablrst_in = "restart_abl" ! suffix of abl restart name (input) 301 cn_ablrst_out = "restart_abl" ! suffix of abl restart name (output) 302 cn_ablrst_indir = "." ! directory to read input abl restarts 303 cn_ablrst_outdir = "." ! directory to write output abl restarts 304 305 ln_hpgls_frc = .false. 306 ln_geos_winds = .false. 307 nn_dyn_restore = 2 ! restoring option for dynamical ABL variables: = 0 no restoring 308 ! = 1 equatorial restoring 309 ! = 2 global restoring 310 rn_ldyn_min = 4.5 ! magnitude of the nudging on ABL dynamics at the bottom of the ABL [hour] 311 rn_ldyn_max = 1.5 ! magnitude of the nudging on ABL dynamics at the top of the ABL [hour] 312 rn_ltra_min = 4.5 ! magnitude of the nudging on ABL tracers at the bottom of the ABL [hour] 313 rn_ltra_max = 1.5 ! magnitude of the nudging on ABL tracers at the top of the ABL [hour] 314 nn_amxl = 0 ! mixing length: = 0 Deardorff 80 length-scale 315 ! = 1 length-scale based on the distance to the PBL height 316 ! = 2 Bougeault & Lacarrere 89 length-scale 317 rn_Cm = 0.0667 ! 0.126 in MesoNH 318 rn_Ct = 0.1667 ! 0.143 in MesoNH 319 rn_Ce = 0.4 ! 0.4 in MesoNH 320 rn_Ceps = 0.7 ! 0.85 in MesoNH 321 rn_Rod = 0.15 ! c0 in RMCA17 mixing length formulation (not yet implemented) 322 rn_Ric = 0.139 ! Critical Richardson number (to compute PBL height and diffusivities) 284 323 / 285 324 !----------------------------------------------------------------------- … … 375 414 nn_chldta = 0 ! RGB : Chl data (=1) or cst value (=0) 376 415 rn_si1 = 23.0 ! 2BD : longest depth of extinction 377 416 378 417 cn_dir = './' ! root directory for the chlorophyl data location 379 418 !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! … … 443 482 / 444 483 !----------------------------------------------------------------------- 445 &namsbc_isf ! Top boundary layer (ISF) (ln_isfcav =T : read (ln_read_cfg=T) 484 &namsbc_isf ! Top boundary layer (ISF) (ln_isfcav =T : read (ln_read_cfg=T) 446 485 !----------------------------------------------------------------------- or set or usr_def_zgr ) 447 ! ! type of top boundary layer 486 ! ! type of top boundary layer 448 487 nn_isf = 1 ! ice shelf melting/freezing 449 ! 1 = presence of ISF ; 2 = bg03 parametrisation 488 ! 1 = presence of ISF ; 2 = bg03 parametrisation 450 489 ! 3 = rnf file for ISF ; 4 = ISF specified freshwater flux 451 490 ! options 1 and 4 need ln_isfcav = .true. (domzgr) … … 470 509 !* nn_isf = 3 case 471 510 sn_rnfisf = 'rnfisf' , -12. ,'sofwfisf' , .false. , .true. , 'yearly' , '' , '' , '' 472 !* nn_isf = 2 and 3 cases 511 !* nn_isf = 2 and 3 cases 473 512 sn_depmax_isf ='rnfisf' , -12. ,'sozisfmax', .false. , .true. , 'yearly' , '' , '' , '' 474 513 sn_depmin_isf ='rnfisf' , -12. ,'sozisfmin', .false. , .true. , 'yearly' , '' , '' , '' … … 477 516 / 478 517 !----------------------------------------------------------------------- 479 &namsbc_iscpl ! land ice / ocean coupling option (ln_isfcav =T : read (ln_read_cfg=T) 518 &namsbc_iscpl ! land ice / ocean coupling option (ln_isfcav =T : read (ln_read_cfg=T) 480 519 !----------------------------------------------------------------------- or set or usr_def_zgr ) 481 520 nn_drown = 10 ! number of iteration of the extrapolation loop (fill the new wet cells) … … 572 611 !----------------------------------------------------------------------- 573 612 ln_tide = .false. ! Activate tides 574 ln_tide_pot = . true.! use tidal potential forcing613 ln_tide_pot = .false. ! use tidal potential forcing 575 614 ln_scal_load = .false. ! Use scalar approximation for 576 615 rn_scal_load = 0.094 ! load potential 577 616 ln_read_load = .false. ! Or read load potential from file 578 617 cn_tide_load = 'tide_LOAD_grid_T.nc' ! filename for load potential 579 ! 618 ! 580 619 ln_tide_ramp = .false. ! Use linear ramp for tides at startup 581 620 rdttideramp = 0. ! ramp duration in days … … 656 695 filtide = 'bdydta/amm12_bdytide_' ! file name root of tidal forcing files 657 696 ln_bdytide_2ddta = .false. ! 658 ln_bdytide_conj = .false. ! 697 ln_bdytide_conj = .false. ! 659 698 / 660 699 … … 683 722 !----------------------------------------------------------------------- 684 723 rn_Cd0 = 1.e-3 ! drag coefficient [-] 685 rn_Uc0 = 0.4 ! ref. velocity [m/s] (linear drag=Cd0*Uc0) 724 rn_Uc0 = 0.4 ! ref. velocity [m/s] (linear drag=Cd0*Uc0) 686 725 rn_Cdmax = 0.1 ! drag value maximum [-] (logarithmic drag) 687 726 rn_ke0 = 2.5e-3 ! background kinetic energy [m2/s2] (non-linear cases) … … 694 733 !----------------------------------------------------------------------- 695 734 rn_Cd0 = 1.e-3 ! drag coefficient [-] 696 rn_Uc0 = 0.4 ! ref. velocity [m/s] (linear drag=Cd0*Uc0) 735 rn_Uc0 = 0.4 ! ref. velocity [m/s] (linear drag=Cd0*Uc0) 697 736 rn_Cdmax = 0.1 ! drag value maximum [-] (logarithmic drag) 698 737 rn_ke0 = 2.5e-3 ! background kinetic energy [m2/s2] (non-linear cases) … … 761 800 nn_cen_v = 4 ! =2/4, vertical 2nd order CEN / 4th order COMPACT 762 801 ln_traadv_fct = .false. ! FCT scheme 763 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 764 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 802 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 803 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 765 804 ln_traadv_mus = .false. ! MUSCL scheme 766 805 ln_mus_ups = .false. ! use upstream scheme near river mouths … … 783 822 ln_traldf_triad = .false. ! iso-neutral (triad operator) 784 823 ! 785 ! ! iso-neutral options: 824 ! ! iso-neutral options: 786 825 ln_traldf_msc = .false. ! Method of Stabilizing Correction (both operators) 787 826 rn_slpmax = 0.01 ! slope limit (both operators) … … 793 832 nn_aht_ijk_t = 0 ! space/time variation of eddy coefficient: 794 833 ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file 795 ! ! = 0 constant 796 ! ! = 10 F(k) =ldf_c1d 797 ! ! = 20 F(i,j) =ldf_c2d 834 ! ! = 0 constant 835 ! ! = 10 F(k) =ldf_c1d 836 ! ! = 20 F(i,j) =ldf_c2d 798 837 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 799 838 ! ! = 30 F(i,j,k) =ldf_c2d * ldf_c1d 800 839 ! ! = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 801 ! ! time invariant coefficients: aht0 = 1/2 Ud*Ld (lap case) 840 ! ! time invariant coefficients: aht0 = 1/2 Ud*Ld (lap case) 802 841 ! ! or = 1/12 Ud*Ld^3 (blp case) 803 842 rn_Ud = 0.01 ! lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30) … … 825 864 nn_aei_ijk_t = 0 ! space/time variation of eddy coefficient: 826 865 ! ! =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file 827 ! ! = 0 constant 828 ! ! = 10 F(k) =ldf_c1d 829 ! ! = 20 F(i,j) =ldf_c2d 866 ! ! = 0 constant 867 ! ! = 10 F(k) =ldf_c1d 868 ! ! = 20 F(i,j) =ldf_c2d 830 869 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 831 870 ! ! = 30 F(i,j,k) =ldf_c2d * ldf_c1d 832 ! ! time invariant coefficients: aei0 = 1/2 Ue*Le 871 ! ! time invariant coefficients: aei0 = 1/2 Ue*Le 833 872 rn_Ue = 0.02 ! lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30) 834 873 rn_Le = 200.e+3 ! lateral diffusive length [m] (nn_aht_ijk_t= 0, 10) … … 870 909 rn_lf_cutoff = 5.0 ! cutoff frequency for low-pass filter [days] 871 910 rn_zdef_max = 0.9 ! maximum fractional e3t deformation 872 ln_vvl_dbg = . true.! debug prints (T/F)911 ln_vvl_dbg = .false. ! debug prints (T/F) 873 912 / 874 913 !----------------------------------------------------------------------- … … 890 929 ln_dynvor_eeT = .false. ! energy conserving scheme (een using e3t) 891 930 ln_dynvor_een = .false. ! energy & enstrophy scheme 892 nn_een_e3f = 0 ! =0 e3f = mi(mj(e3t))/4 931 nn_een_e3f = 0 ! =0 e3f = mi(mj(e3t))/4 893 932 ! ! =1 e3f = mi(mj(e3t))/mi(mj( tmask)) 894 933 ln_dynvor_msk = .false. ! vorticity multiplied by fmask (=T) ==>>> PLEASE DO NOT ACTIVATE … … 935 974 ! ! =-30 read in eddy_viscosity_3D.nc file 936 975 ! ! =-20 read in eddy_viscosity_2D.nc file 937 ! ! = 0 constant 976 ! ! = 0 constant 938 977 ! ! = 10 F(k)=c1d 939 978 ! ! = 20 F(i,j)=F(grid spacing)=c2d … … 941 980 ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) 942 981 ! ! = 32 F(i,j,k)=F(local gridscale and deformation rate) 943 ! ! time invariant coefficients : ahm = 1/2 Uv*Lv (lap case) 982 ! ! time invariant coefficients : ahm = 1/2 Uv*Lv (lap case) 944 983 ! ! or = 1/12 Uv*Lv^3 (blp case) 945 984 rn_Uv = 0.1 ! lateral viscous velocity [m/s] (nn_ahm_ijk_t= 0, 10, 20, 30) … … 1065 1104 ! = 0 constant 10 m length scale 1066 1105 ! = 1 0.5m at the equator to 30m poleward of 40 degrees 1067 rn_eice = 4 ! below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4 1106 rn_eice = 4 ! below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4 1068 1107 / 1069 1108 !----------------------------------------------------------------------- … … 1323 1362 ln_ctl = .FALSE. ! Toggle all report printing on/off (T/F); Ignored if sn_cfctl%l_config is T 1324 1363 sn_cfctl%l_config = .TRUE. ! IF .true. then control which reports are written with the following 1325 sn_cfctl%l_runstat = . FALSE.! switches and which areas produce reports with the proc integer settings.1364 sn_cfctl%l_runstat = .TRUE. ! switches and which areas produce reports with the proc integer settings. 1326 1365 sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure 1327 1366 sn_cfctl%l_oceout = .FALSE. ! that all areas report.
Note: See TracChangeset
for help on using the changeset viewer.