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.  
  7. -INC PPARAM
  8. -INC CCOPTIO
  9. *
  10. REAL*8 TEMP0(0:100),POND(0:100),TR1(0:4)
  11. REAL*8 GAMMA(0:4),TABMOD(0:100,0:10)
  12. REAL*8 A(0:4,0:4),B(0:4)
  13. *
  14. * calcul des coefficients de ponderation
  15. *
  16. NDIM=5
  17. TMP=2.*(TEMP0(NMAX)-TEMP0(0))/FLOAT(NMAX)
  18. DO 10 I=1,NMAX-1
  19. POND(I)=(TEMP0(I+1)-TEMP0(I-1))/TMP
  20. 10 CONTINUE
  21. POND(0)=2.*(TEMP0(1)-TEMP0(0))/TMP
  22. POND(NMAX)=2.*(TEMP0(NMAX)-TEMP0(NMAX-1))/TMP
  23. *
  24. * calcul des temps de relaxation
  25. *
  26. TR1(0)=0.
  27. TR1(1)=1.
  28. TR1(2)=1./10.
  29. TR1(3)=1./100.
  30. TR1(4)=1./1000.
  31. *
  32. * calcul des modules de Maxwell
  33. *
  34. DO 20 I=0,4
  35. DO 21 J=0,4
  36. A(I,J)=.0
  37. 21 CONTINUE
  38. B(I)=.0
  39. 20 CONTINUE
  40. *
  41. *
  42. DO 30 I=0,4
  43. DO 31 J=0,4
  44. DO 32 K=0,NMAX
  45. A(I,J)=A(I,J)+POND(K)*EXP(-TEMP0(K)*TR1(I))
  46. 1 *EXP(-TEMP0(K)*TR1(J))
  47. 32 CONTINUE
  48. 31 CONTINUE
  49. DO 33 K=0,NMAX
  50. B(I)=B(I)+POND(K)*EXP(-TEMP0(K)*TR1(I))
  51. 2 *TABMOD(K,INDICE)
  52. 33 CONTINUE
  53. 30 CONTINUE
  54. *
  55. CALL SODIRI(A,GAMMA,B,5,NDIM,KERRE)
  56. IF(KERRE.EQ.49) CALL ERREUR(49)
  57. *
  58. END
  59.  
  60.  
  61.  
  62.  
  63.  
  64.  

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