C         SUBROUTINE TO CALCULATE PARTIAL SUMS  PSUM
C
          SUBROUTINE PRTSUM(MODULE,ISON)
          COMMON/CENTS/C(3,15)
          COMMON/PROD/PAC(3,3,15)
          COMMON/SUMIT/PSUM(3,16,5)
          COMMON/LINKIT/ISTRT,LINKS(6,15)
          IF(MODULE.EQ.0)GO TO 99
C         PICK OUT PARENT
          IMOD=LINKS(1,MODULE)
          I=1
10        I=I+1
          NUMSON=LINKS(I,MODULE)
          IF(NUMSON.NE.ISON)GO TO 10
          NUMSON=I-1
          IF(IMOD.EQ.0)GO TO 12
          I=1
11        I=I+1
          INSON=LINKS(I,IMOD)
          IF(INSON.NE.MODULE)GO TO 11
          INSON=I-1
          GO TO 13
12        INSON=1
          IMOD=16
13        DO 69 J=1,3
          SUM=0.
          DO 68 I=1,3
68        SUM=PAC(J,I,MODULE)*C(I,ISON)+SUM
69        PSUM(J,MODULE,NUMSON)=SUM
          DO 67 I=1,3
67    PSUM(I,MODULE,NUMSON)=PSUM(I,MODULE,NUMSON)+PSUM(I,
     2IMOD,INSON)
99        RETURN
          END
