/* $Procedure SPKS12 ( S/P Kernel, subset, type 12 ) */
/* Subroutine */ int spks12_(integer *handle, integer *baddr, integer *eaddr,
doublereal *begin, doublereal *end)
{
extern /* Subroutine */ int chkin_(char *, ftnlen), spks08_(integer *,
integer *, integer *, doublereal *, doublereal *), chkout_(char *,
ftnlen);
extern logical return_(void);
/* $ Abstract */
/* Extract a subset of the data in an SPK segment of type 12 */
/* into a new segment. */
/* $ Disclaimer */
/* THIS SOFTWARE AND ANY RELATED MATERIALS WERE CREATED BY THE */
/* CALIFORNIA INSTITUTE OF TECHNOLOGY (CALTECH) UNDER A U.S. */
/* GOVERNMENT CONTRACT WITH THE NATIONAL AERONAUTICS AND SPACE */
/* ADMINISTRATION (NASA). THE SOFTWARE IS TECHNOLOGY AND SOFTWARE */
/* PUBLICLY AVAILABLE UNDER U.S. EXPORT LAWS AND IS PROVIDED "AS-IS" */
/* TO THE RECIPIENT WITHOUT WARRANTY OF ANY KIND, INCLUDING ANY */
/* WARRANTIES OF PERFORMANCE OR MERCHANTABILITY OR FITNESS FOR A */
/* PARTICULAR USE OR PURPOSE (AS SET FORTH IN UNITED STATES UCC */
/* SECTIONS 2312-2313) OR FOR ANY PURPOSE WHATSOEVER, FOR THE */
/* SOFTWARE AND RELATED MATERIALS, HOWEVER USED. */
/* IN NO EVENT SHALL CALTECH, ITS JET PROPULSION LABORATORY, OR NASA */
/* BE LIABLE FOR ANY DAMAGES AND/OR COSTS, INCLUDING, BUT NOT */
/* LIMITED TO, INCIDENTAL OR CONSEQUENTIAL DAMAGES OF ANY KIND, */
/* INCLUDING ECONOMIC DAMAGE OR INJURY TO PROPERTY AND LOST PROFITS, */
/* REGARDLESS OF WHETHER CALTECH, JPL, OR NASA BE ADVISED, HAVE */
/* REASON TO KNOW, OR, IN FACT, SHALL KNOW OF THE POSSIBILITY. */
/* RECIPIENT BEARS ALL RISK RELATING TO QUALITY AND PERFORMANCE OF */
/* THE SOFTWARE AND ANY RELATED MATERIALS, AND AGREES TO INDEMNIFY */
/* CALTECH AND NASA FOR ALL THIRD-PARTY CLAIMS RESULTING FROM THE */
/* ACTIONS OF RECIPIENT IN THE USE OF THE SOFTWARE. */
/* $ Required_Reading */
/* SPK */
/* DAF */
/* $ Keywords */
/* EPHEMERIS */
/* $ Declarations */
/* $ Brief_I/O */
/* Variable I/O Description */
/* -------- --- -------------------------------------------------- */
/* HANDLE I Handle of file containing source segment. */
/* BADDR I Beginning address in file of source segment. */
/* EADDR I Ending address in file of source segment. */
/* BEGIN I Beginning (initial epoch) of subset. */
/* END I End (final epoch) of subset. */
/* $ Detailed_Input */
/* HANDLE, */
/* BADDR, */
/* EADDR are the file handle assigned to an SPK file, and the */
/* beginning and ending addresses of a segment within */
/* that file. Together they determine a complete set of */
/* ephemeris data, from which a subset is to be */
/* extracted. */
/* BEGIN, */
/* END are the initial and final epochs (ephemeris time) */
/* of the subset. */
/* The output segment will be padded to the left of */
/* BEGIN and the right of END with sufficient states to */
/* ensure that the segment yields an ephemeris identical */
/* to that given by the source segment. */
/* $ Detailed_Output */
/* See $Files section. */
/* $ Parameters */
/* None. */
/* $ Exceptions */
/* 1) This routine relies on the caller to ensure that the */
/* interval [BEGIN, END] is contained in the coverage */
/* interval of the segment. */
/* 2) If BEGIN > END, no data is written to the target file. */
/* $ Files */
/* Data is extracted from the file connected to the input */
/* handle, and written to the current DAF open for writing. */
/* The segment descriptor and summary must already have been written */
/* prior to calling this routine. The segment must be ended */
/* external to this routine. */
/* $ Particulars */
/* This routine is intended solely for use as a utility by the */
/* routine SPKSUB. */
/* It transfers a subset of a type 12 SPK data segment to */
/* a properly initialized segment of a second SPK file. */
/* The exact structure of a segment of data type 12 is described */
/* in the section on type 12 in the SPK Required Reading. */
/* $ Examples */
/* This routine is intended only for use as a utility by SPKSUB. */
/* To use this routine successfully, you must: */
/* Open the SPK file from which to extract data. */
/* Locate the segment from which data should be extracted. */
/* Open the SPK file to which this data should be written. */
/* Begin a new segment (array). */
/* Write the summary information for the array. */
/* Call this routine to extract the appropriate data from the */
/* SPK open for read. */
/* End the array to which this routine writes data. */
/* Much of this procedure is carried out by the routine SPKSUB. The */
/* examples of that routine illustrate more fully the process */
/* described above. */
/* $ Restrictions */
/* None. */
/* $ Literature_References */
/* None. */
/* $ Author_and_Institution */
/* N.J. Bachman (JPL) */
/* $ Version */
/* - SPICELIB Version 1.0.0, 25-FEB-2000 (NJB) */
/* -& */
/* $ Index_Entries */
/* subset type_12 spk segment */
/* -& */
/* SPICELIB functions */
/* Standard SPICE error handling. */
if (return_()) {
return 0;
} else {
chkin_("SPKS12", (ftnlen)6);
}
/* The type 8 subsetter knows how to do this job. */
spks08_(handle, baddr, eaddr, begin, end);
chkout_("SPKS12", (ftnlen)6);
return 0;
} /* spks12_ */
/* $Procedure SPKSUB ( S/P Kernel, subset ) */
/* Subroutine */ int spksub_(integer *handle, doublereal *descr, char *ident,
doublereal *begin, doublereal *end, integer *newh, ftnlen ident_len)
{
logical okay;
integer type__, baddr, eaddr;
doublereal alpha, omega;
extern /* Subroutine */ int chkin_(char *, ftnlen), dafps_(integer *,
integer *, doublereal *, integer *, doublereal *), dafus_(
doublereal *, integer *, integer *, doublereal *, integer *);
doublereal ndscr[5];
extern /* Subroutine */ int errdp_(char *, doublereal *, ftnlen), spks01_(
integer *, integer *, integer *, doublereal *, doublereal *),
spks02_(integer *, integer *, integer *, doublereal *, doublereal
*), spks03_(integer *, integer *, integer *, doublereal *,
doublereal *), spks10_(integer *, doublereal *, integer *,
doublereal *, char *, ftnlen), spks05_(integer *, integer *,
integer *, doublereal *, doublereal *), spks12_(integer *,
integer *, integer *, doublereal *, doublereal *), spks13_(
integer *, integer *, integer *, doublereal *, doublereal *),
spks08_(integer *, integer *, integer *, doublereal *, doublereal
*), spks09_(integer *, integer *, integer *, doublereal *,
doublereal *), spks14_(integer *, doublereal *, integer *,
doublereal *, char *, ftnlen), spks15_(integer *, integer *,
integer *, doublereal *, doublereal *), spks17_(integer *,
integer *, integer *, doublereal *, doublereal *), spks18_(
integer *, integer *, integer *, doublereal *, doublereal *),
spks19_(integer *, integer *, integer *, doublereal *, doublereal
*), spks20_(integer *, integer *, integer *, doublereal *,
doublereal *), spks21_(integer *, integer *, integer *,
doublereal *, doublereal *);
doublereal dc[2];
extern /* Subroutine */ int dafbna_(integer *, doublereal *, char *,
ftnlen);
integer ic[6];
extern /* Subroutine */ int dafena_(void), sigerr_(char *, ftnlen),
chkout_(char *, ftnlen), setmsg_(char *, ftnlen), errint_(char *,
integer *, ftnlen);
extern logical return_(void);
/* $ Abstract */
/* Extract a subset of the data in an SPK segment into a */
/* separate segment. */
/* $ Disclaimer */
/* THIS SOFTWARE AND ANY RELATED MATERIALS WERE CREATED BY THE */
/* CALIFORNIA INSTITUTE OF TECHNOLOGY (CALTECH) UNDER A U.S. */
/* GOVERNMENT CONTRACT WITH THE NATIONAL AERONAUTICS AND SPACE */
/* ADMINISTRATION (NASA). THE SOFTWARE IS TECHNOLOGY AND SOFTWARE */
/* PUBLICLY AVAILABLE UNDER U.S. EXPORT LAWS AND IS PROVIDED "AS-IS" */
/* TO THE RECIPIENT WITHOUT WARRANTY OF ANY KIND, INCLUDING ANY */
/* WARRANTIES OF PERFORMANCE OR MERCHANTABILITY OR FITNESS FOR A */
/* PARTICULAR USE OR PURPOSE (AS SET FORTH IN UNITED STATES UCC */
/* SECTIONS 2312-2313) OR FOR ANY PURPOSE WHATSOEVER, FOR THE */
/* SOFTWARE AND RELATED MATERIALS, HOWEVER USED. */
/* IN NO EVENT SHALL CALTECH, ITS JET PROPULSION LABORATORY, OR NASA */
/* BE LIABLE FOR ANY DAMAGES AND/OR COSTS, INCLUDING, BUT NOT */
/* LIMITED TO, INCIDENTAL OR CONSEQUENTIAL DAMAGES OF ANY KIND, */
/* INCLUDING ECONOMIC DAMAGE OR INJURY TO PROPERTY AND LOST PROFITS, */
/* REGARDLESS OF WHETHER CALTECH, JPL, OR NASA BE ADVISED, HAVE */
/* REASON TO KNOW, OR, IN FACT, SHALL KNOW OF THE POSSIBILITY. */
/* RECIPIENT BEARS ALL RISK RELATING TO QUALITY AND PERFORMANCE OF */
/* THE SOFTWARE AND ANY RELATED MATERIALS, AND AGREES TO INDEMNIFY */
/* CALTECH AND NASA FOR ALL THIRD-PARTY CLAIMS RESULTING FROM THE */
/* ACTIONS OF RECIPIENT IN THE USE OF THE SOFTWARE. */
/* $ Required_Reading */
/* SPK */
/* DAF */
/* $ Keywords */
/* EPHEMERIS */
/* $ Declarations */
/* $ Brief_I/O */
/* Variable I/O Description */
/* -------- --- -------------------------------------------------- */
/* HANDLE I Handle of source segment. */
/* DESCR I Descriptor of source segment. */
/* IDENT I Identifier of source segment. */
/* BEGIN I Beginning (initial epoch) of subset. */
/* END I End (final epoch) of subset. */
/* NEWH I Handle of new segment. */
/* $ Detailed_Input */
/* HANDLE, */
/* DESCR, */
/* IDENT are the file handle assigned to a SPK file, the */
/* descriptor for a segment within the file, and the */
/* identifier for that segment. Together they determine */
/* a complete set of ephemeris data, from which a */
/* subset is to be extracted. */
/* BEGIN, */
/* END are the initial and final epochs (ephemeris time) */
/* of the subset. */
/* NEWH is the file handle assigned to the file in which */
/* the new segment is to be written. The file must */
/* be open for write access. NEWH and HANDLE may refer */
/* to the same file. */
/* $ Detailed_Output */
/* See $Files section. */
/* $ Parameters */
/* None. */
/* $ Exceptions */
/* 1) If the condition */
/* ALPHA < BEGIN < END < OMEGA */
/* - - - */
/* is not satisfied (where ALPHA and OMEGA are the initial */
/* and final epochs of the segment respectively), the error */
/* 'SPICE(SPKNOTASUBSET)' is signaled. */
/* 2) If the segment type is not supported by the current */
/* version of SPKSUB, the error 'SPICE(SPKTYPENOTSUPP)' */
/* is signaled. */
/* $ Files */
/* A new segment, which contains a subset of the data in the */
/* segment specified by DESCR and HANDLE, is written to the SPK */
/* file attached to NEWH. */
/* $ Particulars */
/* Sometimes, the segments in official source files---planetary */
/* Developmental Ephemeris (DE) files, archival spacecraft */
/* ephemeris files, and so on---contain more data than is needed */
/* by a particular user. SPKSUB allows a user to extract from a */
/* segment the smallest amount of ephemeris data sufficient to */
/* cover a specific interval. */
/* The new segment is written with the same identifier as the */
/* original segment, and with the same descriptor, with the */
/* following components changed: */
/* 1) ALPHA and OMEGA (DCD(1) and DCD(2)) are assigned the values */
/* specified by BEGIN and END. */
/* 2) The beginning and ending segment addresses (ICD(5) and ICD(6)) */
/* are changed to reflect the location of the new segment. */
/* $ Examples */
/* In the following code fragment, the descriptor for each segment */
/* in a source SPK file is examined. For each segment that covers a */
/* specified time interval, the smallest possible subset of data */
/* from that segment, sufficient to cover the interval, is extracted */
/* into a custom SPK file. */
/* Assume that the source and custom files have been opened, for */
/* read and write access, with handles SRC and CUST respectively. */
/* CALL DAFBFS ( SRC ) */
/* CALL DAFFNA ( FOUND ) */
/* DO WHILE ( FOUND ) */
/* CALL DAFGS ( DESCR ) */
/* CALL DAFUS ( DESCR, 2, 6, DC, IC ) */
/* IF ( DC(1) .LE. BEGIN .AND. END .LE. DC(2) ) THEN */
/* CALL DAFGN ( IDENT ) */
/* CALL SPKSUB ( SRC, DESCR, IDENT, BEGIN, END, CUST ) */
/* END IF */
/* CALL DAFFNA ( FOUND ) */
/* END DO */
/* $ Restrictions */
/* 1) There is no way for SPKSUB to verify that the descriptor and */
/* identifier are the original ones for the segment. Changing */
/* the descriptor can cause the data in the new segment to be */
/* evaluated incorrectly; changing the identifier can destroy */
/* the path from the data back to its original source. */
/* $ Literature_References */
/* NAIF Document 168.0, "S- and P- Kernel (SPK) Specification and */
/* User's Guide" */
/* $ Author_and_Institution */
/* K.R. Gehringer (JPL) */
/* W.L. Taber (JPL) */
/* N.J. Bachman (JPL) */
/* J.M. Lynch (JPL) */
/* R.E. Thurman (JPL) */
/* I.M. Underwood (JPL) */
/* $ Version */
/* - SPICELIB Version 9.0.0, 23-DEC-2013 (NJB) */
/* The routine was updated to handle types 19, 20 and 21. Some */
/* minor changes were made to comments. */
/* - SPICELIB Version 8.0.0, 12-AUG-2002 (NJB) */
/* The routine was updated to handle type 18. */
/* - SPICELIB Version 7.0.0, 06-NOV-1999 (NJB) */
/* The routine was updated to handle types 12 and 13. */
/* - SPICELIB Version 6.0.0, 30-JUN-1997 (WLT) */
/* The routine was updated to handle types 10 and 17. */
/* - SPICELIB Version 5.0.0, 10-MAR-1995 (KRG) */
/* The routine was updated to handle type 14. */
/* - SPICELIB Version 4.0.0, 07-NOV-1994 (WLT) */
/* The routine was updated to handle type 15. */
/* - SPICELIB Version 3.0.0, 05-AUG-1993 (NJB) */
/* The routine was updated to handle types 08 and 09. */
/* - SPICELIB Version 2.0.0, 01-APR-1992 (JML) */
/* 1) The routine was updated to handle type 05. */
/* 2) DESCR was being used as both an input and output */
/* variable when it was only supposed to be used for */
/* input. A new local variable, NDSCR, was added where DESCR */
/* was being altered. */
/* - SPICELIB Version 1.0.1, 22-MAR-1990 (HAN) */
/* Literature references added to the header. */
/* - SPICELIB Version 1.0.0, 31-JAN-1990 (IMU) (RET) */
/* -& */
/* $ Index_Entries */
/* subset of spk file */
/* -& */
/* $ Revisions */
/* - SPICELIB Version 9.0.0, 23-DEC-2013 (NJB) */
/* The routine was updated to handle types 19, 20 and 21. Some */
/* minor changes were made to comments. */
/* - SPICELIB Version 8.0.0, 12-AUG-2002 (NJB) */
/* The routine was updated to handle type 18. */
/* - SPICELIB Version 6.0.0, 30-JUN-1997 (WLT) */
/* The routine was updated to handle types 10 and 17. */
/* - SPICELIB Version 5.0.0, 10-MAR-1995 (KRG) */
/* The routine was updated to handle type 14. */
/* - SPICELIB Version 4.0.0, 07-NOV-1994 (WLT) */
/* The routine was updated to handle type 15. */
/* - SPICELIB Version 3.0.0, 05-AUG-1993 (NJB) */
/* The routine was updated to handle types 08 and 09. */
/* - SPICELIB Version 2.0.0, 01-APR-1992 (JML) */
/* 1) The routine was updated to handle type 05. */
/* 2) DESCR was being used as both an input and output */
/* variable when it was only supposed to be used for */
/* input. A new local variable, NDSCR, was added where DESCR */
/* was being altered. */
/* -& */
/* SPICELIB functions */
/* Local variables */
/* Standard SPICE error handling. */
if (return_()) {
return 0;
} else {
chkin_("SPKSUB", (ftnlen)6);
}
/* Unpack the descriptor. */
dafus_(descr, &c__2, &c__6, dc, ic);
alpha = dc[0];
omega = dc[1];
type__ = ic[3];
baddr = ic[4];
eaddr = ic[5];
/* Make sure the epochs check out. */
okay = alpha <= *begin && *begin <= *end && *end <= omega;
if (! okay) {
setmsg_("Specified interval [#, #] is not a subset of segment interv"
"al [#, #].", (ftnlen)69);
errdp_("#", begin, (ftnlen)1);
errdp_("#", end, (ftnlen)1);
errdp_("#", &alpha, (ftnlen)1);
errdp_("#", &omega, (ftnlen)1);
sigerr_("SPICE(SPKNOTASUBSET)", (ftnlen)20);
chkout_("SPKSUB", (ftnlen)6);
return 0;
}
/* Begin the new segment, with a descriptor containing the subset */
/* epochs. */
dc[0] = *begin;
dc[1] = *end;
dafps_(&c__2, &c__6, dc, ic, ndscr);
/* Let the type-specific (SPKSnn) routines decide what to move. */
if (type__ == 1) {
dafbna_(newh, ndscr, ident, ident_len);
spks01_(handle, &baddr, &eaddr, begin, end);
dafena_();
} else if (type__ == 2) {
dafbna_(newh, ndscr, ident, ident_len);
spks02_(handle, &baddr, &eaddr, begin, end);
dafena_();
} else if (type__ == 3) {
dafbna_(newh, ndscr, ident, ident_len);
spks03_(handle, &baddr, &eaddr, begin, end);
dafena_();
/* Type 04 has not been yet been added to SPICELIB. */
/* ELSE IF ( TYPE .EQ. 04 ) THEN */
/* CALL DAFBNA ( NEWH, NDSCR, IDENT ) */
/* CALL SPKS04 ( HANDLE, BADDR, EADDR, BEGIN, END ) */
/* CALL DAFENA */
} else if (type__ == 5) {
dafbna_(newh, ndscr, ident, ident_len);
spks05_(handle, &baddr, &eaddr, begin, end);
dafena_();
} else if (type__ == 8) {
dafbna_(newh, ndscr, ident, ident_len);
spks08_(handle, &baddr, &eaddr, begin, end);
dafena_();
} else if (type__ == 9) {
dafbna_(newh, ndscr, ident, ident_len);
spks09_(handle, &baddr, &eaddr, begin, end);
dafena_();
} else if (type__ == 10) {
spks10_(handle, descr, newh, ndscr, ident, ident_len);
} else if (type__ == 12) {
dafbna_(newh, ndscr, ident, ident_len);
spks12_(handle, &baddr, &eaddr, begin, end);
dafena_();
} else if (type__ == 13) {
dafbna_(newh, ndscr, ident, ident_len);
spks13_(handle, &baddr, &eaddr, begin, end);
dafena_();
} else if (type__ == 14) {
spks14_(handle, descr, newh, ndscr, ident, ident_len);
} else if (type__ == 15) {
dafbna_(newh, ndscr, ident, ident_len);
spks15_(handle, &baddr, &eaddr, begin, end);
dafena_();
} else if (type__ == 17) {
dafbna_(newh, ndscr, ident, ident_len);
spks17_(handle, &baddr, &eaddr, begin, end);
dafena_();
} else if (type__ == 18) {
dafbna_(newh, ndscr, ident, ident_len);
spks18_(handle, &baddr, &eaddr, begin, end);
dafena_();
} else if (type__ == 19) {
dafbna_(newh, ndscr, ident, ident_len);
spks19_(handle, &baddr, &eaddr, begin, end);
dafena_();
} else if (type__ == 20) {
dafbna_(newh, ndscr, ident, ident_len);
spks20_(handle, &baddr, &eaddr, begin, end);
dafena_();
} else if (type__ == 21) {
dafbna_(newh, ndscr, ident, ident_len);
spks21_(handle, &baddr, &eaddr, begin, end);
dafena_();
} else {
setmsg_("SPK data type # is not supported.", (ftnlen)33);
errint_("#", &type__, (ftnlen)1);
sigerr_("SPICE(SPKTYPENOTSUPP)", (ftnlen)21);
chkout_("SPKSUB", (ftnlen)6);
return 0;
}
chkout_("SPKSUB", (ftnlen)6);
return 0;
} /* spksub_ */
