diff --git a/config/namelist.cvmix_DARS b/config/namelist.cvmix_DARS new file mode 100644 index 000000000..ca8c5a6fe --- /dev/null +++ b/config/namelist.cvmix_DARS @@ -0,0 +1,72 @@ +! namelist for Turbulent Kinetic Energy (TKE) +¶m_tke +tke_c_k = 0.1 +tke_c_eps = 0.7 +tke_alpha = 30.0 +tke_mxl_min = 1.0e-8 +tke_kappaM_min = 0.0 +tke_kappaM_max = 100.0 +tke_cd = 3.75 ! for Dirichlet boundary conditions +!tke_cd = 1.0 ! for Neumann boundary conditions +tke_surf_min = 1.0e-4 +tke_min = 1.0e-6 +! tke_mxl_choice ... Can only be 1 or 2, choice of calculation of mixing +! length; currently only Blanke, B., P. Delecluse option is implemented +tke_mxl_choice = 2 +/ + +! namelist for IDEMIX +¶m_idemix +idemix_tau_v = 172800.0 ! time scale for vertical symmetrisation (sec) +idemix_tau_h = 1296000.0 ! time scale for horizontal symmetrisation +idemix_gamma = 1.570 ! constant of order one derived from the shape of the spectrum in m space (dimensionless) +idemix_jstar = 10.0 ! spectral bandwidth in modes (dimensionless) +idemix_mu0 = 1.33333333 ! dissipation parameter (dimensionless) +idemix_sforcusage = 0.2 +idemix_n_hor_iwe_prop_iter = 5 ! iterations for contribution from horiz. wave propagation +idemix_surforc_file = '/work/ollie/clidyn/forcing/IDEMIX/fourier_smooth_2005_cfsr_inert_rgrid.nc' +idemix_botforc_file = '/albedo/work/projects/p_pool_fesom1/fesom2.0/forcing/idemix/tidal_energy_gx1v6_20090205_rgrid.nc' +/ + +! namelist for PP +¶m_pp +pp_use_fesompp = .true. ! if true use fesom flavor off PP if false use original PP of Pacanowski & Philander +pp_Av0 = 0.01 ! mixing coeff. for PP +pp_alpha = 5.0 ! alpha param eq. (1) in Pacanowski and Philander 1981 +pp_exp = 2.0 ! n param eq. (1) in Pacanowski and Philander 1981 +pp_Avbckg = 1.0e-4 ! const. pp background viscosity +pp_Kvbckg = 1.0e-5 ! const. pp background diffusivity +pp_use_nonconstKvb = .true. ! use lat and depth dependent background diffusivity +/ + +! namelist for KPP +¶m_kpp +kpp_use_fesomkpp = .false. ! if true use cvmix MOM5 like kpp if false use cvmix MOM6 like kpp +kpp_use_enhanceKv = .true. ! If true, add enhanced diffusivity at base of boundary layer +kpp_use_compEkman = .true. ! If true, compute Ekman depth limit for OBLdepth +kpp_use_monob = .true. ! If true, compute Monin-Obukhov limit for OBLdepth +kpp_interptype_ri = "linear" ! Type of interpolation in determining OBL depth: linear,quadratic,cubic +kpp_interptype_atobl= "LMD94" ! Type of interpolation of visc and diff at obl depth: linear,quadratic,cubic,LMD94 +kpp_matchtechc = "ParabolicNonLocal" ! Setting diffusivity and NLT profile functions:SimpleShapes, MatchGradient, MatchBoth, ParabolicNonLocal +kpp_internalmix = "KPP" ! Ri-number dependet mixing scheme below the OBL: 'PP' or 'KPP' +kpp_reduce_tauuice = .false. ! If True, reduce the wind stress (ustar) under sea ice. +kpp_Av0 = 5.0e-3 ! leading coefficient of shear mixing formula, units: m^2/s: default= 5e-3 +kpp_Kv0 = 5.0e-3 +kpp_Ri0 = 0.7 ! critical Richardson number value, units: unitless (0.7 in LMD94) +kpp_pp_Av0 = 0.01 ! Parameter in case of PP mixing below the OBL +kpp_use_nonconstKvb= .true. ! If True use non constant background diffusivity of Qiang from FESOM1.4 +kpp_Avbckg = 1.0e-4 ! Values for const. background viscosity and diffusivity +kpp_Kvbckg = 1.0e-5 +/ + +! namelist for TIDAL +¶m_tidal +tidal_mixscheme = "Simmons" +tidal_efficiency = 0.2 ! units: unitless (fraction); (Gamma in Simmons et al.) +tidal_vert_decayscale= 500.0 ! units: m; zeta in Simmons et al. (used to compute the vertical deposition function) +tidal_max_coeff = 50e-4 ! units: m^2/s; largest acceptable value for diffusivity +tidal_lcl_mixfrac = 0.33 !tidal dissipation efficiency (q in Simmons et al.), i.e. fraction of energy that dissipates locally +tidal_depth_cutoff = 0.0 !depth of the shallowest column where tidal mixing is computed (positive below surface) +tidal_botforc_file = '/albedo/work/projects/p_pool_fesom1/fesom2.0/forcing/idemix/tidal_energy_gx1v6_20090205_rgrid.nc' +/ + diff --git a/config/namelist.dyn_DARS b/config/namelist.dyn_DARS new file mode 100644 index 000000000..2e60ea00d --- /dev/null +++ b/config/namelist.dyn_DARS @@ -0,0 +1,24 @@ +&dynamics_visc +visc_gamma0 = 0.003 ! [m/s], backgroung viscosity= gamma0*len, it should be as small a s possible (keep it < 0.01 m/s). +visc_gamma1 = 0.2 ! [nodim], for computation of the flow aware viscosity +visc_gamma2 = 0.0285 ! [s/m], is only used in easy backscatter option +visc_easybsreturn= 1 + +opt_visc = 7 +! 5=Kinematic (easy) Backscatter +! 6=Biharmonic flow aware (viscosity depends on velocity Laplacian) +! 7=Biharmonic flow aware (viscosity depends on velocity differences) +! 8=Dynamic Backscatter + +use_ivertvisc= .true. +/ + +&dynamics_general +momadv_opt = 2 ! option for momentum advection in moment only =2 +use_freeslip = .false. ! Switch on free slip +use_wsplit = .true. ! Switch for implicite/explicte splitting of vert. velocity +wsplit_maxcfl= 1.0 ! maximum allowed CFL criteria in vertical (0.5 < w_max_cfl < 1.) + ! in older FESOM it used to be w_exp_max=1.e-3 +ldiag_KE=.false. +/ + diff --git a/config/namelist.forcing_DARS b/config/namelist.forcing_DARS new file mode 100644 index 000000000..3e6fab32d --- /dev/null +++ b/config/namelist.forcing_DARS @@ -0,0 +1,72 @@ +! This is the namelist file for forcing + +&forcing_exchange_coeff +Ce_atm_oce=0.00175 ! exchange coeff. of latent heat over open water +Ch_atm_oce=0.00175 ! exchange coeff. of sensible heat over open water +Cd_atm_oce=0.001 ! drag coefficient between atmosphere and water +Ce_atm_ice=0.00175 ! exchange coeff. of latent heat over ice +Ch_atm_ice=0.00175 ! exchange coeff. of sensible heat over ice +Cd_atm_ice=0.0012 ! drag coefficient between atmosphere and ice +Swind =0.0 ! parameterization for coupled current feedback +/ + +&forcing_bulk +AOMIP_drag_coeff=.false. +ncar_bulk_formulae=.false. +ncar_bulk_z_wind=10.0 ! height at which wind forcing is located (CORE, JRA-do: 10m, JRA, NCEP:2m) +ncar_bulk_z_tair=2.0 ! height at which temp forcing is located (CORE, JRA-do: 10m, JRA, NCEP:2m) +ncar_bulk_z_shum=2.0 ! height at which humi forcing is located (CORE, JRA-do: 10m, JRA, NCEP:2m) + +/ + +&land_ice +use_landice_water=.false. +landice_start_mon=5 +landice_end_mon=10 +/ + + +&age_tracer +use_age_tracer=.false. +use_age_mask=.false. +age_tracer_path='./mesh/' +age_start_year=2000 + +/ + +&nam_sbc + nm_xwind_file = '' ! name of file with winds, if nm_sbc=2 + nm_ywind_file = '' ! name of file with winds, if nm_sbc=2 + nm_humi_file = '' ! name of file with humidity + nm_qsr_file = '' ! name of file with solar heat + nm_qlw_file = '' ! name of file with Long wave + nm_tair_file = '' ! name of file with 2m air temperature + nm_prec_file = '' ! name of file with total precipitation + nm_snow_file = '' ! name of file with snow precipitation + nm_mslp_file = '' ! air_pressure_at_sea_level + nm_xwind_var = 'u' ! name of variable in file with wind + nm_ywind_var = 'v' ! name of variable in file with wind + nm_humi_var = 'q' ! name of variable in file with humidity + nm_qsr_var = 'ssrd' ! name of variable in file with solar heat + nm_qlw_var = 'strd' ! name of variable in file with Long wave + nm_tair_var = 't2m' ! name of variable in file with 2m air temperature + nm_prec_var = 'rf' ! name of variable in file with total precipitation + nm_snow_var = 'sf' ! name of variable in file with total precipitation + nm_mslp_var = 'sp' ! name of variable in file with air_pressure_at_sea_level + nm_nc_iyear = 1900 + nm_nc_imm = 1 ! initial month of time axis in netCDF + nm_nc_idd = 1 ! initial day of time axis in netCDF + nm_nc_freq = 1 ! data points per day (i.e. 86400 if the time axis is in seconds) + nm_nc_tmid = 0 ! 1 if the time stamps are given at the mid points of the netcdf file, 0 otherwise (i.e. 1 in CORE1, CORE2; 0 in JRA55) + l_xwind=.false. l_ywind=.false. l_humi=.false. l_qsr=.false. l_qlw=.false. l_tair=.false. l_prec=.false. l_mslp=.false. l_cloud=.false. l_snow=.false. + runoff_data_source ='CORE2' !Dai09, CORE2 + nm_runoff_file ='./forcing/CORE2/runoff.nc' + sss_data_source ='CORE2' + nm_sss_data_file ='./forcing/CORE2/PHC2_salx.nc' + chl_data_source ='Sweeney' !'Sweeney' monthly chlorophyll climatology or 'NONE' for constant chl_const (below). Make use_sw_pene=.TRUE. in namelist.config! + nm_chl_data_file ='./forcing/Sweeney_2005.nc' + chl_const = 0.1 + use_runoff_mapper = .TRUE. + runoff_basins_file = 'runoff_maps_original.nc' + runoff_radius = 500000. +/ diff --git a/config/namelist.tra_DARS b/config/namelist.tra_DARS new file mode 100644 index 000000000..1ed07a2d2 --- /dev/null +++ b/config/namelist.tra_DARS @@ -0,0 +1,57 @@ +&tracer_listsize +num_tracers=100 !number of tracers to allocate. shallbe large or equal to the number of streams in &nml_list +/ + +&tracer_list +nml_tracer_list = +1 , 'MFCT', 'QR4C', 'FCT ', 0., 1., +2 , 'MFCT', 'QR4C', 'FCT ', 0., 1., +!101, 'UPW1', 'UPW1', 'NON ', 0., 0. +/ + +&tracer_init3d ! initial conditions for tracers +n_ic3d = 2 ! number of tracers to initialize +idlist = 2, 1 ! their IDs (0 is temperature, 1 is salinity, etc.). The reading order is defined here! +filelist = 'phc3.0_winter.nc', 'phc3.0_winter.nc' ! list of files in ClimateDataPath to read (one file per tracer), same order as idlist +varlist = 'salt', 'temp' ! variables to read from specified files +t_insitu = .true. ! if T is insitu it will be converted to potential after reading it +/ + +&tracer_init2d ! initial conditions for 2D tracers (sea ice) +n_ic2d = 3 ! number of tracers to initialize +idlist = 1, 2, 3 ! their IDs (0 is a_ice, 1 is m_ice, 3 m_snow). The reading order is defined here! +filelist = 'a_ice.nc', 'm_ice.nc', 'm_snow.nc' ! list of files in ClimateDataPath to read (one file per tracer), same order as idlist +varlist = 'a_ice', 'm_ice', 'm_snow' ! variables to read from specified files +ini_ice_from_file=.false. +/ + +&tracer_general +! bharmonic diffusion for tracers. We recommend to use this option in very high resolution runs (Redi is generally off there). +smooth_bh_tra =.false. ! use biharmonic diffusion (filter implementation) for tracers +gamma0_tra = 0.0005 ! gammaX_tra are analogous to those in the dynamical part +gamma1_tra = 0.0125 +gamma2_tra = 0. +i_vert_diff =.true. +/ + +&tracer_phys +use_momix = .false. ! switch on/off !Monin-Obukhov -> TB04 mixing +momix_lat = -50.0 ! latitidinal treshhold for TB04, =90 --> global +momix_kv = 0.01 ! PP/KPP, mixing coefficient within MO length +use_instabmix = .true. ! enhance convection in case of instable stratification +instabmix_kv = 0.1 +use_windmix = .false. ! enhance mixing trough wind only for PP mixing (for stability) +windmix_kv = 1.e-3 +windmix_nl = 2 +diff_sh_limit=5.0e-3 ! for KPP, max diff due to shear instability +Kv0_const=.true. +double_diffusion=.false. ! for KPP,dd switch +K_ver=1.0e-5 +K_hor=0. +surf_relax_T=0.0 +surf_relax_S=0.0 ! 1.929e-06 ! 50m/300days 6.43e-07! m/s 10./(180.*86400.) +balance_salt_water =.true. ! balance virtual-salt or freshwater flux or not +clim_relax=0.0 ! 1/s, geometrical information has to be supplied +ref_sss_local=.true. +ref_sss=34. +/