Sources Esope | Cast3M

kaam
C KAAM      SOURCE    CHAT      05/01/13    00:50:19     5004
      SUBROUTINE KAAM(FN,HR,PG,PGSQ,AM,K,NEL,NP,NPG,IES,IAXI,RPG)
C************************************************************************
C   CALCUL DE LA MATRICE M ( P*DIV(U) )  --> AM(NEL,NP,IES)
C
C************************************************************************
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
 
C***
      DIMENSION FN(NP,NPG),HR(IES,NP,NPG),PG(NPG),PGSQ(NPG),RPG(NPG)
      DIMENSION AM(NEL,NP,IES)
C?    DIMENSION PORO(NPR,*)
C?    DIMENSION LE(NP),IPADL(*)
C***
 
      IF(IAXI.EQ.0)CALL INITD(RPG,NPG,1.D30)
 
      DO 20 N=1,IES
      IAX=0
      IF((3-IAXI).EQ.N)IAX=1
      IF(IES.EQ.3)IAX=0
 
      IF(IAX.EQ.0)THEN
 
      DO 2 J=1,NP
      SJ=0.D0
C?    NF=IPADL(LE(J))
C?    NFP=IPOR*(NF-1)+1
      DO 1 L=1,NPG
 1    SJ=SJ+HR(N,J,L)*PGSQ(L)
      AM(K,J,N)=SJ
C?    AM(K,J,N)=SJ*PORO(N,NFP)
 2    CONTINUE
 
      ELSE
 
      DO 12 J=1,NP
      SJ=0.D0
C?    NF=IPADL(LE(J))
C?    NFP=IPOR*(NF-1)+1
      DO 11 L=1,NPG
 11   SJ=SJ+(HR(N,J,L)+FN(J,L)/RPG(L))*PGSQ(L)
      AM(K,J,N)=SJ
C?    AM(K,J,N)=SJ*PORO(N,NFP)
 12   CONTINUE
 
      ENDIF
 20   CONTINUE
 
C
C
      RETURN
 1001 FORMAT(20(1X,I5))
 1002 FORMAT(10(1X,1PE11.4))
      END