Esope source files

   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
  48 :       call HEA010(T,TC,HEAS0,HEAF0,DTCS0,ZMP0,TRAV)
  49 :       K0   = DTCS0 * DT
  50 :       T1   = T  + 0.5D0 * DT
  51 :       TCI1 = TC + 0.5D0 * K0
  51 :       call HEA010(T1,TCI1,HEAS1,HEAF1,DTCS1,ZMP1,TRAV)
  52 :       K1   = DTCS1 * DT
  53 :       T2   = T  + 0.5D0 * DT
  54 :       TCI2 = TC + 0.5D0 * K1
  54 :       call HEA010(T2,TCI2,HEAS2,HEAF2,DTCS2,ZMP2,TRAV)
  55 :       K2   = DTCS2 * DT
  56 :       T3   = T  + DT
  57 :       TCI3 = TC + K2
  57 :       call HEA010(T3,TCI3,HEAS3,HEAF3,DTCS3,ZMP3,TRAV)
  58 :       K3   = DTCS3 * DT
  59 : C
  60 :       TC1 = TC + 1.D0/6.D0 * (K0 + 2.D0*K1 + 2.D0*K2 + K3)
  61 :       HEAS = 1.D0/6.D0 * (HEAS0 + 2.D0*HEAS1 + 2.D0*HEAS2 + HEAS3)
  62 :       HEAF = 1.D0/6.D0 * (HEAF0 + 2.D0*HEAF1 + 2.D0*HEAF2 + HEAF3)
  63 : C
  64 :       RETURN
  65 :       END
  66 : 
  67 : 
  68 :  
  69 :  
  70 :  
  71 :  
  72 :  
  73 :