SUBROUTINE ARR_RANG1D( N, ARRAY, DMIN, DMAX, STATUS ) *+ * Name: * ARR_RANG1D * Purpose: * Find range of values in a DOUBLE PRECISION array * Language: * Starlink Fortran * Invocation: * CALL ARR_RANG1D( N, ARRAY, DMIN, DMAX, STATUS ) * Description: * Obtains the minimum and maximum values in a DOUBLE PRECISION array. * All values are consider good. * Arguments: * N = INTEGER (given) * Number of values in array * ARRAY[] = DOUBLE PRECISION (given) * Array of values whose range is to be found * DMIN = DOUBLE PRECISION (returned) * The minimum value in the array * DMAX = DOUBLE PRECISION (returned) * The maximum value in the array * STATUS = INTEGER (given) * The global status. * Examples: * {routine_example_text} * {routine_example_description} * Pitfalls: * {pitfall_description}... * Notes: * {routine_notes}... * Prior Requirements: * {routine_prior_requirements}... * Side Effects: * {routine_side_effects}... * Algorithm: * {algorithm_description}... * Accuracy: * {routine_accuracy} * Timing: * {routine_timing} * External Routines Used: * {name_of_facility_or_package}: * {routine_used}... * Implementation Deficiencies: * {routine_deficiencies}... * References: * ARR Subroutine Guide : http://www.sr.bham.ac.uk/asterix-docs/Programmer/Guides/arr.html * Keywords: * package:arr, usage:public * Copyright: * Copyright (C) University of Birmingham, 1995 * Authors: * DJA: David J. Allan (Jet-X, University of Birmingham) * {enter_new_authors_here} * History: * 7 Mar 1995 (DJA): * Original version. * {enter_changes_here} * Bugs: * {note_any_bugs_here} *- * Type Definitions: IMPLICIT NONE ! No implicit typing * Global Constants: INCLUDE 'SAE_PAR' ! Standard SAE constants INCLUDE 'PRM_PAR' * Arguments Given: INTEGER N ! Number of values DOUBLE PRECISION ARRAY(*) ! Values * Arguments Returned: DOUBLE PRECISION DMIN ! Minimum value DOUBLE PRECISION DMAX ! Maximum value * Status: INTEGER STATUS ! Global status * Local Variables: INTEGER I ! Loop over array *. * Check inherited global status. IF ( STATUS .NE. SAI__OK ) RETURN * Initialise DMIN = VAL__MAXD DMAX = VAL__MIND * Find range DO I = 1, N IF ( ARRAY(I) .LT. DMIN ) DMIN = ARRAY(I) IF ( ARRAY(I) .GT. DMAX ) DMAX = ARRAY(I) END DO END