Commit 8ce3fd2c authored by Lluis Jofre Cruanyes's avatar Lluis Jofre Cruanyes
Browse files

MOVERS scheme modified III

parent 214d085b

src/FlowSolverRHEA.cpp

+3 −3
Original line number Diff line number Diff line
@@ -4551,9 +4551,9 @@ double EckepMoversRHFluxApproximateRiemannSolver::calculateIntercellFlux(const d 
    /// -------------------------------///



    /// Prevent sonic/entropy glitch 

    //if( alpha_S > epsilon ) {

    if( alpha_S > 1.0e-10 ) {						// ... modified for OpenACC

        double theta = 1.0;						/// !! theta needs to be larger than 0.0 !!

    //if( alpha_S < epsilon ) {

    if( alpha_S < 1.0e-10 ) {						// ... modified for OpenACC

        double theta = 0.1;						/// !! theta needs to be larger than 0.0 !!

        alpha_S = ( alpha_S*alpha_S + theta*theta )/( 2.0*theta );

    }