sediment.in: Difference between revisions
No edit summary (change visibility) |
No edit summary (change visibility) |
||
Line 1: | Line 1: | ||
<div class="title">Sediment Model Input Script - <span class="red">sediment.in</span></div> | |||
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. | |||
==Sediment model control switches== | |||
1 | * 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== | |||
* 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> | |||
* Bed load transport rate coefficient, [1:Ngrids] values expected.<div class="box"> [[bedload_coeff|BEDLOAD_COEFF]] == 0.05d0</div> | |||
==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).<div class="box"> [[sd50|MUD_SD50]] == 0.01d0 0.06d0 0.10d0</div> | |||
* Sediment concentration (kg/m3).<div class="box"> [[Csed|MUD_CSED]] == 0.0d0 0.0d0 0.0d0</div> | |||
* Sediment grain density (kg/m3).<div class="box"> [[Srho|MUD_SRHO]] == 2650.0d0 2650.0d0 2650.0d0</div> | |||
* Particle settling velocity (mm/s).<div class="box"> [[Wsed|MUD_WSED]] == 0.1d0 1.0d0 2.0d0</div> | |||
* 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> | |||
* 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> | |||
* 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> | |||
* Vertical diffusion coefficients.<div class="box"> [[Akt_bak|MUD_AKT_BAK]] == 5.0d-6 5.0d-6 5.0d-6 ! m2/s</div> | |||
* Nudging/relaxation time scales, inverse scales will be computed internally.<div class="box"> [[Tnudg|MUD_TNUDG]] == 0.0d0 0.0d0 0.0d0 ! days</div> | |||
* 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> | |||
* 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> | |||
==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).<div class="box"> [[sd50|SAND_SD50]] == 1.0d0</div> | |||
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 | ||
==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
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.
- 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