#include "gclgrid.h" dmatrix ustrans(GCLgrid& g, double lat, double lon);
Sometimes it is useful to work in a local, Cartesian reference frame with a standard convention with x1 positive east, x2 positive north, and x3 positive up. This routine computes a transformation matrix from the standard geographic Cartesian reference frame (origin at Earth's center, x1 pointed at intersection of the prime meridian and the equator, x2 pointed at 90 east and the equator, and x3 pointed at the north pole) to a local coordinate system at (lat,lon).
-lgclgrid
dmatrix(3)
Gary Pavlis Indiana University pavlis@indiana.edu