stepk.f

Go to the documentation of this file.

      SUBROUTINE stepk(DT,DTN)
*
*
*       Block time-steps.
*       -----------------
*
      include 'common6.h'
      DATA  one32 /0.03125/
*     DATA  ONE16,ONE32 /0.0625D0,0.03125/
*
*
*       Truncate any large value to first block-step entry.
      k = 1
      IF (dt.GE.smax) go to 10
*
*       Compare predicted step with discrete values decreasing by 1/32.
      dt1 = dtk(1)
    1 IF (dt1.GT.dt) THEN
*         DT1 = DT1*ONE16
*         K = K + 4
          dt1 = dt1*one32
          k = k + 5
          go to 1
      END IF
*
*       Increase by 2 until predicted step is exceeded (then correct level).
    4 IF (dt1.LT.dt) THEN
          dt1 = 2.0d0*dt1
          k = k - 1
          go to 4
      END IF
      k = k + 1
*
   10 dtn = dtk(k)
*
      RETURN
*
      END