Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

updating namelists from 2.0 to 2.1 (the ones used by mwerner) #311

Open
wants to merge 1 commit into
base: wiso
Choose a base branch
from
Open

updating namelists from 2.0 to 2.1 (the ones used by mwerner) #311

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
69 changes: 23 additions & 46 deletions config/namelist.config
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6,81 +6,58 @@ runid='fesom'

&timestep
step_per_day=32 !96 !96 !72 !72 !45 !72 !96
run_length= 62 !62 !62 !62 !28
run_length= 40 !62 !62 !62 !28
run_length_unit='y' ! y, m, d, s
/

&clockinit ! the model starts at
&clockinit ! the model starts at
timenew=0.0
daynew=1
yearnew=1948
yearnew=1958
/

&paths
MeshPath='/work/ollie/pscholz/mesh_fesom2.0/mesh_CORE2_finaltopo_mean/'
OpbndPath=''
ClimateDataPath='/work/ollie/dsidoren/input/hydrography/phc3.0/'
ForcingDataPath='/work/ollie/dsein/input/forcing/'
TideForcingPath='/work/ollie/dsein/input/tide_forcing/'
ResultPath='../results/linfs/zmeanval/3/'
MeshPath='/work/ollie/projects/clidyn/FESOM2/meshes/core2/'
ClimateDataPath='/work/ollie/projects/clidyn/FESOM2/hydrography/phc3.0/'
ResultPath='/work/ollie/nkolduno/RC12/'
/

&initialization
OceClimaDataName='Winter_PHC3_ts.out' !which T/S data to initial. ocean
!'Winter_PHC3_ts.out'
!'Annual_PHC3_ts.out'
use_prepared_init_ice=.false. !how to init. ice; runid.initial_ice.nc
&restart_log
restart_length=1 !only required for d,h,s cases, y, m take 1
restart_length_unit='y' !output period: y, d, h, s
logfile_outfreq=960 !in logfile info. output frequency, # steps
/

&inout
restartflag='last' !restart from which saved record,'last','#'
restart_length=1 !only required for d,h,s cases, y, m take 1
restart_length_unit='y' !output period: y, d, h, s
restart_offset=64
logfile_outfreq=960 !in logfile info. output frequency, # steps
use_means=.false. !average output, if false prints out snapshots
/

&mesh_def
grid_type=1 !1 z-level, 2 sigma, 3 z+sigma
use_ALE=.true. ! switch on/off ALE
which_ALE='linfs' ! 'linfs','zlevel', 'zstar','zstar-weight', 'ztilde'
use_partial_cell=.false.
min_hnode=0.5
lzstar_lev=4
&ale_def
which_ALE='zstar' ! 'linfs','zlevel', 'zstar'
use_partial_cell=.true.
/

&geometry
cartesian=.false.
fplane=.false.
betaplane=.false.
f_fplane=-1.4e-4 ![1/s]
beta_betaplane=2.0e-11 ![1/s/m]
cyclic_length=360. ![degree]
rotated_grid=.true. !option only valid for coupled model case now
alphaEuler=50. ![degree] Euler angles, convention:
betaEuler=15. ![degree] first around z, then around new x,
gammaEuler=-90. ![degree] then around new z.
cyclic_length=360 ![degree]
rotated_grid=.true. !option only valid for coupled model case now
force_rotation=.false.
alphaEuler=50. ![degree] Euler angles, convention:
betaEuler=15. ![degree] first around z, then around new x,
gammaEuler=-90. ![degree] then around new z.
/

&calendar
include_fleapyear=.false.
include_fleapyear=.true.
/

&run_config
use_ice=.true. ! ocean+ice
use_cavity=.false. !
use_cavity_partial_cell=.false.
use_floatice = .false.
use_sw_pene=.true.
toy_ocean=.false. ! use toy forcing/initialization
lwiso=.false.
/

&machine
! system_arch=1
n_levels=3
!n_part= 2, 4, 36 ! number of partitions on each hierarchy level
n_part= 2, 12, 36 ! 864 number of partitions on each hierarchy level
!n_part= 2, 6, 36 ! 432 number of partitions on each hierarchy level
!n_part= 2, 4, 36 ! 288 number of partitions on each hierarchy level

n_part= 2,4, 36 ! 432 number of partitions on each hierarchy level
/
63 changes: 63 additions & 0 deletions config/namelist.config.linfs
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,63 @@
! This is the namelist file for model general configuration

&modelname
runid='fesom'
/

&timestep
step_per_day=32 !96 !96 !72 !72 !45 !72 !96
run_length= 40 !62 !62 !62 !28
run_length_unit='y' ! y, m, d, s
/

&clockinit ! the model starts at
timenew=0.0
daynew=1
yearnew=1958
/

&paths
MeshPath='/work/ollie/projects/clidyn/FESOM2/meshes/core2/'
ClimateDataPath='/work/ollie/projects/clidyn/FESOM2/hydrography/phc3.0/'
ResultPath='/work/ollie/nkolduno/RC12/'
/

&restart_log
restart_length=1 !only required for d,h,s cases, y, m take 1
restart_length_unit='y' !output period: y, d, h, s
logfile_outfreq=960 !in logfile info. output frequency, # steps
/

&ale_def
which_ALE='linfs' ! 'linfs','zlevel', 'zstar'
use_partial_cell=.true.
/

&geometry
cartesian=.false.
fplane=.false.
cyclic_length=360 ![degree]
rotated_grid=.true. !option only valid for coupled model case now
force_rotation=.false.
alphaEuler=50. ![degree] Euler angles, convention:
betaEuler=15. ![degree] first around z, then around new x,
gammaEuler=-90. ![degree] then around new z.
/

&calendar
include_fleapyear=.true.
/

&run_config
use_ice=.true. ! ocean+ice
use_cavity=.false. !
use_cavity_partial_cell=.false.
use_floatice = .false.
use_sw_pene=.true.
lwiso=.false.
/

&machine
n_levels=3
n_part= 2,4, 36 ! 432 number of partitions on each hierarchy level
/
63 changes: 63 additions & 0 deletions config/namelist.config.zstar
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,63 @@
! This is the namelist file for model general configuration

&modelname
runid='fesom'
/

&timestep
step_per_day=32 !96 !96 !72 !72 !45 !72 !96
run_length= 40 !62 !62 !62 !28
run_length_unit='y' ! y, m, d, s
/

&clockinit ! the model starts at
timenew=0.0
daynew=1
yearnew=1958
/

&paths
MeshPath='/work/ollie/projects/clidyn/FESOM2/meshes/core2/'
ClimateDataPath='/work/ollie/projects/clidyn/FESOM2/hydrography/phc3.0/'
ResultPath='/work/ollie/nkolduno/RC12/'
/

&restart_log
restart_length=1 !only required for d,h,s cases, y, m take 1
restart_length_unit='y' !output period: y, d, h, s
logfile_outfreq=960 !in logfile info. output frequency, # steps
/

&ale_def
which_ALE='zstar' ! 'linfs','zlevel', 'zstar'
use_partial_cell=.true.
/

&geometry
cartesian=.false.
fplane=.false.
cyclic_length=360 ![degree]
rotated_grid=.true. !option only valid for coupled model case now
force_rotation=.false.
alphaEuler=50. ![degree] Euler angles, convention:
betaEuler=15. ![degree] first around z, then around new x,
gammaEuler=-90. ![degree] then around new z.
/

&calendar
include_fleapyear=.true.
/

&run_config
use_ice=.true. ! ocean+ice
use_cavity=.false. !
use_cavity_partial_cell=.false.
use_floatice = .false.
use_sw_pene=.true.
lwiso=.false.
/

&machine
n_levels=3
n_part= 2,4, 36 ! 432 number of partitions on each hierarchy level
/
72 changes: 72 additions & 0 deletions config/namelist.cvmix
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -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 = 86400.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 = '/work/ollie/clidyn/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 = '/work/ollie/pscholz/FORCING/IDEMIX/tidal_energy_gx1v6_20090205_rgrid.nc'
/

Loading