1  Stitle (a80)
ROMS - SED_TOY (consolidation) - Test

2  Sd(1-NST), CSED,  SRHO,   WSED,    ERATE,  TAU_CE, TAU_CD, BED_FRAC(1:NLAY)
      0.0625    0.   2650.    2.     0.0015      0.1     0.1    41*0.25
      0.140     0.   2650.    8.     0.0015      0.1     0.1    41*0.25
      0.004     0.   2650.    0.1    0.0005      0.05      0.1    41*0.25
      0.030     0.   2650.    0.6    0.0005      0.05      0.1    41*0.25

3 BTHK(1:NLAY)
      41*0.001
            
4 BPOR(1:NLAY)
      41*0.6

5 Hrip             
      0.01

6 Lrip
      0.1  

7 bedload_coeff 
      1.

8 morph_fac
      1. 

9 transC
      0.03  
 
10 transN
      0.2

11  tcr_min
      0.030

12  tcr_max
      1.5

13  tcr_slp
      2

14  tcr_off
      3.4d0

15 tcr_tim
      28800.0d0

99 END of sediment input data

!
!=====================================================================
!  Copyright (c) 2004 UCLA/Rutgers
!==================================== Meinte Blaas | Hernan G. Arango =
!
   GLOSSARY
   ========

*** CARD 1: String with a maximum of eighty characters.

  Stitle    Sediment case title.


*** CARD 2: Sediment grain parameters & initial values (NST lines).
            Non-Cohesive sediment first 

  Sd        Diameter of grain size class [mm].

  CSED      Initial concentration (spatially uniform) [mg/l].  
  
  SRHO      Density of sediment material of size class [kg/m3]. 
            Quartz: SRHO=2650 kg/m3  
  
  WSED      Settling velocity of size class [mm/s].   
            Typically WSED  = 1E3*  {visk*(sqrt(10.36^2+1.049*D^3)-10.36)/d50} [mm/s] (Soulsby 97)
            with  D=d50*(g*(srho/rhow-1)/(visk^2))^0.33333
                  rhow=1000 kg/m3; d50=Sd*1.e-3 m; 
                  visk=1.3e-3/rhow m2/s; g=9.8 m/s2  
            
  ERATE     Erosion rate of size class [kg/m2/s].  
            Typically ERATE  =gam0*WSED*1.e-3*SRHO [kg/m2s];
            with gam0=1E-3 - 1E-5 (Smith & McLean, Drake&al)
             
  TAU_CE    Critical shear stress for sediment motion [N/m2]
            (initiation of bedload for coarses, suspension for fines). 
            Typically TAU_CE =tau_cb=thetcr.*(g*(srho-rhow).*d50) [N/m2]; (Shields, bedload)
            with thetcr = 0.3./(1+1.2*D) + 0.055*(1-exp(-0.02*D)) (Soulsby & Whitehouse 97) 
            TAU_CE   = rhow*(0.8 WSED*1E-3)^2 [N/m2] (suspended load)
            
  TAU_CD    Critical shear stress for depostion of cohesive sediments [N/m2]
  
  BED_FRAC  Volume fraction of each size class in each bed layer (NLAY columns) [0<BED_FRAC<1] 


*** CARD 3: Sediment bed thickness, 1st field is top layer ('delt_a')

  BTHK      Initial thicknesses of bed layers [m]
            Bthk(1) active layer thickness, fixed in simulation unless SUM(Bthk(:))<Bthk(1)
  
*** CARD 4: Sediment bed porosity

  BPOR     Initial porosity of bed layers [m]
            used in ana_sediment ifdef ANA_SEDIMENT || not in init.nc*
  
*** CARD 5: Bottom ripple height         

  Hrip     Initial ripple height [m]
            used in ana_sediment ifdef ANA_SEDIMENT || not in init.nc*

*** CARD 6: Bottom ripple length

  Lrip     Initial ripple length [m]
            used in ana_sediment ifdef ANA_SEDIMENT || not in init.nc*
  
*** Card 9 :

  transC   Cohesive transition- Under that value of total mud fraction
            entire bed behaves as a non-cohesive bed

*** Card 10 :

  transN   Noncohesive transition- Over that value of total mud fraction
            entire bed behaves as a cohesive bed

*** Card 11 :

  tcr_min  Minimum shear for erosion

*** Card 12 :

  tcr_max  Maximum shear for erosion

*** Card 13 :

  tcr_slp  Tau_crit profile slope

*** Card 14 :

  tcr_off  Tau_crit profile offset

*** Card 15 :

  tcr_tim  Tau_crit consolidation rate