Esope source files

   1 : C LCPCDI    SOURCE    CHAT      05/01/13    01:13:30     5004           
   2 :        SUBROUTINE LCPCDI(H0,RHOS,E1AN,TMAX,MAILLE,TEMPS
   3 :      1     ,TREL,EVO,YOUNG,JDIM)       
   4 : 
   5 : *
   6 :       IMPLICIT INTEGER(I-N)
   7 :        IMPLICIT REAL*8(A-H,O-Z)      
   8 : -INC CCOPTIO
   9 : *
  10 :        REAL*8 TEMPS(0:1000),REL(0:1000),YOUNG(0:1000)
  11 :        REAL*8 TEMP0(0:100),TR1(0:4),GAMMA(0:4),TR(0:8)   
  12 :        REAL*8 MODULE(0:10),TABMOD(0:100,0:10)
  13 :        REAL*8 TABGAM(0:10,0:4),TREL(0:10)
  14 :        REAL*8 EVO(0:1000,0:10)              
  15 : *
  16 :        TINI=1.
  17 :        TFIN=TMAX      
  18 :        DO 10 I=0,100
  19 :          DO 11 J=0,10
  20 :             TABMOD(I,J)=0.
  21 : 11       CONTINUE
  22 : 10     CONTINUE
  23 :        DO 12 I=0,10
  24 :          DO 13 J=0,4
  25 :             TABGAM(I,J)=0.
  26 : 13       CONTINUE
  27 : 12     CONTINUE
  28 : *
  29 : *     liste des t0 a calculer
  30 : *
  31 :       TEMP0(0)=TINI
  32 :       TMP=LOG((TFIN-TINI)/.25)/.5
  33 :       NMAX=INT(TMP)
  34 :       IF (NMAX.GT.100) THEN
  35 :          NMAX=100
  36 :       ENDIF
  37 :       DO 20 I=1,NMAX
  38 :         TEMP0(I)=.25*EXP(.5*FLOAT(I))+TINI 
  39 : 20    CONTINUE
  40 : *
  41 : *     boucle sur les t0
  42 : *
  43 :       DO 30 L=0,NMAX
  44 :         T0=TEMP0(L)
  45 :         TINI=TEMP0(L)
  45 :         CALL TABTP2(H0,RHOS,E1AN,T0,TINI,TFIN,TEMPS,JDIM)
  45 :         CALL FLURE2(H0,RHOS,E1AN,T0,TEMPS,JDIM,REL) 
  46 : *
  46 :         CALL CALMOD(TEMPS,JDIM,REL,MAILLE,MODULE,TR)
  47 :         IF(IERR.NE.0) RETURN
  48 : *
  49 :         DO 31 I=0,MAILLE
  50 :           TABMOD(L,I)=MODULE(I)
  51 :           TREL(I)=TR(I)       
  52 : 31        CONTINUE  
  53 : 30      CONTINUE 
  54 : *
  55 :       DO 40 INDICE=0,MAILLE
  55 :         CALL CALEVO(TEMP0,NMAX,TABMOD,INDICE,GAMMA,TR1)
  56 :         IF(IERR.NE.0) RETURN
  57 :         DO 41 J=0,4
  58 :           TABGAM(INDICE,J)=GAMMA(J)
  59 : 41      CONTINUE
  60 : 40    CONTINUE
  61 : *    
  62 : *     on recupere la liste des temps dans TEMPS
  63 : * 
  64 :       T0=TEMP0(0)
  65 :       TINI=TEMP0(0)
  65 :       CALL TABTP2(H0,RHOS,E1AN,T0,TINI,TFIN,TEMPS,JDIM)
  66 : *     
  67 : *     on recupere la  liste des gamma(t)
  68 : *
  69 :       DO 50 I=0,MAILLE
  70 :         EVO(0,I)=0.
  71 :         DO 51 J=1,JDIM
  72 :           EVO(J,I)=TABGAM(I,0)
  73 :      1           +TABGAM(I,1)*DEXP(-TR1(1)*TEMPS(J))
  74 :      2           +TABGAM(I,2)*DEXP(-TR1(2)*TEMPS(J))
  75 :      3           +TABGAM(I,3)*DEXP(-TR1(3)*TEMPS(J))
  76 :      4           +TABGAM(I,4)*DEXP(-TR1(4)*TEMPS(J))
  77 : 
  78 : 51      CONTINUE
  79 : 50    CONTINUE
  80 : *
  81 : *     on recupere le module d'young
  82 : *
  83 :       YOUNG(0)=0.
  84 :       DO 60 J=1,JDIM
  85 :         YOUNG(J)=0.
  86 :         DO 61 I=0,MAILLE
  87 :           YOUNG(J)=YOUNG(J)+EVO(J,I)
  88 : 61      CONTINUE
  89 : 60    CONTINUE
  90 : 
  91 :     
  92 :       END
  93 :  
  94 :  
  95 :  
  96 :  
  97 :