Télécharger fgammb.eso

Retour à la liste

Numérotation des lignes :

fgammb
  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.  
  54.  

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