Astronomical Applications Department, U.S. Naval Observatory Sensitivity Integral Page 29
fortran H optimized
=
mode
double
=
precision
double
,
,
,
,
(
=
filename
"d:/FAMEStuff/PhotonSensitivity/error_geom.f" )
For reference, here are the fortran subroutines.
(
)
fortran
,
,
,
F optimized
=
mode
double
=
precision
double
doubleprecision function F(delta,nu1,a0,a1,a2,a3,a4,s,r,g)
doubleprecision delta
doubleprecision nu1
doubleprecision a0
doubleprecision a1
doubleprecision a2
doubleprecision a3
doubleprecision a4
doubleprecision s
doubleprecision r
doubleprecision g
doubleprecision t180
doubleprecision t181
doubleprecision t182
doubleprecision t185
doubleprecision t186
doubleprecision t188
doubleprecision t190
doubleprecision t194
doubleprecision t3
doubleprecision t4
doubleprecision t8
t194 = nu1*(delta+2.D0)
t185 = delta**2
t190 = t185*delta
t186 = nu1**2
t188 = t186**2
t182 = t185**2
t181 = t190**2
t180 = t185*t190
t3 = s**2
t4 = r**2
t8 = (-1.D0+g)**2.D0
F = t3*t4*nu1*t188*delta*t8*(504.D0*a0*(5.D0+10.D0*delta+10.D0*t
#185+5.D0*t190+t182)+(420.D0*a1*(t185+delta+1.D0)*t194+280.D0*a4*(t
#181+6.D0*t180+15.D0*t182+21.D0*t190+18.D0*t185+9.D0*delta+3.D0)*t1
#88)*(t185+3.D0*delta+3.D0)+(315.D0*a3*(t185+2.D0*delta+2.D0)*(t182
#+4.D0*t190+6.D0*t185+4.D0*delta+2.D0)*t194+360.D0*a2*(t181+7.D0*t1
#80+21.D0*t182+35.D0*t190+35.D0*t185+21.D0*delta+7.D0))*t186)/2520.
#D0
return
end
(
)
fortran
,
,
,
G optimized
=
mode
double
=
precision
double
doubleprecision function G(delta,nu1,a0,a1,a2,a3,a4)
doubleprecision delta
doubleprecision nu1
doubleprecision a0
doubleprecision a1
doubleprecision a2
doubleprecision a3
Page 29
