function dist,line,point
;**line=[m,b] for a [lon,lat] rectangular 2D coordinate system
;**point=[lat,lon]
x=(line[0]*point[0]+point[1]-line[0]*line[1])/(line[0]^2+1)
y=line[0]*x+line[1]
return,arclength(point,[y,x])
end