INTEGER*4 FUNCTION IJSYM(IORD,IR,IC)
C
C To find the location of element IR,IC of a symmetric array of order 
C IORD when the array is stored in a column vector in upper triangular form
C
      IF (IR.LE.IC) THEN
        IJSYM= (IR-1)*IORD+IC - (IR-1)*IR/2
      ELSE
        IJSYM= (IC-1)*IORD+IR - (IC-1)*IC/2
      ENDIF
      RETURN
      END