Difference between revisions of "sediment.in"

From WikiROMS
Jump to navigationJump to search
Line 1: Line 1:
'''The text below shows the entries in the default sediment.in file, with links to more detailed explanation where required.'''
<div class="title">Sediment Model Input Script - <span class="red">sediment.in</span></div>


! There are 5 sections to the sediment INPUT File.
The sediment model input script <span class="red">sediment.in</span> is organized in several sections as shown below, with links to more detailed explanation where required.
1) Sediment control switches, for nested grids.
2) Sediment Bed Model parameters.
3) Non-Cohesive sediment paramters.
4) Cohesive sediment paramters.
5) Bed layer and bottom sediment output properties.


!------------------------------------------------------------------------------
==Sediment model control switches==


1) Sediment model control switch.
* Switch is used to control sediment model computation within nested and/or multiple connected grids, [1:Ngrids].<div class="box">    [[Lsediment]] == T</div>


!------------------------------------------------------------------------------
==General sediment bed model controls==


Switch is used to control sediment model computation within nested and/or
* 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.<div class="box">[[newlayer_thick|NEWLAYER_THICK]] == 0.01d0</div>
multiple connected grids, [1:Ngrids].


  Lsediment == T
* Bed load transport rate coefficient, [1:Ngrids] values expected.<div class="box"> [[bedload_coeff|BEDLOAD_COEFF]] == 0.05d0</div>


!------------------------------------------------------------------------------
==Suspended Cohesive Sediment Parameters==


2) General sediment bed model controls.
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).<div class="box">      [[sd50|MUD_SD50]] == 0.01d0  0.06d0  0.10d0</div>


Depositional bed layer thickness criteria to create a new layer (m). If
* Sediment concentration (kg/m3).<div class="box">      [[Csed|MUD_CSED]] == 0.0d0  0.0d0  0.0d0</div>
deposition exceeds this value, then a new layer is created, [1:Ngrids].


    NEWLAYER_THICK == 0.01d0
* Sediment grain density (kg/m3).<div class="box">      [[Srho|MUD_SRHO]] == 2650.0d0  2650.0d0  2650.0d0</div>


Bed load transport rate coefficient. [1:Ngrids].
* Particle settling velocity (mm/s).<div class="box">      [[Wsed|MUD_WSED]] == 0.1d0  1.0d0  2.0d0</div>


     BEDLOAD_COEFF == 0.05d0
* Surface erosion rate (kg/m2/s).<div class="box">     [[Erate|MUD_ERATE]] == 5.0d-5  5.0d-5  5.0d-5</div>


!------------------------------------------------------------------------------
* Critical shear for erosion (N/m2).<div class="box">    [[tau_ce|MUD_TAU_CE]] == 0.05d0  0.05d0  0.05d0</div>


3) Suspended Cohesive Sediment Parameters, [1:NCS,1:Ngrids] values expected.
* Critical stress for deposition (N/m2), not used.<div class="box">    [[tau_cd|MUD_TAU_CD]] == 0.01d0  0.01d0  0.01d0</div>


!------------------------------------------------------------------------------
* Porosity (nondimensional): Vwater/(Vwater+Vsed). Its values range between 0 and 1.<div class="box">    [[poros|MUD_POROS]] == 0.0d0  0.0d0  0.0d0</div>


  Median sediment grain diameter (mm).
* Lateral, harmonic and biharmonic, constant, diffusion coefficient.<div class="box">      [[tnu2|MUD_TNU2]] == 0.0d0  0.0d0  0.0d0            ! m2/s<br />      [[tnu4|MUD_TNU4]] == 0.0d0  0.0d0 0.0d0            ! m4/s</div>


    MUD_SD50 == 0.01d0 0.06d0 0.10d0
* Vertical diffusion coefficients.<div class="box">  [[Akt_bak|MUD_AKT_BAK]] == 5.0d-6 5.0d-6 5.0d-6          ! m2/s</div>


Sediment concentration (kg/m3).
* Nudging/relaxation time scales, inverse scales will be computed internally.<div class="box">     [[Tnudg|MUD_TNUDG]] == 0.0d0 0.0d0 0.0d0               ! days</div>
 
     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
* Morphological time scale factor (greater than or equal to 1.0). A value of 1.0 has no scale effect.<div class="box"> [[morph_fac|MUD_MORPH_FAC]] == 1.0d0 1.0d0 1.0d0              ! nondimensional</div>


  Critical shear for erosion and deposition (N/m2).
* Logical switch (TRUE/FALSE) to activate writing of cohesive sediment into HISTORY output file.<div class="box">  [[Hout|Hout(idmud)]]  == T                              ! suspended concentration<br />  [[Hout|Hout(iMfrac)]] == T                              ! bed layer fraction<br /> [[Hout|Hout(iMmass)]] == T                              ! bed layer mass<br />  [[Hout|Hout(iMUbld)]] == T                              ! bed load at U-points<br />  [[Hout|Hout(iMVbld)]] == T                              ! bed load at V-points</div>


    MUD_TAU_CE == 0.05d0  0.05d0  0.05d0
==Non-cohesive Sediment Parameters==


    Critical stress for deposition is not used.
There are [[NNS]] non-cohesive (sand) sediment tracers per nested grid in ROMS. Therefore, [1:[[NNS]],1:[[Ngrids]]] values are expected for these parameters.
    MUD_TAU_CD == 0.01d0  0.01d0  0.01d0


Porosity (nondimensional: 0.0-1.0):  Vwater/(Vwater+Vsed).
* Median sediment grain diameter (mm).<div class="box">    [[sd50|SAND_SD50]] == 1.0d0</div>
 
    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
 
!------------------------------------------------------------------------------
 
4) Non-cohesive Sediment Parameters, [1:NNS,1:Ngrids] values expected.
 
!------------------------------------------------------------------------------
 
Median sediment grain diameter (mm).
 
  SAND_SD50 == 1.0d0


  Sediment concentration (kg/m3).
  Sediment concentration (kg/m3).
Line 162: Line 102:
Hout(iSVbld) == T                              ! bed load at V-points
Hout(iSVbld) == T                              ! bed load at V-points


!------------------------------------------------------------------------------


5) Bed layer and bottom sediment parameters, [1:Ngrids] values expected.


!------------------------------------------------------------------------------
==Bed layer and bottom sediment parameters==
 
[1:Ngrids] values expected.
 


  Hout(ithck) == T                              ! sediment layer thickness
  Hout(ithck) == T                              ! sediment layer thickness

Revision as of 15:33, 30 June 2007

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

Critical shear for erosion and deposition (N/m2).

SAND_TAU_CE == 0.07d0
SAND_TAU_CD == 0.02d0

Porosity (nondimensional: 0.0-1.0):  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(idsand) == 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
Retrieved from "https://www.myroms.org/wiki/index.php?title=sediment.in&oldid=2356"