sediment.in

From WikiROMS
Jump to navigationJump to search
Sediment Model Input Script - sediment.in

The sediment model input script sediment.in is organized in several sections as shown below, with links to more detailed explanation where required.

Sediment model control switches

  • Switch is used to control sediment model computation within nested and/or multiple connected grids, [1:Ngrids].
    Lsediment == T

General sediment bed model controls

  • Depositional bed layer thickness (m) criteria to create a new layer. If deposition exceeds this value, then a new layer is created, [1:Ngrids] values expected.
    NEWLAYER_THICK == 0.01d0
  • Bed load transport rate coefficient, [1:Ngrids] values expected.
    BEDLOAD_COEFF == 0.05d0

Suspended Cohesive Sediment Parameters

There are NCS cohesive (mud) sediment tracers per nested grid in ROMS. Therefore, [1:NCS,1:Ngrids] values are expected for these parameters.

  • Median sediment grain diameter (mm).
    MUD_SD50 == 0.01d0 0.06d0 0.10d0
  • Sediment concentration (kg/m3).
    MUD_CSED == 0.0d0 0.0d0 0.0d0
  • Sediment grain density (kg/m3).
    MUD_SRHO == 2650.0d0 2650.0d0 2650.0d0
  • Particle settling velocity (mm/s).
    MUD_WSED == 0.1d0 1.0d0 2.0d0
  • Surface erosion rate (kg/m2/s).
    MUD_ERATE == 5.0d-5 5.0d-5 5.0d-5
  • Critical shear for erosion (N/m2).
    MUD_TAU_CE == 0.05d0 0.05d0 0.05d0
  • Critical stress for deposition (N/m2), not used.
    MUD_TAU_CD == 0.01d0 0.01d0 0.01d0
  • Porosity (nondimensional): Vwater/(Vwater+Vsed). Its values range between 0 and 1.
    MUD_POROS == 0.0d0 0.0d0 0.0d0
  • Lateral, harmonic and biharmonic, constant, diffusion coefficient.
    MUD_TNU2 == 0.0d0 0.0d0 0.0d0  ! m2/s
    MUD_TNU4 == 0.0d0 0.0d0 0.0d0  ! m4/s
  • Vertical diffusion coefficients.
    MUD_AKT_BAK == 5.0d-6 5.0d-6 5.0d-6  ! m2/s
  • Nudging/relaxation time scales, inverse scales will be computed internally.
    MUD_TNUDG == 0.0d0 0.0d0 0.0d0  ! days
  • Morphological time scale factor (greater than or equal to 1.0). A value of 1.0 has no scale effect.
    MUD_MORPH_FAC == 1.0d0 1.0d0 1.0d0  ! nondimensional
  • Logical switch (TRUE/FALSE) to activate writing of cohesive sediment into HISTORY output file.
    Hout(idmud) == T  ! suspended concentration
    Hout(iMfrac) == T  ! bed layer fraction
    Hout(iMmass) == T  ! bed layer mass
    Hout(iMUbld) == T  ! bed load at U-points
    Hout(iMVbld) == T  ! bed load at V-points

Non-cohesive Sediment Parameters

There are NNS non-cohesive (sand) sediment tracers per nested grid in ROMS. Therefore, [1:NNS,1:Ngrids] values are expected for these parameters.

  • Median sediment grain diameter (mm).
    SAND_SD50 == 1.0d0
  • Sediment concentration (kg/m3).
    SAND_CSED == 0.0d0
  • Sediment grain density (kg/m3).
    SAND_SRHO == 2650.0d0
  • Particle settling velocity (mm/s).
    SAND_WSED == 1.0d0
  • Surface erosion rate (kg/m2/s).
    SAND_ERATE == 1.0d-5
  • Porosity (nondimensional): Vwater/(Vwater+Vsed).
    SAND_POROS == 0.5d0
  • Lateral, harmonic and biharmonic, constant, diffusion coefficient.
    SAND_TNU2 == 0.0d0  ! m2/s
    SAND_TNU4 == 0.0d0  ! m4/s
  • Vertical diffusion coefficients.
    SAND_AKT_BAK == 5.0d-6  ! m2/s
  • Nudging/relaxation time scales, inverse scales will be computed internally.
    SAND_TNUDG == 0.0d0  ! days
  • Morphological time scale factor (greater than or equal to 1.0). A value of 1.0 has no scale effect.
    SAND_MORPH_FAC == 1.0d0 1.0d0 1.0d0  ! nondimensional
  • Logical switches (TRUE/FALSE) to activate writing of non-cohesive sediment into HISTORY output file.
    Hout(isand) == T  ! suspended concentration
    Hout(iSfrac) == T  ! bed layer fraction
    Hout(iSmass) == T  ! bed layer mass
    Hout(iSUbld) == T  ! bed load at U-points
    Hout(iSVbld) == T  ! bed load at V-points

Bed layer and bottom sediment parameters

[1:Ngrids] values expected.


Hout(ithck) == T                               ! sediment layer thickness
Hout(iaged) == T                               ! sediment layer age
Hout(iporo) == T                               ! sediment layer porosity
Hout(idiff) == F                               ! biodiffusivity