Step 4: determining the energy correction

In step 4 the correction to the energy is determined. The lumped equations have already been solved immediately after the call to mafillsmv1rhs.f and the solution was stored in field v(1..nk,0). Step 4 is only performed for gases (which exhibit a strong coupling of density, temperature and pressure) and for liquids for which initial temperature conditions have been defined. In general the coupling in liquids is rather weak.