verify.f

Go to the documentation of this file.

      SUBROUTINE verify
*
*
*       Input validation.
*       -----------------
*
      include 'common6.h'
*
*
*       Check for unreasonable input parameters (initial & restart).
      IF (n.GE.nmax - 2.OR.nnbmax.GT.lmax - 3) THEN
          WRITE (6,10)  n, nnbmax
   10     FORMAT (/,5x,'FATAL ERROR!   BAD INPUT   N =',i6,
     &                                                  '  NNBMAX =',i4)
          stop
      END IF
*
      IF (etai.GT.0.08.OR.etar.GT.0.16) THEN
          WRITE (6,20)  etai, etar
   20     FORMAT (/,5x,'FATAL ERROR!   BAD INPUT   ETAI =',f6.2,
     &                                                  '  ETAR =',f6.2)
          stop
      END IF
*
      IF (etau.GT.0.5.OR.gmin.GT.0.0001.OR.gmax.GT.0.10) THEN
          WRITE (6,30)  etau, gmin, gmax
   30     FORMAT (/,5x,'FATAL ERROR!   BAD INPUT   ETAU =',f6.2,
     &                                 '  GMIN =',f11.7,'  GMAX =',f7.3)
          stop
      END IF
*
*       Also check for zero or negative values.
      IF (n.LE.0.OR.nnbmax.LE.0.OR.etai.LE.0.0.OR.etar.LE.0.0) THEN
          WRITE (6,40)  n, nnbmax, etai, etar
   40     FORMAT (/,5x,'FATAL ERROR!   BAD INPUT   N =',i5,
     &                  '  NNBMAX =',i4,'  ETAI =',f6.2,'  ETAR =',f6.2)
          stop
      END IF
*
      IF (etau.LE.0.0.OR.gmin.LE.0.0.OR.gmax.LE.0.0) THEN
          WRITE (6,30)  etau, gmin, gmax
          stop
      END IF
*
      IF (dtadj.LE.0.0.OR.deltat.LE.0.0.OR.qe.LE.0.0) THEN
          WRITE (6,50)  dtadj, deltat, qe
   50     FORMAT (/,5x,'FATAL ERROR!   BAD INPUT   DTADJ =',f6.2,
     &                                '  DELTAT =',f6.2,'  QE =',1pe9.1)
          stop
      END IF
*
      RETURN
*
      END