Télécharger fla011.eso

Retour à la liste

Numérotation des lignes :

  1. C FLA011 SOURCE SERRE 12/07/10 21:15:16 7431
  2. SUBROUTINE FLA011(TCPDT,T,DT,TC,TRAV)
  3. C---------------------------------------------------------------------
  4. C Calcul de la température des plaques par la méthode RK4
  5. C---------------------------------------------------------------------
  6. C
  7. C---------------------------
  8. C Parametres Entree/Sortie :
  9. C---------------------------
  10. C
  11. C /S TCPDT : flottant : Température des plaques au temps t+dt (K)
  12. C E/ T : flottant : Temps au début du pas de temps (s)
  13. C E/ DT : flottant : Pas d'intégration (s)
  14. C E/ TC : flottant : Température des plaques au temps t (K)
  15. C
  16. C------------------------------
  17. C Variables de TRAV utilisées :
  18. C------------------------------
  19. C
  20. C Aucune variable n'est utilisée.
  21. C Le pointeur du segment de travail est simplement transmis.
  22. C
  23. C---------------------------------------------------------------------
  24. C
  25. C Langage : ESOPE + FORTRAN 77
  26. C
  27. C Mise en oeuvre : H. Paillère (1997, TTMF)
  28. C
  29. C---------------------------------------------------------------------
  30. IMPLICIT INTEGER(I-N)
  31. IMPLICIT REAL*8 (A-H,O-Z)
  32. REAL*8 K0,K1,K2,K3
  33. segment trav
  34. integer iKALP,iMODEL
  35. real*8 e,L,Lch,Dh,S,sp,Ck
  36. real*8 mc,Cpc
  37. real*8 g,R,deltah
  38. real*8 M(nbesp),cstmod(ncst)
  39. real*8 Cpi(nbesp),al
  40. real*8 eps_mh2,eps_dt,eps_con,u
  41. real*8 XH2MOY,XO2MOY,XN2MOY,XH2OMOY,PRESSION,TEMPENT
  42. real*8 XHEMOY,XCO2MOY,XCOMOY
  43. endsegment
  44. C
  45. K0 = FLA010(T,TC,TRAV) * DT
  46. T1 = T + 0.5D0 * DT
  47. TC1 = TC + 0.5D0 * K0
  48. K1 = FLA010(T1,TC1,TRAV) * DT
  49. T2 = T + 0.5D0 * DT
  50. TC2 = TC + 0.5D0 * K1
  51. K2 = FLA010(T2,TC2,TRAV) * DT
  52. T3 = T + DT
  53. TC3 = TC + K2
  54. K3 = FLA010(T3,TC3,TRAV) * DT
  55. C
  56. TCPDT = TC + 1.D0/6.D0 * (K0 + 2.D0*K1 + 2.D0*K2 +K3)
  57. C
  58. RETURN
  59. END
  60.  
  61.  
  62.  
  63.  
  64.  
  65.  
  66.  
  67.  
  68.  

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