Télécharger la plaquette Cast3M


fr en

Télécharger hea011.eso

Retour à la liste

Numérotation des lignes : 

hea011

  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