gsi-wgs-2.f:after


c     for gsi-wgs program
      character aht
      open(10,file='d:\Myhome\yoreyore\akkeshi.gsi.out',status='old')
      open(15,file='d:\Myhome\yoreyore\akkeshi.txt',status='old')
      open(20,file='d:\Myhome\yoreyore\akkeshi.wgs.data')
    1 continue
      read(10,*,end=99) wlat,wlon,aht
      read(15,*,end=99) zlon,zlat,bht
      lln=int(wlon/10000)
      blon=wlon-lln*10000
      mln=int(blon/100)
      slon=blon-mln*100
      llt=int(wlat/10000)
      blat=wlat-llt*10000
      mlt=int(blat/100)
      slat=blat-mlt*100
c
      aalt=llt+float(mlt)/60+slat/3600
      aaln=lln+float(mln)/60+slon/3600
      write(6,200) aaln,aalt,bht
      write(20,200) aaln,aalt,bht
      go to 1
   99 continue
      close(10)
      close(20)
  200   format(f12.7,f15.7,f10.2)
      stop
      end