def __init__(self, stn_pos, **kwds):
"""
Setup Chapman electron density profile F2 peak height derivative
slant integrator.
"""
z, Nm, Hm, H_O = SYM.symbols('z Nm Hm H_O')
f_sym = chapman_sym_scaled(z, Nm, Hm, H_O)
DHm_sym = SYM.diff(f_sym, Hm)
f = SYM.lambdify((z, Nm, Hm, H_O),
DHm_sym,
modules='numexpr')
wrapper = lambda pos, *args: f(pos.height / 1e3, *args)
super(DHmChapmanSI, self).__init__(wrapper, stn_pos, **kwds)
评论列表
文章目录