Esope source files

   1 : C CALEVO    SOURCE    CHAT      05/01/12    21:45:52     5004           
   2 :        SUBROUTINE CALEVO(TEMP0,NMAX,TABMOD,INDICE,GAMMA,TR1)
   3 : *
   4 :       IMPLICIT INTEGER(I-N)
   5 :        IMPLICIT REAL*8(A-H,O-Z)      
   6 : -INC CCOPTIO
   7 : *
   8 :        REAL*8 TEMP0(0:100),POND(0:100),TR1(0:4)
   9 :        REAL*8 GAMMA(0:4),TABMOD(0:100,0:10) 
  10 :        REAL*8 A(0:4,0:4),B(0:4)
  11 : *
  12 : *      calcul des coefficients de ponderation
  13 : *
  14 :        NDIM=5
  15 :        TMP=2.*(TEMP0(NMAX)-TEMP0(0))/FLOAT(NMAX)
  16 :        DO 10 I=1,NMAX-1
  17 :          POND(I)=(TEMP0(I+1)-TEMP0(I-1))/TMP
  18 : 10     CONTINUE
  19 :        POND(0)=2.*(TEMP0(1)-TEMP0(0))/TMP
  20 :        POND(NMAX)=2.*(TEMP0(NMAX)-TEMP0(NMAX-1))/TMP
  21 : *
  22 : *      calcul des temps de relaxation
  23 : *
  24 :        TR1(0)=0.
  25 :        TR1(1)=1.    
  26 :        TR1(2)=1./10.
  27 :        TR1(3)=1./100.
  28 :        TR1(4)=1./1000.
  29 : *
  30 : *      calcul des modules de Maxwell
  31 : *
  32 :         DO 20 I=0,4
  33 :           DO 21 J=0,4
  34 :             A(I,J)=.0
  35 : 21        CONTINUE
  36 :           B(I)=.0
  37 : 20     CONTINUE
  38 : *
  39 : *
  40 :        DO 30 I=0,4
  41 :          DO 31 J=0,4
  42 :            DO 32 K=0,NMAX
  43 :              A(I,J)=A(I,J)+POND(K)*DEXP(-TEMP0(K)*TR1(I))
  44 :      1                *DEXP(-TEMP0(K)*TR1(J))
  45 : 32         CONTINUE
  46 : 31       CONTINUE
  47 :          DO 33 K=0,NMAX 
  48 :            B(I)=B(I)+POND(K)*DEXP(-TEMP0(K)*TR1(I))
  49 :      2              *TABMOD(K,INDICE)
  50 : 33       CONTINUE
  51 : 30     CONTINUE
  52 : *
  52 :        CALL SODIRI(A,GAMMA,B,5,NDIM,KERRE)
  52 :        IF(KERRE.EQ.49) CALL ERREUR(49)
  53 : *
  54 :        END
  55 :  
  56 :  
  57 :