Assuming calcdist is required, the following result was found.
+ (cos(v_RadLat1) * cos(v_RadLat2) * cos(v_RadTheta)); if (v_Dist > 1) { v_Dist = 1; } v_AcosDist = acos(v_Dist); v_CalcDist = v_AcosDist * 180 / v_MathPi; v_Miles = v_CalcDist * 60 * 1.1515; } return round(v_Miles,2) + " miles"; } For kilometres,...