Télécharger calevo.eso

Retour à la liste

Numérotation des lignes :

  1. C CALEVO SOURCE CB215821 16/04/21 21:15:28 8920
  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)*EXP(-TEMP0(K)*TR1(I))
  44. 1 *EXP(-TEMP0(K)*TR1(J))
  45. 32 CONTINUE
  46. 31 CONTINUE
  47. DO 33 K=0,NMAX
  48. B(I)=B(I)+POND(K)*EXP(-TEMP0(K)*TR1(I))
  49. 2 *TABMOD(K,INDICE)
  50. 33 CONTINUE
  51. 30 CONTINUE
  52. *
  53. CALL SODIRI(A,GAMMA,B,5,NDIM,KERRE)
  54. IF(KERRE.EQ.49) CALL ERREUR(49)
  55. *
  56. END
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  

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