PROCEDURE
[Features]
[Main Page]
[Sample input]
[Rationale]
Column 1: Size fraction range in mm.
The normal list has one (1) fine (silt and clay) and eight (8) coarse (sand and gravel) = nine (9) size fraction ranges.
The long list has two (2) fine (silt and clay) and eight (8) coarse (sand and gravel) = ten (10) size fraction ranges.
Column 2: Geometric mean [of the size fraction range limits] D (mm).
Column 3: Geometric mean [of the size fraction range limits] D (ft).
Column 4: The intensity of shear: ψ = the largest of either 
ψ = 1.65 D_{35} / (SR) or
 ψ = 0.66 D / (SR)
The sampled (measured) suspended load: Q_{sm} = 0.0027 C Q
(in tons/day, or T/D)
K_{s} = D_{65}
Assume x = 1.54
(SR)^{1/2} = u / [32.63 log_{10}(12.27 x d/K_{s})]
The shear velocity: u'_{*} = (gRS)^{1/2}
The kinematic viscosity ν is a function of temperature t (Plate 2, with minor corrections).
The thickness of the laminar sublayer: δ = 11.6 ν / u'_{*}
Given K_{s} / δ, use Plate 3 to find new value of x. Recalculate if different from the assumed value.
P = 2.303 log_{10} (30.2 x d / K_{s})
A' = d_{n} / d_{s}
With P and A', use Plate 4 to find the percentage of flow in sampled zone.
The sediment discharge through the sampled zone: Q'_{s} Total = (percentage of flow in sampled zone) * Q_{sm}
Column 5: Given ψ, use Plate 5 to calculate the intensity of bedload transport
(φ_{*}/2)
Column 6: 1,200 D^{3/2}
Column 7: fraction of bed material i_{b} (input array data)
Column 8: i_{b}q_{b} = 1,200 D^{3/2} i_{b} (φ_{*}/2) =
[Col. 6] [Col 7] [Col.5]
Column 9: 43.2 w
Column 10: i_{B}Q_{B} = (43.2 w) (i_{b}q_{b}) = [Col. 9] [Col. 8]
Column 11: percentage of suspended material Q'_{s} in each size fraction (input array data)
Column 12: [Col. 11] (Q'_{s} Total)
Column 13: Multipliers, or ratios of the 0.7 power of the fall velocities with the reference size.
The reference size is chosen as one of the suspended size fractions containing a substantial amount [the larger percentage] of the sampled load.
By definition, the multiplier for the reference size is 1.
Column 14: First, compute Z' for the reference size by trial and error.
Given Q'_{s} / i_{B}Q_{B}, use Plate 8 to find an initial value of Z'.
LHS = Q'_{s} / i_{B}Q_{B}
Calculate RHS = (I"_{1} / J"_{1}) (P J'_{1} + J'_{2})
With Z' and A" (Col. 15), use Plate 9 to find I"_{1}
With Z' and A" (Col. 15), use Plate 10 to find J"_{1}
With Z' and A', use Plate 10 to find J'_{1}
With Z' and A', use Plate 11 to find J'_{2}
Iterate on Z' until LHS = RHS.
Finally, use multipliers in Col. 13 to calculate Z' for other size fractions, i.e., Col. 14.
Column 15: A" = 2D / d = 2 (Col. 3) / d
Column 16: With Z' and A', use Plate 10 to find J'_{1}
Column 17: With Z' and A', use Plate 11 to find J'_{2}
Column 18: With Z' and A", use Plate 10 to find J"_{1}
Column 19: With Z' and A", use Plate 11 to find J"_{2}
Column 20: [ P (Col. 18)  (Col. 19) ] [ P (Col. 16)  (Col. 17) ]
Column 21: With Z' and A", use Plate 9 to find I"_{1}
Column 22: With Z' and A", use Plate 12 to find I"_{2}
Column 23: [ P (Col. 21)  (Col. 22) + 1 ]
Column 24: Either
a. For finer particle sizes: (Col. 12) (Col. 20).
or
b. For coarser particle sizes: (Col. 10) (Col. 23).
Column 25: [For SI units only] Computed load in Metric Tons/day = (Col. 24) * 0.9072
