From d54fe7c1f704a63824c5bfa0ece65245572e9b27 Mon Sep 17 00:00:00 2001 From: Joseph Hunkeler Date: Wed, 4 Mar 2015 21:21:30 -0500 Subject: Initial commit --- src/slalib/sun67.htx/node228.html | 696 ++++++++++++++++++++++++++++++++++++++ 1 file changed, 696 insertions(+) create mode 100644 src/slalib/sun67.htx/node228.html (limited to 'src/slalib/sun67.htx/node228.html') diff --git a/src/slalib/sun67.htx/node228.html b/src/slalib/sun67.htx/node228.html new file mode 100644 index 0000000..f21ae0e --- /dev/null +++ b/src/slalib/sun67.htx/node228.html @@ -0,0 +1,696 @@ + + + + +SUMMARY OF CALLS + + + + + + + + + + + + +

+ +

+
+ Next: About this document ... +
+Up: SLALIB Positional Astronomy Library +
+ Previous: Numerical Methods +

+SUMMARY OF CALLS +

+The basic trigonometrical and numerical facilities are supplied in both single +and double precision versions. +Most of the more esoteric position and time routines use double precision +arguments only, even in cases where single precision would normally be adequate +in practice. +Certain routines with modest accuracy objectives are supplied in +single precision versions only. +In the calling sequences which follow, no attempt has been made +to distinguish between single and double precision argument names, +and frequently the same name is used on different occasions to +mean different things. +However, none of the routines uses a mixture of single and +double precision arguments; each routine is either wholly +single precision or wholly double precision. +

+In the classified list, below, +subroutine subprograms are those whose names and argument lists +are preceded by `CALL', whereas function subprograms are +those beginning `R=' (when the result is REAL) or `D=' (when +the result is DOUBLE PRECISION). +

+The list is, of course, merely for quick reference; inexperienced +users must refer to the detailed specifications given later. +In particular, don't guess whether arguments are single or +double precision; the result could be a program that happens to +works on one sort of machine but not on another. +

+
String Decoding

CALL sla_INTIN (STRING, NSTRT, IRESLT, JFLAG) +
+Convert free-format string into integer +
CALL sla_FLOTIN (STRING, NSTRT, RESLT, JFLAG) +
+CALL sla_DFLTIN (STRING, NSTRT, DRESLT, JFLAG) +
+Convert free-format string into floating-point number +
CALL sla_AFIN (STRING, NSTRT, RESLT, JFLAG) +
+CALL sla_DAFIN (STRING, NSTRT, DRESLT, JFLAG) +
+Convert free-format string from deg,arcmin,arcsec to radians

+
Sexagesimal Conversions

CALL sla_CTF2D (IHOUR, IMIN, SEC, DAYS, J) +
+CALL sla_DTF2D (IHOUR, IMIN, SEC, DAYS, J) +
+Hours, minutes, seconds to days +
CALL sla_CD2TF (NDP, DAYS, SIGN, IHMSF) +
+CALL sla_DD2TF (NDP, DAYS, SIGN, IHMSF) +
+Days to hours, minutes, seconds +
CALL sla_CTF2R (IHOUR, IMIN, SEC, RAD, J) +
+CALL sla_DTF2R (IHOUR, IMIN, SEC, RAD, J) +
+Hours, minutes, seconds to radians +
CALL sla_CR2TF (NDP, ANGLE, SIGN, IHMSF) +
+CALL sla_DR2TF (NDP, ANGLE, SIGN, IHMSF) +
+Radians to hours, minutes, seconds +
CALL sla_CAF2R (IDEG, IAMIN, ASEC, RAD, J) +
+CALL sla_DAF2R (IDEG, IAMIN, ASEC, RAD, J) +
+Degrees, arcminutes, arcseconds to radians +
CALL sla_CR2AF (NDP, ANGLE, SIGN, IDMSF) +
+CALL sla_DR2AF (NDP, ANGLE, SIGN, IDMSF) +
+Radians to degrees, arcminutes, arcseconds

+
Angles, Vectors and Rotation Matrices

R = sla_RANGE (ANGLE) +
+D = sla_DRANGE (ANGLE) +
+Normalize angle into range $\pm \pi$
R = sla_RANORM (ANGLE) +
+D = sla_DRANRM (ANGLE) +
+Normalize angle into range $0\!-\!2\pi$
CALL sla_CS2C (A, B, V) +
+CALL sla_DCS2C (A, B, V) +
+Spherical coordinates to $[\,x,y,z\,]$
CALL sla_CC2S (V, A, B) +
+CALL sla_DCC2S (V, A, B) +
+ $[\,x,y,z\,]$ to spherical coordinates +
R = sla_VDV (VA, VB) +
+D = sla_DVDV (VA, VB) +
+Scalar product of two 3-vectors +
CALL sla_VXV (VA, VB, VC) +
+CALL sla_DVXV (VA, VB, VC) +
+Vector product of two 3-vectors +
CALL sla_VN (V, UV, VM) +
+CALL sla_DVN (V, UV, VM) +
+Normalize a 3-vector also giving the modulus +
R = sla_SEP (A1, B1, A2, B2) +
+D = sla_DSEP (A1, B1, A2, B2) +
+Angle between two points on a sphere +
R = sla_BEAR (A1, B1, A2, B2) +
+D = sla_DBEAR (A1, B1, A2, B2) +
+Direction of one point on a sphere seen from another +
R = sla_PAV (V1, V2) +
+D = sla_DPAV (V1, V2) +
+Position-angle of one $[\,x,y,z\,]$ with respect to another +
CALL sla_EULER (ORDER, PHI, THETA, PSI, RMAT) +
+CALL sla_DEULER (ORDER, PHI, THETA, PSI, RMAT) +
+Form rotation matrix from three Euler angles +
CALL sla_AV2M (AXVEC, RMAT) +
+CALL sla_DAV2M (AXVEC, RMAT) +
+Form rotation matrix from axial vector +
CALL sla_M2AV (RMAT, AXVEC) +
+CALL sla_DM2AV (RMAT, AXVEC) +
+Determine axial vector from rotation matrix +
CALL sla_MXV (RM, VA, VB) +
+CALL sla_DMXV (DM, VA, VB) +
+Rotate vector forwards +
CALL sla_IMXV (RM, VA, VB) +
+CALL sla_DIMXV (DM, VA, VB) +
+Rotate vector backwards +
CALL sla_MXM (A, B, C) +
+CALL sla_DMXM (A, B, C) +
+Product of two 3x3 matrices +
CALL sla_CS2C6 (A, B, R, AD, BD, RD, V) +
+CALL sla_DS2C6 (A, B, R, AD, BD, RD, V) +
+Conversion of position and velocity in spherical + coordinates to Cartesian coordinates +
CALL sla_CC62S (V, A, B, R, AD, BD, RD) +
+CALL sla_DC62S (V, A, B, R, AD, BD, RD) +
+Conversion of position and velocity in Cartesian + coordinates to spherical coordinates

+
Calendars

CALL sla_CLDJ (IY, IM, ID, DJM, J) +
+Gregorian Calendar to Modified Julian Date +
CALL sla_CALDJ (IY, IM, ID, DJM, J) +
+Gregorian Calendar to Modified Julian Date, + permitting century default +
CALL sla_DJCAL (NDP, DJM, IYMDF, J) +
+Modified Julian Date to Gregorian Calendar, + in a form convenient for formatted output +
CALL sla_DJCL (DJM, IY, IM, ID, FD, J) +
+Modified Julian Date to Gregorian Year, Month, Day, Fraction +
CALL sla_CALYD (IY, IM, ID, NY, ND, J) +
+Calendar to year and day in year, permitting century default +
CALL sla_CLYD (IY, IM, ID, NY, ND, J) +
+Calendar to year and day in year +
D = sla_EPB (DATE) +
+Modified Julian Date to Besselian Epoch +
D = sla_EPB2D (EPB) +
+Besselian Epoch to Modified Julian Date +
D = sla_EPJ (DATE) +
+Modified Julian Date to Julian Epoch +
D = sla_EPJ2D (EPJ) +
+Julian Epoch to Modified Julian Date

+
Timescales

D = sla_GMST (UT1) +
+Conversion from Universal Time to sidereal time +
D = sla_GMSTA (DATE, UT1) +
+Conversion from Universal Time to sidereal time, rounding errors minimized +
D = sla_EQEQX (DATE) +
+Equation of the equinoxes +
D = sla_DAT (DJU) +
+Offset of Atomic Time from Coordinated Universal Time: TAI-UTC +
D = sla_DT (EPOCH) +
+Approximate offset between dynamical time and universal time +
D = sla_DTT (DJU) +
+Offset of Terrestrial Time from Coordinated Universal Time: TT-UTC +
D = sla_RCC (TDB, UT1, WL, U, V) +
+Relativistic clock correction: TDB-TT

+
Precession and Nutation

CALL sla_NUT (DATE, RMATN) +
+Nutation matrix +
CALL sla_NUTC (DATE, DPSI, DEPS, EPS0) +
+Longitude and obliquity components of nutation, and + mean obliquity +
CALL sla_PREC (EP0, EP1, RMATP) +
+Precession matrix (IAU) +
CALL sla_PRECL (EP0, EP1, RMATP) +
+Precession matrix (suitable for long periods) +
CALL sla_PRENUT (EPOCH, DATE, RMATPN) +
+Combined precession/nutation matrix +
CALL sla_PREBN (BEP0, BEP1, RMATP) +
+Precession matrix, old system +
CALL sla_PRECES (SYSTEM, EP0, EP1, RA, DC) +
+Precession, in either the old or the new system

+
Proper Motion

CALL sla_PM (R0, D0, PR, PD, PX, RV, EP0, EP1, R1, D1) +
+Adjust for proper motion

+
FK4/FK5/Hipparcos Conversions

CALL sla_FK425 ( + R1950, D1950, DR1950, DD1950, P1950, V1950, + R2000, D2000, DR2000, DD2000, P2000, V2000) +
+Convert B1950.0 FK4 star data to J2000.0 FK5 +
CALL sla_FK45Z (R1950, D1950, EPOCH, R2000, D2000) +
+Convert B1950.0 FK4 position to J2000.0 FK5 assuming zero + FK5 proper motion and no parallax +
CALL sla_FK524 ( + R2000, D2000, DR2000, DD2000, P2000, V2000, + R1950, D1950, DR1950, DD1950, P1950, V1950) +
+Convert J2000.0 FK5 star data to B1950.0 FK4 +
CALL sla_FK54Z (R2000, D2000, BEPOCH, + R1950, D1950, DR1950, DD1950) +
+Convert J2000.0 FK5 position to B1950.0 FK4 assuming zero + FK5 proper motion and no parallax +
CALL sla_FK52H (R5, D5, DR5, DD5, RH, DH, DRH, DDH) +
+Convert J2000.0 FK5 star data to Hipparcos +
CALL sla_FK5HZ (R5, D5, EPOCH, RH, DH ) +
+Convert J2000.0 FK5 position to Hipparcos assuming zero Hipparcos + proper motion +
CALL sla_H2FK5 (RH, DH, DRH, DDH, R5, D5, DR5, DD5) +
+Convert Hipparcos star data to J2000.0 FK5 +
CALL sla_HFK5Z (RH, DH, EPOCH, R5, D5, DR5, DD5) +
+Convert Hipparcos position to J2000.0 FK5 assuming zero Hipparcos + proper motion +
CALL sla_DBJIN (STRING, NSTRT, DRESLT, J1, J2) +
+Like sla_DFLTIN but with extensions to accept leading `B' and `J' +
CALL sla_KBJ (JB, E, K, J) +
+Select epoch prefix `B' or `J' +
D = sla_EPCO (K0, K, E) +
+Convert an epoch into the appropriate form - `B' or `J'

+
Elliptic Aberration

CALL sla_ETRMS (EP, EV) +
+E-terms +
CALL sla_SUBET (RC, DC, EQ, RM, DM) +
+Remove the E-terms +
CALL sla_ADDET (RM, DM, EQ, RC, DC) +
+Add the E-terms

+
Geographical and Geocentric Coordinates

CALL sla_OBS (NUMBER, ID, NAME, WLONG, PHI, HEIGHT) +
+Interrogate list of observatory parameters +
CALL sla_GEOC (P, H, R, Z) +
+Convert geodetic position to geocentric +
CALL sla_POLMO (ELONGM, PHIM, XP, YP, ELONG, PHI, DAZ) +
+Polar motion +
CALL sla_PVOBS (P, H, STL, PV) +
+Position and velocity of observatory

+
Apparent and Observed Place

CALL sla_MAP (RM, DM, PR, PD, PX, RV, EQ, DATE, RA, DA) +
+Mean place to geocentric apparent place +
CALL sla_MAPPA (EQ, DATE, AMPRMS) +
+Precompute mean to apparent parameters +
CALL sla_MAPQK (RM, DM, PR, PD, PX, RV, AMPRMS, RA, DA) +
+Mean to apparent using precomputed parameters +
CALL sla_MAPQKZ (RM, DM, AMPRMS, RA, DA) +
+Mean to apparent using precomputed parameters, for zero proper + motion, parallax and radial velocity +
CALL sla_AMP (RA, DA, DATE, EQ, RM, DM) +
+Geocentric apparent place to mean place +
CALL sla_AMPQK (RA, DA, AOPRMS, RM, DM) +
+Apparent to mean using precomputed parameters +
CALL sla_AOP ( + RAP, DAP, UTC, DUT, ELONGM, PHIM, HM, XP, YP, + TDK, PMB, RH, WL, TLR, AOB, ZOB, HOB, DOB, ROB) +
+Apparent place to observed place +
CALL sla_AOPPA ( + UTC, DUT, ELONGM, PHIM, HM, XP, YP, + TDK, PMB, RH, WL, TLR, AOPRMS) +
+Precompute apparent to observed parameters +
CALL sla_AOPPAT (UTC, AOPRMS) +
+Update sidereal time in apparent to observed parameters +
CALL sla_AOPQK (RAP, DAP, AOPRMS, AOB, ZOB, HOB, DOB, ROB) +
+Apparent to observed using precomputed parameters +
CALL sla_OAP ( + TYPE, OB1, OB2, UTC, DUT, ELONGM, PHIM, HM, XP, YP, + TDK, PMB, RH, WL, TLR, RAP, DAP) +
+Observed to apparent +
CALL sla_OAPQK (TYPE, OB1, OB2, AOPRMS, RA, DA) +
+Observed to apparent using precomputed parameters

+
Azimuth and Elevation

CALL sla_ALTAZ ( + HA, DEC, PHI, + AZ, AZD, AZDD, EL, ELD, ELDD, PA, PAD, PADD) +
+Positions, velocities etc. for an altazimuth mount +
CALL sla_E2H (HA, DEC, PHI, AZ, EL) +
+CALL sla_DE2H (HA, DEC, PHI, AZ, EL) +
+ $[\,h,\delta\,]$ to $[\,Az,El~]$
CALL sla_H2E (AZ, EL, PHI, HA, DEC) +
+CALL sla_DH2E (AZ, EL, PHI, HA, DEC) +
+ $[\,Az,El~]$ to $[\,h,\delta\,]$
CALL sla_PDA2H (P, D, A, H1, J1, H2, J2) +
+Hour Angle corresponding to a given azimuth +
CALL sla_PDQ2H (P, D, Q, H1, J1, H2, J2) +
+Hour Angle corresponding to a given parallactic angle +
D = sla_PA (HA, DEC, PHI) +
+ $[\,h,\delta\,]$ to parallactic angle +
D = sla_ZD (HA, DEC, PHI) +
+ $[\,h,\delta\,]$ to zenith distance

+
Refraction and Air Mass

CALL sla_REFRO (ZOBS, HM, TDK, PMB, RH, WL, PHI, TLR, EPS, REF) +
+Change in zenith distance due to refraction +
CALL sla_REFCO (HM, TDK, PMB, RH, WL, PHI, TLR, EPS, REFA, REFB) +
+Constants for simple refraction model (accurate) +
CALL sla_REFCOQ (TDK, PMB, RH, WL, REFA, REFB) +
+Constants for simple refraction model (fast) +
CALL sla_ATMDSP ( TDK, PMB, RH, WL1, REFA1, REFB1, WL2, REFA2, REFB2 ) +
+Adjust refraction constants for colour +
CALL sla_REFZ (ZU, REFA, REFB, ZR) +
+Unrefracted to refracted ZD, simple model +
CALL sla_REFV (VU, REFA, REFB, VR) +
+Unrefracted to refracted $[\,Az,El~]$ vector, simple model +
D = sla_AIRMAS (ZD) +
+Air mass

+
Ecliptic Coordinates

CALL sla_ECMAT (DATE, RMAT) +
+Equatorial to ecliptic rotation matrix +
CALL sla_EQECL (DR, DD, DATE, DL, DB) +
+J2000.0 `FK5' to ecliptic coordinates +
CALL sla_ECLEQ (DL, DB, DATE, DR, DD) +
+Ecliptic coordinates to J2000.0 `FK5'

+
Galactic Coordinates

CALL sla_EG50 (DR, DD, DL, DB) +
+B1950.0 `FK4' to galactic +
CALL sla_GE50 (DL, DB, DR, DD) +
+Galactic to B1950.0 `FK4' +
CALL sla_EQGAL (DR, DD, DL, DB) +
+J2000.0 `FK5' to galactic +
CALL sla_GALEQ (DL, DB, DR, DD) +
+Galactic to J2000.0 `FK5'

+
Supergalactic Coordinates

CALL sla_GALSUP (DL, DB, DSL, DSB) +
+Galactic to supergalactic +
CALL sla_SUPGAL (DSL, DSB, DL, DB) +
+Supergalactic to galactic

+
Ephemerides

CALL sla_DMOON (DATE, PV) +
+Approximate geocentric position and velocity of the Moon +
CALL sla_EARTH (IY, ID, FD, PV) +
+Approximate heliocentric position and velocity of the Earth +
CALL sla_EVP (DATE, DEQX, DVB, DPB, DVH, DPH) +
+Barycentric and heliocentric velocity and position of the Earth +
CALL sla_MOON (IY, ID, FD, PV) +
+Approximate geocentric position and velocity of the Moon +
CALL sla_PLANET (DATE, NP, PV, JSTAT) +
+Approximate heliocentric position and velocity of a planet +
CALL sla_RDPLAN (DATE, NP, ELONG, PHI, RA, DEC, DIAM) +
+Approximate topocentric apparent place of a planet +
CALL sla_PLANEL ( + DATE, JFORM, EPOCH, ORBINC, ANODE, PERIH, + AORQ, E, AORL, DM, PV, JSTAT) +
+Heliocentric position and velocity of a planet, asteroid or + comet, starting from orbital elements +
CALL sla_PLANTE ( + DATE, ELONG, PHI, JFORM, EPOCH, ORBINC, ANODE, + PERIH, AORQ, E, AORL, DM, RA, DEC, R, JSTAT) +
+Topocentric apparent place of a Solar-System object whose + heliocentric orbital elements are known +
CALL sla_PV2EL ( + PV, DATE, PMASS, JFORMR, JFORM, EPOCH, ORBINC, + ANODE, PERIH, AORQ, E, AORL, DM, JSTAT) +
+Orbital elements of a planet from instantaneous position and velocity +
CALL sla_PERTEL ( + JFORM, DATE0, DATE1, + EPOCH0, ORBI0, ANODE0, PERIH0, AORQ0, E0, AM0, + EPOCH1, ORBI1, ANODE1, PERIH1, AORQ1, E1, AM1, + JSTAT) +
+Update elements by applying perturbations +
CALL sla_EL2UE ( + DATE, JFORM, EPOCH, ORBINC, ANODE, + PERIH, AORQ, E, AORL, DM, + U, JSTAT) +
+Transform conventional elements to universal elements +
CALL sla_UE2EL ( + U, JFORMR, + JFORM, EPOCH, ORBINC, ANODE, PERIH, + AORQ, E, AORL, DM, JSTAT) +
+Transform universal elements to conventional elements +
CALL sla_PV2UE (PV, DATE, PMASS, U, JSTAT) +
+Package a position and velocity for use as universal elements +
CALL sla_UE2PV (DATE, U, PV, JSTAT) +
+Extract the position and velocity from universal elements +
CALL sla_PERTUE (DATE, U, JSTAT) +
+Update universal elements by applying perturbations +
R = sla_RVEROT (PHI, RA, DA, ST) +
+Velocity component due to rotation of the Earth +
CALL sla_ECOR (RM, DM, IY, ID, FD, RV, TL) +
+Components of velocity and light time due to Earth orbital motion +
R = sla_RVLSRD (R2000, D2000) +
+Velocity component due to solar motion wrt dynamical LSR +
R = sla_RVLSRK (R2000, D2000) +
+Velocity component due to solar motion wrt kinematical LSR +
R = sla_RVGALC (R2000, D2000) +
+Velocity component due to rotation of the Galaxy +
R = sla_RVLG (R2000, D2000) +
+Velocity component due to rotation and translation of the + Galaxy, relative to the mean motion of the local group

+
Astrometry

CALL sla_S2TP (RA, DEC, RAZ, DECZ, XI, ETA, J) +
+CALL sla_DS2TP (RA, DEC, RAZ, DECZ, XI, ETA, J) +
+Transform spherical coordinates into tangent plane +
CALL sla_V2TP (V, V0, XI, ETA, J) +
+CALL sla_DV2TP (V, V0, XI, ETA, J) +
+Transform $[\,x,y,z\,]$ into tangent plane coordinates +
CALL sla_DTP2S (XI, ETA, RAZ, DECZ, RA, DEC) +
+CALL sla_TP2S (XI, ETA, RAZ, DECZ, RA, DEC) +
+Transform tangent plane coordinates into spherical coordinates +
CALL sla_DTP2V (XI, ETA, V0, V) +
+CALL sla_TP2V (XI, ETA, V0, V) +
+Transform tangent plane coordinates into $[\,x,y,z\,]$
CALL sla_DTPS2C (XI, ETA, RA, DEC, RAZ1, DECZ1, RAZ2, DECZ2, N) +
+CALL sla_TPS2C (XI, ETA, RA, DEC, RAZ1, DECZ1, RAZ2, DECZ2, N) +
+Get plate centre from star $[\,\alpha,\delta\,]$ and tangent plane coordinates +
CALL sla_DTPV2C (XI, ETA, V, V01, V02, N) +
+CALL sla_TPV2C (XI, ETA, V, V01, V02, N) +
+Get plate centre from star $[\,x,y,z\,]$ and tangent plane coordinates +
CALL sla_PCD (DISCO, X, Y) +
+Apply pincushion/barrel distortion +
CALL sla_UNPCD (DISCO, X, Y) +
+Remove pincushion/barrel distortion +
CALL sla_FITXY (ITYPE, NP, XYE, XYM, COEFFS, J) +
+Fit a linear model to relate two sets of $[\,x,y\,]$ coordinates +
CALL sla_PXY (NP, XYE, XYM, COEFFS, XYP, XRMS, YRMS, RRMS) +
+Compute predicted coordinates and residuals +
CALL sla_INVF (FWDS, BKWDS, J) +
+Invert a linear model +
CALL sla_XY2XY (X1, Y1, COEFFS, X2, Y2) +
+Transform one $[\,x,y\,]$
CALL sla_DCMPF (COEFFS, XZ, YZ, XS, YS, PERP, ORIENT) +
+Decompose a linear fit into scales etc.

+
Numerical Methods

CALL sla_COMBN (NSEL, NCAND, LIST, J) +
+Next combination (subset from a specified number of items) +
CALL sla_PERMUT (N, ISTATE, IORDER, J) +
+Next permutation of a specified number of items +
CALL sla_SMAT (N, A, Y, D, JF, IW) +
+CALL sla_DMAT (N, A, Y, D, JF, IW) +
+Matrix inversion and solution of simultaneous equations +
CALL sla_SVD (M, N, MP, NP, A, W, V, WORK, JSTAT) +
+Singular value decomposition of a matrix +
CALL sla_SVDSOL (M, N, MP, NP, B, U, W, V, WORK, X) +
+Solution from given vector plus SVD +
CALL sla_SVDCOV (N, NP, NC, W, V, WORK, CVM) +
+Covariance matrix from SVD +
R = sla_RANDOM (SEED) +
+Generate pseudo-random real number in the range $0 \leq x < 1$
R = sla_GRESID (S) +
+Generate pseudo-random normal deviate ( $\equiv$ `Gaussian residual')

+
Real-time

CALL sla_WAIT (DELAY) +
+Interval wait

+ +

+
+ Next: About this document ... +
+Up: SLALIB Positional Astronomy Library +
+ Previous: Numerical Methods +

+SLALIB --- Positional Astronomy Library
Starlink User Note 67
P. T. Wallace
12 October 1999
E-mail:ptw@star.rl.ac.uk +

+ + -- cgit