PRO find_eur,plat,plong,glat,glong,eplat,eplong,ec,er eglat=80.0533 eglong=-86.4161 pd=ACOS(SIN(eplat*!DTOR)*SIN(plat*!DTOR) $ +COS((plong-eplong)*!DTOR)*COS(eplat*!DTOR)*COS(plat*!DTOR)) gd=ACOS(SIN(eglat*!DTOR)*SIN(glat*!DTOR) $ +COS((glong-eglong)*!DTOR)*COS(eglat*!DTOR)*COS(glat*!DTOR)) omin=MIN(SQRT(0.5*(pd^2+gd^2)),oeur) one2two,oeur,plat,ec,er RETURN END