input.f

Go to the documentation of this file.

      SUBROUTINE input
*
*
*       Parameter input.
*       ----------------
*
      include 'common6.h'
      EXTERNAL verify
*
*
*       Make a formal call to define input parameters & counters.
      CALL define
*
*       Read & print the main input parameters.
      READ (5,*)  n, nfix, ncrit, nrand, nnbmax, nrun
      READ (5,*)  etai, etar, rs0, dtadj, deltat, tcrit, qe, rbar, zmbar
      READ (5,*)  (kz(j),j=1,50)
      READ (5,*)  dtmin, rmin, etau, eclose, gmin, gmax
*
      WRITE (6,10)
   10 FORMAT (/////,15x,'N  NFIX  NCRIT  NRAND  NNBMAX  NRUN')
      WRITE (6,12)  n, nfix, ncrit, nrand, nnbmax, nrun
   12 FORMAT (/,i16,i6,2i7,i8,i6)
      WRITE (6,15)
   15 FORMAT (//,12x,'ETAI      ETAR      RS0       DTADJ     DELTAT',
     &                        '    TCRIT     QE        RBAR      ZMBAR')
      WRITE (6,20)  etai, etar, rs0, dtadj, deltat, tcrit, qe, rbar,
     &                                                             zmbar
   20 FORMAT (/,9x,1p,10e10.1)
      WRITE (6,22)
   22 FORMAT (//,12x,'OPTIONS')
      WRITE (6,24)  (j,j=1,50)
   24 FORMAT (/,9x,50i3)
      WRITE (6,26)  (kz(j),j=1,50)
   26 FORMAT (/,9x,50i3)
      WRITE (6,28)
   28 FORMAT (//,12x,'DTMIN     RMIN      ETAU      ECLOSE    GMIN',
     &                                                     '      GMAX')
      WRITE (6,30)  dtmin, rmin, etau, eclose, gmin, gmax
   30 FORMAT (/,9x,1p,6e10.1)
*
*       Perform a simple validation check on main input parameters.
      CALL verify
*
      gprint(1) = 0.0
      deltas = 0.0
      IF (kz(4).GT.0) THEN
*       Read parameters for binary evolution analysis.
          k = kz(4)
          READ (5,*)  deltas, orbits(1), (gprint(j),j=1,k)
          WRITE (6,40)  deltas, orbits(1), (gprint(j),j=1,k)
   40     FORMAT (//,12x,'DELTAS =',f6.2,'  ORBITS(1) =',f6.2,
     &                                            '  GPRINT(J) =',9f7.3)
*       Modify binary output factor by perturbation at different levels.
          DO 50 l = 2,k
              orbits(l) = orbits(1)*(gprint(1)/gprint(l))**0.3333
   50     CONTINUE
      END IF
*
*       Set random number skip for routine DATA.
      idum1 = nrand
*
*       Save square of c.m. approximation parameter (coupled to GMIN).
      cmsep2 = gmin**(-0.666666667)
*
*       Define total particle number & neighbour membership parameters.
      ntot = n
      nzero = n
      nbzero = nnbmax
      znbmin = 0.2*float(nnbmax)
      znbmax = 0.9*float(nnbmax)
*       Save initial ETAI and RMIN**2 (for PHICOR).
      eta0 = etai
      rmin2 = rmin**2
      rsmin = rs0
      rc = rs0
*       Temporary save of initial neighbour sphere for OUTPUT & NBLIST.
*
      RETURN
*
      END