C Program XFAN_LIST, for use with 77-102A-01N,
C                                 78-079A-01O (Solar Wind electron moments).
C                                 and probably other related data sets.

C This program was written in April 1994 by H. K. Hills, using format given by
C the data set documentation. The program reads a logical data record (composed
C of two physical records) and writes the data to a (screen) file using two
C lines, to avoid displaying lines requiring more than 80 characters. 
C Modified October 1994, H. K. Hills (ask input filename, recognize EOF).
C Modified January 26, 1995, H. K. Hills (open as READONLY).
C Modified May 25, 1995, H. K. Hills, to include additional data sets
C that this will handle, and to give user more information at program start.

      Program XFAN_LIST
      real*4 plasma(10)
      character*60 infilename
      write(6,600)
  600 format(/,' For use with the ISEE 1 Solar Wind Ion Moments data',
     1 ' set (77-102A-01N)',/,' from the Cross-Fan (X-Fan) Analyzer,',
     1 ' and with the ISEE 3 Solar Wind Electron',/,
     1 ' Moments data set (78-079A-01O).'/,
     2 ' This program is entered as B44417 in NSSDC''s Technical ',
     3 'Reference File (TRF).',/,
     4 ' It asks user for input filename (uses unit 10), and reads',
     4 ' the pair of long',/,' (85-character) and short records,',
     4 ' then copies them to unit 6, in records short',/,
     4 ' enough (max 73 bytes) for most monitor screens.',/,
     5 ' User input from unit 5; output to unit 6.',//)           

      write(6,*) 'Enter input file name.'
      read(5,401) infilename
  401 format(A60)
      open(unit=10,file=infilename,READONLY, access='sequential',
     1 status='old')

      write(6,*) 'Enter number of records to read and list.'
      read(5,*)irecno
      icnt=0
    1 if(icnt.eq.irecno) then
         write(6,*) irecno,' records read and printed.'
         stop
      endif
      read(10,1000,end=30) iyymmdd,sec,utime,(plasma(i),i=4,10)
 1000 format(1x,I6,1x,2(F8.1,1x),5(1PE11.4,1x),15(/,6(1PE11.4,1x)))

C This format from Gosling seems to be both an input format and an output
C format, and to also be valid for a longer record.
C HKH changed this to a slightly different output format to write on screens
C that don't allow 85 characters on a line. The new format is statement number
C 2000. 

      write(6,2000) iyymmdd,sec,utime,(plasma(i),i=4,10)
 2000 format(1x,I6,1x,2(F8.1,1x),4(1PE11.4,1x),/,26x,3(1PE11.4,1x))
      icnt=icnt+1
      go to 1
   30 continue
      write(6,*) 'EOF encountered.'
      stop
      end