namelist for sediment : parasedim.txt

  • These namelists define sediment behaviour and manage sediment parameters.
  • There are 2 namelist corresponding to the 2 options and 1 namelist for output
  • the option sedim_mudsand is no longer maintained.
  • Definition of each variable is given into the namelist file :

Namelist sedoutput :

variable Observations
name_out_eptot=’EPTOT’
name_out_nbniv=’NBNIV’
name_out_dzs=’DZS FOND SURF’
name_out_tenfon=’TENFON’
name_out_tenfonc=’TENFONC’
name_out_tenfonw=’TENFONW’
riog_valid_min_eptot=0
riog_valid_max_eptot=100
riog_valid_min_nbniv=0
riog_valid_max_nbniv=200
riog_valid_min_dzs=0
riog_valid_max_dzs=100
riog_valid_min_tenfon=0
riog_valid_max_tenfon=50
riog_valid_max_tenfon=50
choice_nivsed_out
nk_nivsed_out
ep_nivsed_out(1 to 5)
epmax_nivsed_out
















5 values of sediment layer thickness (mm)
maximum thickness (mm)

Namelist namsedim : parameters for CPP key_sedim_mixsed only

variable Observations
dzsmin
dzsmaxuni
l_dzsmaxuni

csedmin

l_bathy_smoothing

dt_sed

dtcal_sed

date_start_sed

l_repsed
filrepsed
l_unised
l_initsed_vardiss

l_init_hsed


z0seduni
hseduni
cseduni
ksmiuni
ksmauni
sini_sed
tini_sed
csed_mud_ini


cvolmaxsort
cvolmaxmel
csed_mud_ini
aref_sand



activlayer
ero_option

xexp_ero
frmudcr2
coef_frmudcr1
x1toce_mud, x2toce_mud
n_eros_sand, n_eros_mud
E0_sand_option

E0_sand_para
E0_mud
corfluer1,corfluer2
ero_option




xexp_ero
l_erolat_wet_cell
slopefac
coef_erolat
coef_tenfon_lat
htncrit_eros

l_fricwave
fricwav

cfreshmud
slopefac

l_diffused
choice_flxdiss_diffsed

xdifs1,xdifs2
xdifsi1,xdifsi2
epdifi

l_bioturb
xbioturbmax
xbioturbk
dbiotu0
dbiotum
frmud_db_min
frmud_db_max

l_consolid
cmudcr

xperm1,xperm2

l_morphocoupl
MF

minimum sediment layer thickness (m)
uniforme max thickness for the superficial sediment layer
set to .true. if dzsmax vary with space
user define (in sed_init) dzsmax values
concentration of the upper layer under which there is
fusion with the underlying sediment cell
.true. to apply bathymetric smoothing
(only when l_morphocoupl .TRUE.)
time step for dynamic processes in sediment
(consolidation/diffusion/bioturbation)
sub time step for dynamic processes in sediment
(< or = dt_sed)
starting date for dynamic processes in sediment
initialization
.true. if sed. variables initialized from a previous run
file from which the model is initialized
.true. for a uniform bottom initialization
set to .true. if initialization in sediment of diss.
variables,Temp, Sal will be done with bottom water conc.
.true. to adjust sediment thickness to be coherent with
sediment parameters ==> new hseduni based on cseduni,
cvolmax values, and csed_ini of each sediment
uniform bed roughness (m)
initial uniform sediment thickness(m)
initial sediment concentration
lower grid cell indices in the sediment
upper grid cell indices in the sediment
initial interstitial water uniform salinity
initial interstitial water uniform temperature
mud concentration into initial sediment (if =0.
==> csed_mud_ini=cfreshmud)
mixed sediments properties
maximum volumic concentration of sorted sand
maximum volumic concentration of mixed sediments
mud conc. into initial sed. (if =0. ==> =cfreshmud)
reference height above sediment,used sand deposit
used for sand extrapolation on water column
by default=0.02 correspond to Van Rijn experiments
sedimentary behavior : erosion
active layer thickness (m)
choice of erosion formulation for mixing sand-mud
used only if ero_option=3
critical mud fract. underwhich the behavior is purely sandy
such that critical mud fract. underwhich sandy behaviour
toce = x1toce_mud*(relative mud concentration)**x2toce_mud

choice of formulation for E0_sand evaluation
parametre for erosion flux for san
erosion flux for mud
correction factors for erosion flux computations
choice of erosion formulation for mixing sand-mud
= 0 : pure mud behavior
= 1 : linear interpolation between sand and mud behavior
= 2 : formulation derived from that of J. Vareilles (2013)
= 3 : formulations proposed by B. Mengual (2015)
used only if ero_option=3 adjustment on exp. variation
.true to take into account wet cells lateral erosion
Slope effect multiplicative factor for sliding process
slope effect multiplicative factor for erosion
parameter for lateral stress evaluation
critical water height so as to prevent erosion under a
given threshold
option wave effect in shear stress computation
wave related friction factor
sedimentary behavior : settling
prescribed fresh deposit concentration (kg/m3)
slope effect on deposit (multiplicative factor)
diffusion
.true. if diffusion within sediment is accounted for
choice expression of diss. fluxes at sed-water interface
within sediment,diffusion=xdifs1*porosity**xdifs2
at the water sediment interface
diff thickness in water at sed-water interface
bioturbation
set to .true. si bioturbation diffusion in sediment
max diffusion coeffient by bioturbation Db (in surface)
slope for bioturbation coefficient between Dbmax and 0
max depth beneath the sed. surface below which no bioturb
sed. thickness wherein bioturb coefficient Db is constant
mud fraction limit (min) below which there is no bioturb
mud fraction limit (max) above which Db is max (muddy sed)
consolidation
.true. if sediment consolidation is accounted for
critical relative concentration of the surface layer
above which no mixing is allowed with the underlying sed
permeability=xperm1*d50*d50*voidration**xperm2
morphodynamic
..true. if morphodynamic
morphological factor : multiplication factor for morpho-
-logical evolutions, equivalent to a “time acceleration”

Definition of parameters for particulate-sediment substances in variable.dat

  • particulate variables are defined in variable.dat : Defined simulated substances as :
    • GRAV or SAND or MUDS or PART : constitutive variables (kg/m3)
    • SORB or NoCP : non constitutive variables (X/m3)
  • For GRAV and SAND variables, only diameter of particles, critical stress of deposition and density of particle are given in variable.dat

  • For MUDS, PART and NoCP
    : settling velocity varies with time and space (Settling velocities : Modelling strategy) and 4 lines with 11 parameters are given in variable.dat
    • 1 option for free velocity followed by 6 parameters
      : minimum and maximum setling velocity (m/s) and :
      • If ws_free_opt=0 (Ws constant) ==>> (1)=unused (2)= unused (3)=unused (4)=unused
      • If ws_free_opt=1 (Van Leussen) ==>> (1)=c_factor (2)=c_exponent (3)=VL_a (4)=VL_b
      • If ws_free_opt=2 (Winterwerp) ==>> (1)= Primary Particle Diameter (2)=aggregation factor (3)=breakup factor (4)=fractal dimension
      • If ws_free_opt=3 (Wolanski) ==>> (1)= c_factor (2)=c_exponent (3)=unused (4)=unused
    • 1 option for hindered velocity followed by 2 parameters :
      • If ws_hind_opt=0 (no hindered settling) ==> (1)=unused (2)=unused
      • If ws_hind_opt=1 (Scott) ==> (1)=phi_factor (2)=phi_exponent
      • If ws_hind_opt=2 and ws_free_opt=2 (Winterwerp) ==> (1)=gel concentration (kg/m3) (2)=phi_exponent
      • If ws_hind_opt=3 (Wolanski) ==> (1)=c_constant (2)=c_exponent
    • diameter of particles, critical stress of deposition and density of particle are given after
  • For SORB variables, user must given the name of the associated constitutive particulate variable.

  • Inital concentration in sediment is expressed as :
    • gravel, sand, mud and PART : fraction of total initial sediment concentration (total=1)
    • SORB (or NCSP) : unit/m3 of sediment
    • DISS : unit of m3 of interstitial water
  • For SAND variable only; 2 boolean are given in variable.dat :
    • l_sand2D which is TRUE if this sand variable is treated as 2D variable (used only if key_sand2D)
    • l_outsandrouse which is TRUE if you want to use a ROUSE profil reconstituted for output in water column (used only if ikey_sand2D and l_sand2D is TRUE for this variable)

Special writing of sand in water column and of all variables in sediment (with key_sedim)

Output sand variables in water column (if key_sandbottomcell & key_sedim_mixsed)

If key_sandbottomcell is activated, sand variables with l_sandbottom =. TRUE. (in variable.dat) can be written in two different ways
  • either l_outsandrouse=.FALSE. (in parasedim.txt): only the average value in the water column, determined in 2D by the model, is written at k = 1
  • either l_outsandrouse=.TRUE. : a profile over the entire water column is artificially reconstructed from the calculated average value according to a Rouse profile.

Output variables in sediment

  • There are 4 different ways to write variables in sediment.
_images/sedim_outputsed.jpg
  • User chooses the option with parameters listed in namelist namsedoutput (parasedim.txt)
    • choice_nivsed_out : 1,2 3 or 4
    • nk_nivsed_out : number of saved sediment layers
    • ep_nivsed_out (mm) : 5 values of sediment layer thickness (mm), beginning with surface layer (choice 4). Last layer (6th) is an integrative layer.
    • epmax_nivsed_out (mm) : maximum thickness (mm) for output each layers of sediment (choice 3). Below this thickness, output an integrative layer.