PRO mdmsp,file

eurlat=80.0533
eurlon=-86.4161
inf=qgetfile(file)
nl=N_ELEMENTS(inf)
outf=STRARR(nl)
j=-1
FOR i=0,nl-1 DO BEGIN
    CASE 1 OF
(isnumber(STRMID(inf(i),0,2)) NE 0): $
	BEGIN				; process data
	    j=j+1
	    s=inf(i)
	    dlat=FLOAT(getwrd(s,3))
	    dlon=FLOAT(getwrd(s,4))
	    ll2rb,eurlon,eurlat,dlon,dlat,dist,azi
	    outf(j)=STRING(s,dist*!RADEG,azi,FORMAT='(A,F7.2,F7.1)')
	    END
(STRMID(inf(i),0,2) EQ 'yy'): $
	BEGIN				; augment title
	    j=j+1
	    outf(j)=inf(i)+'   Dist   Azim'
	    END
(STRMID(inf(i),46,5) EQ 'page:'):	; skip page #
ELSE:   BEGIN
	    j=j+1
	    outf(j)=inf(i)
	    END
	ENDCASE
    ENDFOR
putfile,file+'+',outf(0:j)
RETURN
END