Télécharger la plaquette Cast3M


fr en

Télécharger vlh1.eso

Retour à la liste

Numérotation des lignes : 

vlh1

  1. C VLH1      SOURCE    CHAT      05/01/13    04:07:41     5004
  2.       SUBROUTINE VLH1(H1,H1P,X,G)
  3. C************************************************************************
  4. C
  5. C PROJET            :  CASTEM 2000
  6. C
  7. C NOM               :  VLH1 (subroutine du groupe RIECOM.ESO)
  8. C
  9. C DESCRIPTION       :  Fonction H1 de SMOLLER (voir RIEMAN)
  10. C
  11. C LANGAGE           :  FORTRAN 77
  12. C
  13. C AUTEUR            :  P. GALON DRN/DMT/SEMT/TTMF
  14. C
  15. C************************************************************************
  16. C
  17. C APPELES           : aucun
  18. C
  19. C************************************************************************
  20. C
  21. C Entrées : X, G
  22. C
  23. C Sorties :    H1   FONCTION H1
  24. C              H1P  FONCTION DERIVEE DE H1
  25. C
  26. C************************************************************************
  27. C
  28. C HISTORIQUE (Anomalies et modifications éventuelles)
  29. C
  30. C HISTORIQUE :  créé le 1.10.96
  31. C
  32. C HISTORIQUE :  modifié par A. BECCANTINI le 7.1.98 pour
  33. C               ajouter des commentaires et eliminer la
  34. C               partié relative au cas X = 0.
  35. C
  36. C
  37. C************************************************************************
  38. C
  39. C N.B.: Toutes le variables sont DECLAREES
  40. C
  41.       IMPLICIT INTEGER(I-N)
  42.       IMPLICIT REAL*8 (A-H,O-Z)
  43.       REAL*8 H1, H1P, X, G, GM1, UNSGM1, T, ETX, AUX1, AUX2, B, EX
  44. C
  45.       GM1 = G-1.0D0
  46.       UNSGM1 = 1.0D0 / GM1
  47.       T = 0.5D0*GM1/G
  48.       IF(X .GE. 0.D0) THEN
  49. C
  50. C******** X>0
  51. C
  52.          ETX = EXP(-T*X)
  53.          H1 = 2.D0*UNSGM1*(1.D0-ETX)
  54.          H1P = 2.D0*UNSGM1*T*ETX
  55. C
  56. C***** NB Function continue et derivable in X = 0
  57. C
  58. C      X=0
  59. C
  60. C      H1 = 0.0D0
  61. C      H1P = 1.D0/G
  62. C
  63.       ELSE
  64. C
  65. C******** X<0
  66. C
  67.          B = (G+1.0D0)*UNSGM1
  68.          EX = EXP(-X)
  69.          AUX1 = 1.0D0 / SQRT(1.D0+B*EX)
  70.          AUX2 = 2.0D0*SQRT(T)*UNSGM1*AUX1
  71.          H1 = AUX2*(1.0D0-EX)
  72.          H1P = AUX2*EX* ( 1.D0 + 0.5D0*(1.D0-EX)*B*AUX1*AUX1 )
  73.       ENDIF
  74.       RETURN
  75.       END
  76.  
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  

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