C  HEEF - One minute average sorted by local pitch angle
C  Input file = 'd04XXXX.min' ; Output file = 'd04XXXX.lis''

c     05/02/2001: Modified to run on non-SGI UNIX platform by J. F. Cooper

      PROGRAM READ_MIN

      REAL*4 EPHEMA(8),EPHEMB(8),FLUX(19,19),ENERGY(19)
      REAL*4 RS4REV,Y
      INTEGER*2 NOBF(19,19),IDUM
      INTEGER*4 IORBIT,IYEAR,IDAY,IBSEC,IESEC,IBR,IER
      CHARACTER*11 INFILE,OUTFILE
      DATA inFILE/'d04XXXX.min'/

      PRINT*,'ENTER ORBIT NUMBER:'
      READ(*,*)JORBIT
      PRINT*,'ENTER CHANNEL NUMBER:(1-19)'
      READ(*,*)ICH   

      WRITE(INFILE(4:7),'(I4.4)')JORBIT
      WRITE(6,*) INFILE
      OUTFILE=INFILE(1:8)//'lis'
      WRITE(6,*) OUTFILE

c     Use this OPEN statement to run code on non-SGI UNIX platform
      OPEN(1,FILE=INFILE,ACCESS='DIRECT',RECL=558,
     *   CONVERT='BIG_ENDIAN',ERR=100,STATUS='OLD')
c     Use this OPEN statement to run code on SGI UNIX platform
c      OPEN(1,FILE=INFILE,ACCESS='DIRECT',RECL=558,ERR=100,
c     *   STATUS='OLD')
c
      WRITE(6,*) ' OPENED INPUT FILE'
      OPEN(2,FILE=OUTFILE,STATUS='UNKNOWN')
      WRITE(6,*) ' OPENED INPUT AND OUTPUT FILES '
C HEADER FILE
      READ(1,REC=1)IORBIT,IYEAR,IDAY,IBSEC,IESEC,IBR,IER,ENERGY

      WRITE(2,21) IYEAR,IDAY,ENERGY(ICH)
 21   FORMAT(' YEAR=',I4,' DAY=',I4,' ENERGY=',F5.2,' MEV')

      WRITE(2,22)IORBIT,IBSEC,IESEC,IER
 22   FORMAT(' ORBIT=',I4,' BSEC=',I6,' ESEC=',I6,' NO. RECORDS=',
     +         I4)

      DO 10 NREC=ibr,ier
       WRITE(2,*)
       READ(1,REC=NREC,ERR=100)EPHEMB,EPHEMA,FLUX,NOBF,IDUM
       WRITE(2,23)EPHEMB
 23    FORMAT(' BEG EPHEM. ',4F9.1,4F8.3)
       write(2,24)EPHEMA
 24    FORMAT(' AVE EPHEM. ',4F9.1,4F8.3)

       WRITE(2,25)
 25    FORMAT(/,' Pitch Angle   Flux      Obs')

       DO 20 IPA=1,19
         WRITE(2,26)IPA*5-5,FLUX(ICH,IPA),NOBF(ICH,IPA)
 20    CONTINUE

 26    FORMAT(I8,4X,E10.4,i8)

 10   CONTINUE

 100  CONTINUE  

      STOP

      END
c