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.