program READdata

c     pass: orbit number, date:yyyymmdd
c     ut:hhmmss, lt:+-hhmmss (-: LT belongs previous day relative to UT)
c     lat, lon: geodetic latitude, longitude (deg)
c     mlat, mlon: geomagnetic latitude, longitude	(deg)
c     alt: alitude(km), dipag: dip-angle (deg)
c     invlat:invariant latitude (deg), L: L-value
c     sublat,sublon: geodetic latitude,longitude of sub-solar point
c     sza: solar zenith angle (deg)
c     ne:electron density at satellite altitude  (/cm**3)
c     nmf2: maximum density of F2-layer (/cm**3)
c	proton: [H+], helium: [He+], oxygen: [O+](/cm**3)
c     The - sign for ne, nmf2, proton, helium, or oxygen means that the quality
c       of the data is low.  Take absolute value to use that value.
c     The value of 0.0 for the density data measns that it is not available.     
c     maximum record number: 152840

      character fname*30
      integer*4 idata(4)
	integer*4 pass,date,ut,lt
      real*4 rdata(16)
	real*4 proton, helium, oxygen, nmf2, ne
	real*4 lat,lon,alt,mlat,mlon,dipang,invlat,L
	real*4 sublat,sublon,sza

	equivalence (pass,idata(1)),  (date,idata(2))
	equivalence (ut,idata(3)),    (lt,idata(4))
	equivalence (lat,rdata(1)),   (lon,rdata(2)),(alt,rdata(3))
	equivalence (mlat,rdata(4)),  (mlon,rdata(5)),(dipang,rdata(6))
	equivalence (invlat,rdata(7)),(L,rdata(8))
	equivalence (sublat,rdata(9)),(sublon,rdata(10))
      equivalence (sza,rdata(11))
      equivalence (nmf2,rdata(12)), (ne,rdata(13))
	equivalence (proton,rdata(14)),(helium,rdata(15))
	equivalence (oxygen,rdata(16))


      write(*,*) 'Input data-file pass and name, eg., d:\data\ISS-b.dat'
      read(*,10) fname
	open(1,file=fname,status='old')
	

    1	continue
	write(*,*) 'Input Record Number < 152841 ,or O = end '
	read(*,*) nrec
	if(nrec.le.0.or.nrec.gt.152840) go to 99

	do 100 irec=1,nrec
	read(1,40) idata,rdata
	
  100	continue

	write(*,50) nrec,idata,rdata


   10 format(A30)
   40	format(4I10,16F10.1)
   50	format(I10/,4I10/,4(4F10.1/))

   	rewind(1)
	go to 1

   99	continue
	

	close(1)

	pause
      stop
      end