Esope source files

   1 : C FLURE2    SOURCE    CHAT      05/01/13    00:04:43     5004           
   2 :        SUBROUTINE FLURE2(H0,RHOS,E1AN,T0,TABLE,IDIM,REL)
   3 : *
   4 : *      inversion de la fonction de fluage de LCPC - Granger
   5 : *
   6 : *      entree 
   7 : *            H0,RHOS,E1AN,T0 : parametres de la fonction 
   8 : *                              de fluage de LCPC-Granger
   9 : *            TABLE : tableau des temps de calcul
  10 : *            IDIM : dimension de TABLE
  11 : *            
  12 : *      sortie
  13 : *            REL : tableau des valeurs de la fonction de relaxation 
  14 : *                  calculee aux temps de TABLE
  15 : *  
  16 :       IMPLICIT INTEGER(I-N)
  17 :        IMPLICIT REAL*8(A-H,O-Z) 
  18 : *
  19 :        REAL*8 TABLE(0:1000),REL(0:1000)
  20 : *
  21 : *      initialisation de REL a 0
  22 : *
  23 :        DO 5 I=0,1000
  24 :          REL(I)=0.
  25 : 5      CONTINUE
  26 : *
  27 : *      calcul de REL(0)=rel(t0,t0)
  28 : *
  28 :        CALL LCPC(H0,RHOS,E1AN,T0,T0,FLU)
  29 :        REL(0)=1./FLU
  30 :        T1=TABLE(1)
  31 : *
  32 : *      calcul de REL(K)=rel(t0,table(k))
  33 : *
  34 :        DO 10 K=1,IDIM
  35 :          TK=TABLE(K)
  36 :          TKM=TABLE(K-1)
  36 :          CALL LCPC(H0,RHOS,E1AN,T0,TK,FLU0)
  36 :          CALL LCPC(H0,RHOS,E1AN,T1,TK,FLU1)
  37 :          SOM=2.+(FLU1-FLU0)*REL(0)
  38 :          IF (K.GE.2) THEN
  39 :            DO 20 I=1,K-1
  40 :              TIM=TABLE(I-1)
  41 :              TIP=TABLE(I+1)
  41 :              CALL LCPC(H0,RHOS,E1AN,TIP,TK,FLUP)
  41 :              CALL LCPC(H0,RHOS,E1AN,TIM,TK,FLUM)
  42 :              SOM=SOM+(FLUP-FLUM)*REL(I)
  43 : 20         CONTINUE
  44 :          ENDIF
  44 :          CALL LCPC(H0,RHOS,E1AN,TK,TK,FLUK)
  44 :          CALL LCPC(H0,RHOS,E1AN,TKM,TK,FLUKM)
  45 :          REL(K)=SOM/(FLUK+FLUKM)
  46 : 10     CONTINUE       
  47 :        END
  48 :  
  49 :