reordi |
Table of contents
ProcedureREORDI ( Reorder an integer array ) SUBROUTINE REORDI ( IORDER, NDIM, ARRAY ) AbstractReorder the elements of an integer array according to a given order vector. Required_ReadingNone. KeywordsARRAY SORT DeclarationsIMPLICIT NONE INTEGER IORDER ( * ) INTEGER NDIM INTEGER ARRAY ( * ) Brief_I/OVARIABLE I/O DESCRIPTION -------- --- -------------------------------------------------- IORDER I Order vector to be used to re-order ARRAY. NDIM I Dimension of ARRAY. ARRAY I-O Array to be re-ordered. Detailed_InputIORDER is the order vector to be used to re-order the input array. The first element of IORDER is the index of the first item of the re-ordered array, and so on. Note that the order imposed by REORDI is not the same order that would be imposed by a sorting routine. In general, the order vector will have been created (by one of the ORDER routines) for a related array, as illustrated in the example below. NDIM is the number of elements in the input array. ARRAY on input, is an array containing some number of elements in unspecified order. Detailed_OutputARRAY on output, is the same array, with the elements re-ordered as specified by IORDER. ParametersNone. ExceptionsError free. 1) If NDIM < 2, this routine executes a no-op. This case is not an error. FilesNone. ParticularsREORDI uses a cyclical algorithm to re-order the elements of the array in place. After re-ordering, element IORDER(1) of the input array is the first element of the output array, element IORDER(2) is the input array is the second element of the output array, and so on. The order vector used by REORDI is typically created for a related array by one of the ORDER routines, as shown in the example below. 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) Sort four related arrays (containing the names, masses, integer ID codes, and flags indicating whether they have a ring system, for a group of planets). Example code begins here. PROGRAM REORDI_EX1 IMPLICIT NONE C C Local constants. C INTEGER NDIM PARAMETER ( NDIM = 8 ) INTEGER STRLEN PARAMETER ( STRLEN = 7 ) C C Local variables. C CHARACTER*(STRLEN) NAMES ( NDIM ) DOUBLE PRECISION MASSES ( NDIM ) INTEGER CODES ( NDIM ) INTEGER I INTEGER IORDER ( NDIM ) LOGICAL RINGS ( NDIM ) C C Set the arrays containing the names, masses (given as C ratios to of Solar GM to barycenter GM), integer ID C codes, and flags indicating whether they have a ring C system. C DATA NAMES / . 'MERCURY', 'VENUS', 'EARTH', 'MARS', . 'JUPITER', 'SATURN', 'URANUS', 'NEPTUNE' / DATA MASSES / . 22032.080D0, 324858.599D0, . 398600.436D0, 42828.314D0, . 126712767.881D0, 37940626.068D0, . 5794559.128D0, 6836534.065D0 / DATA CODES / 199, 299, 399, 499, . 599, 699, 799, 899 / DATA RINGS / . .FALSE., .FALSE., .FALSE., .FALSE., . .TRUE., .TRUE., .TRUE., .TRUE. / C C Sort the object arrays by name. C CALL ORDERC ( NAMES, NDIM, IORDER ) CALL REORDC ( IORDER, NDIM, NAMES ) CALL REORDD ( IORDER, NDIM, MASSES ) CALL REORDI ( IORDER, NDIM, CODES ) CALL REORDL ( IORDER, NDIM, RINGS ) C C Output the resulting table. C WRITE(*,'(A)') ' Planet Mass(GMS/GM) ID Code Rings?' WRITE(*,'(A)') '------- ------------- ------- ------' DO I = 1, NDIM WRITE(*,'(A,F15.3,I9,L5)') NAMES(I), MASSES(I), . CODES(I), RINGS(I) END DO END When this program was executed on a Mac/Intel/gfortran/64-bit platform, the output was: Planet Mass(GMS/GM) ID Code Rings? ------- ------------- ------- ------ EARTH 398600.436 399 F JUPITER 126712767.881 599 T MARS 42828.314 499 F MERCURY 22032.080 199 F NEPTUNE 6836534.065 899 T SATURN 37940626.068 699 T URANUS 5794559.128 799 T VENUS 324858.599 299 F RestrictionsNone. Literature_ReferencesNone. Author_and_InstitutionJ. Diaz del Rio (ODC Space) W.L. Taber (JPL) I.M. Underwood (JPL) VersionSPICELIB Version 1.1.0, 27-AUG-2021 (JDR) Added IMPLICIT NONE statement. Edited the header to comply with NAIF standard. Added complete code example. Added entry #1 in $Exceptions section. 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:36:42 2021