Program: MKMIDPOINT
mkmidpoint........compute middle of great circle arc between two points
call: [midlon,midlat]=mkmidpoint(lon1,lat1,lon2,lat2);
lon1: (number) [deg]
longitude of 1st point
lat1: (number) [deg]
latitude of 1st point
lon2: (number) [deg]
longitude of 2nd point
lat2: (number) [deg]
latitude of 2nd point
result: midlon: (number) [deg]
longitude of the point in the middle
midlat: (number) [deg]
latitude of the point in the middle
The mid-point computation is needed for the refinement of triangular
grids, where the edges of triangles are divided in the middle to form the
edges of new, smaller triangles.
The computation is a quite simple vector addition procedure.
Martin Knapmeyer, 15.09.2004
Read M-File Source Code
eof.