[AG1994]

Hello,
i'm not sure if I am in the right chat to ask my question, but I have problems with programming an low mach scheme for the euler equation without sourceterm. I am writing a python code to approximate the euler equations in all mach systems. The approximate Riemann solver I am using is the suliciu relaxation model with pressure splitting. My problem lies with in the scaling of the variables, I don't know how to scale the pressure terms. I am using the dimensionless form of the euler equation for the flux computation and use the initial data of the sod shock tube problem. I only get an ok solution if i multiply my flux function with the mach number during the update of the state vector for the next time step, but I know that is not the right way to scale the variables. Does anybody have an idea how i can implement the scaling of the variables into my finite volume method?
thank you in advance