modify.f

Go to the documentation of this file.

      SUBROUTINE modify(KSTART)
*
*
*       Parameter modification at restart.
*       ----------------------------------
*
      include 'common6.h'
      EXTERNAL verify
*
*
*       Read first, second or both lines (KSTART = 3, 4, 5).
      IF (kstart.EQ.4) go to 10
*
*       Read new DTADJ, DELTAT, TADJ, TNEXT, TCRIT, QE & KZ(J) (if > 0).
      READ (5,*)  dta, dt, ta, tn, tc, qe1, j, k
*
*       Copy old parameters if corresponding input is zero.
      IF (dta.LE.0.0) THEN
          dta = dtadj
      END IF
*
      IF (dt.LE.0.0) THEN
          dt = deltat
      END IF
*
      IF (ta.LE.0.0) THEN
          tadj = max(tadj - dtadj + dta,time)
      ELSE
          tadj = max(ta-toff,time)
      END IF
*
      IF (tn.LE.0.0) THEN
          tnext = max(tnext - deltat + dt,time)
      ELSE
          tnext = max(tn-toff,time)
      END IF
*
      dtadj = dta
      deltat = dt
      IF (tc.GT.0.0) tcrit = tc
      IF (qe1.GT.0.0) qe = qe1
*
*       See whether any options should be changed.
      IF (j.GT.0) kz(j) = k
*
      WRITE (6,5)  dtadj, deltat, tcrit, qe, j, k
    5 FORMAT (///,7x,'RESTART PARAMETERS:   DTADJ =',f7.3,'  DELTAT =',
     &                            f7.3,'  TCRIT =',f7.1,'  QE =',1pe9.1,
     &                                            '  KZ(',i2,') =',i2,/)
*
*       Read new ETAI, ETAR, ETAU, DTMIN, RMIN, NCRIT (if > 0 & KSTART >= 4).
   10 IF (kstart.GE.4) THEN
          READ (5,*)  eta1, eta2, eta3, dtm, rm, newcr
*
*       Check modification of integration parameters.
          IF (eta1.GT.0.0) etai = eta1
          IF (eta2.GT.0.0) etar = eta2
          IF (eta3.GT.0.0) etau = eta3
          IF (dtm.GT.0.0) THEN
              dtmin = dtm
              smin = 2.0*dtm
          END IF
          IF (rm.GT.0.0) THEN
              rmin = rm
              rmin2 = rm**2
              rmin22 = 4.0*rmin2
          END IF
          IF (newcr.GT.0) ncrit = newcr
*
          WRITE (6,15)  etai, etar, etau, dtmin, rmin, ncrit
   15     FORMAT (/,7x,'RESTART PARAMETERS:   ETAI =',f6.3,'  ETAR =',
     &                        f6.3,'  ETAU =',f6.3,'  DTMIN =',1p,e8.1,
     &                             '  RMIN =',e8.1,'  NCRIT =',0p,i5,/)
      END IF
*
*       Perform a simple validation check on main input parameters.
      CALL verify
*
*       Save the new parameters on tape/disc unit #1 just in case.
      CALL mydump(1,1)
*
      RETURN
*
      END