Numérotation des lignes :

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)CC************************************************************************      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.D0C?    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)=SJC?    AM(K,J,N)=SJ*PORO(N,NFP) 2    CONTINUE       ELSE       DO 12 J=1,NP      SJ=0.D0C?    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)=SJC?    AM(K,J,N)=SJ*PORO(N,NFP) 12   CONTINUE       ENDIF 20   CONTINUE CC      RETURN 1001 FORMAT(20(1X,I5)) 1002 FORMAT(10(1X,1PE11.4))      END  

© Cast3M 2003 - Tous droits réservés.
Mentions légales