DARS mesh namelist

config/namelist.cvmix_DARS

@@ -0,0 +1,72 @@
+! namelist for Turbulent Kinetic Energy (TKE) 
+&param_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
+&param_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
+&param_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
+&param_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
+&param_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'
+/
+

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.
+/
+

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.
+/

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.
+/