Télécharger fbg2.eso

Retour à la liste

Numérotation des lignes :

  1. C FBG2 SOURCE CHAT 05/01/12 23:57:00 5004
  2.  
  3.  
  4. FUNCTION FBG2(T)
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8 (A-H,O-Z)
  7.  
  8. t1 = t**2
  9. t2 = t1**2
  10. t4 = t2**2
  11. t5 = t4**2
  12. t6 = t5**2
  13. t7 = t6**2
  14. t9 = 0.2126494641379152D-89*t7*t2*t1+1.D0
  15. t10 = 1.D0/t9
  16. t11 = t2*t
  17. t12 = 0.1316872427983539D0*t11+1.D0
  18. t13 = 1.D0/t12
  19. t14 = -0.156710172162D4-3020.D0*t10-0.595089827838D4*t13-0.1146775
  20. #430405D2*t1+0.458710172162D3*t
  21. t16 = t9**2
  22. t19 = 1.D0/t16*t7*t11
  23. t20 = t12**2
  24. t22 = 1.D0/t20*t2
  25. t23 = 0.4495409671875527D-84*t19+0.3918286932266667D4*t22-0.229355
  26. #086081D2*t+0.458710172162D3
  27. t25 = 1.D0/t14
  28. t26 = sign(1.D0,t14)*t23*t25
  29. t27 = abs(t14)
  30. t28 = t14**2
  31. t31 = t27/t28*t23
  32. t33 = 0.458710172162D4+3020.D0*t10+0.595089827838D4*t13+0.11467754
  33. #30405D2*t1-0.458710172162D3*t
  34. t35 = t27*t25
  35. t36 = 1.D0/2.D0+t35/2.D0
  36. t37 = -0.4495409671875527D-84*t19-0.3918286932266667D4*t22+0.22935
  37. #5086081D2*t-0.458710172162D3
  38. t40 = 0.156710172162D4+3020.D0*t10+0.595089827838D4*t13+0.11467754
  39. #30405D2*t1-0.458710172162D3*t
  40. t42 = 1.D0/2.D0-t35/2.D0
  41. s1 = 31536000.D0
  42. s3 = -0.8765881258141834D0*t26+0.8765881258141834D0*t31+0.58052193
  43. #7625287D-3*(t26/2.D0-t31/2.D0)*t33+0.580521937625287D-3*t36*t37+0.
  44. #3152398523064585D-3*(-t26/2.D0+t31/2.D0)*t40+0.3152398523064585D-3
  45. #*t42*t37
  46. s4 = exp(-0.2514499507768026D2-0.8765881258141834D0*t35+0.58052193
  47. #7625287D-3*t36*t33+0.3152398523064585D-3*t42*t40)
  48. s2 = s3*s4
  49. t50 = s1*s2
  50. FBG2=T50
  51. RETURN
  52. END
  53.  
  54.  
  55.  

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