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