! nbs 1/24/22 compute headcut parameters compute travel time coefficients - delete when finished with flood plain check IF bottom width (b) is < 0 compute travel time coefficients - delete when finished with flood plain
compute rating curve set Muskingum parameters compute storage discharge for Muskingum at 0.1 and 1.0 times bankfull depth c^2=a^2+b^2 - a=dep; a/b=slope; b^2=a^2/slope^2 0.1bankfull storage discharge coef bankfull storage discharge coef Compute storage time constant for reach (msk_co1 + msk_co2 = 1.) Muskingum numerical stability -Jaehak Jeong, 2011 Discretize time interval to meet the stability criterion initial inflow-outflow compute coefficients c1+c2+c3 must equal 1 initialize rating curves - inflow and outflow at current time step initial volume is frac of flow depth - fracmmkm*1000. = m3
convert concentration to mass partition water between channel and flood plain depth below bankfull depth above bankfull save initial water if calibrating and rerunning calculate volume of active river bed sediment layer - m3 calculate mixing velocity using molecular weight and porosity