Télécharger darcytra.notice

Retour à la liste thématique

Afficher cette notice en

Numérotation des lignes :
   1 : $$$$ DARCYTRA NOTICE  CHAT      11/09/12    21:15:43     7124           
   2 :                                              DATE     11/09/12
   3 :                                              
   4 :     Procedure DARCYTRA                       Voir aussi : PECHE
   5 :     ------------------                                    DARCYSAT  
   6 :     DARCYTRA TAB1 ;
   7 : 
   8 :              TAB1.'SOUSTYPE'.'MODELE'.'DOMAINE'.
   9 :                   'CARACTERISTIQUES'.'EMMAGASINEMENT'.'CONVECTION'.
  10 :                   'TEMPS'.'TRACE_CHARGE'.'CHARGE'.'FLUX'.
  11 :                   'BLOCAGE'.'TRACE_IMPOSE'.'FLUX_IMPOSE'.'SOURCE'.
  12 :                   'TEMPS_CALCULES'.'TEMPS_SAUVES'.
  13 :                   'THETA'.'THETA_CONVECTION'
  14 : 
  15 :          ou  TAB1.'SOUSTYPE'.'MODELE'.'DOMAINE'.'ORIENTATION'.
  16 :                   'CARACTERISTIQUES'.'POROSITE'.'DECROISSANCE'.
  17 :                   'COEF_RETARD'.'LANGMUIR'.'FREUNDLICH'.
  18 :                   'LIMITE_SOLUBILITE'.'COEF_DISSOLUTION'.
  19 :                   'CONVECTION'.'TEMPS'.'TRACE_CONC'.'CONCENTRATION'.
  20 :                   'FLUX'.'PRECIPITE'.'DISSOLUTION'.
  21 :                   'BLOCAGE'.'TRACE_IMPOSE'.'FLUX_IMPOSE'.
  22 :                   'DISSOLUTION_IMPOSEE'.'SOURCE'.'TEMPS_CALCULES'.
  23 :                   'TEMPS_SAUVES'.'THETA_DIFF'.'THETA_CONVECTION'.
  24 :                   'THETA_DEC'.'THETA_DISS'.'PENALISATION'.
  25 :                   'EPSI_LIM'.'ITMAX_LIM'.'EPSI_RET'.'EPSI_COR'.
  26 :                   'ITMAX_RET'
  27 : 
  28 :                                     
  29 : 
  30 :     Objet :
  31 :     _______
  32 : 
  33 :     Cette procedure a deux fonctions.
  34 : 
  35 :  1) En presence de l'indice 'CHARGE', on resoud les equations
  36 : de DARCY en transitoire pour l'ecoulement par une methode d'elements
  37 : finis mixtes hybrides (EFMH).
  38 :     Les inconnues du probleme sont la charge ('H'),
  39 : la trace de charge ('TH') et le flux diffusif ('FLUX').
  40 : 
  41 :  2) En l'absence de l'indice 'CHARGE', resoud l'equation
  42 : de transport par diffusion-convection d'un champ scalaire actif
  43 : par un fluide dont la vitesse est connue. L'espece peut se trouver
  44 : sous trois formes : solute, adsorbat et precipite, dont les lois
  45 : d'echange doivent etre specifiees. A chaque loi correspondent un ou
  46 : plusieurs algorithmes auxquels des parametres numeriques doivent etre
  47 : fournis. On utilise la modelisation Darcy EFMH.
  48 :     Les inconnues du probleme sont la concentration ('H'),
  49 : la trace de concentration ('TH') et le flux diffusif ('FLUX').
  50 : 
  51 : 
  52 :     Commentaire :
  53 :     _____________
  54 : 
  55 :     En entree, TAB1 sert a definir les options et les parametres du
  56 : calcul. Les indices de la table TAB1 sont des mots (a coder tel quel)
  57 : dont voici la description :
  58 : 
  59 :   ___________________________________________________________________
  60 :   |                                                                 |
  61 :   | Indice      Contenu                                             |
  62 :   |                                                                 |
  63 :   -------------------------------------------------------------------
  64 :   |                                                                 |
  65 :   |------------------------------------------------                 |
  66 :   |Donnees physiques, geometriques et materielles :                 |
  67 :   |------------------------------------------------                 |
  68 :   |                                                                 |
  69 :   |   ------ Indices communs a l'ecoulement et au transport ------  |
  70 :   |   ------------------------------------------------------------  |
  71 :   |                                                                 |
  72 :   |'SOUSTYPE'   'DARCY' (type MOT)                                  |
  73 :   |                                                                 |
  74 :   |'MODELE'      Objet modele (MMODEL cree par MODE)                |
  75 :   |                                                                 |
  76 :   |'DOMAINE'     References geometriques (TABLE creee par DOMA)     |
  77 :   |                                                                 |
  78 :   |   ------ 1ere possibilite : Resolution de l'ecoulement ------   |
  79 :   |   -----------------------------------------------------------   |
  80 :   |                                                                 |
  81 :   |'CARACTERISTIQUES' Donnees physiques et materielles :            |
  82 :   |              conductivite hydraulique (CHAMELEM cree par MATE)  |
  83 :   |                                                                 |
  84 :   |'EMMAGASINEMENT' Valeur du coefficient d'emmagasinement          |
  85 :   |              (Type CHPO Centre, Comp 'CK', ou FLOTTANT)         |
  86 :   |              - Defaut 1.                                        |
  87 :   |                                                                 |
  88 :   |   ------ 2eme possibilite : Resolution du transport     ------  |
  89 :   |   ------------------------------------------------------------  |
  90 :   |                                                                 |
  91 :   |'CARACTERISTIQUES' Donnees physiques et materielles :            |
  92 :   |              diffusivite effective (CHAMELEM cree par MATE)     |
  93 :   |                                                                 |
  94 :   |'POROSITE'    Valeur de la porosite (Type CHPO Centre, Comp      |
  95 :   |              'CK', ou FLOTTANT) - Defaut 1.                     |
  96 :   |                                                                 |
  97 :   |'DECROISSANCE' Valeur du terme de decroissance (Type FLOTTANT)   |
  98 :   |              Tel que dC/dt = - Lambda * C     - Defaut 0.       |
  99 :   |                                                                 |
 100 :   |'COEF_RETARD' Coefficient de retard lineaire dans le cas simple, |
 101 :   |              ou  Pente a l'origine de la fonction F(C) dans le  |
 102 :   |                  cas d'isotherme non lineaire de Langmuir       |
 103 :   |              ou  Coefficient K de l'isotherme de Freundlich     |
 104 :   |              (Type CHPO Centre 'SCAL', ou FLOTTANT)             |
 105 :   |                                                                 |
 106 :   |'LANGMUIR'    Quantite maximale adsorbee sur le solide           |
 107 :   |              rapportee a l'unite de volume du fluide et exprimee|
 108 :   |              dans la meme unite que le solute.                  |
 109 :   |              (Type CHPO Centre 'SCAL', ou FLOTTANT).            |
 110 :   |              F = (R-1) C / [1 + ((R-1) C / Fsat)]               |
 111 :   |              Si cet indice et le suivant sont absents,          |
 112 :   |              l'equilibre d'adsorption est lineaire. Cet indice a|
 113 :   |              priorite sur l'indice FREUNDLICH.                  |
 114 :   |                                                                 |
 115 :   |'FREUNDLICH'  Exposant de la loi de Freundlich F = K (C ^ 1/n)   |
 116 :   |              (Type FLOTTANT).                                   |
 117 :   |              Dans ce cas (et si l'indice LANGMUIR n'existe pas),|
 118 :   |              l'indice 'RETARD' contient le coefficient          |
 119 :   |              K ramene a une unite de volume de fluide.          |
 120 :   |              - Non disponible pour l'instant -                  |
 121 :   |                                                                 |
 122 :   |'LIMITE_SOLUBILITE' Limite de solubilite (Type MCHAML, Comp 'H') |
 123 :   |                                                                 |
 124 :   |'COEF_DISSOLUTION' Coef. de dissolution (Type CHPO Centre, Comp  |
 125 :   |              'SCAL'). Tel que dC/dt = Coef * (Csat - C) - Par   |
 126 :   |              defaut, la dissolution est instantanee             |
 127 :   |                                                                 |
 128 :   |'CONVECTION'  Debit integre de la vitesse convective a travers   |
 129 :   |              chaque face des elements (Type CHPO Face, comp.    |
 130 :   |              'FLUX')                                            |
 131 :   |                                                                 |
 132 :   |----------------------                                           |
 133 :   |Conditions initiales :                                           |
 134 :   |----------------------                                           |
 135 :   |                                                                 |
 136 :   |   ------ Indices communs a l'ecoulement et au transport ------  |
 137 :   |   ------------------------------------------------------------  |
 138 :   |                                                                 |
 139 :   |'TEMPS'       TABLE contenant a l'indice 0 la valeur du temps    |
 140 :   |              initial (FLOTTANT)                                 |
 141 :   |                                                                 |
 142 :   |                                                                 |
 143 :   |   ------ 1ere possibilite : Resolution de l'ecoulement ------   |
 144 :   |   -----------------------------------------------------------   |
 145 :   |                                                                 |
 146 :   |'CHARGE'      TABLE contenant a l'indice 0 la charge hydraulique |
 147 :   |              (quantite d'element par unite de volume d'eau)     |
 148 :   |              (Type CHPO Centre, Comp 'H')                       |
 149 :   |                                                                 |
 150 :   |'TRACE_CHARGE' TABLE contenant a l'indice 0 la trace de          |
 151 :   |              charge initiale (CHPO, 'TH')                       |
 152 :   |                                                                 |
 153 :   |'FLUX'        TABLE contenant a l'indice 0 le flux hydraulique   |
 154 :   |              initial integre sur chaque face (Type CHPO Face,   |
 155 :   |              comp. 'FLUX')                                      |
 156 :   |                                                                 |
 157 :   |   ------ 2eme possibilite : Resolution du transport     ------  |
 158 :   |   ------------------------------------------------------------  |
 159 :   |                                                                 |
 160 :   |'CONCENTRATION' TABLE contenant a l'indice 0 la concentration    |
 161 :   |              (quantite d'element par unite de volume d'eau)     |
 162 :   |              (Type CHPO Centre, Comp 'H')                       |
 163 :   |                                                                 |
 164 :   |'TRACE_CONC'  TABLE contenant a l'indice 0 la trace de           |
 165 :   |              concentration initiale (Type CHPO Face, Comp 'TH') |
 166 :   |                                                                 |
 167 :   |'FLUX'        TABLE contenant a l'indice 0 le flux diffusif      |
 168 :   |              initial integre sur chaque face (Type CHPO Face,   |
 169 :   |              comp. 'FLUX')                                      |
 170 :   |                                                                 |
 171 :   |'PRECIPITE'   TABLE contenant a l'indice 0 la quantite initiale  |
 172 :   |              de precipite par unite de volume de milieu         |
 173 :   |              (Type CHPO Centre, Comp 'H')                       |
 174 :   |                                                                 |
 175 :   |'DISSOLUTION' TABLE contenant a l'indice 0 la quantite initiale  |
 176 :   |              pour estimer la dissolution au premier pas de temps|
 177 :   |              (Type CHPO, Comp 'H'), voir plus loin.             |
 178 :   |                                                                 |
 179 :   |--------------------------------------                           |
 180 :   |Conditions aux limites / chargements :                           |
 181 :   |--------------------------------------                           |
 182 :   |                                                                 |
 183 :   |   ------ Indices communs a l'ecoulement et au transport ------  |
 184 :   |   ------------------------------------------------------------  |
 185 :   |                                                                 |
 186 :   |'BLOCAGE'     Contient les matrices de blocage (RIGIDITE)        |
 187 :   |                                                                 |
 188 :   |'TRACE_IMPOSE' Valeurs des traces imposees (charge ou concentra- |
 189 :   |              -tion) (CHARGEMENT 'TH' - Obligatoire si BLOCAGE)  |
 190 :   |                                                                 |
 191 :   |'FLUX_IMPOSE' Valeurs des flux imposes integres par face         |
 192 :   |              (Type CHARGEMENT de CHPO Face, comp. 'FLUX'-       |
 193 :   |              defaut 0.)                                         |
 194 :   |                                                                 |
 195 :   |'SOURCE'      Valeurs du terme source par maille et par unite de |
 196 :   |              temps (ex : puits, filiation)                      |
 197 :   |              Les valeurs a l'indice i sont les valeurs entre    |
 198 :   |              les temps i-1 et i.                                |
 199 :   |              (CHARGEMENT de CHPO Centre, comp 'SOUR'- defaut 0.)|
 200 :   |                                                                 |
 201 :   |   ------ 2eme possibilite : Resolution du transport     ------  |
 202 :   |   ------------------------------------------------------------  |
 203 :   |                                                                 |
 204 :   |'DISSOLUTION_IMPOSEE' Valeurs des dissolutions imposees par unite|
 205 :   |              de temps et par maille. (Type CHARGEMENT de CHPO,  |
 206 :   |              Comp 'H'). Les valeurs a l'indice i sont les       |
 207 :   |              valeurs moyennes de dissolution par unite de temps |
 208 :   |              entre les temps i-1 et i.                          |
 209 :   |              Priorite de la dissolution imposee sur les         |
 210 :   |              cinetiques.                                        |
 211 :   |                                                                 |
 212 :   |--------------------                                             |
 213 :   |Donnees numeriques :                                             |
 214 :   |--------------------                                             |
 215 :   |                                                                 |
 216 :   |   ------ Indices communs a l'ecoulement et au transport ------  |
 217 :   |   ------------------------------------------------------------  |
 218 :   |                                                                 |
 219 :   |'TEMPS_CALCULES' Valeur des temps calcules (LISTREEL)            |
 220 :   |              Contient obligatoirement le temps final.           |
 221 :   |                                                                 |
 222 :   |'TEMPS_SAUVES' Valeur des temps sauvegardes (LISTREEL - defaut : |
 223 :   |              on sauve tous les pas de temps)                    |
 224 :   |                                                                 |
 225 :   |   ------ 1ere possibilite : Resolution de l'ecoulement ------   |
 226 :   |   -----------------------------------------------------------   |
 227 :   |                                                                 |
 228 :   |'THETA'       Coefficient de relaxation compris entre 0. et 1.   |
 229 :   |              (theta-methode diffusion) (FLOTTANT - defaut 1.)   |
 230 :   |              Possibilite de non-convergence lorsque theta<1/2   |
 231 :   |              Valeurs de theta generalement utilisees :          |
 232 :   |                   Schema de Euler explicite : 0.                |
 233 :   |                   Schema de Crank-Nicholson : 1/2               |
 234 :   |                   Schema de Galerkin        : 2/3               |
 235 :   |                   Schema de Euler implicite : 1.                |
 236 :   |                                                                 |
 237 :   |   ------ 2eme possibilite : Resolution du transport     ------  |
 238 :   |   ------------------------------------------------------------  |
 239 :   |                                                                 |
 240 :   |'THETA_DIFF'  Coefficient de relaxation compris entre 0. et 1.   |
 241 :   |              (theta-methode diffusion) ('FLOTTANT' - defaut 1.) |
 242 :   |                                                                 |
 243 :   |'THETA_CONVECTION' Idem pour la convection                       |
 244 :   |              ('FLOTTANT', Defaut = THETA_DIFF)                  |
 245 :   |                                                                 |
 246 :   |'THETA_DEC'   Idem mais pour la decroissance                     |
 247 :   |              ('FLOTTANT' - defaut 1/2)                          |
 248 :   |                                                                 |
 249 :   |'THETA_DISS'   Idem mais pour la dissolution                     |
 250 :   |              ('FLOTTANT' - defaut 1.)                           |
 251 :   |                                                                 |
 252 :   |'PENALISATION' Coefficient de penalisation pour la prise en      |
 253 :   |              compte de la limite de solubilite. La presence de  |
 254 :   |              cet indice ou du suivant indique quel schema a ete |
 255 :   |              choisi.                                            |
 256 :   |              (Type 'FLOTTANT') - Valeur conseillee 1.D7         |
 257 :   |                                                                 |
 258 :   |'EPSI_LIM'    Precision relative d'arret pour le shema limite de |
 259 :   |              solubilite predicteur-correcteur iteratif          |
 260 :   |              (Type FLOTTANT) - Valeur conseillee 5.D-3          |
 261 :   |                                                                 |
 262 :   |'ITMAX_LIM'   Nombre maxi d'iterations correspondant aux modules |
 263 :   |              de dissolution avant d'abandonner                  |
 264 :   |              (Type 'ENTIER') - Defaut 50                        |
 265 :   |                                                                 |
 266 :   |'EPSI_RET'    Precision relative d'arret pour la resolution      |
 267 :   |              iterative (Picard) de l'adsorption non lineaire    |
 268 :   |              (Type FLOTTANT) - Defaut 1.D-4                     |
 269 :   |                                                                 |
 270 :   |'EPSI_COR'    Petit saut de concentration pour calculer le coef. |
 271 :   |              de retard par la methode de la corde lorsque le    |
 272 :   |              retard est non-lineaire.                           |
 273 :   |              (Type FLOTTANT) - Defaut 1.D-4                     |
 274 :   |                                                                 |
 275 :   |'ITMAX_RET'   Nombre maxi d'iterations correspondant au retard   |
 276 :   |              non lineaire avant d'abandonner.                   |
 277 :   |              (Type 'ENTIER') - Defaut 20                        |
 278 :   |_________________________________________________________________|
 279 : 
 280 : 
 281 : 
 282 : ---------------------------------
 283 :  Resultats (stockes dans TABLE) :
 284 : ---------------------------------
 285 : 
 286 :   ___________________________________________________________________
 287 :   |                                                                 |
 288 :   | Indice      Contenu                                             |
 289 :   |                                                                 |
 290 :   -------------------------------------------------------------------
 291 :   |                                                                 |
 292 :   |   ------ Indices communs a l'ecoulement et au transport ------  |
 293 :   |   ------------------------------------------------------------  |
 294 :   |                                                                 |
 295 :   |'TEMPS'       TABLE contenant les temps sauvegardes (FLOTTANT)   |
 296 :   |                                                                 |
 297 :   |   ------ 1ere possibilite : Resolution de l'ecoulement  ------  |
 298 :   |   ------------------------------------------------------------  |
 299 :   |                                                                 |
 300 :   |'CHARGE'      TABLE contenant les charges                        |
 301 :   |              (Type CHPO Centre, Comp 'H')                       |
 302 :   |                                                                 |
 303 :   |'TRACE_CHARGE' TABLE contenant les traces de charge              |
 304 :   |              (Type CHPO Face, Comp 'TH')                        |
 305 :   |                                                                 |
 306 :   |'FLUX'        TABLE contenant les debits hydrauliques integres   |
 307 :   |              par face (Type CHPO Face, comp. 'FLUX')            |
 308 :   |                                                                 |
 309 :   |   ------ 2eme possibilite : Resolution du transport     ------  |
 310 :   |   ------------------------------------------------------------  |
 311 :   |                                                                 |
 312 :   |'CONCENTRATION' TABLE contenant les concentrations               |
 313 :   |              (Type CHPO Centre, Comp 'H')                       |
 314 :   |                                                                 |
 315 :   |'TRACE_CONC'  TABLE contenant les traces de concentration        |
 316 :   |              (Type CHPO Face, Comp 'TH')                        |
 317 :   |                                                                 |
 318 :   |'FLUX'        TABLE contenant les debits diffusifs integres      |
 319 :   |              par face (Type CHPO Face, comp. 'FLUX')            |
 320 :   |                                                                 |
 321 :   |'PRECIPITE'   TABLE contenant la quantite de precipite par maille|
 322 :   |              (Type CHPO Centre, Comp 'H')                       |
 323 :   |                                                                 |
 324 :   |'DISSOLUTION' TABLE contenant la quantite de precipite dissoute  |
 325 :   |              entre deux pas de temps par unite de volume et par |
 326 :   |              unite de temps. La valeur stockee a l'indice i,    |
 327 :   |              est valable entre les temps i-1 et i               |
 328 :   |              (Type CHPO, Comp 'H').                             |
 329 :   |              ATTENTION, les valeurs de cette table resultat     |
 330 :   |              n'ont aucun sens lorsque les temps sauvegardes ne  |
 331 :   |              sont pas les memes que les temps calcules. Toute   |
 332 :   |              tentative d'exploitation donnera alors des         |
 333 :   |              resultats incoherents (erreurs de bilan)           |
 334 :   |                                                                 |
 335 :   |'RETARD'      Si cet indice a ete prealablement defini comme une |
 336 :   |              TABLE, alors il contient les valeurs du coefficient|
 337 :   |              de retard (Type 'CHPO' centre, Comp 'SCAL'). Sinon,|
 338 :   |              les valeurs du coefficient de retard ne sont pas   |
 339 :   |              sauvegardees.                                      |
 340 :   |_________________________________________________________________|
 341 : 
 342 :   
 343 :   ___________________________________________________________________
 344 :   |                                                                 |
 345 :   | Indice       Contenu                                            |
 346 :   |                                                                 |
 347 :   -------------------------------------------------------------------
 348 :   |                                                                 |
 349 :   |'TEMPS'       TABLE contenant les temps sauvegardes (FLOTTANT)   |
 350 :   |                                                                 |
 351 :   |'TRACE_CONC'  TABLE contenant les traces de concentration        |
 352 :   |              (Type 'CHPO' Centre, Comp 'TH')                    |
 353 :   |                                                                 |
 354 :   |'CONCENTRATION' TABLE contenant les concentrations               |
 355 :   |              (Type 'CHPO' Centre, Comp 'H')                     |
 356 :   |                                                                 |
 357 :   |'FLUX'        TABLE contenant les debits diffusifs integres par  |
 358 :   |              face (Type 'CHPO' Face, Comp 'FLUX')               |
 359 :   |                                                                 |
 360 :   |'PRECIPITE'   TABLE contenant la quantite de precipite par maille|
 361 :   |              (Type 'CHPO' Centre, Comp 'H')                     |
 362 :   |                                                                 |
 363 :   |'DISSOLUTION' TABLE contenant la quantite de precipite dissoute  |
 364 :   |              entre deux pas de temps par unite de volume et par |
 365 :   |              unite de temps. La valeur stockee a l'indice i,    |
 366 :   |              est valable entre les temps i-1 et i               |
 367 :   |              (Type 'CHPO' Centre, Comp 'H')                     |
 368 :   |              ATTENTION, les valeurs de cette table resultat     |
 369 :   |              n'ont aucun sens lorsque les temps sauvegardes ne  |
 370 :   |              sont pas les memes que les temps calcules. Toute   |
 371 :   |              tentative d'exploitation donnera alors des         |
 372 :   |              resultats incoherents (erreurs de bilan)           |
 373 :   |                                                                 |
 374 :   |'RETARD'      Si cet indice a ete prealablement defini comme une |
 375 :   |              TABLE, alors il contient les valeurs du coefficient|
 376 :   |              de retard (Type 'CHPO' centre, Comp 'SCAL'). Sinon,|
 377 :   |              les valeurs du coefficient de retard ne sont pas   |
 378 :   |              sauvegardees.                                      |
 379 :   |_________________________________________________________________|
 380 : 
 381 :   ___________________________________________________________________
 382 :   |                                                                 |
 383 :   | Les tables resultats sont indicees par des entiers variant de 0 |
 384 :   | a N .                                                           |
 385 :   | A l'indice 0 on stocke les valeurs initiales, aux indices       |
 386 :   | suivants les champs correspondant au temps de sortie TEMPS.I .  |
 387 :   | Les champs servant en cas de reprise sont ceux correpondant au  |
 388 :   | dernier indice.                                                 |
 389 :   |_________________________________________________________________|
 390 : 
 391 :        Exemple : pour lister le CHPOINT de charge calcule pour la
 392 :                  valeur du parametre d'evolution 2.5, il faudra coder :
 393 :                            LIST ( PECHE TAB1 CHARGE 2.5 ) ;
 394 :                  ou
 395 : 
 396 :                  si on sait que l'indice i de la table TEMPS contient
 397 :                  la valeur du parametre d'evolution 2.5, on peut coder
 398 :                            LIST ( TAB1 . CHARGE . i ) ;
 399 : 
 400 : 
 401 :     Remarques :
 402 :     ___________
 403 : 
 404 :     Les resultats etant reperes dans TAB1 il n'y a pas d'objets nommes
 405 : crees par cette procedure.
 406 : 
 407 :     Une fois sortie de DARCYTRA on peut y re-entrer en definissant de
 408 : nouvelles valeurs pour les LISTREELS TEMPS_CALCULES et TEMPS_SAUVES
 409 : et en invoquant de nouveau DARCYTRA avec les memes operandes que
 410 : lors du premier appel.
 411 : 
 412 :     Les FLUX imposes et les flux
 413 : calcules sont les flux diffusifs. Les flux convectifs s'obtiennent en
 414 : post-traitement en faisant le produit du flux de vitesse (stocke a
 415 : l'indice CONVECTION) par la trace de charge (stockee a l'indice
 416 : TRACE_CHARGE).
 417 : 

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