pckfrm |
Table of contents
ProcedurePCKFRM ( PCK, get reference frame class ID set ) SUBROUTINE PCKFRM ( PCKFNM, IDS ) AbstractFind the set of reference frame class ID codes of all frames in a specified binary PCK file. Required_ReadingCELLS DAF SETS PCK KeywordsORIENTATION UTILITY DeclarationsIMPLICIT NONE INTEGER LBCELL PARAMETER ( LBCELL = -5 ) CHARACTER*(*) PCKFNM INTEGER IDS ( LBCELL : * ) Brief_I/OVARIABLE I/O DESCRIPTION -------- --- -------------------------------------------------- PCKFNM I Name of PCK file. IDS I-O Set of frame class ID codes of frames in PCK file. Detailed_InputPCKFNM is the name of a binary PCK file. IDS is an initialized SPICE set data structure. IDS optionally may contain a set of ID codes on input; on output, the data already present in IDS will be combined with ID code set found for the file PCKFNM. If IDS contains no data on input, its size and cardinality still must be initialized. Detailed_OutputIDS is a SPICE set data structure which contains the union of its contents upon input with the set of reference frame class ID codes of each frame for which data are present in the indicated PCK file. The elements of SPICE sets are unique; hence each ID code in IDS appears only once, even if the PCK file contains multiple segments for that ID code. See the $Examples section below for a complete example program showing how to retrieve the ID codes from IDS. ParametersNone. Exceptions1) If the input file has transfer format, the error SPICE(INVALIDFORMAT) is signaled. 2) If the input file is not a transfer file but has architecture other than DAF, the error SPICE(INVALIDARCHTYPE) is signaled. 3) If the input file is a binary DAF file of type other than PCK, the error SPICE(INVALIDFILETYPE) is signaled. 4) If the PCK file cannot be opened or read, an error is signaled by a routine in the call tree of this routine. 5) If the size of the output set argument IDS is insufficient to contain the actual number of ID codes of frames covered by the indicated PCK file, an error is signaled by a routine in the call tree of this routine. FilesNone. ParticularsThis routine provides an API via which applications can determine the set of reference frames for which there are data in a specified PCK file. ExamplesThe numerical results shown for this example may differ across platforms. The results depend on the SPICE kernels used as input, the compiler and supporting libraries, and the machine specific arithmetic implementation. 1) Display the coverage for each frame in a specified PCK file. Find the set of frames in the file. Loop over the contents of the ID code set: find the coverage for each item in the set and display the coverage. Example code begins here. PROGRAM PCKFRM_EX1 IMPLICIT NONE C C SPICELIB functions C INTEGER WNCARD INTEGER CARDI C C Local parameters C C C Declare the coverage window. Make enough room C for MAXIV intervals. C INTEGER FILSIZ PARAMETER ( FILSIZ = 255 ) INTEGER LBCELL PARAMETER ( LBCELL = -5 ) INTEGER MAXIV PARAMETER ( MAXIV = 1000 ) INTEGER WINSIZ PARAMETER ( WINSIZ = 2 * MAXIV ) INTEGER TIMLEN PARAMETER ( TIMLEN = 50 ) INTEGER MAXFRM PARAMETER ( MAXFRM = 1000 ) C C Local variables C CHARACTER*(FILSIZ) LSK CHARACTER*(FILSIZ) PCKFNM CHARACTER*(TIMLEN) TIMSTR DOUBLE PRECISION B DOUBLE PRECISION COVER ( LBCELL : WINSIZ ) DOUBLE PRECISION E INTEGER I INTEGER IDS ( LBCELL : MAXFRM ) INTEGER J INTEGER NIV C C Load a leapseconds kernel for output time conversion. C PCKCOV itself does not require a leapseconds kernel. C CALL PROMPT ( 'Name of leapseconds kernel > ', LSK ) CALL FURNSH ( LSK ) C C Get name of PCK file. C CALL PROMPT ( 'Name of PCK file > ', PCKFNM ) C C Initialize the set IDS. C CALL SSIZEI ( MAXFRM, IDS ) C C Initialize the window COVER. C CALL SSIZED ( WINSIZ, COVER ) C C Find the set of frames in the PCK file. C CALL PCKFRM ( PCKFNM, IDS ) C C We want to display the coverage for each frame. Loop C over the contents of the ID code set, find the coverage C for each item in the set, and display the coverage. C DO I = 1, CARDI( IDS ) C C Find the coverage window for the current frame. C Empty the coverage window each time so C we don't include data for the previous frame. C CALL SCARDD ( 0, COVER ) CALL PCKCOV ( PCKFNM, IDS(I), COVER ) C C Get the number of intervals in the coverage C window. C NIV = WNCARD( COVER ) C C Display a simple banner. C WRITE (*,*) '========================================' WRITE (*,*) 'Coverage for reference frame ', IDS(I) C C Convert the coverage interval start and stop C times to TDB calendar strings. C DO J = 1, NIV C C Get the endpoints of the Jth interval. C CALL WNFETD ( COVER, J, B, E ) C C Convert the endpoints to TDB calendar C format time strings and display them. C CALL TIMOUT ( B, . 'YYYY MON DD HR:MN:SC.### ' // . '(TDB) ::TDB', . TIMSTR ) WRITE (*,*) ' ' WRITE (*,*) 'Interval: ', J WRITE (*,*) 'Start: ', TIMSTR CALL TIMOUT ( E, . 'YYYY MON DD HR:MN:SC.### ' // . '(TDB) ::TDB', . TIMSTR ) WRITE (*,*) 'Stop: ', TIMSTR WRITE (*,*) ' ' END DO WRITE (*,*) '========================================' END DO END When this program was executed on a Mac/Intel/gfortran/64-bit platform, using the LSK file named naif0012.tls, and the PCK file named earth_720101_070426.bpc, the output was: Name of leapseconds kernel > naif0012.tls Name of PCK file > earth_720101_070426.bpc ======================================== Coverage for reference frame 3000 Interval: 1 Start: 1962 JAN 20 00:00:41.184 (TDB) Stop: 2007 APR 26 00:01:05.185 (TDB) ======================================== Restrictions1) If an error occurs while this routine is updating the set IDS, the set may be corrupted. Literature_ReferencesNone. Author_and_InstitutionN.J. Bachman (JPL) J. Diaz del Rio (ODC Space) E.D. Wright (JPL) VersionSPICELIB Version 1.1.0, 08-OCT-2021 (JDR) (NJB) Changed input argument name "PCK" to "PCKFNM" for consistency with other routines. Bug fix: added call to FAILED after call to GETFAT. Edited the header to comply with NAIF standard. Added example's solution. Corrected short error message in entries #2 and #3 in $Exceptions section. SPICELIB Version 1.0.1, 03-JAN-2014 (EDW) Minor edits to $Procedure; clean trailing whitespace. SPICELIB Version 1.0.0, 01-DEC-2007 (NJB) |
Fri Dec 31 18:36:37 2021