//Font settings
mathcolor="Black";
mathfontsize="1.1em";
mathfontfamily = "Arial";  
var mathfontweight = "bold";  

//General
define("!<=","\u2270");
define("!>=","\u2271");
define("undefined","O/");
define("dgr","\u00BA");
define("tns","(t ons)");

//Mathematics
define("DDX","d/(dx)");
define("DDY","d/(dy)");
define("DDZ","d/(dz)");
define("DDX2","d^2/(dx^2)");
define("DDY2","d^2/(dy^2)");
define("DDZ2","d^2/(dz^2)");
define("ddx","del/(delx)");
define("ddy","del/(dely)");
define("ddz","del/(delz)");
define("ddx2","del^2/(delx^2)");
define("ddy2","del^2/(dely^2)");
define("ddz2","del^2/(delz^2)");


//Dredging general
define("concentration","(rho_m-rho_w)/(rho_q-rho_w)");
define("Cv","(rho_m-rho_w)/(rho_q-rho_w)");


//Dredging cutting
define("Fhnc","(c_1*rho_w*g*h_i^2*v_c*w*v_d)/k_m");
define("Fvnc","(c_2*rho_w*g*h_i^2*v_c*w*v_d)/k_m");
define("Fhc","d_1*rho_w*g*(z+10)*h_i*w");
define("Fvc","d_2*rho_w*g*(z+10)*h_i*w");

//TSHD sedimentation
define("TSE1","eta_t=eta_g^0*(1-.184*eta_g^(+.885-.20*eta_g)*(1-tanh(eta_g^(-.13-.80*eta_g)*(log(v_s/s_o)-.2614-.5*log(lambda)+eta_g^(-.33-.94*eta_g)))))");
define("TSE2","eta_t=eta_g^(-1)*(1-.184*eta_g^(-.69-.38*eta_g)*(1-tanh(eta_g^(+.77-.08*eta_g)*(log(v_s/s_o)-.2614-.5*log(lambda)+eta_g^(+1.01-.18*eta_g)))))");
define("HLP","v_o=(Q/(L*W))");
define("FlowV","s_o=((Q)/(W*H)");
define("Stokes","v_s=424*R_d*d^2");
define("Budryck","v_s=8.925*((sqrt(1+95*R_d*d^3)-1)/(d))");
define("Rittinger","v_s=87*sqrt(R_d*d)");
define("GrainE","eta_g=(v_s/v_o)");
define("SwameeJain","lambda=0.25/(log(d/(3.7*D_h)+5.75/Re^0.9))^2");
define("Reynolds","Re=(v*D_h)/nu");
define("Dhr","D_h=(4*A)/O=(4*H*W)/(2*H+W)~~4*H"); 
define("Dhp","D_h=(4*A)/O=(4*(pi/4*D^2))/(pi*D)=D"); 

//Turbulence Settling Efficiency
function TSE(vsvo,vsso,labda) {
  var etat;
  if (vsvo<=1) {
    etat = (1-0.184*pow(vsvo,0.885-0.20*vsvo)*(1-tanh(pow(vsvo,-0.13-0.88*vsvo)*(log(vsso)-0.2614-0.5*log(labda)+pow(vsvo,-0.33-0.94*vsvo)))))*vsvo;
  } else {
    etat = (1-0.184*pow(vsvo,-0.690-0.38*vsvo)*(1-tanh(pow(vsvo,+0.77-0.08*vsvo)*(log(vsso)-0.2614-0.5*log(labda)+pow(vsvo,+1.01-0.18*vsvo)))));
  }
  return etat;
}
var rho_w = 1.025;
var rho_q = 2.650;
var Rd = (rho_q-rho_w)/rho_w; 
//Settling velocity according to Stokes
function Stokes(d) {
  return 424*Rd*pow(d,2);
}
//Settling velocity according to Budryck
function Budryck(d) {
  return 8.925*(sqrt(1+95*Rd*pow(d,3))-1)/d;
}
//Settling velocity according to Rittinger
function Rittinger(d) {
  return 87*sqrt(Rd*d);
}
//Swamee Jain equation for the Moody diagram, alert(SwameeJain(0.000001,1,10,0.000001));
function SwameeJain(d_r,D_h,Re) {
  return 0.25/sq(log(d_r/(3.7*D_h)+5.75/pow(Re,0.9)));
}