Télécharger la plaquette Cast3M


fr en

Télécharger fla011.eso

Retour à la liste

Numérotation des lignes : 

fla011

  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