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 453 and 455 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.