diff options
Diffstat (limited to 'math/slalib/doc/slalib.hlp.sav')
| -rw-r--r-- | math/slalib/doc/slalib.hlp.sav | 591 |
1 files changed, 591 insertions, 0 deletions
diff --git a/math/slalib/doc/slalib.hlp.sav b/math/slalib/doc/slalib.hlp.sav new file mode 100644 index 00000000..8da465cf --- /dev/null +++ b/math/slalib/doc/slalib.hlp.sav @@ -0,0 +1,591 @@ +.help slalib Nov95 "Immatch Package" +.ih +NAME +slalib -- Starlink library of positional astronomy routines + +.ih +DESCRIPTION +SLALIB is a library of Fortran 77 routines intended to make accurate +and reliable positional-astronomy applications easier to write. Most +SLALIB library routines are concerned with astronomical position and time, +but a number have wider trigonometrical, numerical or general applications. +SLALIB contains routines covering the following topics: 1) string +decoding and sexagesimal conversions, 2) angles, vectors and rotation +matrices, 3) calendars and timescales, 4) precession and nutation, 5) +proper motion, 6) FK4/5 and elliptic aberration, 7) geocentric coordinates, +8) apparent and observed place, 9) azimuth and elevation, 10) refraction +and air mass, 11) ecliptic, galactic, and supergalactic coordinates, +12) ephemerides, 13) astrometry, and 14) numerical methods. + +The labels and calling sequences of the SLALIB are listed below grouped +by function. To get more detailed help on any individual routine type +the help command followed by the label, e.g the command "help flotin" +will give detailed help on the subroutine slrfli. + +.ih +STRING DECODING +.nf + intin -- call subroutine slinti (string, nstrt, ireslt, jflag) + Convert free-format string into integer + +flotin -- call subroutine slrfli (string, nstrt, reslt, jflag) +dfltin -- call subroutine sldfli (string, nstrt, dreslt, jflag) + Convert free-format string into floating point number + + afin -- call subroutine slafin (string, iptr, a, j) + dafin -- call subroutine sldafn (string, iptr, a, j) + Convert free-format string from deg, armin, arcsec to radians +.fi +.ih +SEXAGESIMAL CONVERSIONS +.nf + ctf2d -- call subroutine slctfd (ihour, imin, sec, days, j) + dtf2d -- call subroutine sldtfd (ihour, imin, sec, days, j) + Hours, minutes, seconds to days + + cd2tf -- call subroutine slcdtf (ndp, days, sign, ihmsf) + dd2tf -- call subroutine slddtf (ndp, days, sign, ihmsf) + Days to hours, minutes, and seconds + + ctf2r -- call subroutine slctfr (ihour, imin, sec, rad, j) + dtf2r -- call subroutine sldtfr (ihour, imin, sec, rad, j) + Hours, minutes, seconds to radians + + cr2tf -- call subroutine slcrtf (ndp, angle, sign, ihmsf) + dr2tf -- call subroutine sldrtf (ndp, angle, sign, ihmsf) + Radians to hours, minutes, seconds + + caf2r -- call subroutine slcafr (ideg, iamin, asec, rad, j) + daf2r -- call subroutine sldafr (ideg, iamin, asec, rad, j) + Degrees, arcminutes, arcseconds to radians + + cr2af -- call subroutine slcraf (ndp, angle, sign, idmsf) + dr2af -- call subroutine sldraf (ndp, angle, sign, idmsf) + Radians to degrees, arcminutes, arcseconds +.fi +.ih +ANGLES, VECTORS AND ROTATION MATRICES +.nf + range -- r = slra1p (angle) +drange -- d = slda1p (angle) + Normalize angle into range [-pi,pi] + +ranorm -- r = slra2p (angle) +dranrm -- d = slda2p (angle) + Normalize angle into range [0,2pi] + + cs2c -- call subroutine slcs2c (a, b, v) + dcs2c -- call subroutine slds2c (a, b, v) + Spherical coordinates to [x,y,z] + + cc2s -- call subroutine slcc2s (v, a, b) + dcc2s -- call subroutine sldc2s (v, a, b) + [x,y,z] to spherical coordinates + + vdv -- r = slvdv (va, vb) + dvdv -- d = sldvdv (va, vb) + Scalar product of two 3-vectors + + vxv -- call subroutine slvxv (va, vb, vc) + dvxv -- call subroutine sldvxv (va, vb, vc) + Vector product of two 3-vectors + + vn -- call subroutine slvn (v, uv, vm) + dvn -- call subroutine sldvn (v, uv, vm) + Normalize a 3-vector also giving the modulus + + sep -- s = slsep (a1, b1, a2, b2) + dsep -- d = sldsep (a1, b1, a2, b2) + Angle between two points on a sphere + + bear -- s = slbear (a1, b1, a2, b2) + dbear -- d = sldber (a1, b1, a2, b2) + pav -- s = slpav (v1, v2) + dpav -- d = sldpav (v1, v2) + Direction of one point on a sphere seen from another + + euler -- call subroutine sleulr (order, phi, theta, psi, rmat) +deuler -- call subroutine sldeul (order, phi, theta, psi, rmat) + Form form rotation matrix from three Euler angles + + + av2m -- call subroutine slav2m (axvec, rmat) + dav2m -- call subroutine sldavm (axvec, rmat) + Form rotation matrix from axial vector + + m2av -- call subroutine slm2av (rmat, axvec) + dm2av -- call subroutine sldmav (rmat, axvec) + Determine axial vector from rotation matrix + + dmxv -- call subroutine sldmxv (dm, va, vb) + mxv -- call subroutine slmxv (rm, va, vb) + Rotate vector forwards + + imxv -- call subroutine slimxv (rm, va, vb) + dimxv -- call subroutine sldimv (dm, va, vb) + Rotate vector backwards + + dmxm -- call subroutine sldmxm (a, b, c) + mxm -- call subroutine slmxm (a, b, c) + Product of two 3X3 matrices + + cs2c6 -- call subroutine sls2c6 (a, b, r, ad, bd, rd, v) + ds2c6 -- call subroutine sldsc6 (a, b, r, ad, bd, rd, v) + Conversion of position/velocity from spherical to Cartesian + coordinates + + cc62s -- call subroutine slc62s (v, a, b, r, ad, bd, rd) + dc62s -- call subroutine sldc6s (v, a, b, r, ad, bd, rd) + Conversion of position/velocity from Cartesian to spherical + coordinates +.fi +.ih +CALENDARS +.nf + cldj -- call subroutine slcadj (iy, im, id, djm, j) + Gregorian calendar to Modified Julian Date + + caldj -- call subroutine slcadj (iy, im, id, djm, j) + Gregorian calendar to Modified Julian Date, permitting century by + default + + djcal -- call subroutine sldjca (ndp, djm, iymdf, j) + Modified Julian Date to Gregorian calendar, in a from convenient + for formatted output + + djcl -- call subroutine sldjcl (djm, iy, im, id, fd, j) + Modified Julian Date to Gregorian Year, Month, Day, Fraction + + calyd -- call subroutine slcayd (iy, im, id, ny, nd, j) + Calendar to year and day in year, permitting century default + + clyd -- call subroutine slclyd (iy, im, id, ny, nd, jstat) + Calendar to year and day in year + + epb -- d = slepb (date) + Modified Julian Date to Besselian Epoch + + epb2d -- d = sleb2d (epb) + Besselian epoch to Modified Julian Date + + epj -- d = slepj (date) + Modified Julian Date to Julian Epoch + + epj2d -- d = slej2d (epj) + Julian epoch to Modified Julian Date +.fi +.fi +.ih +TIMESCALES +.nf + gmst -- d = slgmst (ut1) + Conversion from Universal Time to siderial time + + gmsta -- d = slgmsa (date, ut) + Conversion from Universal Time to siderial time, rounding errors + minimized + + eqeqx -- d = sleqex (date) + Equation of the equinoxes + + dat -- d = sldat (utc) + Offset of Atomic Time from Coordinated Universal Time: + TAI - UTC + + dt -- d = sldt (epoch) + Approximate offset between dynamical time and universal time + + dtt -- d = sldtt (utc) + Offset of Terrestrial Time from Coordinated Universal Time: + TT - UTC + + rcc -- d = slrcc (tdb, ut1, wl, u, v) + Relativistic clock correction: TDB - TT +.fi +.ih +PRECESSION AND NUTATION +.nf + nut -- call subroutine slnut (date, rmatn) + Nutation matrix + + nutc -- call subroutine slnutc (date, dpsi, deps, eps0) + Longitude and obliquity components of nutation, mean obliquity + + prec -- call subroutine slprec (ep0, ep1, rmatp) + Precession matrix (IAU) + + precl -- call subroutine slprel (ep0, ep1, rmatp) + Precession matrix (suitable for long periods) + +prenut -- call subroutine slprnu (epoch, date, rmatpn) + Combined precession/nutation matrix + + prebn -- call subroutine slprbn (bep0, bep1, rmatp) + Precession matrix (old system) + +preces -- call subroutine slprce (system, ep0, ep1, ra, dc) + Precession in either the old or new system, character string + ep0 and ep1 + +precss -- call subroutine slprcs (system, ep0, ep1, ra, dc) + Precession in either the old or new system, integer ep0 and ep1 +.fi +.fi +.ih +PROPER MOTION +.nf + pm -- call subroutine slpm (r0, d0, pr, pd, px, rv, ep0, ep1, r1, d1) + Adjust for proper motion +.fi +.ih +FK4/5/ICRS CONVERSIONS +.nf + fk425 -- call subroutine slfk45 (r1950, d1950, dr1950, dd1950, + p1950, v1950, r2000, d2000, dr2000, dd2000, p2000, v2000) + Convert B1950.0 FK4 star data to J2000.0 FK5 + + fk45z -- call subroutine slf45z (r1950, d1950, bepoch, r2000, d2000) + Convert B1950.0 FK4 position to J2000.0 FK5 assuming zero proper + motion in an inertial frame and no parallax + + fk524 -- call subroutine slfk54 (r2000, d2000, dr2000, dd2000, + p2000, v2000, r1950, d1950, dr1950, dd1950, p1950, v1950) + Convert J2000.0 FK5 star data to B1950.0 FK4 + + fk54z -- call subroutine slf54z (r2000, d2000, bepoch, r1950, d1950, + dr1950, dd1950) + Convert J2000.0 FK5 star data to B1950.0 FK4 assuming zero proper + motion in an inertial frame and no parallax + + fk52h -- call subroutine slfk5h (r5, d5, dr5, dd5, rh, dh, drh, ddh) + Convert J2000.0 FK5 star data to ICRS J2000.0 data + + fk5hz -- call subroutine slf5hz (r5, d5, epoch, rh, dh) + Convert J2000.0 FK5 star data to ICRS J2000.0 data assuming + no Hipparcos proper motion. + + h2fk5 -- call subroutine slhfk5 (rh, dh, drh, ddh, r5, d5, dr5, dd5) + Convert ICRS J2000.0 data to J2000.0 Fk5 star data. + + hfk5z -- call subroutine slhf5z (rh, dh, epoch, r5, d5) + Convert ICRS J2000.0 data to J2000.0 Fk5 star data assuming no + Hipparchos proper motion. + + dbjin -- call subroutine sldbji (string, nstrt, dreslt, j1, j2) + Like dfltin but with extensions to accept leading 'B' and 'J' + + kbj -- call subroutine slkbj (jb, e, k, j) + Select epoch prefix 'B' or 'J' + + epco -- d = slepco (k0, k, e) + Convert an epoch into the appropriate form 'B' or 'J' +.fi +.ih +ELLIPTIC ABERRATIONS +.nf + etrms -- call subroutine sletrm (ep, ev) + E-terms + + subet -- call subroutine slsuet (rc, dc, eq, rm, dm) + Remove the E-terms + + addet -- call subroutine sladet (rm, dm, eq, rc, dc) + Add the E-terms +.fi +.ih +GEOCENTRIC COORDINATES +.nf + obs -- call subroutine slobs (n, c, name, w, p, h) + Interrogate list of observatory parameters + + geoc -- call subroutine slgeoc (p, h, r, z) + Convert geodetic position to geocentric + + pvobs -- call subroutine slpvob (p, h, stl, pv) + Position and velocity of observatory +.fi +.ih +APPARENT AND OBSERVED PLACE +.nf + map -- call subroutine slmap (rm, dm, pr, pd, px, rv, eq, date, ra, da) + Mean place to geocentric apparent place + + mappa -- call subroutine slmapa (eq, date, amprms) + Precompute mean to apparent parameters + + mapqk -- call subroutine slmapq (rm, dm, pr, pd, px, rv, amprms, ra, da) + Mean to apparent place using precomputed parameters + +mapqkz -- call subroutine slmapz (rm, dm, amprms, ra, da) + Mean to apparent place using precomputed parameters, for zero + proper motion, parallax, and radial velocity + + amp -- call subroutine slamp (ra, da, date, eq, rm, dm) + Geocentric apparent place to mean place + + ampqk -- call subroutine slampq (ra, da, amprms, rm, dm) + Apparent to mean place using precomputed parameters + + aop -- call subroutine slaop (rap, dap, date, dut, elongm, phim, hm, + xp, yp, tdk, pmb, rh, wl, tlr, aob, zob, hob, dob, rob) + Apparent place to observed place + + aoppa -- call subroutine slaopa (date, dut, elongm, phim, hm, xy, yp, + tdk, pmb, rh, wl, tlr, aoprms) + Precompute apparent to observed parameters + +aoppat -- call subroutine slaopt (date, aoprms) + Update siderial time in apparent to observed parameters + + aopqk -- call subroutine slaopq (rap, dap, aoprms, aob, zob, hob, dob, rob) + Apparent to observed using precomputed parameters + + oap -- call subroutine sloap (type, ob1, ob2, date, dut, elongm, phim, + xp, yp, tdk, pmb, rh, wl, tlr, rap, dap) + Observed to apparent + + oapqk -- call subroutine sloapq (type, ob1, ob2, aoprms, rap, dap) + Observed to apparent using precomputed parameters + + polmo -- call subroutine slplmo (elongim, phim, xp, yp, elong, phi, daz) + Correct site longitude and latitude for polar motion +.fi +.ih +AZIMUTH AND ELEVATION +.nf + altaz -- call subroutine slalaz (ha, dec, phi, + Positions, velocities, etc. for an altazimuth mount + + e2h -- call subroutine sle2h (ha, dec, phi, az, el) + de2h -- call subroutine slde2h (ha, dec, phi, az, el) + Hour angle and declination to azimuth and elevation + + h2e -- call subroutine slh2e (az, el, phi, ha, dec) + dh2e -- call subroutine sldh2e (az, el, phi, ha, dec) + Azimuth and elevation to hour angle and declination + + pda2h -- call subroutine slpdah (p, d, a, h1, j1, h2, j2) + Hour angle corresponding to a given azimuth + + pdq2h -- call subroutine slpdqh (p, d, q, h1, j1, h2, j2) + Hour angle corresponding to a given parallactic angle + + pa -- d = slpa (ha, dec, phi) + Hour angle and declination to parallactic angle + + zd -- d = slzd (ha, dec, phi) + Hour angle and declination to zenith distance +.fi +.ih +REFRACTION AND AIR MASS +.nf + refro -- call subroutine slrfro (zobs, hm, tdk, pmb, rh, wl, phi, tlr, + eps, ref) + Change in zenith distance due to refraction + + refco -- call subroutine slrfco (hm, tdk, pmb, rh, wl, phi, tlr, eps, + refa, refb) + Constants for simple refraction model + +refcoq -- call subroutine slrfcq (tdk, pmb, rl, wl, refa, refb) + Constants for simple refraction model (quick version) + +atmdsp -- call subroutine slatmd (tdk, pmb, rh, wl1, a1, b1, wl2, a2, b2) + Adjust refraction constants for color + + refz -- call subroutine slrefz (zu, refa, refb, zr) + Unrefracted to refracted zenith distance, simple model + + refv -- call subroutine slrefv (vu, refa, refb, vr) + Unrefracted to refracted azimuth and elevation, simple model + +airmas -- d = slarms (zd) + Air mass +.fi +.ih +ECLIPTIC COORDINATES +.nf + ecmat -- call subroutine slecma (date, rmat) + Equatorial to ecliptic rotation matrix + + eqecl -- call subroutine sleqec (dr, dd, date, dl, db) + J2000.0 FK5 to ecliptic coordinates + + ecleq -- call subroutine sleceq (dl, db, date, dr, dd) + Ecliptic to J2000.0 FK5 coordinates +.fi +.ih +GALACTIC COORDINATES +.nf + eg50 -- call subroutine sleg50 (dr, dd, dl, db) + B1950.0 FK4 to galactic coordinates + + ge50 -- call subroutine slge50 (dl, db, dr, dd) + Galactic to B1950.0 FK4 coordinates + + eqgal -- call subroutine sleqga (dr, dd, dl, db) + J2000.0 FK5 to galactic coordinates + + galeq -- call subroutine slgaeq (dl, db, dr, dd) + Galactic to J2000.0 FK5 coordinates +.fi +.ih +SUPERGALACTIC COORDINATES +.nf +galsup -- call subroutine slgasu (dl, db, dsl, dsb) + Galactic to supergalactic coordinates + +supgal -- call subroutine slsuga (dsl, dsb, dl, db) + Supergalactic to galactic coordinates +.fi +.ih +EPHEMERIDES +.nf + dmoon -- call subroutine sldmon (date, pv) + Approximate geocentric position and velocity of moon + + earth -- call subroutine slerth (iy, id, fd, pv) + Approximate heliocentric position and velocity of earth + + evp -- call subroutine slevp (date, deqx, dvb, dpb, dvh, dph) + Barycentric and heliocentric velocity and position of earth + + moon -- call subroutine slmoon (iy, id, fd, pv) + Approximate geocentric position and velocity of moon + +planet -- call subroutine slplnt (date, np, pv, jstat) + Approximate heliocentric position and velocity of planet + +rdplan -- call subroutine slrdpl (date, np, elong, phi, ra, dec, diam) + Approximate topocentric apparent place of a planet + +planel -- call subroutine slplnl (date, jform, epoch, orbinc, anode, + perih, aorg, e, aorl, dm, pv, jstat) + Approximate heliocentric position and velocity of planet + +plante -- call subroutine slplte (date, elong, phi, jform, epoch, orbinc, + anode, perih, aorq, e, aorl, dm, ra, dec, r, jstat) + Approximate topocentric apparent place of a planet + + pv2el -- call subroutine slpvel (pv, date, pmass, jformr, jform, epoch, + orbinc, anode, perih, aorg, e, aorl, dm, jstat) + Convert J2000 position and velocity to equivalent osculating elements + + el2ue -- call subroutine slelue (date, jform, epoch, orbinc, anode, perih, + aorq, e, aorl, dm, u, jstat) + Convert conventional osculating orbital elements into universal + form. + + ue2el -- call subroutine slueel (u, jformr, jform, epoch, orbinc, anode, + perih, aorq, e, aorl, dm, jstat) + Convert universal elements into conventional heliocentric osculating + form. + + pv2ue -- call subroutine slpvue (pv, date, pmass, u, jstat) + Construct a universal element set based on instantaneous position + and velocity. + + ue2pv -- call subroutine sluepv (date, u, pv, jstat) + Compute heliocentric position and velocity of a planet, asteroid, or + comet, starting from orbital elements in the "universal variables" + form. + +pertel -- call subroutine slprtl (jform, date0, date1, epoch0, epoch1, + orbi0, anode0, perih0, aorq0, e0, am0, epoch1, orbi1, anode1, + perih1, aorq1m e1, am1, jstat) + Update the osculating elements of a comet or asteroid by applying + planetary perturbations. + +pertue -- call subroutine slprue (date, u, jstat) + Update universal elements of a comet or asteroid by applying planetary + perturbations. +.fi +.ih +RADIAL VELOCITIES +rverot -- s = slrver (phi, ra, da, st) + Velocity component due to rotation of the earth + + ecor -- call subroutine slecor (rm, dm, iy, id, fd, rv, tl) + Components of velocity and light time due to earth orbital motion + +rvlsrd -- r = slrvld (r2000, d2000) + Velocity component due to solar motion wrt dynamical LSR + +rvlsrk -- r = slrvlk (r2000, d2000) + Velocity component due to solar motion wrt kinematical LSR + +rvgalc -- r = slrvga (r2000, d2000) + Velocity component due to rotation of the Galaxy + + rvlg -- r = slrvlg (r2000, d2000) + Velocity component due to rotation and translation of the Galaxy, + relative to the mean motion of the local group +.fi +.ih +ASTROMETRY +.nf + s2tp -- call subroutine sls2tp (ra, dec, raz, decz, xi, eta, j) + ds2tp -- call subroutine sldstp (ra, dec, raz, decz, xi, eta, j) + Transform spherical into tangent plane coordinates + + v2tp -- call subroutine slv2tp (v, v0, xi, eta, j) + dv2tp -- call subroutine sldvtp (v, v0, xi, eta, j) + Transform [x,y,z] into tangent plane coordinates + + tp2s -- call subroutine sltp2s (xi, eta, raz, decz, ra, dec) + dtp2s -- call subroutine sldtps (xi, eta, raz, decz, ra, dec) + Transform tangent plane into spherical coordinates + + tp2v -- call subroutine sltp2v (xi, eta, v0, v) + dtp2v -- call subroutine sldtpv (xi, eta, v0, v) + Transform tangent plane coordinates into [x,y,z] + + tps2c -- call subroutine sltpsc (xi, eta, ra, dec, raz1, decz1, + raz2, decz2, n) +dtps2c -- call subroutine sldpsc (xi, eta, ra, dec, raz1, decz1, + raz2, decz2, n) + Get plate center from tangent plane and spherical coordinates + + tpv2c -- call subroutine sltpvc (xi, eta, v, v01, v02, n) +dtpv2c -- call subroutine sldpvc (xi, eta, v, v01, v02, n) + Get plate center from [x,y,x] and tangent plane coordinates + + pcd -- call subroutine slpcd (disco, x, y) + Apply pincushion/barrel distortion + + unpcd -- call subroutine slupcd (disco, x, y) + Remove pincushion/barrel distortion + + fitxy -- call subroutine slftxy (itype, np, xye, xym, coeffs, j) + Fit a linear model to relate two sets of [x,y] coordinates + + pxy -- call subroutine slpxy (np, xye, xym, coeffs, xyp, + xrms, yrms, rrms) + Compute predicted coordinates and residuals + + invf -- call subroutine slinvf (fwds, bkwds, j) + Invert a linear model + + xy2xy -- call subroutine slxyxy (x1, y1, coeffs, x2, y2) + Transform one set of [x,y] coordinates + + dcmpf -- call subroutine sldcmf (coeffs, xz, yz, xs, ys, perp, orient) + Decompose a linear fit into geometric parameters +.fi +.ih +NUMERICAL METHODS +.nf + smat -- call subroutine slsmat (n, a, y, d, jf, iw) + dmat -- call subroutine sldmat (n, a, y, d, jf, iw) + Matrix inversion and solution of simultaneous equations + + svd -- call subroutine slsvd (m, n, mp, np, a, w, v, work, jstat) + Singular value decomposition of a matrix + +svdsol -- call subroutine slsvds (m, n, mp, np, b, u, w, v, work, x) + Solution from a given vector plus SVD + +svdcov -- call subroutine slsvdc (n, np, nc, w, v, work, cvm) + Covariance matrix from SVD +.fi +.endhelp |