program convrt
cc
c This program converts the World2 data base files from
c   the raw ASCII form to direct access binary and it also
c   builds the corresponding index files.
c  Center for Seismic Studies - Louis Huszar - Mar 84
cc
	character*1 lad,lod
c	open(3,file='cil.dat',form='formatted',
c	open(3,file='bdy.dat',form='formatted',
c	open(3,file='riv.dat',form='formatted',
	open(3,file='pby.dat',form='formatted',
&	status='old',access='sequential',iostat=istat)
	if(istat.ne.0) stop 1
c	open(4,file='cil.unf',form='unformatted',
c	open(4,file='bdy.unf',form='unformatted',
c	open(4,file='riv.unf',form='unformatted',
	open(4,file='pby.unf',form='unformatted',
&	status='new',access='direct',recl=8)
	inxct=0
	istrt=1
	irec=0
3	read(3,end=20,fmt='(i7,i2,i6)') lsn,irnk,nop
c read in line segment number, rank, no of points
c	write(6,'(1x,i7,i3,i6)') lsn,irnk,nop
	do 1 i=1,nop
	read(3,'(3i2,a1,i3,2i2,a1,i5)')ladg,lam,las,lad,
&	lodg,lom,los,lod,ipn
	if(i.ne.ipn) stop 2
	alt=float(ladg)+float(lam)/60.+float(las)/3600.
	if(lad.eq.'S') alt=-alt
	aln=float(lodg)+float(lom)/60.+float(los)/3600.
	if(lod.eq.'W') aln=-aln
c	if(i.eq.1.or.i.eq.nop)
c&	write(6,'(i6,2f10.4,1x,3i2,a1,i3,2i2,a1)') ipn,
c&	alt,aln,ladg,lam,las,lad,lodg,lom,los,lod
	irec=irec+1
	write(4,rec=irec) alt,aln
	if(i.gt.1) go to 2
	altmn=alt
	altmx=alt
	alnmn=aln
	alnmx=aln
	go to 1
2	if(alt.lt.altmn) altmn=alt
	if(alt.gt.altmx) altmx=alt
	if(aln.lt.alnmn) alnmn=aln
	if(aln.gt.alnmx) alnmx=aln
1	continue
	inxct=inxct+1
	write(6,'(1x,i6,i8,i3,i8,i6,4f10.5)') inxct,lsn,irnk,
&	istrt,nop,altmn,alnmn,altmx,alnmx
	istrt=istrt+nop
	go to 3
20	stop
	end