Sources Esope | Cast3M

fbg2
C FBG2      SOURCE    CHAT      05/01/12    23:57:00     5004
 
 
             FUNCTION FBG2(T)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
 
      t1 = t**2
      t2 = t1**2
      t4 = t2**2
      t5 = t4**2
      t6 = t5**2
      t7 = t6**2
      t9 = 0.2126494641379152D-89*t7*t2*t1+1.D0
      t10 = 1.D0/t9
      t11 = t2*t
      t12 = 0.1316872427983539D0*t11+1.D0
      t13 = 1.D0/t12
      t14 = -0.156710172162D4-3020.D0*t10-0.595089827838D4*t13-0.1146775
     #430405D2*t1+0.458710172162D3*t
      t16 = t9**2
      t19 = 1.D0/t16*t7*t11
      t20 = t12**2
      t22 = 1.D0/t20*t2
      t23 = 0.4495409671875527D-84*t19+0.3918286932266667D4*t22-0.229355
     #086081D2*t+0.458710172162D3
      t25 = 1.D0/t14
      t26 = sign(1.D0,t14)*t23*t25
      t27 = abs(t14)
      t28 = t14**2
      t31 = t27/t28*t23
      t33 = 0.458710172162D4+3020.D0*t10+0.595089827838D4*t13+0.11467754
     #30405D2*t1-0.458710172162D3*t
      t35 = t27*t25
      t36 = 1.D0/2.D0+t35/2.D0
      t37 = -0.4495409671875527D-84*t19-0.3918286932266667D4*t22+0.22935
     #5086081D2*t-0.458710172162D3
      t40 = 0.156710172162D4+3020.D0*t10+0.595089827838D4*t13+0.11467754
     #30405D2*t1-0.458710172162D3*t
      t42 = 1.D0/2.D0-t35/2.D0
      s1 = 31536000.D0
      s3 = -0.8765881258141834D0*t26+0.8765881258141834D0*t31+0.58052193
     #7625287D-3*(t26/2.D0-t31/2.D0)*t33+0.580521937625287D-3*t36*t37+0.
     #3152398523064585D-3*(-t26/2.D0+t31/2.D0)*t40+0.3152398523064585D-3
     #*t42*t37
      s4 = exp(-0.2514499507768026D2-0.8765881258141834D0*t35+0.58052193
     #7625287D-3*t36*t33+0.3152398523064585D-3*t42*t40)
      s2 = s3*s4
      t50 = s1*s2
      FBG2=T50
      RETURN
      END