wncard |
Table of contents
ProcedureWNCARD ( Cardinality of a double precision window ) INTEGER FUNCTION WNCARD ( WINDOW ) AbstractReturn the cardinality (number of intervals) of a double precision window. Required_ReadingWINDOWS KeywordsWINDOWS DeclarationsIMPLICIT NONE INTEGER LBCELL PARAMETER ( LBCELL = -5 ) DOUBLE PRECISION WINDOW ( LBCELL:* ) Brief_I/OVARIABLE I/O DESCRIPTION -------- --- -------------------------------------------------- WINDOW I Input window. The function returns the cardinality of the input window. Detailed_InputWINDOW is a window containing zero or more intervals. Detailed_OutputThe function returns the cardinality of (number of intervals in) the input window. ParametersNone. Exceptions1) If the number of elements in WINDOW is not even, the error SPICE(INVALIDSIZE) is signaled. FilesNone. ParticularsThe window cardinality (WNCARD) function simply wraps a CARD call then divides the result by 2. A common error when using the SPICE windows function is to use the CARDD value as the number of window intervals rather than the CARDD/2 value. 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) The following code example demonstrates how to insert an interval into an existing double precision SPICE window, and how to loop over all its intervals to extract their left and right points. Example code begins here. PROGRAM WNCARD_EX1 IMPLICIT NONE C C SPICELIB functions. C INTEGER WNCARD C C Local parameters. C INTEGER LBCELL PARAMETER ( LBCELL = -5 ) INTEGER WNSIZE PARAMETER ( WNSIZE = 10 ) C C Local variables. C DOUBLE PRECISION WINDOW ( LBCELL:WNSIZE ) DOUBLE PRECISION LEFT DOUBLE PRECISION RIGHT INTEGER I C C Validate the window with size WNSIZE and zero elements. C CALL WNVALD( WNSIZE, 0, WINDOW ) C C Insert the intervals C C [ 1, 3 ] [ 7, 11 ] [ 23, 27 ] C C into WINDOW. C CALL WNINSD( 1.D0, 3.D0, WINDOW ) CALL WNINSD( 7.D0, 11.D0, WINDOW ) CALL WNINSD( 23.D0, 27.D0, WINDOW ) C C Loop over the number of intervals in WINDOW, output C the LEFT and RIGHT endpoints for each interval. C DO I=1, WNCARD(WINDOW) CALL WNFETD( WINDOW, I, LEFT, RIGHT ) WRITE(*,'(A,I2,2(A,F8.3),A)') 'Interval', I, . ' [', LEFT,',', RIGHT, ']' END DO END When this program was executed on a Mac/Intel/gfortran/64-bit platform, the output was: Interval 1 [ 1.000, 3.000] Interval 2 [ 7.000, 11.000] Interval 3 [ 23.000, 27.000] RestrictionsNone. Literature_ReferencesNone. Author_and_InstitutionJ. Diaz del Rio (ODC Space) E.D. Wright (JPL) VersionSPICELIB Version 1.1.0, 06-JUL-2021 (JDR) Added IMPLICIT NONE statement. Updated to remove unnecessary lines of code in the Standard SPICE error handling CHKIN statements. Edited the header to comply to NAIF standard. Created complete code example from code fragment and added example's problem statement. SPICELIB Version 1.0.1, 24-APR-2010 (EDW) Minor edit to code comments eliminating typo. SPICELIB Version 1.0.0, 10-AUG-2007 (EDW) |
Fri Dec 31 18:37:07 2021