Esope source files

   1 : C FGAMMB    SOURCE    CHAT      05/01/12    23:59:20     5004           
   2 : C***********************************************************************
   3 : C*                                                                     *
   4 : C* PROJET      : Opérateur CLMI                                        *
   5 : C* NOM         : FGAMMB                                                *
   6 : C* DESCRIPTION : Dans ce fichier, les deux fonctions permettent de     *
   7 : C*               calculer la valeur de la relation f(gamma) et de sa   *
   8 : C*               dérivée en X (Relations de fermeture de Cousteix)     *
   9 : C* LANGAGE     : Esope                                                 *
  10 : C* AUTEUR      : Guillaume VENCO - DRN/DMT/SEMT/LTMF                   *
  11 : C*                                                                     *
  12 : C***********************************************************************
  13 : C*                                                                     *
  14 : C* APPELES     : Aucun                                                 *
  15 : C*                                                                     *
  16 : C***********************************************************************
  17 : C*                                                                     *
  18 : C* ENTREES     : D1  : -Epaisseur de déplacement                       *
  19 : C*               UE  : -Champ de vitesse extérieure                    *
  20 : C*               HN  : -Facteur de forme                               *
  21 : C*               X   : -Point de calcul                                *
  22 : C*                                                                     * 
  23 : C* SORTIES     : fgamma: -Valeur de la fonction f en X                 *
  24 : C*               dfga  : -Valeur de la dérivée en X                    *
  25 : C*                                                                     *
  26 : C***********************************************************************
  27 : C*                                                                     *
  28 : C* VERSION     : 21/04/2000                                            *
  29 : C* CREATION    :                                                       *
  30 : C*                                                                     *
  31 : C***********************************************************************
  32 : 
  33 : 
  34 : C*       *********************************
  35 : C*       ********** Fonction f  **********
  36 : C*       *********************************
  37 : 
  38 :       FUNCTION fgammb(UE,D1,HN,X)
  39 :       IMPLICIT INTEGER(I-N)
  40 :       IMPLICIT REAL*8 (A-H,O-Z)
  41 :       REAL*8 UE,D1,HN,X
  42 :       REAL*8 a1,a2,a3
  43 : 
  44 :       a1 = (2.12+(1./0.41)*LOG((UE*D1)/1.D-6))
  45 :       a2 = -4.25*(((HN-1)/HN)**(0.5))
  46 :       a3 = (HN-2)/HN
  47 :       
  48 :       fgammb = a1*X+a2*(X**0.5)+a3
  49 : 
  50 :       RETURN
  51 :       END
  52 :  
  53 :