C     PROGRAM READHEL
C     Version 3/19/93:
C     3/19/93: corrected proton flow angle indices from 23,24 to 26,27.
      REAL*4 PBLOCK(3000)
      CHARACTER*70 INFILE,OUTFIL
C
      WRITE(6,*) ' Enter Helios input file name:'
      READ(5,10) INFILE
10    FORMAT(A70)
      WRITE(6,*) ' Enter Helios output file name:'
      READ(5,10) OUTFIL
      OPEN(UNIT=10,FILE=INFILE,BLOCKSIZE=3000,FORM='UNFORMATTED',
     *   RECORDTYPE='VARIABLE',STATUS='OLD')
      OPEN(UNIT=20,FILE=OUTFIL,STATUS='NEW')
      WRITE(6,*) ' Enter year of output file (YY): '
      READ(5,*) MY
100   READ(10) PBLOCK
      DO I=1,100
         J0= (I-1)*30
         ICODE= PBLOCK(J0+2) ! 1=PLASMA 2=MAG 3=BOTH
         NY= PBLOCK(J0+3)
         ND= PBLOCK(J0+4)
         NH= PBLOCK(J0+5)
         FD1972= PBLOCK(J0+6)
         SCRAU= PBLOCK(J0+7)   ! RADIAL DISTANCE OF S/C FROM SUN IN AU
         NCARROT= PBLOCK(J0+8) ! CARRINGTON ROTATION NUMBER
         CARLON= PBLOCK(J0+9)  ! CARRINGTON LONGITUDE OF S/C
         HELLAT= PBLOCK(J0+10)  ! HELIOGRAPHIC LATITUDE OF S/C
         ESSCSANG= PBLOCK(J0+11) ! EARTH-SUN-S/C SEPARATION ANGLE
         BXSE= PBLOCK(J0+12) ! IMF COMPONENT ALONG X-SOLAR ECLIPTIC AXIS -> SUN
         BYSE= PBLOCK(J0+13) ! IMF COMPONENT ALONG Y-SOLAR ECLIPTIC AXIS (NT)
         BZSE= PBLOCK(J0+14) ! IMF COMPONENT ALONG Z-SOLAR ECLIPTIC AXIS (NT)
         BMAG= PBLOCK(J0+15) ! IMF VECTOR MAGNITUDE (NT)
         SBXSE= PBLOCK(J0+16) ! STAND. DEV. OF BXSE (NT)
         SBYSE= PBLOCK(J0+17) !  STAND. DEV. OF BYSE (NT)
         SBZSE= PBLOCK(J0+18) ! STAND. DEV. OF BZSE (NT)
         SBMAG= PBLOCK(J0+19) ! STAND. DEV. OF BMAG (NT)
         BTHETA= PBLOCK(J0+28) ! ELEVATION ANGLE THETA OF VECTOR B IN SE SYSTEM
         BPHI= PBLOCK(J0+29)  ! AZIMUTH ANGLE PHI OF VECTOR B IN SE SYSTEM
         SWSPD= PBLOCK(J0+20) ! SOLAR WIND BULK FLOW SPEED (KM/SEC)
         TPRADKK= PBLOCK(J0+21) ! RADIAL COMPONENT OF PROTON TEMPERATURE 
C                                 (10*3 DEGREES K.)
         TPRADK= TPRADKK*1000.0 ! RADIAL P TEMPERATURE (DEG. K)
         SWPDEN= PBLOCK(J0+22) ! SOLAR WIND NUMBER DENSITY (PROTONS/CC)
C        The following two lines were corrected on 3/19/93. The item fields
C         were corrected to 26 and 27 from 23 and 24, respect. JFC
         PFELEV= PBLOCK(J0+26) ! PROTON FLOW ELEVATION ANGLE (DEG, + = NORTH)
         PFAZIM= PBLOCK(J0+27) ! PROTON FLOW AZIMUTH ANGLE (DEG, + = EASTWARD
C         WRITE(6,9876) (PBLOCK(J0+N),N=26,27)
C9876     FORMAT(1X,2E12.4)
C                                OR CCW FROM SUNWARD DIRECTION)
         IF(MY.EQ.NY) THEN
C             WRITE(20,140) NY,ND,NH,ICODE,SCRAU,BXSE,SWSPD
C140          FORMAT(1X,I2,1X,I3,1X,I2,1X,I1,1X,3E12.4)
C            WRITE(20,110) NY,ND,NH,NCARROT,SCRAU,CARLON,HELLAT,
C    *         ESSCSANG,ICODE ! WRITE STATEMENT FOR TRAJECTORY DATA
C             IF(ICODE.NE.1) THEN   ! WRITE MAGNETOMETER DATA
C                WRITE(20,120) NY,ND,NH,BXSE,BYSE,BZSE,BMAG,BTHETA,BPHI
C             ENDIF
              IF(ICODE.NE.2) THEN   ! WRITE PLASMA DATA
               WRITE(20,130) NY,ND,NH,SWSPD,TPRADK,SWPDEN,PFELEV,PFAZIM
              ENDIF
         ELSE IF(NY.GT.MY) THEN
            GOTO 999
         ENDIF
110      FORMAT(1X,I2,1X,I3,1X,I2,1X,I4,1X,F6.4,1X,F7.3,1X,F6.3,1X,
     *      F7.3,1X,I1)
120      FORMAT(1X,I2,1X,I3,1X,I2,3F8.2,F6.1,2F7.2)
130      FORMAT(1X,I2,1X,I3,1X,I2,F8.2,2E10.3,2F7.2)
      ENDDO
      GOTO 100
999   STOP
      END