C FLUPR1    SOURCE    AM        09/12/08    21:15:18     6582           
      SUBROUTINE FLUPR1(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
     $                  AN,BN,CN,KERRE)
C
C ROUTINE APPELE DANS CAMECA
C_______________________________________________________________________
C
C SOUS PROGRAMME QUI CALCUL LES MATRICES DE DEFORMATION 
C  DE FLUAGE PROPRE 
C
C_______________________________________________________________________
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION VARF(NVARI),CMAT(NMAT)
      DIMENSION AN(NSTRS),BN(NSTRS,NSTRS),CN(NSTRS,NSTRS)
      DIMENSION A1(10),B1(10,10),C1(10,10)
      DIMENSION A2(10),B2(10,10),C2(10,10)
      DIMENSION A3(10),B3(10,10),C3(10,10)
      DIMENSION A4(10),B4(10,10),C4(10,10)
      DIMENSION A5(10),B5(10,10),C5(10,10)
      DIMENSION A6(10),B6(10,10),C6(10,10)
      DIMENSION AA(10),BA(10,10),CA(10,10)
      DIMENSION AD(10),BD(10,10),CD(10,10)

C  TEST SUR LA DIMENSION

      NDI = 10
      IF(NSTRS.GT.NDI) THEN
         PRINT *, ' FLUD - DIMENSION INSUFFISANTE'
         KERRE=99
         RETURN
      ENDIF
C
C CALCUL DES DEFORMATIONS DE FLUAGE
C   POUR LES 6 CHAINES DE KELVIN VOIGT
C
C   Chaine 1
C
      NN=1
      CALL FLUD(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
     $                       A1,B1,C1,NN,NDI)
C
C   Chaine 2
C
      NN=2
      CALL FLUD(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
     $                       A2,B2,C2,NN,NDI)
C
C   Chaine 3
C
      NN=3
      CALL FLUD(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
     $                       A3,B3,C3,NN,NDI)
C
C   Chaine 4
C
      NN=4
      CALL FLUD(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
     $                       A4,B4,C4,NN,NDI)
C
C   Chaine 5
C
      NN=5
      CALL FLUD(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
     $                       A5,B5,C5,NN,NDI)
C
C   Chaine 6
C
      NN=6
      CALL FLUD(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
     $                       A6,B6,C6,NN,NDI)
C
C   Amortisseur
C
      CALL FLUA(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
     $                       AA,BA,CA,NDI)
C
C   fluage dessiccation
C
      CALL FLUDE1(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
     $                       AD,BD,CD,NDI)
C
C
C   ASSEMBLAGE DES MATRICES
C
      DO 10 I=1,NSTRS
        AN(I) = A1(I) + A2(I) + A3(I) + A4(I) + A5(I)
     $          + A6(I) + AA(I) + AD(I)
        DO 10 J=1,NSTRS
          BN(I,J) = B1(I,J) + B2(I,J) + B3(I,J)
     $              + B4(I,J) + B5(I,J) + B6(I,J)
     $              + BA(I,J) + BD(I,J)
          CN(I,J) = C1(I,J) + C2(I,J) + C3(I,J)
     $              + C4(I,J) + C5(I,J) + C6(I,J)
     $              + CA(I,J) + CD(I,J)
   10   CONTINUE 
C
      END