| dsksrf |
|
Table of contents
Procedure
DSKSRF ( DSK, get surface IDs for body )
SUBROUTINE DSKSRF ( DSKFNM, BODYID, SRFIDS )
Abstract
Find the set of surface ID codes for all surfaces associated with
a given body in a specified DSK file.
Required_Reading
CELLS
DAS
DSK
NAIF_IDS
SETS
Keywords
COVERAGE
SURFACE
TOPOGRAPHY
Declarations
IMPLICIT NONE
INCLUDE 'dla.inc'
INCLUDE 'dskdsc.inc'
INTEGER LBCELL
PARAMETER ( LBCELL = -5 )
CHARACTER*(*) DSKFNM
INTEGER BODYID
INTEGER SRFIDS ( LBCELL : * )
Brief_I/O
VARIABLE I/O DESCRIPTION
-------- --- --------------------------------------------------
DSKFNM I Name of DSK file.
BODYID I Integer body ID code.
SRFIDS I-O Set of ID codes of surfaces in DSK file.
Detailed_Input
DSKFNM is the name of a DSK file. This file will be opened for
read access by this routine.
BODYID is the integer ID code of a body for which topographic
data are present in the specified DSK file.
SRFIDS is an initialized SPICE set data structure.
SRFIDS optionally may contain a set of surface ID codes
on input; on output, the ID codes already present in
SRFIDS will be combined with surface ID code set found
for the body designated by BODYID in the file DSKFNM.
If SRFIDS contains no data on input, its size and
cardinality still must be initialized.
Detailed_Output
SRFIDS is a SPICE set data structure that contains the union
of its contents upon input with the set of ID codes of
the surfaces associated with the body designated by
BODYID, for which segments were found in the indicated
DSK file.
The elements of SPICE sets are unique; each ID code in
SRFIDS appears only once, even if the DSK file contains
multiple segments for that ID code.
See the $Examples section below for a complete example
program showing how to retrieve body and surface ID codes
from a DSK file.
Parameters
None.
Exceptions
1) 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 DAS, the error SPICE(INVALIDARCHTYPE) is signaled.
3) If the input file is a binary DAS file of type other than DSK,
the error SPICE(INVALIDFILETYPE) is signaled.
4) If the DSK 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 SRFIDS is insufficient
to contain the actual number of ID codes of surfaces covered
by the indicated DSK file, the error SPICE(CELLTOOSMALL) is
signaled.
Files
See the description of the argument DSKFNM above.
Particulars
This routine provides an API via which applications can determine
the set of surfaces associated with a given body in a specified
DSK file. This routine is normally used together with DSKOBJ.
Examples
The 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) Examine a DSK file and identify the set of central bodies
associated with the segments in the file. For each body, find
the set of surfaces associated with that body.
Example code begins here.
C
C Examine a DSK file and identify the set of
C central bodies associated with the segments
C in the file. For each body, find the
C set of surfaces associated with that body.
C
PROGRAM DSKSRF_EX1
IMPLICIT NONE
C
C SPICELIB functions
C
INTEGER CARDI
C
C Local parameters
C
INTEGER LBCELL
PARAMETER ( LBCELL = -5 )
INTEGER FILSIZ
PARAMETER ( FILSIZ = 255 )
INTEGER MAXID
PARAMETER ( MAXID = 10000 )
C
C Local variables
C
CHARACTER*(FILSIZ) DSKFNM
INTEGER BODIDS ( LBCELL : MAXID )
INTEGER I
INTEGER J
INTEGER SRFIDS ( LBCELL : MAXID )
C
C Initialize body ID and surface ID cells.
C
CALL SSIZEI ( MAXID, BODIDS )
CALL SSIZEI ( MAXID, SRFIDS )
C
C Prompt for the name of a DSK file.
C
CALL PROMPT ( 'Enter name of DSK file > ', DSKFNM )
C
C Obtain body ID set for the DSK.
C
CALL DSKOBJ ( DSKFNM, BODIDS )
DO I = 1, CARDI( BODIDS )
WRITE (*,*) ' '
WRITE (*,*) 'Body ID: ', BODIDS(I)
C
C Get the surface IDs for the Ith body.
C
CALL DSKSRF ( DSKFNM, BODIDS(I), SRFIDS )
DO J = 1, CARDI( SRFIDS )
WRITE (*,*) ' Surface ID: ', SRFIDS(J)
END DO
END DO
END
When this program was executed on a Mac/Intel/gfortran/64-bit
platform, using the DSK file named phobos512.bds, the output
was:
Enter name of DSK file > phobos512.bds
Body ID: 401
Surface ID: 401
Restrictions
1) If an error occurs while this routine is updating the set
SRFIDS, the set may be corrupted.
Literature_References
None.
Author_and_Institution
N.J. Bachman (JPL)
J. Diaz del Rio (ODC Space)
Version
SPICELIB Version 1.1.0, 08-OCT-2021 (JDR) (NJB)
Changed input argument name "DSK" to "DSKFNM" for consistency
with other routines.
Bug fix: added call to FAILED after call to GETFAT.
Edited the header to comply with NAIF standard. Update problem
statement in $Examples section. Corrected short error message
in entries #2, #3 and #5 in $Exceptions section.
SPICELIB Version 1.0.0, 22-AUG-2016 (NJB)
|
Fri Dec 31 18:36:15 2021