PROGRAM READIT
         CHARACTER CJUNK
         PARAMETER (PI=3.141592654)
         real ZO,ZMAX,DZ,angle
         real year, month,day,hour,lat,long,height
         real lat0,long0,mlat
         lat0=78.5
         long0=-291.
         open(10,file="msis2.dat")
         read(10,*) cjunk,year,cjunk,month,cjunk,day,cjunk
         read(10,*) cjunk,lat,cjunk,long,cjunk,height
         read(10,*)cjunk,cjunk,cjunk,Z0,cjunk,ZMAX,cjunk,dz
         do 10 i=1,9
  10        read(10,*)cjunk
         read(10,*) xx
         close(10)
         lat=torad(lat)
         lat0=torad(lat0)
         long=torad(long)
         long0=torad(long0)
        mlat=asin(sin(lat)*sin(lat0)+cos(lat)*cos(lat0)*cos(long-long0))
         angle=todeg(atan(2.*tan(mlat)))
         write(6,*)year,month,day
         write(6,*) lat,long,height
         write(6,*) z0,zmax,dz
         write(6,*)angle,xx
         end

         real function torad(x)
           real x
           PARAMETER (PI=3.141592654)
           torad=x*pi/180.
           return
           end

         real function todeg(x)
           real x
           PARAMETER (PI=3.141592654)
           todeg=x*180./pi
           return
           end