Numérotation des lignes :

calevo
C CALEVO    SOURCE    CB215821  16/04/21    21:15:28     8920       SUBROUTINE CALEVO(TEMP0,NMAX,TABMOD,INDICE,GAMMA,TR1)*      IMPLICIT INTEGER(I-N)       IMPLICIT REAL*8(A-H,O-Z) -INC PPARAM-INC CCOPTIO*       REAL*8 TEMP0(0:100),POND(0:100),TR1(0:4)       REAL*8 GAMMA(0:4),TABMOD(0:100,0:10)       REAL*8 A(0:4,0:4),B(0:4)**      calcul des coefficients de ponderation*       NDIM=5       TMP=2.*(TEMP0(NMAX)-TEMP0(0))/FLOAT(NMAX)       DO 10 I=1,NMAX-1         POND(I)=(TEMP0(I+1)-TEMP0(I-1))/TMP10     CONTINUE       POND(0)=2.*(TEMP0(1)-TEMP0(0))/TMP       POND(NMAX)=2.*(TEMP0(NMAX)-TEMP0(NMAX-1))/TMP**      calcul des temps de relaxation*       TR1(0)=0.       TR1(1)=1.       TR1(2)=1./10.       TR1(3)=1./100.       TR1(4)=1./1000.**      calcul des modules de Maxwell*        DO 20 I=0,4          DO 21 J=0,4            A(I,J)=.021        CONTINUE          B(I)=.020     CONTINUE**       DO 30 I=0,4         DO 31 J=0,4           DO 32 K=0,NMAX             A(I,J)=A(I,J)+POND(K)*EXP(-TEMP0(K)*TR1(I))     1                *EXP(-TEMP0(K)*TR1(J))32         CONTINUE31       CONTINUE         DO 33 K=0,NMAX           B(I)=B(I)+POND(K)*EXP(-TEMP0(K)*TR1(I))     2              *TABMOD(K,INDICE)33       CONTINUE30     CONTINUE*       CALL SODIRI(A,GAMMA,B,5,NDIM,KERRE)       IF(KERRE.EQ.49) CALL ERREUR(49)*       END      

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