Télécharger flurel.eso

Retour à la liste

Numérotation des lignes :

flurel
  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. *
  28. CALL EURO(HO,ROH,FCM,S,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)
  37. CALL EURO(HO,ROH,FCM,S,T0,TK,FLU0)
  38. CALL EURO(HO,ROH,FCM,S,T1,TK,FLU1)
  39. SOM=2.+(FLU1-FLU0)*REL(0)
  40. IF (K.GE.2) THEN
  41. DO 20 I=1,K-1
  42. TIM=TABLE(I-1)
  43. TIP=TABLE(I+1)
  44. CALL EURO(HO,ROH,FCM,S,TIP,TK,FLUP)
  45. CALL EURO(HO,ROH,FCM,S,TIM,TK,FLUM)
  46. SOM=SOM+(FLUP-FLUM)*REL(I)
  47. 20 CONTINUE
  48. ENDIF
  49. CALL EURO(HO,ROH,FCM,S,TK,TK,FLUK)
  50. CALL EURO(HO,ROH,FCM,S,TKM,TK,FLUKM)
  51. REL(K)=SOM/(FLUK+FLUKM)
  52. 10 CONTINUE
  53. END
  54.  
  55.  
  56.  

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