pp=ddread('~steele/PACE/pace.pole',rows=INDGEN(25))
alt=pp(0,*)
glat=pp(1,*)
glong=(pp(2,*)+360) MOD 360

FOR i=0,N_ELEMENTS(alt)-1 DO BEGIN
    h=alt(i)
    glt=glat(i)
    gln=glong(i)
    mconvert,1995,1,1,0,0,0,glt,gln,h,1,mlat,mlong,mlt
    del=90.-mlat
    REPEAT BEGIN
        rb2ll,gln,glt,del*!DTOR,FINDGEN(8)*45,dln,dlt
	dmlat=FLTARR(8)
	dmlong=dmlat
        FOR j=0,7 DO BEGIN
	    mconvert,1995,1,1,0,0,0,dlt(j),dln(j),h,1,xx,yy,zz
	    dmlat(j)=xx
	    dmlong(j)=yy
	    ENDFOR
	PRINT,dlt,dln,dmlat,dmlong
        mxmlat=MAX(dmlat,jx)
	glt=dlt(jx)
	gln=dln(jx)
	del=90.-mxmlat
	PRINT,jx,glt,gln,del $
	     ,FORMAT='(i1,": (",F7.4,",",F8.4,") del = ",F6.4," deg")'
	ENDREP UNTIL (del LT 0.002)
    glat(i)=glt
    glong(i)=gln
    PRINT,alt(i),pp(1,i),pp(2,i),glt,gln $
	 ,FORMAT='(F4.0,": (",F7.4,",",F8.4,") -> (",F7.4,",",F8.4,")")'
    ENDFOR

OPENW,lun,'~/pace/pace.pole',/GET_LUN
FOR i=0,N_ELEMENTS(alt)-1 DO PRINTF,lun,alt(i),glat(i),glong(i)
FREE_LUN,lun
END