Mercator Sailing: Course and Distance

08252018, 06:08 AM
(This post was last modified: 08252018 09:47 PM by Dieter.)
Post: #5




RE: Mercator Sailing: Course and Distance
(08242018 07:21 PM)Gene222 Wrote: I modified Eddie's program where South latitudes and West longitudes must be entered as negative numbers in degrees, but I ran into a problem when the two points have the same latitude, such as I wrote a Mercator Saling program for the 35s and encountered the same problem. However, a bit of calculus yields the following result: Code: 1  e²·sin²(lat) The course c then is 0° (exactly East or West). Or 90° / 270° true course. Here "e" is the eccentricity of the assumed ellipsoid (around 0,082). This is not explicitely considered in the approximate formula given in the book (as used in Namir's program). But we can do better: calculate the meridional parts (M) as follows using hyperbolic functions: Code: M = 60·180/pi · [artanh(sin(lat)) – e·artanh(e·sin(lat))] You can also do it without hyperbolics: Code: M = 60·180/pi · ln[tan(45° + lat/2) · ((1–e·sin(lat))/(1+e·sin(lat))^(e/2)] BTW, the constant 7915,7045 in the approximate formula is 60·180/pi · ln10 due to the use of the base 10 logarithm instead of the natural log. Both formulas are more accurate and allow implementing different ellipsoids. I prefer the one with the two artanhs since it is more compact and ...elegant. ;) For "e" you may use these values: Clarke 1880 ellipsoid: e = 0,082483399 WGS84 / GRS80 ellipsoid: e = 0,081819191 For your example the distance is 230,7 nmi (for both ellipsoids the results differ by less than 0,02 nmi). Check: use your program "as is" and make the latitudes slightly different, e.g. 40° and 40,0001°. This should give (nearly) the same result. Remark: instead of calculating the arctan you can determine the angle with a R→P conversion (automatically selects the correct quadrant) followed by mod 360° to get a positive angle. This way the four IF cases are not required (BTW, are you sure the last two are both correct?). Dieter 

« Next Oldest  Next Newest »

User(s) browsing this thread: 1 Guest(s)