def __init__(self, stn_pos, **kwds):
"""
Setup Chapman electron density profile plasma scale 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)
DH_O_sym = SYM.diff(f_sym, H_O)
f = SYM.lambdify((z, Nm, Hm, H_O),
DH_O_sym,
modules='numexpr')
wrapper = lambda pos, *args: f(pos.height / 1e3, *args)
super(DH_OChapmanSI, self).__init__(wrapper, stn_pos, **kwds)
评论列表
文章目录