Esope source files

   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 :