SUM:  PROC(IN_ARRAY) RETURNS(FLOAT BINARY(53));

/* ROUTINE SUMS A DOUBLE PRECISION FLOATING POINT ARRAY 
   AUTHOR:  AL JAWORSKI
   DATE:    MARCH 5, 1982 
   WESTAT FUNCTION LIBRARY */

DCL (IN_ARRAY(*),SUM_VALUE) FLOAT BINARY(53);
DCL I FIXED BINARY(31);

SUM_VALUE=0;

DO I=LBOUND(IN_ARRAY,1) TO HBOUND(IN_ARRAY,1);
  SUM_VALUE=SUM_VALUE+IN_ARRAY(I);
END;

RETURN(SUM_VALUE);

END SUM;