Accueil | Cast3M

   1 : $$$$ HT_PRO   NOTICE  FANDEUR   13/03/15    21:15:00     7712           
   2 :                                              DATE     13/03/15
   3 : 
   4 :   Voir aussi : DARCYSAT
   5 :                                                       KR_PRO  
   6 :     ----------------
   7 : 
   8 :     CHPO2 CHPO3 CHPO4 = HT_PRO  TAB1 (CHPO1) ;
   9 : 
  10 : 
  11 :     Objet :
  12 :     _______
  13 : 
  14 :     Cette procedure calcule la saturation, la teneur en eau et la 
  15 :     capacite en fonction de la pression d'eau P.
  16 :     Cette procedure est utilisee a partir de la procedure DARCYSAT
  17 :     Elle correspond a une loi definie par le soustype de TAB1 :
  18 :     - VAN_GENUCHTEN
  19 :        Saturation reduite dans [0,1] : S = [1 + (-beta.Pw)**n]**(-m)
  20 :     - EXPONENTIELLE
  21 :        S = C / ([exp(-beta.Pw)]**N + C - 1)
  22 :     - LOGARITHMIQUE
  23 :        (-beta.Pw) tronque a 1 quand inferieur a 1.
  24 :        S = C / ([ln(-beta.Pw)]**N + C)
  25 : 
  26 :     avec S = (TH2O - teta_r) / (poros - teta_r)
  27 : 
  28 : 
  29 :     Commentaires :
  30 :     --------------
  31 : 
  32 :     TAB1  : table contenant les caracteristiques physiques, ayant
  33 :             pour indices
  34 :             BHETA, NEXP, MEXP : coefficients de la loi VAN_GENUCHTEN   
  35 :                                 (attention a l'unite de beta)
  36 :             BHETA, COEF_N, COEF_C : coef pour la loi log ou exp.
  37 :             PORO              : porosite (s. d.), poros
  38 :                                 (FLOTTANT ou CHAMP-POINT centre,
  39 :                                  comp 'SCAL')
  40 :             TERESIDU          : teneur en eau residuelle (s. d.), teta_r
  41 :                                 (FLOTTANT ou CHAMP-POINT centre,
  42 :                                  comp 'SCAL')
  43 : 
  44 :     CHPO1 : Pression d'eau (negative en non sature)
  45 :             seulement pour les lois exponentielle et logarithme.
  46 : 
  47 :     CHPO2 : saturation reduite, S
  48 :             (CHAMP-POINT support de Pw)
  49 : 
  50 :     CHPO3 : teneur en eau, TH2O
  51 :             (CHAMP-POINT support de Pw)
  52 : 
  53 :     CHPO4 : capacite calculee analytiquement
  54 :             (CHAMP-POINT support de Pw)
  55 : 
  56 : 
  57 :     Remarque :
  58 :     ----------
  59 :      Si le champ CHPO1 est absent, on teste les arguments de la procedure 
  60 :      et on ne fait pas de calculs.
  61 :