wnfetd |
Table of contents
ProcedureWNFETD ( Fetch an interval from a DP window ) SUBROUTINE WNFETD ( WINDOW, N, LEFT, RIGHT ) AbstractFetch a particular interval from a double precision window. Required_ReadingWINDOWS KeywordsWINDOWS DeclarationsIMPLICIT NONE INTEGER LBCELL PARAMETER ( LBCELL = -5 ) DOUBLE PRECISION WINDOW ( LBCELL:* ) INTEGER N DOUBLE PRECISION LEFT DOUBLE PRECISION RIGHT Brief_I/OVARIABLE I/O DESCRIPTION -------- --- -------------------------------------------------- WINDOW I Input window. N I Index of interval to be fetched. LEFT, RIGHT O Left, right endpoints of the Nth interval. Detailed_InputWINDOW is a window containing zero or more intervals. N is the index of a particular interval within the window. Indices range from 1 to NINT, where NINT is the number of intervals in the window (CARDD(WINDOW)/2). Detailed_OutputLEFT, RIGHT are the left and right endpoints of the N'th interval in the input window. If the interval is not found, LEFT and RIGHT are not defined. ParametersNone. Exceptions1) If N is less than one, the error SPICE(NOINTERVAL) is signaled. 2) If the interval does not exist, i.e. N > CARDD(WINDOW)/2, the error SPICE(NOINTERVAL) is signaled. 3) The cardinality of the input WINDOW must be even. Left endpoints of stored intervals must be strictly greater than preceding right endpoints. Right endpoints must be greater than or equal to corresponding left endpoints. Invalid window data are not diagnosed by this routine and may lead to unpredictable results. FilesNone. ParticularsNone. 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 WNFETD_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_InstitutionN.J. Bachman (JPL) J. Diaz del Rio (ODC Space) W.L. Taber (JPL) I.M. Underwood (JPL) E.D. Wright (JPL) VersionSPICELIB Version 1.1.0, 15-MAR-2021 (JDR) (NJB) Added IMPLICIT NONE statement. Edited the header to comply with NAIF standard. Added complete code example, problem statement and solution. Added entry #3 in $Exceptions section. Improved description of argument N in $Detailed_Input. SPICELIB Version 1.0.2, 30-JUL-2007 (EDW) Removed erroneous description in the $Examples section indicating "Undefined" as a return state after an error event caused by an invalid value of N. SPICELIB Version 1.0.1, 10-MAR-1992 (WLT) Comment section for permuted index source lines was added following the header. SPICELIB Version 1.0.0, 31-JAN-1990 (WLT) (IMU) |
Fri Dec 31 18:37:07 2021