The first action within the major loop is the determination of the time increment in subroutine compdt. The formulas for doing so are Formulas 176 and 178 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.