.sigma - cos.u1 * cos.sigma * cos.alpha1
lat2 <- atan2(sin.u1 * cos.sigma + cos.u1 * sin.sigma * cos.alpha1,
(1 - ellips["f"]) * sqrt(sin.alpha ^ 2 + tmp ^ 2))
lambda <- atan2(sin.sigma * sin.alpha1, cos.u1 * cos.sigma - sin.u1 *
sin.sigma * cos.alpha1)
cap.C <- ellips["f
ry to roll my own. Thanks for any pointers.
Cheers,
eric
--
Eric Archer, Ph.D.
NOAA-SWFSC
8604 La Jolla Shores Dr.
La Jolla, CA 92037
858-546-7121,7003(FAX)
[EMAIL PROTECTED]
http://swfsc.noaa.gov/prd-etp.aspx
"Innocence about Science is the worst crime today."
- Sir Cha