Determining the time increment

The first action within the major loop is the determination of the time increment in subroutine compdt. The formulas for doing so are Formulas (475) and (477) for liquids and gases, respectively. Notice that the solution influences the time increment, so the increment has to be recalculated at the start of the major loop. For steady state calculations the time step is only recalculated in iterations $2^n$, $n \in \mathbb{N}$, since the change in solution decreases the closer steady state. This accelerates the steady state calculations.