Esope source files

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