eulerstep_stag_com.mtotenthalpy = etotstar + pstar; totenthL = extrap(totenthalpy, limtype);rhoL = extrap(rhostar, limtype);etotstar(1) = etotold(1)...
eulerstep_stag_com.mtotenthalpy = etotstar + pstar; totenthL = extrap(totenthalpy, limtype);rhoL = extrap(rhostar, limtype);etotstar(1) = etotold(1)...