spke19 |
Table of contents
ProcedureSPKE19 ( SPK, evaluate record, type 19 ) SUBROUTINE SPKE19 ( ET, RECORD, STATE ) AbstractEvaluate a single data record from a type 19 SPK segment. Required_ReadingSPK KeywordsEPHEMERIS SPK DeclarationsIMPLICIT NONE INCLUDE 'spk19.inc' INCLUDE 'spkrec.inc' DOUBLE PRECISION ET DOUBLE PRECISION RECORD ( * ) DOUBLE PRECISION STATE ( 6 ) Brief_I/OVARIABLE I/O DESCRIPTION -------- --- -------------------------------------------------- MAXREC P Maximum size of SPK record. ET I Epoch for which a state is desired. RECORD I Record from a type 19 SPK segment valid for ET. STATE O State (position and velocity) at epoch ET. Detailed_InputET is the epoch for which a state vector is desired, expressed as seconds past J2000 TDB. RECORD is a record from a type 19 SPK segment which, when evaluated at epoch ET, will give the state (position and velocity) of some body, relative to some center, in the reference frame associated with the data packets. Usually the body, center and frame are those of the SPK segment from which the packets were read. The structure of the record is as follows: +----------------------+ | subtype code | +----------------------+ | number of packets (n)| +----------------------+ | packet 1 | +----------------------+ . . . +----------------------+ | packet n | +----------------------+ | epoch 1 | +----------------------+ . . . +----------------------+ | epoch n | +----------------------+ Detailed_OutputSTATE is the state vector at epoch ET. Its contents are, in order, X, Y, Z, X', Y', and Z'. Units are km and km/sec. ParametersMAXREC is the maximum size of SPK record. See the SPICELIB routine SPKPVN for details. Exceptions1) Most types of errors in the input record cannot be diagnosed by this routine. This routine assumes that the input record is valid. 2) If the subtype code in the input record is invalid, an error is signaled by a routine in the call tree of this routine. FilesNone. ParticularsThe exact format and structure of type 19 (ESOC/DDID piecewise interpolation) SPK segments is described in the SPK Required Reading. ExamplesThe SPKEnn routines are almost always used in conjunction with the corresponding SPKRnn routines, which read the records from SPK files. The data returned by the SPKRnn routine is in a raw form, taken directly from the segment. As such, it will be not be directly useful to a user unless they have a complete understanding of the structure of the data type. Given that understanding, however, the SPKRnn routines could be used to "dump" and check segment data for a particular epoch before evaluating the record to obtain a state vector, as in the example that follows. C C Get a segment applicable to a specified body and epoch. C CALL SPKSFS ( BODY, ET, HANDLE, DESCR, IDENT, FOUND ) C C Look at parts of the descriptor. C CALL DAFUS ( DESCR, 2, 6, DCD, ICD ) CENTER = ICD( 2 ) REF = ICD( 3 ) TYPE = ICD( 4 ) IF ( TYPE .EQ. 19 ) THEN CALL SPKR19 ( HANDLE, DESCR, ET, RECORD ) . . Look at the RECORD data. . CALL SPKE19 ( ET, RECORD, STATE ) . . Check out the evaluated state. . END IF Restrictions1) This routine assumes that the input record is valid. Any checking of the input data is assumed to have been performed when the source SPK file was created. Literature_ReferencesNone. Author_and_InstitutionN.J. Bachman (JPL) J. Diaz del Rio (ODC Space) B.V. Semenov (JPL) VersionSPICELIB Version 2.0.1, 12-AUG-2021 (JDR) Edited the header to comply with NAIF standard. SPICELIB Version 2.0.0, 11-MAY-2015 (NJB) Updated to support subtype 2. Now performs computations in-line, rather than calling SPKE18. SPICELIB Version 1.0.0, 14-MAR-2014 (NJB) (BVS) |
Fri Dec 31 18:36:52 2021