C PROGRAM TO CALCULATE NMR FREQUENCY VS. PARTICLE ENERGY TABLES
C
C THE ORIGINS OF THIS CODE ARE UNKNOWN; IT WAS FOUND AT THE BACK
C OF THE PREVIOUSLY USED TABLES.  THE PREVIOUSLY USED TABLES WERE
C WORN OUT, AND HAD TO BE REPLACED.  THE NUMBERS ARE IDENTICAL.
C
C THE OUTPUT FORMAT WAS ALTERED SLIGHTLY TO ALLOW FOR USE ON THE
C 6-LPI PRINTER; THE ORIGINAL VERSION WAS DESIGNED FOR AN 8-LPI
C OUTPUT.  THE HEADER AND TRAILER LINES WERE THE ONLY THINGS CHANGED.
C
C NOVEMBER 1987; D.E. CARTER AND J.E. O'DONNELL
C
      REAL*8 M,MZERO,K,NU(10),E,EC,Z
      INTEGER ZIN
      DIMENSION NAME(20)
      DATA MZERO/1.8765592D03/
C     Header
   1  FORMAT(1x,'E(MEV) VS. NMR (MHZ) FOR',1X,
     >20A1,2X,'MASS = ',F9.6,2X,'CHG STATE = +',I1,2X,'K = ',1PE10.3
     >/9x,'E',10X,'.000',9x,'.001',9x,'.002',9x,'.003',9x,
     > '.004',9x,'.005',9x,'.006',9x,'.007',9x,'.008',9x,'.009')

   2  FORMAT(5X,F6.2,10(3X,F8.4))
   3  FORMAT(1X)
   4  FORMAT(20A1,F10.4,4X,I1,5X,F10.4,5X,F5.0,5X,F5.0)
C     Footer
   5  FORMAT(1H ,48X,20A1,2X,'CHARGE STATE = +',I1,5X,
     >  'E = ',F7.3,' TO ',F7.3,' MEV ')
   6  FORMAT(/' ')
      OPEN (FILE='NMRTABLE.INP',UNIT=1)
      OPEN (FILE='NMRTABLE.OUT',UNIT=3,status="UNKNOWN")
1001  READ (1,4) NAME,M,ZIN,K,EIN,EOUT
      IF (M.EQ.0.00) STOP
      IF (ZIN.EQ.0) ZIN = 1
      Z = ZIN
      E = EIN
C     L = 2.0D0*(EOUT-EIN)
      L = 3.0D0*(EOUT-EIN)
      DO 1005 IIII=1,L
c         print *,iiii
C        WRITE (3,6)
         WRITE (3,1) NAME,M,ZIN,K
         DO 1004 III=1,3
            DO 1003 II=1,10
               DO 1002 I=1,10

C    EC = particle energy step by 0.001 MeV

               EC = E+0.001D0*DFLOAT(I-1)
1002           NU(I) = DSQRT(EC*M/K*(1.0D0+EC/(M*MZERO)))/Z

C    Write out the frequency

               WRITE (3,2) E,NU

C   step up by 0.010 MeV

1003        E = E + 0.01D0

1004     WRITE (3,3)
         ET = E - 0.5D0

1005  WRITE (3,5) NAME,ZIN,ET,EC
      GOTO 1001
      END