def nrgaiatot_3D_deriv(rad):
lrho = np.log(rhogaiatot_3D(rad))
lr = np.log(rad)
import gl_helper as gh
return -gh.derivcoarse(lrho, lr)
def get_kzpars():
Kz_zthd = -2.0 * F * zth # [with pc]
Sigz_zth = abs(Kz_zthd) / (2.0*np.pi*gp.G1)
denth = gh.derivcoarse(Sigz_zth,zth)
kzpars = gh.ipol(zth,abs(denth),gp.xipol)*(2.0*np.pi*gp.G1)
return kzpars