We consider the isothermal Euler model with non-vacuum initial data. We extract the Riemann invariants of the isothermal Euler model, which admits vital applications. We also design the modified Rusanov (mR) scheme to solve the isothermal Euler model. This scheme consists of two steps, the first step of the scheme depends on a local parameter allowing to control diffusion. The second stage recovers conservation equation. This technique is a straightforward to implement and precise. We compare this scheme with the Rusanov scheme via three numerical examples. This numerical study verifies the efficiency of the mR scheme. Finally, the mR scheme can be used to solve many other models in applied science.