Télécharger hea011.eso

Retour à la liste

Numérotation des lignes :

  1. C HEA011 SOURCE SERRE 12/07/10 21:15:19 7431
  2. SUBROUTINE HEA011(TC1,HEAS,HEAF,T,DT,TC,TRAV)
  3. C---------------------------------------------------------------------
  4. C Calcul de la température des plaques par la méthode RK4 et des
  5. C puissances transmises aux plaques et au fluide
  6. C---------------------------------------------------------------------
  7. C
  8. C---------------------------
  9. C Parametres Entree/Sortie :
  10. C---------------------------
  11. C
  12. C /S TC1 : flottant : Température des plaques au temps t+dt (K)
  13. C /S HEAS : flottant : Puissance générée dans le solide (J/s)
  14. C /S HEAF : flottant : Puissance transmise au fluide (J/s)
  15. C E/ T : flottant : Temps au début du pas de temps (s)
  16. C E/ DT : flottant : Pas d'intégration (s)
  17. C E/ TC : flottant : Température des plaques au temps t (K)
  18. C E/ TRAV : entier : Pointeur du segment de travail
  19. C
  20. C------------------------------
  21. C Variables de TRAV utilisées :
  22. C------------------------------
  23. C
  24. C Aucune variable n'est utilisée.
  25. C Le pointeur du segment de travail est simplement transmis.
  26. C
  27. C---------------------------------------------------------------------
  28. C
  29. C Langage : ESOPE + FORTRAN 77
  30. C
  31. C Mise en oeuvre : H. Paillère (1997, TTMF)
  32. C
  33. C---------------------------------------------------------------------
  34. IMPLICIT INTEGER(I-N)
  35. IMPLICIT REAL*8 (A-H,O-Z)
  36. REAL*8 K0,K1,K2,K3
  37. segment trav
  38. integer iKALP,iMODEL
  39. real*8 e,L,Lch,Dh,S,sp,Ck
  40. real*8 mc,Cpc
  41. real*8 g,R,deltah
  42. real*8 M(nbesp),cstmod(ncst)
  43. real*8 Cpi(nbesp),al
  44. real*8 eps_mh2,eps_dt,eps_con,u
  45. real*8 XH2MOY,XO2MOY,XN2MOY,XH2OMOY,PRESSION,TEMPENT
  46. real*8 XHEMOY,XCO2MOY,XCOMOY
  47. endsegment
  48. C
  49. call HEA010(T,TC,HEAS0,HEAF0,DTCS0,ZMP0,TRAV)
  50. K0 = DTCS0 * DT
  51. T1 = T + 0.5D0 * DT
  52. TCI1 = TC + 0.5D0 * K0
  53. call HEA010(T1,TCI1,HEAS1,HEAF1,DTCS1,ZMP1,TRAV)
  54. K1 = DTCS1 * DT
  55. T2 = T + 0.5D0 * DT
  56. TCI2 = TC + 0.5D0 * K1
  57. call HEA010(T2,TCI2,HEAS2,HEAF2,DTCS2,ZMP2,TRAV)
  58. K2 = DTCS2 * DT
  59. T3 = T + DT
  60. TCI3 = TC + K2
  61. call HEA010(T3,TCI3,HEAS3,HEAF3,DTCS3,ZMP3,TRAV)
  62. K3 = DTCS3 * DT
  63. C
  64. TC1 = TC + 1.D0/6.D0 * (K0 + 2.D0*K1 + 2.D0*K2 + K3)
  65. HEAS = 1.D0/6.D0 * (HEAS0 + 2.D0*HEAS1 + 2.D0*HEAS2 + HEAS3)
  66. HEAF = 1.D0/6.D0 * (HEAF0 + 2.D0*HEAF1 + 2.D0*HEAF2 + HEAF3)
  67. C
  68. RETURN
  69. END
  70.  
  71.  
  72.  
  73.  
  74.  
  75.  
  76.  
  77.  
  78.  

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