def _overpressureatscaledtime(r, y, alt, t):
sgr = r / y**(1.0/3)
shob = alt / y**(1.0/3)
x_m = _scaledmachstemformationrange(y, alt)
v = _slantrangescalingfactor(r, y, alt)
r1 = _scale_slant_range(r, y, alt) / v
ta_air = _scaledfreeairblastwavetoa(r1) * v
dp = _scaledopposphasedur(r, y, alt)
return _opatscaledtime(r, y, alt, sgr, shob, x_m, ta_air, dp, t)
# In lieu of the 20-point Gauss-Legendre quadrature used in the original
# BLAST.EXE, this fuction uses scipy.integrate.quad to call the venerable FORTRAN
# library QUADPACK and perform a Gauss-Kronod quadrature. This appears
# to be more accurate than the BLAST.EXE help file claims for the original
# approach, which is not surprising as it uses an adaptive algorithm that
# attempts to reduce error to within a particular tolerance.
# IPTOTAL
评论列表
文章目录