1 | 1 Stitle (a80) |
---|
2 | ROMS - Sediment - Test |
---|
3 | |
---|
4 | 2 Sd(1-NST), CSED, SRHO, WSED, ERATE, TAU_CE, TAU_CD, BED_FRAC(1:NLAY) |
---|
5 | 0.125 9.9 2650. 9.4 25.0e-5 0.05 0.14 0.4 0.4 |
---|
6 | 0.050 0.0 2650. 1.6 4.0e-5 0.01 0.14 0.6 0.6 |
---|
7 | |
---|
8 | 3 BTHK(1:NLAY) |
---|
9 | 0.3E-2 10. |
---|
10 | |
---|
11 | 4 BPOR(1:NLAY) |
---|
12 | 0.41 0.42 |
---|
13 | |
---|
14 | 5 Hrip |
---|
15 | 0.03 |
---|
16 | |
---|
17 | 6 Lrip |
---|
18 | 0.14 |
---|
19 | |
---|
20 | 99 END of sediment input data |
---|
21 | |
---|
22 | ! |
---|
23 | !===================================================================== |
---|
24 | ! Copyright (c) 2004 UCLA/Rutgers |
---|
25 | !==================================== Meinte Blaas | Hernan G. Arango = |
---|
26 | ! |
---|
27 | GLOSSARY |
---|
28 | ======== |
---|
29 | |
---|
30 | *** CARD 1: String with a maximum of eighty characters. |
---|
31 | |
---|
32 | Stitle Sediment case title. |
---|
33 | |
---|
34 | |
---|
35 | *** CARD 2: Sediment grain parameters & initial values (NST lines). |
---|
36 | |
---|
37 | Sd Diameter of grain size class [mm]. |
---|
38 | |
---|
39 | CSED Initial concentration (spatially uniform) [mg/l]. |
---|
40 | |
---|
41 | SRHO Density of sediment material of size class [kg/m3]. |
---|
42 | Quartz: SRHO=2650 kg/m3 |
---|
43 | |
---|
44 | WSED Settling velocity of size class [mm/s]. |
---|
45 | Typically WSED = 1E3* {visk*(sqrt(10.36^2+1.049*D^3)-10.36)/d50} [mm/s] (Soulsby 97) |
---|
46 | with D=d50*(g*(srho/rhow-1)/(visk^2))^0.33333 |
---|
47 | rhow=1000 kg/m3; d50=Sd*1.e-3 m; |
---|
48 | visk=1.3e-3/rhow m2/s; g=9.8 m/s2 |
---|
49 | |
---|
50 | ERATE Erosion rate of size class [kg/m2/s]. |
---|
51 | Typically ERATE =gam0*WSED*1.e-3*SRHO [kg/m2s]; |
---|
52 | with gam0=1E-3 - 1E-5 (Smith & McLean, Drake&al) |
---|
53 | |
---|
54 | TAU_CE Critical shear stress for sediment motion [N/m2] |
---|
55 | (initiation of bedload for coarses, suspension for fines). |
---|
56 | Typically TAU_CE =tau_cb=thetcr.*(g*(srho-rhow).*d50) [N/m2]; (Shields, bedload) |
---|
57 | with thetcr = 0.3./(1+1.2*D) + 0.055*(1-exp(-0.02*D)) (Soulsby & Whitehouse 97) |
---|
58 | TAU_CE = rhow*(0.8 WSED*1E-3)^2 [N/m2] (suspended load) |
---|
59 | |
---|
60 | TAU_CD Critical shear stress for depostion of cohesive sediments [N/m2] |
---|
61 | |
---|
62 | BED_FRAC Volume fraction of each size class in each bed layer (NLAY columns) [0<BED_FRAC<1] |
---|
63 | |
---|
64 | |
---|
65 | *** CARD 3: Sediment bed thickness, 1st field is top layer ('delt_a') |
---|
66 | |
---|
67 | BTHK Initial thicknesses of bed layers [m] |
---|
68 | Bthk(1) active layer thickness, fixed in simulation unless SUM(Bthk(:))<Bthk(1) |
---|
69 | |
---|
70 | *** CARD 4: Sediment bed porosity |
---|
71 | |
---|
72 | BPOR Initial porosity of bed layers [m] |
---|
73 | used in ana_sediment ifdef ANA_SEDIMENT || not in init.nc* |
---|
74 | |
---|
75 | *** CARD 5: Bottom ripple height |
---|
76 | |
---|
77 | Hrip Initial ripple height [m] |
---|
78 | used in ana_sediment ifdef ANA_SEDIMENT || not in init.nc* |
---|
79 | |
---|
80 | *** CARD 6: Bottom ripple length |
---|
81 | |
---|
82 | Lrip Initial ripple length [m] |
---|
83 | used in ana_sediment ifdef ANA_SEDIMENT || not in init.nc* |
---|
84 | |
---|