c
      subroutine rddat
c
c Read satellite data from file satfile, defining /dat/:
c
      include "tgcmsat.h"
      character*80 chrec
c
      lenflnm = lenstr(satfile)
      open(lusat,file=satfile(1:lenflnm),status='OLD',err=900,
     +  iostat=ios)
      write(6,"('Opened file ',a,' (lusat=',i2,')')")
     +  satfile(1:lenflnm),lusat
      npts = 0
 100  continue
      read(lusat,"(a)",end=101) chrec
      if (chrec(1:1).eq.';') then
c       write(6,"(a)") chrec(1:lenstr(chrec))
        goto 100
      endif
      read(chrec,"(i2,i4,f7.2,2f9.2,f7.2)")
     +  ihr,min,sec,glat,glon,slt
      npts = npts+1
      satut(npts) = float(ihr)+float(min)/60.+sec/3600.
      satlat(npts) = glat
      satlon(npts) = glon
      if (satlon(npts).ge.180.) satlon(npts) = satlon(npts) - 360.
      satslt(npts) = slt
c     write(6,"('npts=',i4,' satut=',f6.3,' satlat,lon=',f7.2,
c    +  f8.2,' satslt=',f7.2)") npts,satut(npts),satlat(npts),
c    +  satlon(npts),satslt(npts)
      goto 100
 101  write(6,"('EOF encountered on file ',a,' at npts=',i4)")
     +  satfile(1:lenflnm),npts
      return
 900  continue
      write(6,"('>>> rddat: error opening file ',a,' ios=',i5)")
     +  satfile(1:lenflnm),ios
      stop 'satfile'
      return
      end