Esope source files

   1 : C FLA016    SOURCE    TTMF3     12/07/05    21:15:38     7425           
   2 :       FUNCTION FLA016(X,ROE,A,B,C,DD,EE,FF)
   3 : C------------------------------------------------------------------------
   4 : C Calcul de f(x) dans le cas où le débit n'est pas limité par les 
   5 : C conditions d'entrée : recherche de la racine par dichotomie
   6 : C (voir fla015.eso)
   7 : C------------------------------------------------------------------------
   8 : C
   9 : C---------------------------
  10 : C Parametres Entree/Sortie :
  11 : C---------------------------
  12 : C
  13 : C E/  X      : flottant : Debit au travers du PAR (kg/s) 
  14 : C E/  ROE    : flottant : Densité du mélange à l'entrée du PAR (kg/m3)
  15 : C E/  A      : flottant : Coefficient devant X^2
  16 : C E/  B      : flottant : Coefficient devant X
  17 : C E/  DD     : flottant : Coefficient
  18 : C E/  EE     : flottant : Coefficient 
  19 : C E/  FF     : flottant : Coefficient
  20 : C  /S FLA016 : flottant : Valeur de la fonction f(X)
  21 : C
  22 : C------------------------------------------------------------------------
  23 : C
  24 : C Langage : FORTRAN 77
  25 : C
  26 : C Mise en oeuvre : H. Paillère (1997, TTMF)
  27 : C
  28 : C------------------------------------------------------------------------
  29 :       IMPLICIT INTEGER(I-N)
  30 :       IMPLICIT REAL*8 (A-H,O-Z)
  31 : C
  32 : C Polynome de degré 6 tirée de Bernouilli dont la plus grande racine au
  33 : C dela des poles (EE-DD) et FF, si elle existe, est le débit cherché
  34 :       P61 = A*X*X + B*X + C*ROE*ROE
  35 :       P62 = (X-EE+DD)*(X-EE+DD)*(X-FF)*(X-FF)
  36 :       P63 = C*ROE*ROE*(X-EE)*(X-EE)*X*X
  37 : C
  38 :       FLA016 = P61*P62 - P63
  39 : C
  40 :       RETURN
  41 :       END
  42 : 
  43 :  
  44 :  
  45 :  
  46 :  
  47 :