Calculation of total sediment discharge by the Modified Einstein Procedure (USBR method), Victor Miguel Ponce, San Diego State University

Calculation of total sediment discharge by the Modified Einstein Procedure (USBR method), Victor Miguel Ponce, San Diego State University =$x[$n]) { $k = $n-1; $yl = (($y[$n]-$y[$k])*$xl + ($y[$k]*$x[$n] - $y[$n]*$x[$k]))/ ($x[$n] - $x[$k]); } else { while($j<$n) { $k = $j+1; if($xl>=$x[$j] and $xl<$x[$k]) { $m = $k; } $j++; } $l = $m - 1.; $yl = (($y[$m]-$y[$l])*$xl + ($y[$l]*$x[$m] - $y[$m]*$x[$l]))/ ($x[$m] - $x[$l]); } return($yl); } //end function pwlin0//////////////////////////////////////////////////////////////// //beginning of function pwlin0logx/////////////////////////////////////////////////// //developed 070514 //this function defines two paired arrays, x and y, of size n [0 to n], and a //logarithmically interpolated location xl in the x axis, and returns a linearly //interpolated location yl in the y axis. //the first value of x cannot be zero (undefined) function pwlin0logx($x,$y,$n,$xl,$yl) { $vl = log10($xl); $i= 0; while($i<=$n) { $v[$i] = log10($x[$i]); $i++; } $j= 0; if($vl<$v[$j]) { $k = $j + 1; $yl = (($y[$k]-$y[$j])*$vl + ($y[$j]*$v[$k] - $y[$k]*$v[$j]))/ ($v[$k] - $v[$j]); } elseif($vl>=$v[$n]) { $k = $n - 1; $yl = (($y[$n]-$y[$k])*$vl + ($y[$k]*$v[$n] - $y[$n]*$v[$k]))/ ($v[$n] - $v[$k]); } else { while($j<$n) { $k = $j + 1; if($vl>=$v[$j] and $vl<$v[$k]) { $m = $k; } $j++; } $l = $m - 1; $yl = (($y[$m]-$y[$l])*$vl + ($y[$l]*$v[$m] - $y[$m]*$v[$l]))/ ($v[$m] - $v[$l]); } return($yl); } //end of function pwlin0logx/////////////////////////////////////////////////////////// //beginning of function pwlin0logy/////////////////////////////////////////////////// //developed 070514 //this function defines two paired arrays, x and y, of size n [0 to n], and a //linearly interpolated location xl in the x axis, and returns a logarithmically //interpolated location yl in the y axis. //the first value of x cannot be zero (undefined) function pwlin0logy($x,$y,$n,$xl,$yl) { $i= 0; while($i<=$n) { $w[$i] = log10($y[$i]); $i++; } $j= 0; if($xl<$x[$j]) { $k = $j + 1; $wl = (($w[$k]-$w[$j])*$xl + ($w[$j]*$x[$k] - $w[$k]*$x[$j]))/ ($x[$k] - $x[$j]); } elseif($xl>=$x[$n]) { $k = $n - 1; $wl = (($w[$n]-$w[$k])*$xl + ($w[$k]*$x[$n] - $w[$n]*$x[$k]))/ ($x[$n] - $x[$k]); } else { while($j<$n) { $k = $j + 1; if($xl>=$x[$j] and $xl<$x[$k]) { $m = $k; } $j++; } $l = $m - 1; $wl = (($w[$m]-$w[$l])*$xl + ($w[$l]*$x[$m] - $w[$m]*$x[$l]))/ ($x[$m] - $x[$l]); } $yl = pow(10,$wl); return($yl); } //end of function pwlin0logy/////////////////////////////////////////////////////////// //beginning of function pwlin0logxlogy///////////////////////////////////////////////// //developed 070514 //this function defines two paired arrays, x and y, of size n [0 to n], and a //logarithmically interpolated location xl in the x axis, and returns a //logarithmically interpolated location yl in the y axis. //the first value of x cannot be zero (undefined) function pwlin0logxlogy($x,$y,$n,$xl,$yl) { $vl = log10($xl); $wl = log10($yl); $i= 0; while($i<=$n) { $v[$i] = log10($x[$i]); $w[$i] = log10($y[$i]); $i++; } $j= 0; if($vl<$v[$j]) { $k = $j + 1; $wl = (($w[$k]-$w[$j])*$vl + ($w[$j]*$v[$k] - $w[$k]*$v[$j]))/ ($v[$k] - $v[$j]); } elseif($vl>=$v[$n]) { $k = $n - 1; $wl = (($w[$n]-$w[$k])*$vl + ($w[$k]*$v[$n] - $w[$n]*$v[$k]))/ ($v[$n] - $v[$k]); } else { while($j<$n) { $k = $j + 1; if($vl>=$v[$j] and $vl<$v[$k]) { $m = $k; } $j++; } $l = $m - 1; $wl = (($w[$m]-$w[$l])*$vl + ($w[$l]*$v[$m] - $w[$m]*$v[$l]))/ ($v[$m] - $v[$l]); } $yl = pow(10,$wl); return($yl); } //end of function pwlin0logxlogy////////////////////////////////////////////////////// //beginning of function plate09/////////////////////////////////////////////////////// function modeinsplate09($x,$z,$y) { $ncplate09 = 11; $naplate09 = 7; $curplate09[0] = 0.; $curplate09[1] = 0.2; $curplate09[2] = 0.4; $curplate09[3] = 0.6; $curplate09[4] = 0.8; $curplate09[5] = 1.0; $curplate09[6] = 1.2; $curplate09[7] = 1.5; $curplate09[8] = 2.0; $curplate09[9] = 3.0; $curplate09[10] = 4.0; $curplate09[11] = 5.0; $absplate09[0] = 0.00003; $absplate09[1] = 0.0001; $absplate09[2] = 0.001; $absplate09[3] = 0.01; $absplate09[4] = 0.02; $absplate09[5] = 0.04; $absplate09[6] = 0.07; $absplate09[7] = 0.1; $ord2plate09[0][0] = 7400.; $ord2plate09[1][0] = 970.; $ord2plate09[2][0] = 148.; $ord2plate09[3][0] = 27.5; $ord2plate09[4][0] = 6.4; $ord2plate09[5][0] = 2.04; $ord2plate09[6][0] = 0.91; $ord2plate09[7][0] = 0.43; $ord2plate09[8][0] = 0.22; $ord2plate09[9][0] = 0.11; $ord2plate09[10][0] = 0.072; $ord2plate09[11][0] = 0.054; $ord2plate09[0][1] = 2200.; $ord2plate09[1][1] = 370.; $ord2plate09[2][1] = 72.; $ord2plate09[3][1] = 17.; $ord2plate09[4][1] = 4.8; $ord2plate09[5][1] = 1.8; $ord2plate09[6][1] = 0.87; $ord2plate09[7][1] = 0.42; $ord2plate09[8][1] = 0.22; $ord2plate09[9][1] = 0.11; $ord2plate09[10][1] = 0.071; $ord2plate09[11][1] = 0.054; $ord2plate09[0][2] = 220.; $ord2plate09[1][2] = 58.; $ord2plate09[2][2] = 18.; $ord2plate09[3][2] = 6.3; $ord2plate09[4][2] = 2.6; $ord2plate09[5][2] = 1.3; $ord2plate09[6][2] = 0.73; $ord2plate09[7][2] = 0.4; $ord2plate09[8][2] = 0.215; $ord2plate09[9][2] = 0.11; $ord2plate09[10][2] = 0.071 ; $ord2plate09[11][2] = 0.054; $ord2plate09[0][3] = 24.; $ord2plate09[1][3] = 9.; $ord2plate09[2][3] = 4.2; $ord2plate09[3][3] = 2.2; $ord2plate09[4][3] = 1.26; $ord2plate09[5][3] = 0.79; $ord2plate09[6][3] = 0.542; $ord2plate09[7][3] = 0.348; $ord2plate09[8][3] = 0.202; $ord2plate09[9][3] = 0.11; $ord2plate09[10][3] = 0.07 ; $ord2plate09[11][3] = 0.0535; $ord2plate09[0][4] = 10.8; $ord2plate09[1][4] = 5.1; $ord2plate09[2][4] = 2.7; $ord2plate09[3][4] = 1.55; $ord2plate09[4][4] = 0.97; $ord2plate09[5][4] = 0.65; $ord2plate09[6][4] = 0.47; $ord2plate09[7][4] = 0.318; $ord2plate09[8][4] = 0.19; $ord2plate09[9][4] = 0.103; $ord2plate09[10][4] = 0.069 ; $ord2plate09[11][4] = 0.0525; $ord2plate09[0][5] = 5.2; $ord2plate09[1][5] = 2.8; $ord2plate09[2][5] = 1.66; $ord2plate09[3][5] = 1.05; $ord2plate09[4][5] = 0.71; $ord2plate09[5][5] = 0.52; $ord2plate09[6][5] = 0.39; $ord2plate09[7][5] = 0.275; $ord2plate09[8][5] = 0.174; $ord2plate09[9][5] = 0.098; $ord2plate09[10][5] = 0.0655; $ord2plate09[11][5] = 0.05; $ord2plate09[0][6] = 2.9; $ord2plate09[1][6] = 1.7; $ord2plate09[2][6] = 1.1; $ord2plate09[3][6] = 0.75; $ord2plate09[4][6] = 0.535; $ord2plate09[5][6] = 0.41; $ord2plate09[6][6] = 0.32; $ord2plate09[7][6] = 0.24; $ord2plate09[8][6] = 0.156; $ord2plate09[9][6] = 0.09; $ord2plate09[10][6] = 0.062; $ord2plate09[11][6] = 0.0477; $ord2plate09[0][7] = 1.95; $ord2plate09[1][7] = 1.2; $ord2plate09[2][7] = 0.84; $ord2plate09[3][7] = 0.58; $ord2plate09[4][7] = 0.43; $ord2plate09[5][7] = 0.34; $ord2plate09[6][7] = 0.27; $ord2plate09[7][7] = 0.21; $ord2plate09[8][7] = 0.14; $ord2plate09[9][7] = 0.085; $ord2plate09[10][7] = 0.058; $ord2plate09[11][7] = 0.0457; $iflag = 0; //$x = 0.00118; //$z = 0.8; if($z == 0) { $zn = 0; $iflag = 1; } elseif($z > 0 && $z < 0.2) { $znt = 0; $znb = 1; $iflag = 2; } elseif($z == 0.2) { $zn = 1; $iflag = 1; } elseif($z > 0.2 && $z < 0.4) { $znt = 1; $znb = 2; $iflag = 2; } elseif($z == 0.4) { $zn = 2; $iflag = 1; } elseif($z > 0.4 && $z < 0.6) { $znt = 2; $znb = 3; $iflag = 2; } elseif($z == 0.6) { $zn = 3; $iflag = 1; } elseif($z > 0.6 && $z < 0.8) { $znt = 3; $znb = 4; $iflag = 2; } elseif($z == 0.8) { $zn = 4; $iflag = 1; } elseif($z > 0.8 && $z < 1.0) { $znt = 4; $znb = 5; $iflag = 2; } elseif($z == 1.0) { $zn = 5; $iflag = 1; } elseif($z > 1.0 && $z < 1.2) { $znt = 5; $znb = 6; $iflag = 2; } elseif($z == 1.2) { $zn = 6; $iflag = 1; } elseif($z > 1.2 && $z < 1.5) { $znt = 6; $znb = 7; $iflag = 2; } elseif($z == 1.5) { $zn = 7; $iflag = 1; } elseif($z > 1.5 && $z < 2.0) { $znt = 7; $znb = 8; $iflag = 2; } elseif($z == 2.0) { $zn = 8; $iflag = 1; } elseif($z > 2.0 && $z < 3.0) { $znt = 8; $znb = 9; $iflag = 2; } elseif($z == 3.0) { $zn = 9; $iflag = 1; } elseif($z > 3.0 && $z < 4.0) { $znt = 9; $znb = 10; $iflag = 2; } elseif($z == 4.0) { $zn = 10; $iflag = 1; } elseif($z > 4.0 && $z < 5.0) { $znt = 10; $znb = 11; $iflag = 2; } elseif($z == 5.0) { $zn = 11; $iflag = 1; } while($iflag == 1) { $i = 0.; while($i <= $naplate09) { $ordplate09[$i] = $ord2plate09[$zn][$i]; $i++; } $y = pwlin0logxlogy($absplate09,$ordplate09,$naplate09,$x,$y); //echo $x . "sofarsogood 5 " . $y . "
"; $iflag = 0; } while($iflag == 2) { $i = 0.; while($i <= $naplate09) { $ordtplate09[$i] = $ord2plate09[$znt][$i]; $i++; } $i = 0.; while($i <= $naplate09) { $ordbplate09[$i] = $ord2plate09[$znb][$i]; $i++; } $yt = pwlin0logxlogy($absplate09,$ordtplate09,$naplate09,$x,$yt); $yb = pwlin0logxlogy($absplate09,$ordbplate09,$naplate09,$x,$yb); //echo $x . "sofarsogood 7 " . " " . $yt . " " . $yb . "
"; $tempcurvex[0] = $curplate09[$znt]; $tempcurvex[1] = $curplate09[$znb]; //echo $tempcurvex[0] . "sofarsogood 9 " . " " . $tempcurvex[1] . "
"; $tempcurvey[0] = $yt; $tempcurvey[1] = $yb; $y = pwlin0logy($tempcurvex,$tempcurvey,1,$z,$y); //echo $x . "sofarsogood 7 " . " " . $yt . " " . $yb . " " . $y . "
"; $iflag = 0; } //$i1doubleprime = $y; //echo $i1doubleprime . "
"; //echo "so far so good end plate 09 i1doubleprime function modeinsplate09" . "
"; return($y); } //end of function plate09///////////////////////////////////////////////////////////// //beginning of function plate10/////////////////////////////////////////////////////// function modeinsplate10($x,$z,$y) { $ncplate10 = 9; $naplate10 = 9; $curplate10[0] = 0.; $curplate10[1] = 0.2; $curplate10[2] = 0.4; $curplate10[3] = 0.6; $curplate10[4] = 0.8; $curplate10[5] = 1.0; $curplate10[6] = 1.2; $curplate10[7] = 1.5; $curplate10[8] = 1.75; $curplate10[9] = 2.0; $absplate10[0] = 0.00001; $absplate10[1] = 0.0001; $absplate10[2] = 0.001; $absplate10[3] = 0.01; $absplate10[4] = 0.04; $absplate10[5] = 0.1; $absplate10[6] = 0.2; $absplate10[7] = 0.3; $absplate10[8] = 0.4; $absplate10[9] = 0.5; $ord2plate10[0][0] = 1.; $ord2plate10[1][0] = 1.1; $ord2plate10[2][0] = 1.36; $ord2plate10[3][0] = 1.96; $ord2plate10[4][0] = 3.80; $ord2plate10[5][0] = 11.; $ord2plate10[6][0] = 44.; $ord2plate10[7][0] = 630.; $ord2plate10[8][0] = 7000.; $ord2plate10[9][0] = 100000.; $ord2plate10[0][1] = 1.; $ord2plate10[1][1] = 1.08; $ord2plate10[2][1] = 1.32; $ord2plate10[3][1] = 1.92; $ord2plate10[4][1] = 3.5; $ord2plate10[5][1] = 8.2; $ord2plate10[6][1] = 25.5; $ord2plate10[7][1] = 95.; $ord2plate10[8][1] = 1300.; $ord2plate10[9][1] = 10000.; $ord2plate10[0][2] = 0.99; $ord2plate10[1][2] = 1.05; $ord2plate10[2][2] = 1.3; $ord2plate10[3][2] = 1.8; $ord2plate10[4][2] = 3.04; $ord2plate10[5][2] = 5.95; $ord2plate10[6][2] = 13.7; $ord2plate10[7][2] = 59.; $ord2plate10[8][2] = 230.; $ord2plate10[9][2] = 1000.; $ord2plate10[0][3] = 0.98; $ord2plate10[1][3] = 1.02; $ord2plate10[2][3] = 1.2; $ord2plate10[3][3] = 1.6; $ord2plate10[4][3] = 2.3; $ord2plate10[5][3] = 3.6; $ord2plate10[6][3] = 6.2; $ord2plate10[7][3] = 15.6; $ord2plate10[8][3] = 36.8; $ord2plate10[9][3] = 92.; $ord2plate10[0][4] = 0.96; $ord2plate10[1][4] = 0.975; $ord2plate10[2][4] = 1.1; $ord2plate10[3][4] = 1.3; $ord2plate10[4][4] = 1.68; $ord2plate10[5][4] = 2.26; $ord2plate10[6][4] = 3.2; $ord2plate10[7][4] = 5.9; $ord2plate10[8][4] = 10.5; $ord2plate10[9][4] = 18.7; $ord2plate10[0][5] = 0.865; $ord2plate10[1][5] = 0.87; $ord2plate10[2][5] = 0.92; $ord2plate10[3][5] = 1.; $ord2plate10[4][5] = 1.16; $ord2plate10[5][5] = 1.4; $ord2plate10[6][5] = 1.76; $ord2plate10[7][5] = 2.55; $ord2plate10[8][5] = 3.75; $ord2plate10[9][5] = 5.4; $ord2plate10[0][6] = 0.8; $ord2plate10[1][6] = 0.72; $ord2plate10[2][6] = 0.71; $ord2plate10[3][6] = 0.72; $ord2plate10[4][6] = 0.74; $ord2plate10[5][6] = 0.8; $ord2plate10[6][6] = 0.91; $ord2plate10[7][6] = 1.1; $ord2plate10[8][6] = 1.26; $ord2plate10[9][6] = 1.56; $ord2plate10[0][7] = 0.7; $ord2plate10[1][7] = 0.6; $ord2plate10[2][7] = 0.55; $ord2plate10[3][7] = 0.55; $ord2plate10[4][7] = 0.55; $ord2plate10[5][7] = 0.58; $ord2plate10[6][7] = 0.60; $ord2plate10[7][7] = 0.60; $ord2plate10[8][7] = 0.61; $ord2plate10[9][7] = 0.62; $ord2plate10[0][8] = 0.6; $ord2plate10[1][8] = 0.48; $ord2plate10[2][8] = 0.4; $ord2plate10[3][8] = 0.37; $ord2plate10[4][8] = 0.34; $ord2plate10[5][8] = 0.31; $ord2plate10[6][8] = 0.29; $ord2plate10[7][8] = 0.28; $ord2plate10[8][8] = 0.27; $ord2plate10[9][8] = 0.26; $ord2plate10[0][9] = 0.5; $ord2plate10[1][9] = 0.38; $ord2plate10[2][9] = 0.31; $ord2plate10[3][9] = 0.26; $ord2plate10[4][9] = 0.22; $ord2plate10[5][9] = 0.19; $ord2plate10[6][9] = 0.165; $ord2plate10[7][9] = 0.14; $ord2plate10[8][9] = 0.128; $ord2plate10[9][9] = 0.115; $iflag = 0; if($z == 0) { $zn = 0; $iflag = 1; } elseif($z > 0 && $z < 0.2) { $znb = 0; $znt = 1; $iflag = 2; } elseif($z == 0.2) { $zn = 1; $iflag = 1; } elseif($z > 0.2 && $z < 0.4) { $znb = 1; $znt = 2; $iflag = 2; } elseif($z == 0.4) { $zn = 2; $iflag = 1; } elseif($z > 0.4 && $z < 0.6) { $znb = 2; $znt = 3; $iflag = 2; } elseif($z == 0.6) { $zn = 3; $iflag = 1; } elseif($z > 0.6 && $z < 0.8) { $znb = 3; $znt = 4; $iflag = 2; } elseif($z == 0.8) { $zn = 4; $iflag = 1; } elseif($z > 0.8 && $z < 1.0) { $znb = 4; $znt = 5; $iflag = 2; } elseif($z == 1.0) { $zn = 5; $iflag = 1; } elseif($z > 1.0 && $z < 1.2) { $znb = 5; $znt = 6; $iflag = 2; } elseif($z == 1.2) { $zn = 6; $iflag = 1; } elseif($z > 1.2 && $z < 1.5) { $znb = 6; $znt = 7; $iflag = 2; } elseif($z == 1.5) { $zn = 7; $iflag = 1; } elseif($z > 1.5 && $z < 1.75) { $znb = 7; $znt = 8; $iflag = 2; } elseif($z == 1.75) { $zn = 8; $iflag = 1; } elseif($z > 1.75 && $z < 2.0) { $znb = 8; $znt = 9; $iflag = 2; } elseif($z == 2.0) { $zn = 9; $iflag = 1; } while($iflag == 1) { $i = 0.; while($i <= $naplate10) { $ordplate10[$i] = $ord2plate10[$zn][$i]; $i++; } $y = pwlin0logxlogy($absplate10,$ordplate10,$naplate10,$x,$y); //echo $x . " sofarsogood 11 " . $y . "
"; $iflag = 0; } while($iflag == 2) { $i = 0.; while($i <= $naplate10) { $ordbplate10[$i] = $ord2plate10[$znb][$i]; $i++; } $i = 0.; while($i <= $naplate10) { $ordtplate10[$i] = $ord2plate10[$znt][$i]; $i++; } $yb = pwlin0logxlogy($absplate10,$ordbplate10,$naplate10,$x,$yb); $yt = pwlin0logxlogy($absplate10,$ordtplate10,$naplate10,$x,$yt); //echo $x . " sofarsogood 12 " . " " . $yb . " " . $yt . "
"; $tempcurvex[0] = $curplate10[$znb]; $tempcurvex[1] = $curplate10[$znt]; //echo $tempcurvex[0] . "sofarsogood 13 " . " " . $tempcurvex[1] . "
"; $tempcurvey[0] = $yb; $tempcurvey[1] = $yt; $y = pwlin0logy($tempcurvex,$tempcurvey,1,$z,$y); //echo $z . " sofarsogood 13 " . " " . $yb . " " . $yt . " " . $y . "
"; $iflag = 0; } $j1doubleprime = $y; //echo $j1doubleprime . "
"; //echo "so far so good end plate 10 j1doubleprime" . "
"; //echo "so far so good end plate 10 j1doubleprime function modeinsplate10" . "
"; return($y); } //end of function plate10///////////////////////////////////////////////////////////// //beginning of function plate11/////////////////////////////////////////////////////// function modeinsplate11($x,$z,$y) { $ncplate11 = 9; $naplate11 = 9; $curplate11[0] = 0.; $curplate11[1] = 0.2; $curplate11[2] = 0.4; $curplate11[3] = 0.6; $curplate11[4] = 0.8; $curplate11[5] = 1.0; $curplate11[6] = 1.2; $curplate11[7] = 1.5; $curplate11[8] = 1.75; $curplate11[9] = 2.0; $absplate11[0] = 0.00001; $absplate11[1] = 0.0001; $absplate11[2] = 0.001; $absplate11[3] = 0.01; $absplate11[4] = 0.04; $absplate11[5] = 0.1; $absplate11[6] = 0.2; $absplate11[7] = 0.3; $absplate11[8] = 0.4; $absplate11[9] = 0.5; $ord2plate11[0][0] = 1.; $ord2plate11[1][0] = 1.5; $ord2plate11[2][0] = 2.6; $ord2plate11[3][0] = 5.6; $ord2plate11[4][0] = 16.; $ord2plate11[5][0] = 66.; $ord2plate11[6][0] = 360.; $ord2plate11[7][0] = 6200.; $ord2plate11[8][0] = 72000.; $ord2plate11[9][0] = 1000000.; $ord2plate11[0][1] = 1.; $ord2plate11[1][1] = 1.5; $ord2plate11[2][1] = 2.5; $ord2plate11[3][1] = 5.2; $ord2plate11[4][1] = 13.2; $ord2plate11[5][1] = 42.; $ord2plate11[6][1] = 154.; $ord2plate11[7][1] = 1450.; $ord2plate11[8][1] = 10000.; $ord2plate11[9][1] = 84000.; $ord2plate11[0][2] = 0.995; $ord2plate11[1][2] = 1.46; $ord2plate11[2][2] = 2.35; $ord2plate11[3][2] = 4.3; $ord2plate11[4][2] = 9.6; $ord2plate11[5][2] = 23.; $ord2plate11[6][2] = 62.; $ord2plate11[7][2] = 310.; $ord2plate11[8][2] = 1300.; $ord2plate11[9][2] = 5900.; $ord2plate11[0][3] = 0.94; $ord2plate11[1][3] = 1.3; $ord2plate11[2][3] = 1.93; $ord2plate11[3][3] = 3.2; $ord2plate11[4][3] = 5.4; $ord2plate11[5][3] = 9.7; $ord2plate11[6][3] = 18.8; $ord2plate11[7][3] = 53.; $ord2plate11[8][3] = 130.; $ord2plate11[9][3] = 340.; $ord2plate11[0][4] = 0.82; $ord2plate11[1][4] = 1.1; $ord2plate11[2][4] = 1.4; $ord2plate11[3][4] = 2.; $ord2plate11[4][4] = 2.9; $ord2plate11[5][4] = 4.35; $ord2plate11[6][4] = 6.8; $ord2plate11[7][4] = 13.8; $ord2plate11[8][4] = 24.5; $ord2plate11[9][4] = 48.; $ord2plate11[0][5] = 0.67; $ord2plate11[1][5] = 0.78; $ord2plate11[2][5] = 0.95; $ord2plate11[3][5] = 1.18; $ord2plate11[4][5] = 1.52; $ord2plate11[5][5] = 2.; $ord2plate11[6][5] = 2.65; $ord2plate11[7][5] = 4.2; $ord2plate11[8][5] = 6.3; $ord2plate11[9][5] = 9.5; $ord2plate11[0][6] = 0.475; $ord2plate11[1][6] = 0.51; $ord2plate11[2][6] = 0.56; $ord2plate11[3][6] = 0.63; $ord2plate11[4][6] = 0.71; $ord2plate11[5][6] = 0.81; $ord2plate11[6][6] = 0.94; $ord2plate11[7][6] = 1.125; $ord2plate11[8][6] = 1.5; $ord2plate11[9][6] = 1.9; $ord2plate11[0][7] = 0.33; $ord2plate11[1][7] = 0.34; $ord2plate11[2][7] = 0.35; $ord2plate11[3][7] = 0.37; $ord2plate11[4][7] = 0.38; $ord2plate11[5][7] = 0.39; $ord2plate11[6][7] = 0.41; $ord2plate11[7][7] = 0.45; $ord2plate11[8][7] = 0.49; $ord2plate11[9][7] = 0.57; $ord2plate11[0][8] = 0.23; $ord2plate11[1][8] = 0.215; $ord2plate11[2][8] = 0.21; $ord2plate11[3][8] = 0.205; $ord2plate11[4][8] = 0.20; $ord2plate11[5][8] = 0.195; $ord2plate11[6][8] = 0.19; $ord2plate11[7][8] = 0.187; $ord2plate11[8][8] = 0.183; $ord2plate11[9][8] = 0.18; $ord2plate11[0][9] = 0.16; $ord2plate11[1][9] = 0.135; $ord2plate11[2][9] = 0.125; $ord2plate11[3][9] = 0.12; $ord2plate11[4][9] = 0.115; $ord2plate11[5][9] = 0.11; $ord2plate11[6][9] = 0.10; $ord2plate11[7][9] = 0.09; $ord2plate11[8][9] = 0.08; $ord2plate11[9][9] = 0.07; $iflag = 0; if($z == 0) { $zn = 0; $iflag = 1; } elseif($z > 0 && $z < 0.2) { $znb = 0; $znt = 1; $iflag = 2; } elseif($z == 0.2) { $zn = 1; $iflag = 1; } elseif($z > 0.2 && $z < 0.4) { $znb = 1; $znt = 2; $iflag = 2; } elseif($z == 0.4) { $zn = 2; $iflag = 1; } elseif($z > 0.4 && $z < 0.6) { $znb = 2; $znt = 3; $iflag = 2; } elseif($z == 0.6) { $zn = 3; $iflag = 1; } elseif($z > 0.6 && $z < 0.8) { $znb = 3; $znt = 4; $iflag = 2; } elseif($z == 0.8) { $zn = 4; $iflag = 1; } elseif($z > 0.8 && $z < 1.0) { $znb = 4; $znt = 5; $iflag = 2; } elseif($z == 1.0) { $zn = 5; $iflag = 1; } elseif($z > 1.0 && $z < 1.2) { $znb = 5; $znt = 6; $iflag = 2; } elseif($z == 1.2) { $zn = 6; $iflag = 1; } elseif($z > 1.2 && $z < 1.5) { $znb = 6; $znt = 7; $iflag = 2; } elseif($z == 1.5) { $zn = 7; $iflag = 1; } elseif($z > 1.5 && $z < 1.75) { $znb = 7; $znt = 8; $iflag = 2; } elseif($z == 1.75) { $zn = 8; $iflag = 1; } elseif($z > 1.75 && $z < 2.0) { $znb = 8; $znt = 9; $iflag = 2; } elseif($z == 2.0) { $zn = 9; $iflag = 1; } while($iflag == 1) { $i = 0.; while($i <= $naplate11) { $ordplate11[$i] = $ord2plate11[$zn][$i]; $i++; } $y = pwlin0logxlogy($absplate11,$ordplate11,$naplate11,$x,$y); //echo $x . " sofarsogood 15 " . $y . "
"; $iflag = 0; } while($iflag == 2) { $i = 0.; while($i <= $naplate11) { $ordbplate11[$i] = $ord2plate11[$znb][$i]; $i++; } $i = 0.; while($i <= $naplate11) { $ordtplate11[$i] = $ord2plate11[$znt][$i]; $i++; } $yb = pwlin0logxlogy($absplate11,$ordbplate11,$naplate11,$x,$yb); $yt = pwlin0logxlogy($absplate11,$ordtplate11,$naplate11,$x,$yt); //echo $x . " sofarsogood 12 " . " " . $yb . " " . $yt . "
"; $tempcurvex[0] = $curplate11[$znb]; $tempcurvex[1] = $curplate11[$znt]; //echo $tempcurvex[0] . "sofarsogood 13 " . " " . $tempcurvex[1] . "
"; $tempcurvey[0] = $yb; $tempcurvey[1] = $yt; $y = pwlin0logy($tempcurvex,$tempcurvey,1,$z,$y); //echo $z . " sofarsogood 16 " . " " . $yb . " " . $yt . " " . $y . "
"; $iflag = 0; } $j2prime = $y; //echo $j2prime . "
"; //echo "so far so good end plate 11 j2prime" . "
"; //echo "so far so good end plate 11 j2prime function modeinsplate11" . "
"; return($y); } //end of function plate11///////////////////////////////////////////////////////////// //beginning of function plate12/////////////////////////////////////////////////////// function modeinsplate12($x,$z,$y) { $ncplate12 = 11; $naplate12 = 7; $curplate12[0] = 0.; $curplate12[1] = 0.2; $curplate12[2] = 0.4; $curplate12[3] = 0.6; $curplate12[4] = 0.8; $curplate12[5] = 1.0; $curplate12[6] = 1.2; $curplate12[7] = 1.5; $curplate12[8] = 2.0; $curplate12[9] = 3.0; $curplate12[10] = 4.0; $curplate12[11] = 5.0; $absplate12[0] = 0.00003; $absplate12[1] = 0.0001; $absplate12[2] = 0.001; $absplate12[3] = 0.01; $absplate12[4] = 0.02; $absplate12[5] = 0.04; $absplate12[6] = 0.07; $absplate12[7] = 0.1; $ord2plate12[0][0] = 7400.; $ord2plate12[1][0] = 1330.; $ord2plate12[2][0] = 290.; $ord2plate12[3][0] = 75.; $ord2plate12[4][0] = 25.6; $ord2plate12[5][0] = 11.6; $ord2plate12[6][0] = 6.45; $ord2plate12[7][0] = 3.65; $ord2plate12[8][0] = 2.05; $ord2plate12[9][0] = 1.08; $ord2plate12[10][0] = 0.735; $ord2plate12[11][0] = 0.555; $ord2plate12[0][1] = 2170.; $ord2plate12[1][1] = 500.; $ord2plate12[2][1] = 138.; $ord2plate12[3][1] = 44.; $ord2plate12[4][1] = 18; $ord2plate12[5][1] = 9.; $ord2plate12[6][1] = 5.3; $ord2plate12[7][1] = 3.15; $ord2plate12[8][1] = 1.8; $ord2plate12[9][1] = 0.95; $ord2plate12[10][1] = 0.645; $ord2plate12[11][1] = 0.49; $ord2plate12[0][2] = 218.; $ord2plate12[1][2] = 77.5; $ord2plate12[2][2] = 32.5; $ord2plate12[3][2] = 15.3; $ord2plate12[4][2] = 8.1; $ord2plate12[5][2] = 5.; $ord2plate12[6][2] = 3.29; $ord2plate12[7][2] = 2.09; $ord2plate12[8][2] = 1.028; $ord2plate12[9][2] = 0.69; $ord2plate12[10][2] = 0.475; $ord2plate12[11][2] = 0.36; $ord2plate12[0][3] = 20.5; $ord2plate12[1][3] = 11.; $ord2plate12[2][3] = 6.6; $ord2plate12[3][3] = 4.2; $ord2plate12[4][3] = 2.9; $ord2plate12[5][3] = 2.1; $ord2plate12[6][3] = 1.6; $ord2plate12[7][3] = 1.15; $ord2plate12[8][3] = 0.75; $ord2plate12[9][3] = 0.44; $ord2plate12[10][3] = 0.303; $ord2plate12[11][3] = 0.232; $ord2plate12[0][4] = 9.6; $ord2plate12[1][4] = 5.9; $ord2plate12[2][4] = 3.9; $ord2plate12[3][4] = 2.7; $ord2plate12[4][4] = 1.97; $ord2plate12[5][4] = 1.5; $ord2plate12[6][4] = 1.175; $ord2plate12[7][4] = 0.86; $ord2plate12[8][4] = 0.59; $ord2plate12[9][4] = 0.35; $ord2plate12[10][4] = 0.25; $ord2plate12[11][4] = 0.19; $ord2plate12[0][5] = 4.4; $ord2plate12[1][5] = 3.0; $ord2plate12[2][5] = 2.15; $ord2plate12[3][5] = 1.6; $ord2plate12[4][5] = 1.24; $ord2plate12[5][5] = 0.98; $ord2plate12[6][5] = 0.805; $ord2plate12[7][5] = 0.61; $ord2plate12[8][5] = 0.435; $ord2plate12[9][5] = 0.268; $ord2plate12[10][5] = 0.19; $ord2plate12[11][5] = 0.149; $ord2plate12[0][6] = 2.25; $ord2plate12[1][6] = 1.63; $ord2plate12[2][6] = 1.25; $ord2plate12[3][6] = 0.98; $ord2plate12[4][6] = 0.8; $ord2plate12[5][6] = 0.65; $ord2plate12[6][6] = 0.55; $ord2plate12[7][6] = 0.435; $ord2plate12[8][6] = 0.34; $ord2plate12[9][6] = 0.208; $ord2plate12[10][6] = 0.145; $ord2plate12[11][6] = 0.115; $ord2plate12[0][7] = 1.48; $ord2plate12[1][7] = 1.08; $ord2plate12[2][7] = 0.85; $ord2plate12[3][7] = 0.7; $ord2plate12[4][7] = 0.56; $ord2plate12[5][7] = 0.475; $ord2plate12[6][7] = 0.412; $ord2plate12[7][7] = 0.337; $ord2plate12[8][7] = 0.25; $ord2plate12[9][7] = 0.165; $ord2plate12[10][7] = 0.117; $ord2plate12[11][7] = 0.094; $iflag = 0; if($z == 0) { $zn = 0; $iflag = 1; } elseif($z > 0 && $z < 0.2) { $znt = 0; $znb = 1; $iflag = 2; } elseif($z == 0.2) { $zn = 1; $iflag = 1; } elseif($z > 0.2 && $z < 0.4) { $znt = 1; $znb = 2; $iflag = 2; } elseif($z == 0.4) { $zn = 2; $iflag = 1; } elseif($z > 0.4 && $z < 0.6) { $znt = 2; $znb = 3; $iflag = 2; } elseif($z == 0.6) { $zn = 3; $iflag = 1; } elseif($z > 0.6 && $z < 0.8) { $znt = 3; $znb = 4; $iflag = 2; } elseif($z == 0.8) { $zn = 4; $iflag = 1; } elseif($z > 0.8 && $z < 1.0) { $znt = 4; $znb = 5; $iflag = 2; } elseif($z == 1.0) { $zn = 5; $iflag = 1; } elseif($z > 1.0 && $z < 1.2) { $znt = 5; $znb = 6; $iflag = 2; } elseif($z == 1.2) { $zn = 6; $iflag = 1; } elseif($z > 1.2 && $z < 1.5) { $znt = 6; $znb = 7; $iflag = 2; } elseif($z == 1.5) { $zn = 7; $iflag = 1; } elseif($z > 1.5 && $z < 2.0) { $znt = 7; $znb = 8; $iflag = 2; } elseif($z == 2.0) { $zn = 8; $iflag = 1; } elseif($z > 2.0 && $z < 3.0) { $znt = 8; $znb = 9; $iflag = 2; } elseif($z == 3.0) { $zn = 9; $iflag = 1; } elseif($z > 3.0 && $z < 4.0) { $znt = 9; $znb = 10; $iflag = 2; } elseif($z == 4.0) { $zn = 10; $iflag = 1; } elseif($z > 4.0 && $z < 5.0) { $znt = 10; $znb = 11; $iflag = 2; } elseif($z == 5.0) { $zn = 11; $iflag = 1; } while($iflag == 1) { $i = 0.; while($i <= $naplate12) { $ordplate12[$i] = $ord2plate12[$zn][$i]; $i++; } $y = pwlin0logxlogy($absplate12,$ordplate12,$naplate12,$x,$y); //echo $x . "sofarsogood 5 " . $y . "
"; $iflag = 0; } while($iflag == 2) { $i = 0.; while($i <= $naplate12) { $ordtplate12[$i] = $ord2plate12[$znt][$i]; $i++; } $i = 0.; while($i <= $naplate12) { $ordbplate12[$i] = $ord2plate12[$znb][$i]; $i++; } $yt = pwlin0logxlogy($absplate12,$ordtplate12,$naplate12,$x,$yt); $yb = pwlin0logxlogy($absplate12,$ordbplate12,$naplate12,$x,$yb); //echo $x . "sofarsogood 7 " . " " . $yt . " " . $yb . "
"; $tempcurvex[0] = $curplate12[$znt]; $tempcurvex[1] = $curplate12[$znb]; //echo $tempcurvex[0] . "sofarsogood 9 " . " " . $tempcurvex[1] . "
"; $tempcurvey[0] = $yt; $tempcurvey[1] = $yb; $y = pwlin0logy($tempcurvex,$tempcurvey,1,$z,$y); //echo $x . "sofarsogood 7 " . " " . $yt . " " . $yb . " " . $y . "
"; $iflag = 0; } //$i2doubleprime = $y; //echo $i2doubleprime . "
"; //echo "so far so good end plate 12 i2doubleprime function modeinsplate12" . "
"; return($y); } //end of function plate12///////////////////////////////////////////////////////////// ?>

"; echo " online_modified_einstein:
*** Error No. 1F (Fatal)
*** Please enter only positive number for discharge Q.
*** Execution stops here. "; echo "
"; die; } if($_POST["vein"] <= 0) { if($_POST["Submit"] == 'Calculate') echo "

"; echo " online_modified_einstein:
*** Error No. 2F (Fatal)
*** Please enter only positive number for mean velocity u.
*** Execution stops here. "; echo "
"; die; } if($_POST["dein"] <= 0) { if($_POST["Submit"] == 'Calculate') echo "

"; echo " online_modified_einstein:
*** Error No. 3F (Fatal)
*** Please enter only positive number for average depth d.
*** Execution stops here. "; echo "
"; die; } if($_POST["dein"] <= 0) { if($_POST["Submit"] == 'Calculate') echo "

"; echo " online_modified_einstein:
*** Error No. 4F (Fatal)
*** Please enter only positive number for average depth d.
*** Execution stops here. "; echo "
"; die; } if($_POST["wdin"] <= 0) { if($_POST["Submit"] == 'Calculate') echo "

"; echo " online_modified_einstein:
*** Error No. 5F (Fatal)
*** Please enter only positive number for width d.
*** Execution stops here. "; echo "
"; die; } if($_POST["d65in"] <= 0) { if($_POST["Submit"] == 'Calculate') echo "

"; echo " online_modified_einstein:
*** Error No. 6F (Fatal)
*** Please enter only positive number for D₆₅.
*** Execution stops here. "; echo "
"; die; } if($_POST["d35in"] <= 0) { if($_POST["Submit"] == 'Calculate') echo "

"; echo " online_modified_einstein:
*** Error No. 7F (Fatal)
*** Please enter only positive number for D₃₅.
*** Execution stops here. "; echo "
"; die; } if($_POST["coin"] <= 0) { if($_POST["Submit"] == 'Calculate') echo "

"; echo " online_modified_einstein:
*** Error No. 8F (Fatal)
*** Please enter only positive number for measured concentration C.
*** Execution stops here. "; echo "
"; die; } if($_POST["dnin"] <= 0) { if($_POST["Submit"] == 'Calculate') echo "

"; echo " online_modified_einstein:
*** Error No. 9F (Fatal)
*** Please enter only positive number for vertical distance not sampled d_n.
*** Execution stops here. "; echo "
"; die; } if($_POST["dsin"] <= 0) { if($_POST["Submit"] == 'Calculate') echo "

"; echo " online_modified_einstein:
*** Error No. 10F (Fatal)
*** Please enter only positive number for mean water depth at sampled verticals d_s .
*** Execution stops here. "; echo "
"; die; } if($unitsystems == 0) { if($_POST["tein"] < 32 || $_POST["tein"] > 104) { if($_POST["Submit"] == 'Calculate') echo "

"; echo " online_modified_einstein:
*** Error No. 11F (Fatal)
*** Please enter only values of temperature within allowable range: 32 ^oF <= temperature <= 104 ^oF.
*** Execution stops here. "; echo "
"; die; } } elseif($unitsystems == 1) { if($_POST["tein"] < 0 || $_POST["tein"] > 40) { if($_POST["Submit"] == 'Calculate') echo "

"; echo " online_modified_einstein:
*** Error No. 12F (Fatal)
*** Please enter only values of temperature within allowable range: 0 ^oC <= temperature <= 40 ^oC.
*** Execution stops here. "; echo "
"; die; } } //////////////////////////////////////////////////////////////////// //conversion from si to us units if($unitsystems == 1) { $qdin = $qdin/0.0283168; //echo "qdin = " . $qdin . "
"; $vein = $vein/0.3048; //echo "vein = " . $vein . "
"; $dein = $dein/0.3048; //echo "dein = " . $dein . "
"; $wdin = $wdin/0.3048; //echo "wdin = " . $wdin . "
"; $tein = $tein*1.8 + 32; //echo "tein = " . $tein . "
"; $d65in = $d65in*0.0032808; //echo "d65in = " . $d65in . "
"; $d35in = $d35in*0.0032808; //echo "d35in = " . $d35in . "
"; $dnin = $dnin/0.3048; //echo "dnin = " . $dnin . "
"; $dsin = $dsin/0.3048; //echo "dsin = " . $dsin . "
"; } /////////////////////////////////////////////////////////////////// //calculation of viscosity from temperature $numbertempvisc = 13; $temperaturetable[0] = 32.; $temperaturetable[1] = 40.; $temperaturetable[2] = 50.; $temperaturetable[3] = 60.; $temperaturetable[4] = 70.; $temperaturetable[5] = 80.; $temperaturetable[6] = 90.; $temperaturetable[7] = 100.; $temperaturetable[8] = 120.; $temperaturetable[9] = 140.; $temperaturetable[10] = 160.; $temperaturetable[11] = 180.; $temperaturetable[12] = 200.; $temperaturetable[13] = 212.; $viscositytable[0] = 0.00001931; $viscositytable[1] = 0.00001664; $viscositytable[2] = 0.00001410; $viscositytable[3] = 0.00001217; //linsley //$viscositytable[3] = 0.00001210; //rouse $viscositytable[4] = 0.00001058; //linsley //$viscositytable[4] = 0.00001050; //rouse $viscositytable[5] = 0.00000930; $viscositytable[6] = 0.00000826; $viscositytable[7] = 0.00000739; $viscositytable[8] = 0.00000609; $viscositytable[9] = 0.00000514; $viscositytable[10] = 0.00000442; $viscositytable[11] = 0.00000386; $viscositytable[12] = 0.00000341; $viscositytable[13] = 0.00000319; $viscosity = pwlin0($temperaturetable,$viscositytable,$numbertempvisc,$tein,$viscosity); //////////////////////////////////////////////////////////////////////////////////////// //plate 3 $numberksdx = 21; $ksovdeltatable[0] = 0.2; $ksovdeltatable[1] = 0.3; $ksovdeltatable[2] = 0.4; $ksovdeltatable[3] = 0.5; $ksovdeltatable[4] = 0.6; $ksovdeltatable[5] = 0.7; $ksovdeltatable[6] = 0.8; $ksovdeltatable[7] = 0.9; $ksovdeltatable[8] = 1.0; $ksovdeltatable[9] = 1.3; $ksovdeltatable[10] = 1.5; $ksovdeltatable[11] = 2.; $ksovdeltatable[12] = 2.5; $ksovdeltatable[13] = 3.; $ksovdeltatable[14] = 3.5; $ksovdeltatable[15] = 4.; $ksovdeltatable[16] = 5.; $ksovdeltatable[17] = 6.; $ksovdeltatable[18] = 7.; $ksovdeltatable[19] = 8.; $ksovdeltatable[20] = 9.; $ksovdeltatable[21] = 10.; $xtable[0] = 0.71; $xtable[1] = 1.02; $xtable[2] = 1.24; $xtable[3] = 1.39; $xtable[4] = 1.49; $xtable[5] = 1.55; $xtable[6] = 1.59; $xtable[7] = 1.608; $xtable[8] = 1.614; $xtable[9] = 1.58; $xtable[10] = 1.52; $xtable[11] = 1.38; $xtable[12] = 1.27; $xtable[13] = 1.18; $xtable[14] = 1.13; $xtable[15] = 1.092; $xtable[16] = 1.05; $xtable[17] = 1.025; $xtable[18] = 1.012; $xtable[19] = 1.004; $xtable[20] = 1.001; $xtable[21] = 1.000; ////////////////////////////////////////////////////////////////////////////////// //plate 4 $numberplate4 = 18; $aprimetable[0] = 0.01; $aprimetable[1] = 0.02; $aprimetable[2] = 0.03; $aprimetable[3] = 0.04; $aprimetable[4] = 0.05; $aprimetable[5] = 0.06; $aprimetable[6] = 0.07; $aprimetable[7] = 0.08; $aprimetable[8] = 0.09; $aprimetable[9] = 0.10; $aprimetable[10] = 0.20; $aprimetable[11] = 0.30; $aprimetable[12] = 0.40; $aprimetable[13] = 0.50; $aprimetable[14] = 0.60; $aprimetable[15] = 0.70; $aprimetable[16] = 0.80; $aprimetable[17] = 0.90; $aprimetable[18] = 1.00; $percnotable14[0] = 99.8; $percnotable14[1] = 98.9; $percnotable14[2] = 98.1; $percnotable14[3] = 97.4; //ok $percnotable14[4] = 96.5; $percnotable14[5] = 95.6; $percnotable14[6] = 94.8; //ok $percnotable14[7] = 94.0; $percnotable14[8] = 93.1; $percnotable14[9] = 92.2; //ok $percnotable14[10] = 83.7; $percnotable14[11] = 73.7; $percnotable14[12] = 63.8; $percnotable14[13] = 53.0; $percnotable14[14] = 43.4; $percnotable14[15] = 32.9; $percnotable14[16] = 20.9; $percnotable14[17] = 10.9; $percnotable14[18] = 0.; $percnotable11[0] = 99.8; $percnotable11[1] = 99.0; $percnotable11[2] = 98.2; $percnotable11[3] = 97.5; //ok $percnotable11[4] = 96.6; $percnotable11[5] = 95.8; $percnotable11[6] = 95.0; //ok $percnotable11[7] = 94.2; $percnotable11[8] = 93.3; $percnotable11[9] = 92.5; //ok $percnotable11[10] = 84.5; $percnotable11[11] = 74.5; $percnotable11[12] = 64.5; $percnotable11[13] = 54.0; $percnotable11[14] = 43.5; $percnotable11[15] = 33.0; $percnotable11[16] = 21.0; $percnotable11[17] = 11.0; $percnotable11[18] = 0.; $percnotable08[0] = 99.9; $percnotable08[1] = 99.2; $percnotable08[2] = 98.4; $percnotable08[3] = 97.7; //ok $percnotable08[4] = 97.2; $percnotable08[5] = 97.0; $percnotable08[6] = 96.0; //ok $percnotable08[7] = 95.2; $percnotable08[8] = 94.5; $percnotable08[9] = 94.0; //ok $percnotable08[10] = 85.5; $percnotable08[11] = 76.5; $percnotable08[12] = 66.0; $percnotable08[13] = 56.0; $percnotable08[14] = 45.2; $percnotable08[15] = 33.5; $percnotable08[16] = 21.5; $percnotable08[17] = 11.2; $percnotable08[18] = 0.; $percnotable04[0] = 100.; $percnotable04[1] = 100.; $percnotable04[2] = 100.; $percnotable04[3] = 100.; //ok $percnotable04[4] = 100.; $percnotable04[5] = 100.; $percnotable04[6] = 100.; //ok $percnotable04[7] = 100.; $percnotable04[8] = 100.; $percnotable04[9] = 99.5; //ok $percnotable04[10] = 93.8; $percnotable04[11] = 85.5; $percnotable04[12] = 75.0; $percnotable04[13] = 64.0; $percnotable04[14] = 51.8; $percnotable04[15] = 39.5; $percnotable04[16] = 26.0; $percnotable04[17] = 13.0; $percnotable04[18] = 0.; ////////////////////////////////////////////////////////////////////////////////////////////// //setting nsize $nsizenorm = 8.; $nsizelong = 9.; //$nsizelong2 = 10.; //echo "sofarsogood111 " . "$sizefractype" . "
"; if($sizefractype == 0) { $nsizeapp = $nsizenorm + 1; } elseif($sizefractype == 1) { $nsizeapp = $nsizelong + 1; } //echo "sofarsogood112 " . "$nsizeapp" . "
"; ///////////////////////////////////////////////////////////////////////////////////////////////// $sizelownormal[0] = 0.002; $sizelownormal[1] = 0.0625; $sizelownormal[2] = 0.125; $sizelownormal[3] = 0.25; $sizelownormal[4] = 0.5; $sizelownormal[5] = 1.; $sizelownormal[6] = 2.; $sizelownormal[7] = 4.; $sizelownormal[8] = 8.; $sizehighnormal[0] = 0.0625; $sizehighnormal[1] = 0.125; $sizehighnormal[2] = 0.25; $sizehighnormal[3] = 0.5; $sizehighnormal[4] = 1.; $sizehighnormal[5] = 2.; $sizehighnormal[6] = 4.; $sizehighnormal[7] = 8.; $sizehighnormal[8] = 16.; $sizelowlong[0] = 0.002; $sizelowlong[1] = 0.0156; $sizelowlong[2] = 0.0625; $sizelowlong[3] = 0.125; $sizelowlong[4] = 0.25; $sizelowlong[5] = 0.5; $sizelowlong[6] = 1.; $sizelowlong[7] = 2.; $sizelowlong[8] = 4.; $sizelowlong[9] = 8.; $sizehighlong[0] = 0.0156; $sizehighlong[1] = 0.0625; $sizehighlong[2] = 0.125; $sizehighlong[3] = 0.25; $sizehighlong[4] = 0.5; $sizehighlong[5] = 1.; $sizehighlong[6] = 2.; $sizehighlong[7] = 4.; $sizehighlong[8] = 8.; $sizehighlong[9] = 16.; $const = 304.8; if($sizefractype == 0) { $i = 0; while($i < $nsizeapp) { $sizediamnormalmm[$i] = pow(($sizelownormal[$i]*$sizehighnormal[$i]),0.5); $sizediamnormalft[$i] = $sizediamnormalmm[$i]/$const; $sizediamappft[$i] = $sizediamnormalft[$i]; $sizediamappmm[$i] = $sizediamnormalmm[$i]; $i++; } } elseif($sizefractype == 1) { $i = 0; while($i < $nsizeapp) { $sizediamlongmm[$i] = pow(($sizelowlong[$i]*$sizehighlong[$i]),0.5); $sizediamlongft[$i] = $sizediamlongmm[$i]/$const; $sizediamappft[$i] = $sizediamlongft[$i]; $sizediamappmm[$i] = $sizediamlongmm[$i]; $i++; } } ///////////////////////////////////////////////////////////////// //plate 5 $nplate5 = 24.; $psitable[0] = 0.08; $psitable[1] = 0.09; $psitable[2] = 0.1; $psitable[3] = 0.2; $psitable[4] = 0.3; $psitable[5] = 0.4; $psitable[6] = 0.5; $psitable[7] = 0.6; $psitable[8] = 0.7; $psitable[9] = 0.8; $psitable[10] = 0.9; $psitable[11] = 1.; $psitable[12] = 2.; $psitable[13] = 3.; $psitable[14] = 4.; $psitable[15] = 5.; $psitable[16] = 6.; $psitable[17] = 7.; $psitable[18] = 8.; $psitable[19] = 9.; $psitable[20] = 10.; $psitable[21] = 15.; $psitable[22] = 20.; $psitable[23] = 25.; $psitable[24] = 27.; $phitable[0] = 98.; $phitable[1] = 88.; $phitable[2] = 79.; $phitable[3] = 40.; $phitable[4] = 26.5; $phitable[5] = 19.6; $phitable[6] = 15.5; $phitable[7] = 12.8; $phitable[8] = 10.9; $phitable[9] = 9.5; $phitable[10] = 8.4; $phitable[11] = 7.4; $phitable[12] = 3.25; $phitable[13] = 1.8; $phitable[14] = 1.05; $phitable[15] = 0.65; $phitable[16] = 0.41; $phitable[17] = 0.27; $phitable[18] = 0.18; $phitable[19] = 0.125; $phitable[20] = 0.088; $phitable[21] = 0.017; $phitable[22] = 0.0032; $phitable[23] = 0.00033; $phitable[23] = 0.0001; ////////////////////////////////////////////////////////////////////////////////////// //plate 8 $nplate08 = 1; $absplate08[0] = 400.; $ordplate08[0] = 0.21; $absplate08[1] = 0.03; $ordplate08[1] = 1.6; ////////////////////////////////////////////////////////////////////////////////////// //units output if ($unitsystems == 0) { $dischargeunits = "cfs"; $velocityunits = "fps"; $depthunits = "ft"; $areaunits = "ft²"; $temperatureunits = "^oF"; $viscosityunits = "ft²/s"; $diameterunits = "ft"; $concentrationunits = "ppm"; $loadunits = "tons/day"; } elseif($unitsystems == 1) { $dischargeunits = "m³/s"; $velocityunits = "m/s"; $depthunits = "m"; $areaunits = "m²"; $temperatureunits = "^oC"; $viscosityunits = "m²/s"; $diameterunits = "mm"; $concentrationunits = "ppm"; $loadunits = "tons/day"; $loadunitsout = "M.Tons/day"; } ////////////////////////////////////////////////////////////////////////////////////// ////////main computation $area = $wdin*$dein; $samplesuspload = 0.0027*$qdin*$coin; if($unitsystems == 1) { $areametric = $area*0.092903; $samplesusploadmetric = 0.9072*$samplesuspload; $viscositymetric = $viscosity*0.092903; } //calculation of x by trial and error $ks = $d65in; $flag = 0; $xin = 1.54; while ($flag == 0) { //echo "xin " . $xin . "
"; $temp1 = 12.27*$xin*$dein/$ks; //echo "temp1 " . $temp1 . "
"; $temp2 = log10($temp1); //echo "temp2 " . $temp2 . "
"; $temp3 = 32.63*$temp2; //echo "temp3 " . $temp3 . "
"; $sqsr = $vein/$temp3; //echo "sqsr " . $sqsr . "
"; $sr = pow($sqsr,2); //echo "sr " . $sr . "
"; $gravsr = 5.68; //$gravsr = 5.672; //exact value $ups = $gravsr*$sqsr; //echo "ups " . $ups . "
"; $delta = 11.6*$viscosity/$ups; //echo "delta " . $delta . "
"; $ksovdelta = $ks/$delta; //echo "ksovdelta " . $ksovdelta . "
"; $xout = pwlin0logx($ksovdeltatable,$xtable,$numberksdx,$ksovdelta,$xout); //echo "xout " . $xout . "
"; $diff = $xout - $xin; //echo "diff " . $diff . "
"; if(abs($diff) < 0.005) { $flag = 1; } else { $xin= $xout; } } $xxx = $xout; //echo "sofarsogood_die"; $ppp = 2.303*log10(30.2*$xxx*$dein/$ks); //echo "ppp " . $ppp . "
"; $aprime= $dnin/$dsin; //echo "aprime " . $aprime . "
"; if($ppp <= 4) { $pfsampled = pwlin0logx($aprimetable,$percnotable04,$numberplate4,$aprime,$pfsampled); //echo "pfsampled " . $pfsampled . "
"; } elseif($ppp < 8) { $pfsampled04 = pwlin0logx($aprimetable,$percnotable04,$numberplate4,$aprime,$pfsampled04); $pfsampled08 = pwlin0logx($aprimetable,$percnotable08,$numberplate4,$aprime,$pfsampled08); $numberpo = 1.; $abscissa[0] = 4.; $abscissa[1] = 8.; $ordinate[0] = $pfsampled04; $ordinate[1] = $pfsampled08; $pfsampled = pwlin0logx($abscissa,$ordinate,$numberpo,$ppp,$pfsampled); //echo "pfsampled " . $pfsampled . "
"; } elseif($ppp == 8) { $pfsampled = pwlin0logx($aprimetable,$percnotable08,$numberplate4,$aprime,$pfsampled); //echo "pfsampled " . $pfsampled . "
"; } elseif($ppp < 11) { $pfsampled08 = pwlin0logx($aprimetable,$percnotable08,$numberplate4,$aprime,$pfsampled08); $pfsampled11 = pwlin0logx($aprimetable,$percnotable11,$numberplate4,$aprime,$pfsampled11); $numberpo = 1.; $abscissa[0] = 8.; $abscissa[1] = 11.; $ordinate[0] = $pfsampled08; $ordinate[1] = $pfsampled11; $pfsampled = pwlin0logx($abscissa,$ordinate,$numberpo,$ppp,$pfsampled); //echo "pfsampled " . $pfsampled . "
"; } elseif($ppp == 11) { $pfsampled = pwlin0logx($aprimetable,$percnotable11,$numberplate4,$aprime,$pfsampled); //echo "pfsampled " . $pfsampled . "
"; } elseif($ppp < 14) { $pfsampled11 = pwlin0logx($aprimetable,$percnotable11,$numberplate4,$aprime,$pfsampled11); $pfsampled14 = pwlin0logx($aprimetable,$percnotable14,$numberplate4,$aprime,$pfsampled14); $numberpo = 1.; $abscissa[0] = 11.; $abscissa[1] = 14.; $ordinate[0] = $pfsampled11; $ordinate[1] = $pfsampled14; $pfsampled = pwlin0logx($abscissa,$ordinate,$numberpo,$ppp,$pfsampled); //echo "pfsampled " . $pfsampled . "
"; } elseif($ppp >= 14) { $pfsampled = pwlin0logx($aprimetable,$percnotable14,$numberplate4,$aprime,$pfsampled); //echo "pfsampled " . $pfsampled . "
"; } $pfsampledw = round($pfsampled,0); //echo "pfsampledw " . $pfsampledw . "
"; $samplesuspload = 0.0027*$qdin*$coin; //echo "samplesuspload " . $samplesuspload . "
"; $qsprime = $pfsampledw*$samplesuspload/100.; //$qsprime = 0.03; $qsprimew = round($qsprime,0); $qsprimew2 = number_format($qsprime,0); if($qsprimew == 0.) { $qsprimew = round($qsprime,2); } else { $qsprimew = round($qsprime,0); } //echo "qsprimew " . $qsprimew . "
"; //echo "qsprime " . $qsprime . "
"; $aprime = $dnin/$dsin; /////////////////////////////////////////////////////////////////////////////////////////// //explosion of arrays $ibvar = explode(",",$ibvarin); //echo "ibvar[0]" . $ibvar[0] . "
"; //echo "ibvar[1]" . $ibvar[1] . "
"; //echo "ibvar[2]" . $ibvar[2] . "
"; //echo "ibvar[3]" . $ibvar[3] . "
"; //echo "ibvar[4]" . $ibvar[4] . "
"; //echo "ibvar[5]" . $ibvar[5] . "
"; //echo "ibvar[6]" . $ibvar[6] . "
"; //echo "ibvar[7]" . $ibvar[7] . "
"; //echo "ibvar[8]" . $ibvar[8] . "
"; if($sizefractype == 1) { //echo "ibvar[9]" . $ibvar[9] . "
"; } $qspvar = explode(",",$qspvarin); //echo "qspvar[0]" . $qspvar[0] . "
"; //echo "qspvar[1]" . $qspvar[1] . "
"; //echo "qspvar[2]" . $qspvar[2] . "
"; //echo "qspvar[3]" . $qspvar[3] . "
"; //echo "qspvar[4]" . $qspvar[4] . "
"; //echo "qspvar[5]" . $qspvar[5] . "
"; //echo "qspvar[6]" . $qspvar[6] . "
"; //echo "qspvar[7]" . $qspvar[7] . "
"; //echo "qspvar[8]" . $qspvar[8] . "
"; if($sizefractype == 1) { //echo "qspvar[9]" . $qspvar[9] . "
"; } ///////////////////////////////////////////////////////////////////////////////////////////// //testing of array size read in if($sizefractype == 0) { $i = 0; while($i < $nsizeapp) { //echo "i " . $i . " " . $ibvar[$i] . "
"; if($ibvar[$i] < 0) { if($_POST["Submit"] == 'Calculate') //echo "so far so good" . "
"; echo "

"; echo " online_modified_einstein:
*** Error No. 13F (Fatal)
*** Please enter only zero or positive number for fraction of bed material i_b, size fraction " . number_format($i+1,0) . "
*** Execution stops here. "; echo "
"; die; } $i++; } } elseif($sizefractype == 1) { $i = 0; while($i < $nsizeapp) { if($ibvar[$i] < 0) { if($_POST["Submit"] == 'Calculate') echo "

"; echo " online_modified_einstein:
*** Error No. 13F (Fatal)
*** Please enter only zero or positive number for fraction of bed material i_b. size fraction " . number_format($i+1,0) . "
*** Execution stops here. "; echo "
"; die; } $i++; } } if($sizefractype == 0) { $i = 0; while($i < $nsizeapp) { //echo "i " . $i . " " . $qspvar[$i] . "
"; if($qspvar[$i] < 0) { if($_POST["Submit"] == 'Calculate') //echo "so far so good" . "
"; echo "

"; echo " online_modified_einstein:
*** Error No. 14F (Fatal)
*** Please enter only zero or positive number for percentage of suspended material Q'_s, size fraction " . number_format($i+1,0) . "
*** Execution stops here. "; echo "
"; die; } $i++; } } elseif($sizefractype == 1) { $i = 0; while($i < $nsizeapp) { if($qspvar[$i] < 0) { if($_POST["Submit"] == 'Calculate') echo "

"; echo " online_modified_einstein:
*** Error No. 14F (Fatal)
*** Please enter only zero or positive number for percentage of suspended material Q'_s. size fraction " . number_format($i+1,0) . "
*** Execution stops here. "; echo "
"; die; } $i++; } } /////////////////////////////////////////////////////////////////////////////////////////// //testing for common size fractions in bed material and suspended material 150228 $i = 0; while($i < $nsizeapp) { //echo "i = " . $i . " ibvar = " . $ibvar[$i] . "
"; $i++; } $i = 0; while($i < $nsizeapp) { //echo "i = " . $i . " qspvar = " . $qspvar[$i] . "
"; $i++; } $i = 0; $qspvarmax = 0.; $iqspvarmax = 0; while($i < $nsizeapp) { if($qspvar[$i] > $qspvarmax) { $qspvarmax = $qspvar[$i]; $iqspvarmax = $i; } //echo "i = " . $i . " qspvarmax = " . $qspvarmax . " iqspvarmax = " . $iqspvarmax . "
"; $i++; } if($ibvar[$iqspvarmax] < 0.01) { //echo "so far so good " . "
"; //die; if($_POST["Submit"] == 'Calculate') { echo "

"; echo ""; echo " online_modified_einstein:" . "
"; echo "
*** Error No. 15F (Fatal)" . "
"; echo "
*** The reference size fraction, which is the maximum SUSPENDED MATERIAL fraction (Q_s'_max)" . "
" . "must have a nonzero entry (at least 0.01) in the corresponding bed material size fraction i_b." . "
"; echo "
*** Execution stops here. "; echo "
"; die; } } //die; //////////////////////////////////////////////////////////////////////////////////////////// //calculation of psi (Column 4) $psi01 = 1.65*$d35in/$sr; $i = 0; while($i < $nsizeapp) { $psi02[$i] = 0.66*$sizediamappft[$i] /$sr; $i++; } $i = 0; while($i < $nsizeapp) { if($psi02[$i] > $psi01) { $psi[$i] = $psi02[$i]; } else { $psi[$i] = $psi01; } $i++; } /////////////////////////////////////////////////////////////////////////////////////////////// //column 5 $i = 0.; while ($i < $nsizeapp) { if($psi[$i] > 27.) { $phistar[$i] = 0.; $phistar2[$i] = 0.; } else { $phistar[$i] = pwlin0logxlogy($psitable,$phitable,$nplate5,$psi[$i],$phistar[$i]); $phistar2[$i] = 0.5*$phistar[$i]; } //echo "phistar2[i] " . $i, " " . $phistar2[$i] . "
"; $i++; } /////////////////////////////////////////////////////////////////////////////////////////////// //column 9 $w432 = 43.2*$wdin; //echo "w432 " . " " . $w432 . "
"; /////////////////////////////////////////////////////////////////////////////////////////////// //column 6, 7, 8, 10 $i = 0.; while ($i < $nsizeapp) { $d1200[$i] = 1200.*pow($sizediamappft[$i],1.5); //echo "d1200[i] " . $i, " " . $sizediamappft[$i] . " " . $d1200[$i] . "
"; $ibqb[$i] = $d1200[$i]*$ibvar[$i]*$phistar2[$i]; //echo "ibqb[i] " . $i, " " . $ibqb[$i] . "
"; $ibqbcap[$i] = $ibqb[$i]*$w432; //echo "ibqbcap[i] " . $i, " " . $ibqbcap[$i] . "
"; $i++; } /////////////////////////////////////////////////////////////////////////////////////////////// //column 12 $i = 0.; while ($i < $nsizeapp) { $qsprimei[$i] = $qspvar[$i]*$qsprimew/100.; //echo "qsprimei[i] " . $i, " " . $qsprimei[$i] . "
"; $i++; } $i = 0.; $qsprimelocmax = $i; $qsprimeimax = $qsprimei[$i]; while ($i < $nsizeapp) { //echo "qsprimei[i] " . $i, " " . $qsprimei[$i] . " " . $qsprimeimax . " " . $qsprimelocmax . "
"; $i++; if($qsprimei[$i] > $qsprimeimax) { $qsprimeimax = $qsprimei[$i]; $qsprimelocmax = $i; } } //echo "qsprimelocmax" . " " . $qsprimelocmax . "
"; $consfall = 0.6666667*1.65*32.17; $i = 0.; while ($i < $nsizeapp) { //echo $i. " " . $viscosity . " " . $sizediamappft[$i] . "
"; $tempnu = $viscosity/$sizediamappft[$i]; //echo $i. " " . $tempnu . "
"; $tempnu2 = $tempnu*$tempnu; $fallvelocity[$i] = pow((36.*$tempnu2 + $consfall*$sizediamappft[$i]),0.5) - 6.*$tempnu; //echo "fallvelocity[i]" . $i . " " . $fallvelocity[$i] . "
"; $i++; } /////////////////////////////////////////////////////////////////////////////////////////////// //column 13, 15 $i = 0.; $loc = $qsprimelocmax; while ($i < $nsizeapp) { $multiplier[$i] = pow(($fallvelocity[$i]/$fallvelocity[$loc]),0.7); //echo "multiplier[i]" . " " . $i . " " . $multiplier[$i] . "
"; $adoubleprime[$i] = 2.*$sizediamappft[$i]/$dein; //echo "adoubleprime[i]" . " " . $i . " " . $adoubleprime[$i] . "
"; $i++; } //echo "sofarsogood 1"; ////////////////////////////////////////////////////////////////////////////////////////////// //iteration to calculate zprime for reference size fraction $k = 0; $qsovibqb0 = $qsprimei[$qsprimelocmax+$k]/$ibqbcap[$qsprimelocmax+$k]; //echo "
" . "k = " . $k . " qsovibqb0 " . $qsovibqb0 . "
"; $ztrial = pwlin0logx($absplate08,$ordplate08,$nplate08,$qsovibqb0,$zprimetrial); //echo "ztrial " . $ztrial . "
"; $i = 0; //echo "i " . $i; $deltaz = 0.01; $z = $ztrial; $zstore[$i] = $z; //echo "
" . "z " . $z . "
"; $x = $adoubleprime[$qsprimelocmax+$k]; $y = modeinsplate09($x,$z,$y);"
"; $i1doubleprime = $y; //echo "i1doubleprime " . $i1doubleprime . "
"; $y = modeinsplate10($x,$z,$y);"
"; $j1doubleprime = $y; //echo "j1doubleprime " . $j1doubleprime . "
"; $x = $aprime; $y = modeinsplate10($x,$z,$y);"
"; $j1prime = $y; //echo "j1prime " . $j1prime . "
"; $y = modeinsplate11($x,$z,$y);"
"; $j2prime = $y; //echo "j2prime " . $j2prime . "
"; $righthandside = ($i1doubleprime/$j1doubleprime)*($ppp*$j1prime - $j2prime); //echo "righthandside " . $righthandside . "
"; $difference = $righthandside - $qsovibqb0; $differencestore[$i] = $difference; //echo "difference " . $difference . "
"; if($differencestore[0] < 0) { while($difference < 0) { $i++; //echo "i " . $i; $z = $ztrial - $i*$deltaz; $zstore[$i] = $z; //echo "
" . "z " . $z . "
"; $x = $adoubleprime[$qsprimelocmax+$k]; $y = modeinsplate09($x,$z,$y);"
"; $i1doubleprime = $y; //echo "i1doubleprime " . $i1doubleprime . "
"; $y = modeinsplate10($x,$z,$y);"
"; $j1doubleprime = $y; //echo "j1doubleprime " . $j1doubleprime . "
"; $x = $aprime; $y = modeinsplate10($x,$z,$y);"
"; $j1prime = $y; //echo "j1prime " . $j1prime . "
"; $y = modeinsplate11($x,$z,$y);"
"; $j2prime = $y; //echo "j2prime " . $j2prime . "
"; $righthandside = ($i1doubleprime/$j1doubleprime)*($ppp*$j1prime - $j2prime); //echo "righthandside " . $righthandside . "
"; $difference = $righthandside - $qsovibqb0; $differencestore[$i] = $difference; //echo "difference " . $difference . "
"; } //echo "i " . $i; $zfinal = $zstore[$i] + ($zstore[$i-1] - $zstore[$i])*$differencestore[$i]/($differencestore[$i] - $differencestore[$i-1]); //echo "
" . "zfinal " . $zfinal; $zfinalw = number_format($zfinal,2); //echo "
" . "zfinalw " . $zfinalw ."
"; } elseif($differencestore[0] > 0) { while($difference > 0) { $i++; //echo "i " . $i; $z = $ztrial + $i*$deltaz; $zstore[$i] = $z; //echo "
" . "z " . $z . "
"; $x = $adoubleprime[$qsprimelocmax+$k]; $y = modeinsplate09($x,$z,$y);"
"; $i1doubleprime = $y; //echo "i1doubleprime " . $i1doubleprime . "
"; $y = modeinsplate10($x,$z,$y);"
"; $j1doubleprime = $y; //echo "j1doubleprime " . $j1doubleprime . "
"; $x = $aprime; $y = modeinsplate10($x,$z,$y);"
"; $j1prime = $y; //echo "j1prime " . $j1prime . "
"; $y = modeinsplate11($x,$z,$y);"
"; $j2prime = $y; //echo "j2prime " . $j2prime . "
"; $righthandside = ($i1doubleprime/$j1doubleprime)*($ppp*$j1prime - $j2prime); //echo "righthandside " . $righthandside . "
"; $difference = $righthandside - $qsovibqb0; $differencestore[$i] = $difference; //echo "difference " . $difference . "
"; } //echo "i " . $i; $zfinal = $zstore[$i] + ($zstore[$i-1] - $zstore[$i])*$differencestore[$i]/($differencestore[$i] - $differencestore[$i-1]); //echo "
" . "zfinal " . $zfinal; $zfinalw = number_format($zfinal,2); //echo "
" . "zfinalw " . $zfinalw ."
"; } $zprimeout[$qsprimelocmax+$k] = $zfinal; //echo "sofarsogood 2"; ////////////////////////////////////////////////////////////////////////////////////////////// //iteration to calculate zprime for reference size fraction + 1 if($zprimemode == 1) { $k = 1; $qsovibqb0 = $qsprimei[$qsprimelocmax+$k]/$ibqbcap[$qsprimelocmax+$k]; //echo "
" . "k = " . $k . " qsovibqb0 " . $qsovibqb0 . "
"; $ztrial = pwlin0logx($absplate08,$ordplate08,$nplate08,$qsovibqb0,$zprimetrial); //echo "ztrial " . $ztrial . "
"; $i = 0; //echo "i " . $i; $deltaz = 0.01; $z = $ztrial; $zstore[$i] = $z; //echo "
" . "z " . $z . "
"; $x = $adoubleprime[$qsprimelocmax+$k]; $y = modeinsplate09($x,$z,$y);"
"; $i1doubleprime = $y; //echo "i1doubleprime " . $i1doubleprime . "
"; $y = modeinsplate10($x,$z,$y);"
"; $j1doubleprime = $y; //echo "j1doubleprime " . $j1doubleprime . "
"; $x = $aprime; $y = modeinsplate10($x,$z,$y);"
"; $j1prime = $y; //echo "j1prime " . $j1prime . "
"; $y = modeinsplate11($x,$z,$y);"
"; $j2prime = $y; //echo "j2prime " . $j2prime . "
"; $righthandside = ($i1doubleprime/$j1doubleprime)*($ppp*$j1prime - $j2prime); //echo "righthandside " . $righthandside . "
"; $difference = $righthandside - $qsovibqb0; $differencestore[$i] = $difference; //echo "difference " . $difference . "
"; if($differencestore[0] < 0) { while($difference < 0) { $i++; //echo "i " . $i; $z = $ztrial - $i*$deltaz; $zstore[$i] = $z; //echo "
" . "z " . $z . "
"; $x = $adoubleprime[$qsprimelocmax+$k]; $y = modeinsplate09($x,$z,$y);"
"; $i1doubleprime = $y; //echo "i1doubleprime " . $i1doubleprime . "
"; $y = modeinsplate10($x,$z,$y);"
"; $j1doubleprime = $y; //echo "j1doubleprime " . $j1doubleprime . "
"; $x = $aprime; $y = modeinsplate10($x,$z,$y);"
"; $j1prime = $y; //echo "j1prime " . $j1prime . "
"; $y = modeinsplate11($x,$z,$y);"
"; $j2prime = $y; //echo "j2prime " . $j2prime . "
"; $righthandside = ($i1doubleprime/$j1doubleprime)*($ppp*$j1prime - $j2prime); //echo "righthandside " . $righthandside . "
"; $difference = $righthandside - $qsovibqb0; $differencestore[$i] = $difference; //echo "difference " . $difference . "
"; } //echo "i " . $i; $zfinal = $zstore[$i] + ($zstore[$i-1] - $zstore[$i])*$differencestore[$i]/($differencestore[$i] - $differencestore[$i-1]); //echo "
" . "zfinal " . $zfinal; $zfinalw = number_format($zfinal,2); //echo "
" . "zfinalw " . $zfinalw ."
"; } elseif($differencestore[0] > 0) { while($difference > 0) { $i++; //echo "i " . $i; $z = $ztrial + $i*$deltaz; $zstore[$i] = $z; //echo "
" . "z " . $z . "
"; $x = $adoubleprime[$qsprimelocmax+$k]; $y = modeinsplate09($x,$z,$y);"
"; $i1doubleprime = $y; //echo "i1doubleprime " . $i1doubleprime . "
"; $y = modeinsplate10($x,$z,$y);"
"; $j1doubleprime = $y; //echo "j1doubleprime " . $j1doubleprime . "
"; $x = $aprime; $y = modeinsplate10($x,$z,$y);"
"; $j1prime = $y; //echo "j1prime " . $j1prime . "
"; $y = modeinsplate11($x,$z,$y);"
"; $j2prime = $y; //echo "j2prime " . $j2prime . "
"; $righthandside = ($i1doubleprime/$j1doubleprime)*($ppp*$j1prime - $j2prime); //echo "righthandside " . $righthandside . "
"; $difference = $righthandside - $qsovibqb0; $differencestore[$i] = $difference; //echo "difference " . $difference . "
"; } //echo "i " . $i; $zfinal = $zstore[$i] + ($zstore[$i-1] - $zstore[$i])*$differencestore[$i]/($differencestore[$i] - $differencestore[$i-1]); //echo "
" . "zfinal " . $zfinal; $zfinalw = number_format($zfinal,2); //echo "
" . "zfinalw " . $zfinalw ."
"; } $zprimeout[$qsprimelocmax+$k] = $zfinal; //echo "zprimeout(loc+0) = " . $zprimeout[$qsprimelocmax+0] . " fall velocity = " . $fallvelocity[$qsprimelocmax+0] . "
"; //echo "zprimeout(loc+1) = " . $zprimeout[$qsprimelocmax+1] . " fall velocity = " . $fallvelocity[$qsprimelocmax+1] . "
"; //calculation of lara's 1966 zprime $xfallveloc[0] = $fallvelocity[$qsprimelocmax+0]; $xfallveloc[1] = $fallvelocity[$qsprimelocmax+1]; $yzprimeout[0] = $zprimeout[$qsprimelocmax+0]; $yzprimeout[1] = $zprimeout[$qsprimelocmax+1]; $i = 0.; while ($i < $nsizeapp) { $zprimeoutlara[$i] = pwlin0logxlogy($xfallveloc,$yzprimeout,1,$fallvelocity[$i],$zprimeoutlara[$i]); //echo "zprimeoutlara[i]" . " " . $i . " " . $zprimeoutlara[$i] . "
"; $i++; } } //echo "sofarsogood 3"; /////////////////////////////////////////////////////////////////////////////////////////////// //calculation of usbr 1955 z prime if($zprimemode == 0) { $i = 0.; while ($i < $nsizeapp) { $zprimeout[$i] = $multiplier[$i]*$zprimeout[$qsprimelocmax]; //echo "zprimeout[i]" . " " . $i . " " . $zprimeout[$i] . "
"; $i++; } } elseif($zprimemode == 1) { $i = 0.; while ($i < $nsizeapp) { $zprimeout[$i] = $zprimeoutlara[$i]; //echo "zprimeout[i]" . " " . $i . " " . $zprimeout[$i] . "
"; $i++; } } /////////////////////////////////////////////////////////////////////////////////////////////// //columns 16, 17, 18, 19, 20, 24 (finer) $i = 0.; while($i < $loc) { $j1primecol16[$i] = modeinsplate10($aprime,$zprimeout[$i],$j1primecol16[$i]); //echo "
" . "j1primecol16 " . $i . " " . $j1primecol16[$i]; $j2primecol17[$i] = modeinsplate11($aprime,$zprimeout[$i],$j2primecol17[$i]); //echo "
" . "j2primecol17 " . $i . " " . $j2primecol17[$i]; $j1doubleprimecol18[$i] = modeinsplate10($adoubleprime[$i],$zprimeout[$i],$j1doubleprimecol18[$i]); //echo "
" . "j1doubleprimecol18 " . $i . " " . $j1doubleprimecol18[$i]; $j2doubleprimecol19[$i] = modeinsplate11($adoubleprime[$i],$zprimeout[$i],$j2doubleprimecol19[$i]); //echo "
" . "j2doubleprimecol19 " . $i . " " . $j2doubleprimecol19[$i]; $column20[$i] = ($ppp*$j1doubleprimecol18[$i] - $j2doubleprimecol19[$i])/($ppp*$j1primecol16[$i] - $j2doubleprimecol17[$i]); //echo "
" . "column20 " . $i . " " . $column20[$i]; $column24[$i] = $qsprimei[$i]*$column20[$i]; //echo "
" . "column24 " . $i . " " . $column24[$i]; $i++; } /////////////////////////////////////////////////////////////////////////////////////////////// //columns 21, 22, 23, 24 (coarser) while($i < $nsizeapp) { //echo "
" . "$i"; if($zprimeout[$i] > 5.) { $zprimeloc[$i] = 5.; } else { $zprimeloc[$i] = $zprimeout[$i]; } $i1doubleprimecol21[$i] = modeinsplate09($adoubleprime[$i],$zprimeloc[$i],$i1doubleprimecol21[$i]); //echo "
" . "i1doubleprimecol21 " . $i . " " . $i1doubleprimecol21[$i]; $i2doubleprimecol22[$i] = modeinsplate12($adoubleprime[$i],$zprimeloc[$i],$i2doubleprimecol22[$i]); //echo "
" . "i2doubleprimecol22 " . $i . " " . $i2doubleprimecol22[$i]; if($zprimeout[$i] > 5.) { $temporateabs[0] = 4.; $temporateabs[1] = 5.; $temporateord[0] = modeinsplate09($adoubleprime[$i],4.,$temporateord[0]); $temporateord[1] = modeinsplate09($adoubleprime[$i],5.,$temporateord[1]); $i1doubleprimecol21[$i] = pwlin0logxlogy($temporateabs,$temporateord,1,$zprimeout[$i],$i1doubleprimecol21[$i]); //echo "
" . "i1doubleprimecol21 alt " . $i . " " . $i1doubleprimecol21[$i]; $temporateord[0] = modeinsplate12($adoubleprime[$i],4.,$temporateord[0]); $temporateord[1] = modeinsplate12($adoubleprime[$i],5.,$temporateord[1]); $i2doubleprimecol22[$i] = pwlin0logxlogy($temporateabs,$temporateord,1,$zprimeout[$i],$i2doubleprimecol22[$i]); //echo "
" . "i2doubleprimecol22 alt " . $i . " " . $i2doubleprimecol22[$i]; } $column23[$i] = ($ppp*$i1doubleprimecol21[$i] - $i2doubleprimecol22[$i] + 1.); //echo "
" . "column23 " . $i . " " . $column23[$i]; $column24[$i] = $ibqbcap[$i]*$column23[$i]; //echo "
" . "column24 " . $i . " " . $column24[$i]; $i++; } /////////////////////////////////////////////////////////////////////////////////////////////// //sum of column 24 $i = 0; $sumload = 0.; while($i < $nsizeapp) { $sumload = $sumload + $column24[$i]; //echo "
" . "sumload " . $i . " " . $sumload; $i++; } //echo "
" . "sum load " . $i . " " . $sumload; if($unitsystems == 1) { $sumload25 = 0.; $i = 0; while($i < $nsizeapp) { $column25[$i] = 0.9072*$column24[$i]; $sumload25 = $sumload25 + $column25[$i]; //echo "
" . "sumload " . $i . " " . $sumload; $i++; } //echo "
" . "sum load " . $i . " " . $sumload; } } //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// ?>

    online_modified_einstein: Total sediment discharge by the Modified Einstein Procedure    (Version 2.0 Beta 150228)

References

Einstein, H. A., 1950. The bed-load function for sediment transportation in open channel flows. USDA Soil Conservation Service, Technical Bulletin No. 1026, Washington, D.C., September.

Colby, B. R. and C. H. Hembree, 1955. Computations of total sediment discharge, Niobrara river, near Cody, Nebraska. U.S. Geological Survey Water Supply Paper 1357, Washington, D.C.

U.S. Bureau of Reclamation, 1955. Step method for computing total sediment load by the Modified Einstein Procedure. Sedimentation Section, Hydrology Branch, Project Investigations Division, July.

U.S. Bureau of Reclamation, 1966. Computation of "Z's" for use in the Modified Einstein Procedure. Sedimentation Section, Hydrology Branch, Project Investigations Division, June.

INPUT DATA:

[Sample input]     [Rationale]     [Procedure]     [Features]
Stream or river (optional):
Select:        U.S. Customary or SI units [Default: U.S. Customary Units].

1. Discharge Q [cfs] (m³/s):      " . $dischargeunits . ""; } elseif($unitsystems == 0) { echo "" . $dischargeunits . ""; } ?>
2. Mean velocity u [fps] (m/s):      " . $velocityunits . ""; } elseif($unitsystems == 0) { echo "" . $velocityunits . ""; } ?>
3. Average depth d [ft] (m):      " . $depthunits . ""; } elseif($unitsystems == 0) { echo "" . $depthunits . ""; } ?>
4. Width w [ft] (m):      " . $depthunits . ""; } elseif($unitsystems == 0) { echo "" . $depthunits . ""; } ?> " . "5. Flow area A = " . round($areametric,1) . "  "; echo "" . $areaunits . ""; } elseif($unitsystems == 0) { echo "
" . "5. Flow area A = " . round($area,0) . "  "; echo "" . $areaunits . ""; } ?>
6. Water temperature t [^oF] (^oC):      " . $temperatureunits . ""; } elseif($unitsystems == 0) { echo "" . $temperatureunits . ""; } ?>
" . "7. Kinematic viscosity ν = " . number_format($viscositymetric,9) . " "; echo "" . $viscosityunits . ""; } elseif($unitsystems == 0) { echo "
" . "7. Kinematic viscosity ν = " . number_format($viscosity,8) . " "; echo "" . $viscosityunits . ""; } ?>

8. Diameter D₆₅ [ft] (mm):      " . $diameterunits . ""; } elseif($unitsystems == 0) { echo "" . $diameterunits . ""; } ?>
9. Diameter D₃₅ [ft] (mm):      " . $diameterunits . ""; } elseif($unitsystems == 0) { echo "" . $diameterunits . ""; } ?>
10. Measured concentration C [ppm] (ppm):      " . $concentrationunits . ""; } elseif($unitsystems == 0) { echo "" . $concentrationunits . ""; } ?>
" . "11. Sampled suspended load Q_sm = " . round($samplesusploadmetric,0) . "  "; echo "" . $loadunitsout . ""; } elseif($unitsystems == 0) { echo "
" . "11. Sampled suspended load Q_sm = " . round($samplesuspload,0) . "  "; echo "" . $loadunits . ""; } ?>
12. Vertical distance not sampled d_n [ft] (m):      " . $depthunits . ""; } elseif($unitsystems == 0) { echo "" . $depthunits . ""; } ?>
13. Mean water depth at sampled verticals d_s [ft] (m):      " . $depthunits . ""; } elseif($unitsystems == 0) { echo "" . $depthunits . ""; } ?>

14. Number of size fractions:     "; ?>   Normal has 9 size fractions (1 fine and 8 coarse); Long has 10 size fractions (2 fine and 8 coarse) [Default: Normal].
15. Enter fraction of bed material in each size range (make sure it is 9 values for Normal; 10 values for Long) (separate each value with a comma; maximum character length is 100 characters).
i_b:
16. Enter percentage of suspended material in each size range (make sure it is 9 values for Normal; 10 values for Long) (separate each value with a comma; maximum character length is 100 characters).
Q'_s:
17. Z' mode of computation:       See No. 6 in Features [Default: USBR (1955)].

OUTPUT:
"; echo "
"; echo "NORMAL"; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo "
0"; echo " 1"; echo " 2"; echo " 3"; echo " 4"; echo " 5"; echo " 6"; echo " 7"; echo " 8"; echo " 9"; echo " 10"; echo " 11"; echo " 12"; echo " 13"; echo " 14"; echo " 15"; echo " 16"; echo " 17"; echo " 18"; echo " 19"; echo " 20"; echo " 21"; echo " 22"; echo " 23"; echo " 24"; echo "
i"; echo " Size
fraction
(mm)"; echo " D
(mm)"; echo " D
(ft)"; echo " ψ"; echo " φ_*/2"; echo " 1200 D^3/2"; echo " i_b"; echo " i_bq_b"; echo " 43.2 w"; echo " i_BQ_B
(T/D)"; echo "
% Size
Q'_s"; echo " Q'_s
(T/D)"; echo "
Mult."; echo " Z'"; echo " A\""; echo " J'₁"; echo " -J'₂"; echo " J\"₁"; echo " -J\"₂"; echo " PJ\"₁+J\"₂
________
PJ'₁+J'₂ "; echo "
I\"₁"; echo " -I\"₂"; echo " PI\"₁+I\"₂+1"; echo " Computed
load
(T/D)"; echo "
1"; echo " $sizelownormal[0]-$sizehighnormal[0]"; echo " " . number_format($sizediamnormalmm[0],4); echo " " . number_format($sizediamnormalft[0],6); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo "
2"; echo " $sizelownormal[1]-$sizehighnormal[1]"; echo " " . number_format($sizediamnormalmm[1],4); echo " " . number_format($sizediamnormalft[1],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo "
3"; echo " $sizelownormal[2]-$sizehighnormal[2]"; echo " " . number_format($sizediamnormalmm[2],3); echo " " . number_format($sizediamnormalft[2],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo "
4"; echo " $sizelownormal[3]-$sizehighnormal[3]"; echo " " . number_format($sizediamnormalmm[3],3); echo " " . number_format($sizediamnormalft[3],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo "
5"; echo " $sizelownormal[4]-$sizehighnormal[4]"; echo " " . number_format($sizediamnormalmm[4],3); echo " " . number_format($sizediamnormalft[4],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo "
6"; echo " $sizelownormal[5]-$sizehighnormal[5]"; echo " " . number_format($sizediamnormalmm[5],2); echo " " . number_format($sizediamnormalft[5],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo "
7"; echo " $sizelownormal[6]-$sizehighnormal[6]"; echo " " . number_format($sizediamnormalmm[6],2); echo " " . number_format($sizediamnormalft[6],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo "
8"; echo " $sizelownormal[7]-$sizehighnormal[7]"; echo " " . number_format($sizediamnormalmm[7],2); echo " " . number_format($sizediamnormalft[7],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo "
9"; echo " $sizelownormal[8]-$sizehighnormal[8]"; echo " " . number_format($sizediamnormalmm[8],2); echo " " . number_format($sizediamnormalft[8],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo "
-"; echo " Total"; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo "
"; echo "
LONG"; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo "
0"; echo " 1"; echo " 2"; echo " 3"; echo " 4"; echo " 5"; echo " 6"; echo " 7"; echo " 8"; echo " 9"; echo " 10"; echo " 11"; echo " 12"; echo " 13"; echo " 14"; echo " 15"; echo " 16"; echo " 17"; echo " 18"; echo " 19"; echo " 20"; echo " 21"; echo " 22"; echo " 23"; echo " 24"; echo "
i"; echo " Size
fraction
(mm)"; echo " D
(mm)"; echo " D
(ft)"; echo " ψ"; echo " φ_*/2"; echo " 1200 D^3/2"; echo " i_b"; echo " i_bq_b"; echo " 43.2 w"; echo " i_BQ_B
(T/D)"; echo "
% Size
Q'_s"; echo " Q'_s
(T/D)"; echo "
Mult."; echo " Z'"; echo " A\""; echo " J'₁"; echo " -J'₂"; echo " J\"₁"; echo " -J\"₂"; echo " PJ\"₁+J\"₂
________
PJ'₁+J'₂ "; echo "
I\"₁"; echo " -I\"₂"; echo " PI\"₁+I\"₂+1"; echo " Computed
load
(T/D)"; echo "
1"; echo " $sizelowlong[0]-$sizehighlong[0]"; echo " " . number_format($sizediamlongmm[0],4); echo " " . number_format($sizediamlongft[0],6); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo "
2"; echo " $sizelowlong[1]-$sizehighlong[1]"; echo " " . number_format($sizediamlongmm[1],4); echo " " . number_format($sizediamlongft[1],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo "
3"; echo " $sizelowlong[2]-$sizehighlong[2]"; echo " " . number_format($sizediamlongmm[2],4); echo " " . number_format($sizediamlongft[2],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo "
4"; echo " $sizelowlong[3]-$sizehighlong[3]"; echo " " . number_format($sizediamlongmm[3],3); echo " " . number_format($sizediamlongft[3],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo "
5"; echo " $sizelowlong[4]-$sizehighlong[4]"; echo " " . number_format($sizediamlongmm[4],3); echo " " . number_format($sizediamlongft[4],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo "
6"; echo " $sizelowlong[5]-$sizehighlong[5]"; echo " " . number_format($sizediamlongmm[5],3); echo " " . number_format($sizediamlongft[5],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo "
7"; echo " $sizelowlong[6]-$sizehighlong[6]"; echo " " . number_format($sizediamlongmm[6],2); echo " " . number_format($sizediamlongft[6],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo "
8"; echo " $sizelowlong[7]-$sizehighlong[7]"; echo " " . number_format($sizediamlongmm[7],2); echo " " . number_format($sizediamlongft[7],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo "
9"; echo " $sizelowlong[8]-$sizehighlong[8]"; echo " " . number_format($sizediamlongmm[8],2); echo " " . number_format($sizediamlongft[8],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo "
10"; echo " $sizelowlong[9]-$sizehighlong[9]"; echo " " . number_format($sizediamlongmm[9],2); echo " " . number_format($sizediamlongft[9],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo "
-"; echo " Total"; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo "
"; echo "
"; } elseif($_POST["Submit"] == 'Calculate') { //echo "sizefractype = " . $sizefractype . "
"; echo "
"; if($sizefractype == 0) { echo ""; echo "
"; echo "Stream or river: " . "  " . $alpha . " " . "
" ."
"; if($unitsystems == 0) { echo "
" . "•   " . "U.S. Customary Units." . "   •"; } elseif($unitsystems == 1) { echo "
" . "•   " . "SI (Metric) Units." . "   •"; } echo "   " . "Number of size fractions: Normal (9)" . "   •"; if($zprimemode == 0) { echo "   " . "UBSR (1955) Z' mode of computation" . "   •"; } if($zprimemode == 1) { echo "   " . "USBR (1966) Z' mode of computation " . "   •"; } echo "   " . "Computed load (Column 24) in T/D (tons/day)" . "   •"; if($unitsystems == 1) { echo "   " . "Computed load (Column 25) in Metric tons/day" . "   •"; } echo "
" . "•   " . "RS = " . " " . round($sr,6) . "   "; echo "•   " . "u'_* = " . " " . round($ups,4) . "   "; echo "•   " . "δ = " . " " . round($delta,5) . "   "; echo "•   " . "k_s/δ = " . " " . round($ksovdelta,2) . "   "; echo "•   " . "x = " . " " . round($xxx,2) . "   "; echo "•   " . "P = " . " " . round($ppp,1) . "   "; echo "•   " . "A' = " . " " . round($aprime,3) . "   •" . "
"; echo "
" . "•   " . "Percentage of flow in sampled zone = " . " " . $pfsampledw . "%" . "   "; if($unitsystems == 0) { echo "•   " . "Sediment discharge through the sampled zone Q'_s = " . " " . $qsprimew2 . " " . $loadunits . "   •" . "
"; } elseif($unitsystems == 1) { echo "•   " . "Sediment discharge through the sampled zone Q'_s = " . " " . $qsprimew2 . " " . $loadunits . " = " . number_format(($qsprime*0.9072),0) . " " . $loadunitsout . "   •" . "
"; } echo "•   " . "Reference size fraction D = " . number_format($sizediamappft[$qsprimelocmax],5) . " ft.    •   "; if($unitsystems == 1) { echo "Reference size fraction D = " . number_format($sizediamappmm[$qsprimelocmax],3) . " mm.    •   "; } echo "Reference size fraction Z' = " . number_format($zprimeout[$qsprimelocmax],2) . "   •" . "
"; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo "
0"; echo " 1"; echo " 2"; echo " 3"; echo " 4"; echo " 5"; echo " 6"; echo " 7"; echo " 8"; echo " 9"; echo " 10"; echo " 11"; echo " 12"; echo " 13"; echo " 14"; echo " 15"; echo " 16"; echo " 17"; echo " 18"; echo " 19"; echo " 20"; echo " 21"; echo " 22"; echo " 23"; echo " 24"; echo " 25"; echo "
i"; echo " Size
fraction
(mm)"; echo " D
(mm)"; echo " D
(ft)"; echo " ψ"; echo " φ_*/2"; echo " 1200 D^3/2"; echo " i_b"; echo " i_bq_b"; echo " 43.2 w"; echo " i_BQ_B
(T/D)"; echo "
% Size
Q'_s"; echo " Q'_s
(T/D)"; echo "
Mult."; echo " Z'"; echo " A\""; echo " J'₁"; echo " -J'₂"; echo " J\"₁"; echo " -J\"₂"; echo " PJ\"₁+J\"₂
________
PJ'₁+J'₂ "; echo "
I\"₁"; echo " -I\"₂"; echo " PI\"₁+I\"₂
+1"; echo "
Computed
load
(T/D)"; echo "
Computed
load
(Metric tons/day)"; echo "
1"; echo " $sizelownormal[0]-$sizehighnormal[0]"; echo " " . number_format($sizediamnormalmm[0],4); echo " " . number_format($sizediamnormalft[0],6); echo " " . number_format($psi[0],2) . ""; echo " " . number_format($phistar2[0],4) . ""; echo " " . number_format($d1200[0],5) . ""; echo " $ibvar[0]"; echo " " . number_format($ibqb[0],6); echo ""; echo " " . number_format($w432,0); echo ""; echo " " . number_format($ibqbcap[0],2); echo ""; echo " " . number_format($qspvar[0],1); echo " "; echo " " . number_format($qsprimei[0],1); echo " "; echo " " . number_format($multiplier[0],3); echo " "; echo " " . number_format($zprimeout[0],3); echo " "; echo " " . number_format($adoubleprime[0],6); echo " "; echo " " . number_format($j1primecol16[0],2); echo " "; echo " " . number_format($j2primecol17[0],2); echo " "; echo " " . number_format($j1doubleprimecol18[0],2); echo " "; echo " " . number_format($j2doubleprimecol19[0],2); echo " "; echo " " . number_format($column20[0],2); echo " "; echo " "; echo " "; echo " "; echo " " . number_format($column24[0],1); echo " "; echo " " . number_format($column25[0],1); echo " "; echo "
2"; echo " $sizelownormal[1]-$sizehighnormal[1]"; echo " " . number_format($sizediamnormalmm[1],4); echo " " . number_format($sizediamnormalft[1],5); echo " " . number_format($psi[1],2) . ""; echo " " . number_format($phistar2[1],4) . ""; echo " " . number_format($d1200[1],5) . ""; echo " $ibvar[1]"; echo " " . number_format($ibqb[1],6); echo ""; echo " " . number_format($w432,0); echo ""; echo " " . number_format($ibqbcap[1],2); echo ""; echo " " . number_format($qspvar[1],1); echo " "; echo " " . number_format($qsprimei[1],1); echo " "; echo " " . number_format($multiplier[1],3); echo " "; echo " " . number_format($zprimeout[1],2); echo " "; echo " " . number_format($adoubleprime[1],5); echo " "; echo " " . number_format($j1primecol16[1],2); echo " "; echo " " . number_format($j2primecol17[1],2); echo " "; echo " " . number_format($j1doubleprimecol18[1],2); echo " "; echo " " . number_format($j2doubleprimecol19[1],2); echo " "; echo " " . number_format($column20[1],2); echo " "; echo " "; echo " "; echo " "; echo " " . number_format($column24[1],1); echo " "; echo " " . number_format($column25[1],1); echo " "; echo "
3"; echo " $sizelownormal[2]-$sizehighnormal[2]"; echo " " . number_format($sizediamnormalmm[2],3); echo " " . number_format($sizediamnormalft[2],5); echo " " . number_format($psi[2],2) . ""; echo " " . number_format($phistar2[2],4) . ""; echo " " . number_format($d1200[2],4) . ""; echo " $ibvar[2]"; echo " " . number_format($ibqb[2],6); echo ""; echo " " . number_format($w432,0); echo ""; echo " " . number_format($ibqbcap[2],2); echo ""; echo " " . number_format($qspvar[2],1); echo " "; echo " " . number_format($qsprimei[2],1); echo " "; echo " " . number_format($multiplier[2],3); echo " "; echo " " . number_format($zprimeout[2],2); echo " "; echo " " . number_format($adoubleprime[2],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " " . number_format($i1doubleprimecol21[2],2); echo " "; echo " " . number_format($i2doubleprimecol22[2],2); echo " "; echo " " . number_format($column23[2],2); echo " "; echo " " . number_format($column24[2],1); echo " "; echo " " . number_format($column25[2],1); echo " "; echo "
4"; echo " $sizelownormal[3]-$sizehighnormal[3]"; echo " " . number_format($sizediamnormalmm[3],3); echo " " . number_format($sizediamnormalft[3],5); echo " " . number_format($psi[3],2) . ""; echo " " . number_format($phistar2[3],4) . ""; echo " " . number_format($d1200[3],4) . ""; echo " $ibvar[3]"; echo " " . number_format($ibqb[3],6); echo ""; echo " " . number_format($w432,0); echo ""; echo " " . number_format($ibqbcap[3],2); echo ""; echo " " . number_format($qspvar[3],1); echo " "; echo " " . number_format($qsprimei[3],1); echo " "; echo " " . number_format($multiplier[3],2); echo " "; echo " " . number_format($zprimeout[3],2); echo " "; echo " " . number_format($adoubleprime[3],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " " . number_format($i1doubleprimecol21[3],2); echo " "; echo " " . number_format($i2doubleprimecol22[3],2); echo " "; echo " " . number_format($column23[3],2); echo " "; echo " " . number_format($column24[3],1); echo " "; echo " " . number_format($column25[3],1); echo " "; echo "
5"; echo " $sizelownormal[4]-$sizehighnormal[4]"; echo " " . number_format($sizediamnormalmm[4],3); echo " " . number_format($sizediamnormalft[4],5); echo " " . number_format($psi[4],2) . ""; echo " " . number_format($phistar2[4],4) . ""; echo " " . number_format($d1200[4],4) . ""; echo " $ibvar[4]"; echo " " . number_format($ibqb[4],6); echo ""; echo " " . number_format($w432,0); echo ""; echo " " . number_format($ibqbcap[4],2); echo ""; echo " " . number_format($qspvar[4],1); echo " "; echo " " . number_format($qsprimei[4],1); echo " "; echo " " . number_format($multiplier[4],2); echo " "; echo " " . number_format($zprimeout[4],2); echo " "; echo " " . number_format($adoubleprime[4],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " " . number_format($i1doubleprimecol21[4],2); echo " "; echo " " . number_format($i2doubleprimecol22[4],2); echo " "; echo " " . number_format($column23[4],2); echo " "; echo " " . number_format($column24[4],1); echo " "; echo " " . number_format($column25[4],1); echo " "; echo "
6"; echo " $sizelownormal[5]-$sizehighnormal[5]"; echo " " . number_format($sizediamnormalmm[5],2); echo " " . number_format($sizediamnormalft[5],5); echo " " . number_format($psi[5],2) . ""; echo " " . number_format($phistar2[5],4) . ""; echo " " . number_format($d1200[5],4) . ""; echo " $ibvar[5]"; echo " " . number_format($ibqb[5],6); echo ""; echo " " . number_format($w432,0); echo ""; echo " " . number_format($ibqbcap[5],2); echo ""; echo " " . number_format($qspvar[5],1); echo " "; echo " " . number_format($qsprimei[5],1); echo " "; echo " " . number_format($multiplier[5],2); echo " "; echo " " . number_format($zprimeout[5],2); echo " "; echo " " . number_format($adoubleprime[5],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " " . number_format($i1doubleprimecol21[5],2); echo " "; echo " " . number_format($i2doubleprimecol22[5],2); echo " "; echo " " . number_format($column23[5],2); echo " "; echo " " . number_format($column24[5],1); echo " "; echo " " . number_format($column25[5],1); echo " "; echo "
7"; echo " $sizelownormal[6]-$sizehighnormal[6]"; echo " " . number_format($sizediamnormalmm[6],2); echo " " . number_format($sizediamnormalft[6],5); echo " " . number_format($psi[6],2) . ""; echo " " . number_format($phistar2[6],4) . ""; echo " " . number_format($d1200[6],4) . ""; echo " $ibvar[6]"; echo " " . number_format($ibqb[6],6); echo ""; echo " " . number_format($w432,0); echo ""; echo " " . number_format($ibqbcap[6],2); echo ""; echo " " . number_format($qspvar[6],1); echo " "; echo " " . number_format($qsprimei[6],1); echo " "; echo " " . number_format($multiplier[6],2); echo " "; echo " " . number_format($zprimeout[6],2); echo " "; echo " " . number_format($adoubleprime[6],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " " . number_format($i1doubleprimecol21[6],2); echo " "; echo " " . number_format($i2doubleprimecol22[6],2); echo " "; echo " " . number_format($column23[6],2); echo " "; echo " " . number_format($column24[6],1); echo " "; echo " " . number_format($column25[6],1); echo " "; echo "
8"; echo " $sizelownormal[7]-$sizehighnormal[7]"; echo " " . number_format($sizediamnormalmm[7],2); echo " " . number_format($sizediamnormalft[7],5); echo " " . number_format($psi[7],2) . ""; echo " " . number_format($phistar2[7],4) . ""; echo " " . number_format($d1200[7],4) . ""; echo " $ibvar[7]"; echo " " . number_format($ibqb[7],6); echo ""; echo " " . number_format($w432,0); echo ""; echo " " . number_format($ibqbcap[7],2); echo ""; echo " " . number_format($qspvar[7],1); echo " "; echo " " . number_format($qsprimei[7],1); echo " "; echo " " . number_format($multiplier[7],2); echo " "; echo " " . number_format($zprimeout[7],2); echo " "; echo " " . number_format($adoubleprime[7],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " " . number_format($i1doubleprimecol21[7],2); echo " "; echo " " . number_format($i2doubleprimecol22[7],2); echo " "; echo " " . number_format($column23[7],2); echo " "; echo " " . number_format($column24[7],1); echo " "; echo " " . number_format($column25[7],1); echo " "; echo "
9"; echo " $sizelownormal[8]-$sizehighnormal[8]"; echo " " . number_format($sizediamnormalmm[8],2); echo " " . number_format($sizediamnormalft[8],5); echo " " . number_format($psi[8],2) . ""; echo " " . number_format($phistar2[8],4) . ""; echo " " . number_format($d1200[8],4) . ""; echo " $ibvar[8]"; echo " " . number_format($ibqb[8],6); echo ""; echo " " . number_format($w432,0); echo ""; echo " " . number_format($ibqbcap[8],2); echo ""; echo " " . number_format($qspvar[8],1); echo " "; echo " " . number_format($qsprimei[8],1); echo " "; echo " " . number_format($multiplier[8],2); echo " "; echo " " . number_format($zprimeout[8],2); echo " "; echo " " . number_format($adoubleprime[8],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " " . number_format($i1doubleprimecol21[8],2); echo " "; echo " " . number_format($i2doubleprimecol22[8],2); echo " "; echo " " . number_format($column23[8],2); echo " "; echo " " . number_format($column24[8],1); echo " "; echo " " . number_format($column25[8],1); echo " "; echo "
-"; echo " Total"; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " 100"; echo " " . number_format($qsprime,0) .""; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " " . number_format($sumload,1); echo " "; echo " " . number_format($sumload25,1); echo " "; echo "
"; echo "
"; } elseif($sizefractype == 1) { echo ""; echo "
"; echo "Stream or river: " . "  " . $alpha . " " . "
" ."
"; if($unitsystems == 0) { echo "
" . "•   " . "U.S. Customary Units." . "   •"; } elseif($unitsystems == 1) { echo "
" . "•   " . "SI (Metric) Units." . "   •"; } echo "   " . "Number of size fractions: Long (10)" . "   •"; if($zprimemode == 0) { echo "   " . "UBSR (1955) Z' mode of computation" . "   •"; } if($zprimemode == 1) { echo "   " . "USBR (1966) Z' mode of computation " . "   •"; } echo "   " . "Computed load (Column 24) in T/D (tons/day)" . "   •"; if($unitsystems == 1) { echo "   " . "Computed load (Column 25) in Metric tons/day" . "   •"; } echo "
" . "•   " . "RS = " . " " . round($sr,6) . "   "; echo "•   " . "u'_* = " . " " . round($ups,4) . "   "; echo "•   " . "δ = " . " " . round($delta,5) . "   "; echo "•   " . "k_s/δ = " . " " . round($ksovdelta,2) . "   "; echo "•   " . "x = " . " " . round($xxx,2) . "   "; echo "•   " . "P = " . " " . round($ppp,1) . "   "; echo "•   " . "A' = " . " " . round($aprime,3) . "   •" . "
"; echo "
" . "•   " . "Percentage of flow in sampled zone = " . " " . $pfsampledw . "%" . "   "; if($unitsystems == 0) { echo "•   " . "Sediment discharge through the sampled zone Q'_s = " . " " . $qsprimew2 . " " . $loadunits . "   •" . "
"; } elseif($unitsystems == 1) { echo "•   " . "Sediment discharge through the sampled zone Q'_s = " . " " . $qsprimew2 . " " . $loadunits . " = " . number_format(($qsprime*0.9072),0) . " " . $loadunitsout . "   •" . "
"; } echo "•   " . "Reference size fraction D = " . number_format($sizediamappft[$qsprimelocmax],5) . " ft.    •   "; if($unitsystems == 1) { echo "Reference size fraction D = " . number_format($sizediamappmm[$qsprimelocmax],3) . " mm.    •   "; } echo "Reference size fraction Z' = " . number_format($zprimeout[$qsprimelocmax],2) . "   •" . "
"; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; echo ""; if($unitsystems == 1) { echo ""; } echo ""; echo "
0"; echo " 1"; echo " 2"; echo " 3"; echo " 4"; echo " 5"; echo " 6"; echo " 7"; echo " 8"; echo " 9"; echo " 10"; echo " 11"; echo " 12"; echo " 13"; echo " 14"; echo " 15"; echo " 16"; echo " 17"; echo " 18"; echo " 19"; echo " 20"; echo " 21"; echo " 22"; echo " 23"; echo " 24"; echo " 25"; echo "
i"; echo " Size
fraction
(mm)"; echo " D
(mm)"; echo " D
(ft)"; echo " ψ"; echo " φ_*/2"; echo " 1200 D^3/2"; echo " i_b"; echo " i_bq_b"; echo " 43.2 w"; echo " i_BQ_B
(T/D)"; echo "
% Size
Q'_s"; echo " Q'_s
(T/D)"; echo "
Mult."; echo " Z'"; echo " A\""; echo " J'₁"; echo " -J'₂"; echo " J\"₁"; echo " -J\"₂"; echo " PJ\"₁+J\"₂
________
PJ'₁+J'₂ "; echo "
I\"₁"; echo " -I\"₂"; echo " PI\"₁+I\"₂
+1"; echo "
Computed
load
(T/D)"; echo "
Computed
load
(Metric tons/day)"; echo "
1"; echo " $sizelowlong[0]-$sizehighlong[0]"; echo " " . number_format($sizediamlongmm[0],4); echo " " . number_format($sizediamlongft[0],6); echo " " . number_format($psi[0],2) . ""; echo " " . number_format($phistar2[0],4) . ""; echo " " . number_format($d1200[0],5) . ""; echo " $ibvar[0]"; echo " " . number_format($ibqb[0],6); echo ""; echo " " . number_format($w432,0); echo ""; echo " " . number_format($ibqbcap[0],2); echo ""; echo " " . number_format($qspvar[0],1); echo " "; echo " " . number_format($qsprimei[0],1); echo " "; echo " " . number_format($multiplier[0],3); echo " "; echo " " . number_format($zprimeout[0],3); echo " "; echo " " . number_format($adoubleprime[0],6); echo " "; echo " " . number_format($j1primecol16[0],2); echo " "; echo " " . number_format($j2primecol17[0],2); echo " "; echo " " . number_format($j1doubleprimecol18[0],2); echo " "; echo " " . number_format($j2doubleprimecol19[0],2); echo " "; echo " " . number_format($column20[0],2); echo " "; echo " "; echo " "; echo " "; echo " " . number_format($column24[0],1); echo " "; echo " " . number_format($column25[0],1); echo " "; echo "
2"; echo " $sizelowlong[1]-$sizehighlong[1]"; echo " " . number_format($sizediamlongmm[1],4); echo " " . number_format($sizediamlongft[1],5); echo " " . number_format($psi[1],2) . ""; echo " " . number_format($phistar2[1],4) . ""; echo " " . number_format($d1200[1],5) . ""; echo " $ibvar[1]"; echo " " . number_format($ibqb[1],6); echo ""; echo " " . number_format($w432,0); echo ""; echo " " . number_format($ibqbcap[1],2); echo ""; echo " " . number_format($qspvar[1],1); echo " "; echo " " . number_format($qsprimei[1],1); echo " "; echo " " . number_format($multiplier[1],3); echo " "; echo " " . number_format($zprimeout[1],2); echo " "; echo " " . number_format($adoubleprime[1],5); echo " "; echo " " . number_format($j1primecol16[1],2); echo " "; echo " " . number_format($j2primecol17[1],2); echo " "; echo " " . number_format($j1doubleprimecol18[1],2); echo " "; echo " " . number_format($j2doubleprimecol19[1],2); echo " "; echo " " . number_format($column20[1],2); echo " "; echo " "; echo " "; echo " "; echo " " . number_format($column24[1],1); echo " "; echo " " . number_format($column25[1],1); echo " "; echo "
3"; echo " $sizelowlong[2]-$sizehighlong[2]"; echo " " . number_format($sizediamlongmm[2],3); echo " " . number_format($sizediamlongft[2],5); echo " " . number_format($psi[2],2) . ""; echo " " . number_format($phistar2[2],4) . ""; echo " " . number_format($d1200[2],4) . ""; echo " $ibvar[2]"; echo " " . number_format($ibqb[2],6); echo ""; echo " " . number_format($w432,0); echo ""; echo " " . number_format($ibqbcap[2],2); echo ""; echo " " . number_format($qspvar[2],1); echo " "; echo " " . number_format($qsprimei[2],1); echo " "; echo " " . number_format($multiplier[2],3); echo " "; echo " " . number_format($zprimeout[2],2); echo " "; echo " " . number_format($adoubleprime[2],5); echo " "; echo " " . number_format($j1primecol16[2],2); echo " "; echo " " . number_format($j2primecol17[2],2); echo " "; echo " " . number_format($j1doubleprimecol18[2],2); echo " "; echo " " . number_format($j2doubleprimecol19[2],2); echo " "; echo " " . number_format($column20[2],2); echo " "; echo " "; echo " "; echo " "; echo " " . number_format($column24[2],1); echo " "; echo " " . number_format($column25[2],1); echo " "; echo "
4"; echo " $sizelowlong[3]-$sizehighlong[3]"; echo " " . number_format($sizediamlongmm[3],3); echo " " . number_format($sizediamlongft[3],5); echo " " . number_format($psi[3],2) . ""; echo " " . number_format($phistar2[3],4) . ""; echo " " . number_format($d1200[3],4) . ""; echo " $ibvar[3]"; echo " " . number_format($ibqb[3],6); echo ""; echo " " . number_format($w432,0); echo ""; echo " " . number_format($ibqbcap[3],2); echo ""; echo " " . number_format($qspvar[3],1); echo " "; echo " " . number_format($qsprimei[3],1); echo " "; echo " " . number_format($multiplier[3],2); echo " "; echo " " . number_format($zprimeout[3],2); echo " "; echo " " . number_format($adoubleprime[3],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " " . number_format($i1doubleprimecol21[3],2); echo " "; echo " " . number_format($i2doubleprimecol22[3],2); echo " "; echo " " . number_format($column23[3],2); echo " "; echo " " . number_format($column24[3],1); echo " "; echo " " . number_format($column25[3],1); echo " "; echo "
5"; echo " $sizelowlong[4]-$sizehighlong[4]"; echo " " . number_format($sizediamlongmm[4],3); echo " " . number_format($sizediamlongft[4],5); echo " " . number_format($psi[4],2) . ""; echo " " . number_format($phistar2[4],4) . ""; echo " " . number_format($d1200[4],4) . ""; echo " $ibvar[4]"; echo " " . number_format($ibqb[4],6); echo ""; echo " " . number_format($w432,0); echo ""; echo " " . number_format($ibqbcap[4],2); echo ""; echo " " . number_format($qspvar[4],1); echo " "; echo " " . number_format($qsprimei[4],1); echo " "; echo " " . number_format($multiplier[4],2); echo " "; echo " " . number_format($zprimeout[4],2); echo " "; echo " " . number_format($adoubleprime[4],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " " . number_format($i1doubleprimecol21[4],2); echo " "; echo " " . number_format($i2doubleprimecol22[4],2); echo " "; echo " " . number_format($column23[4],2); echo " "; echo " " . number_format($column24[4],1); echo " "; echo " " . number_format($column25[4],1); echo " "; echo "
6"; echo " $sizelowlong[5]-$sizehighlong[5]"; echo " " . number_format($sizediamlongmm[5],2); echo " " . number_format($sizediamlongft[5],5); echo " " . number_format($psi[5],2) . ""; echo " " . number_format($phistar2[5],4) . ""; echo " " . number_format($d1200[5],4) . ""; echo " $ibvar[5]"; echo " " . number_format($ibqb[5],6); echo ""; echo " " . number_format($w432,0); echo ""; echo " " . number_format($ibqbcap[5],2); echo ""; echo " " . number_format($qspvar[5],1); echo " "; echo " " . number_format($qsprimei[5],1); echo " "; echo " " . number_format($multiplier[5],2); echo " "; echo " " . number_format($zprimeout[5],2); echo " "; echo " " . number_format($adoubleprime[5],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " " . number_format($i1doubleprimecol21[5],2); echo " "; echo " " . number_format($i2doubleprimecol22[5],2); echo " "; echo " " . number_format($column23[5],2); echo " "; echo " " . number_format($column24[5],1); echo " "; echo " " . number_format($column25[5],1); echo " "; echo "
7"; echo " $sizelowlong[6]-$sizehighlong[6]"; echo " " . number_format($sizediamlongmm[6],2); echo " " . number_format($sizediamlongft[6],5); echo " " . number_format($psi[6],2) . ""; echo " " . number_format($phistar2[6],4) . ""; echo " " . number_format($d1200[6],4) . ""; echo " $ibvar[6]"; echo " " . number_format($ibqb[6],6); echo ""; echo " " . number_format($w432,0); echo ""; echo " " . number_format($ibqbcap[6],2); echo ""; echo " " . number_format($qspvar[6],1); echo " "; echo " " . number_format($qsprimei[6],1); echo " "; echo " " . number_format($multiplier[6],2); echo " "; echo " " . number_format($zprimeout[6],2); echo " "; echo " " . number_format($adoubleprime[6],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " " . number_format($i1doubleprimecol21[6],2); echo " "; echo " " . number_format($i2doubleprimecol22[6],2); echo " "; echo " " . number_format($column23[6],2); echo " "; echo " " . number_format($column24[6],1); echo " "; echo " " . number_format($column25[6],1); echo " "; echo "
8"; echo " $sizelowlong[7]-$sizehighlong[7]"; echo " " . number_format($sizediamlongmm[7],2); echo " " . number_format($sizediamlongft[7],5); echo " " . number_format($psi[7],2) . ""; echo " " . number_format($phistar2[7],4) . ""; echo " " . number_format($d1200[7],4) . ""; echo " $ibvar[7]"; echo " " . number_format($ibqb[7],6); echo ""; echo " " . number_format($w432,0); echo ""; echo " " . number_format($ibqbcap[7],2); echo ""; echo " " . number_format($qspvar[7],1); echo " "; echo " " . number_format($qsprimei[7],1); echo " "; echo " " . number_format($multiplier[7],2); echo " "; echo " " . number_format($zprimeout[7],2); echo " "; echo " " . number_format($adoubleprime[7],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " " . number_format($i1doubleprimecol21[7],2); echo " "; echo " " . number_format($i2doubleprimecol22[7],2); echo " "; echo " " . number_format($column23[7],2); echo " "; echo " " . number_format($column24[7],1); echo " "; echo " " . number_format($column25[7],1); echo " "; echo "
9"; echo " $sizelowlong[8]-$sizehighlong[8]"; echo " " . number_format($sizediamlongmm[8],2); echo " " . number_format($sizediamlongft[8],5); echo " " . number_format($psi[8],2) . ""; echo " " . number_format($phistar2[8],4) . ""; echo " " . number_format($d1200[8],4) . ""; echo " $ibvar[8]"; echo " " . number_format($ibqb[8],6); echo ""; echo " " . number_format($w432,0); echo ""; echo " " . number_format($ibqbcap[8],2); echo ""; echo " " . number_format($qspvar[8],1); echo " "; echo " " . number_format($qsprimei[8],1); echo " "; echo " " . number_format($multiplier[8],2); echo " "; echo " " . number_format($zprimeout[8],2); echo " "; echo " " . number_format($adoubleprime[8],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " " . number_format($i1doubleprimecol21[8],2); echo " "; echo " " . number_format($i2doubleprimecol22[8],2); echo " "; echo " " . number_format($column23[8],2); echo " "; echo " " . number_format($column24[8],1); echo " "; echo " " . number_format($column25[8],1); echo " "; echo "
10"; echo " $sizelowlong[9]-$sizehighlong[9]"; echo " " . number_format($sizediamlongmm[9],2); echo " " . number_format($sizediamlongft[9],5); echo " " . number_format($psi[9],2) . ""; echo " " . number_format($phistar2[9],4) . ""; echo " " . number_format($d1200[9],4) . ""; echo " $ibvar[9]"; echo " " . number_format($ibqb[9],6); echo ""; echo " " . number_format($w432,0); echo ""; echo " " . number_format($ibqbcap[9],2); echo ""; echo " " . number_format($qspvar[9],1); echo " "; echo " " . number_format($qsprimei[9],1); echo " "; echo " " . number_format($multiplier[9],2); echo " "; echo " " . number_format($zprimeout[9],2); echo " "; echo " " . number_format($adoubleprime[9],5); echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " " . number_format($i1doubleprimecol21[9],2); echo " "; echo " " . number_format($i2doubleprimecol22[9],2); echo " "; echo " " . number_format($column23[9],2); echo " "; echo " " . number_format($column24[9],1); echo " "; echo " " . number_format($column25[9],1); echo " "; echo "
-"; echo " Total"; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " 100"; echo " " . number_format($qsprime,0) . ""; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " "; echo " " . number_format($sumload,1); echo " "; echo " " . number_format($sumload25,1); echo " "; echo "
"; echo "
"; } echo "
"; } ?>



Your request was processed at "; echo " "; echo date('h:i:s a \o\n F jS, Y'); ?>    [ ].

Thank you for running online_modified_einstein. [080808]

online calc

normal depth critical depth normal and critical depth discharge in channel critical slope

normal depth in culvert critical depth in culvert discharge in culvert discharge sluice discharge weir

M1 wsprofile M2 wsprofile M3 wsprofile S1 wsprofile S2 wsprofile S3 wsprofile

C1 wsprofile H2 wsprofile A2 wsprofile C3 wsprofile H3 wsprofile A3 wsprofile

sequent depth HJ energy loss HJ initial sequent HJ efficiency HJ critical width constriction

ogee spillway Hazen-Williams parallel pipes three reservoirs tractive force

V-notch weir V-notch partially contracted Cipolletti weir Rectangular weir Standard contracted rectangular weir Standard suppressed rectangular weir

Froude number Vedernikov number Limiting contracting ratio

Creager rational slope-area linear reservoir storage indication 1 storage indication 2

Muskingum Muskingum-Cunge time-area Clark UH Cascade of linear reservoirs

USGS Methods for magnitude of floods in California Kinematic wave applicability Diffusion wave applicability Clark's unit hydrography compared to Ponce's version Correlation coefficient of a joint probability distribution Storage volume of a detention basin

Blaney-Criddle Penman Penman-Monteith reference crop Thornthwaite Priestley-Taylor Penman-Monteith ecosystems
Gumbel Gumbel 2 Log Pearson Log Pearson 2 TR-55 graphical curve number

Overland flow using the diffusion wave method Dynamic hydraulic diffusivity convolution S-hydrograph time of concentration water balance

UH cascade dimensionless UH cascade general UH cascade series UH cascade all series UH cascade

one-predictor linear one-predictor nonlinear two-predictor linear two-predictor nonlinear hyperbolic regression

fall velocity Lane & Koelzer USLE USLE2 Dendy-Bolton Shields

Duboys Meyer-Peter Colby 1957 Colby reservoir design life Equilibrium channel top width using the Lane et al. theory

Modified Einstein Procedure bridge scour using Melville equation

DO sag DO sag analysis Oxygenation Salinity (EC to TDS)