从TLE的B*提取航天器SI弹道系数：p₀的正确值与单位

我正在编写Python脚本，通过分析近地轨道（LEO）航天器半长轴平方根随时间的变化，计算大气阻力的影响。用到的控制方程如下：

np.sqrt(a_final) - np.sqrt(a_initial) = (- np.sqrt(mu) * C_d * A * p * delta_time) / (2 * m)

其中各参数定义：

• a_initial：初始半长轴
• a_final：最终半长轴
• mu：引力参数（单位：m³/s²）
• C_d：阻力系数（无量纲）
• A：横截面积（单位：m²）
• m：航天器质量（单位：kg）
• p：轨道高度处的大气密度（单位：kg/m³）
• delta_time：时间间隔（单位：秒）

相关的弹道系数与B*的计算公式如下：

Ballistic_coefficient = (C_d * A) / m  # 单位：m²/kg

B* = (Ballistic_coefficient * p_0) / 2 # 单位：1/m

其中p₀为参考大气密度。

我的困惑：
Celestrak仅将p₀定义为参考大气密度，但定义不够明确；而Spacetrack报告中给出的p₀单位为kg/m²，这会导致单位不匹配，在从B*反推弹道系数时产生混淆。

请问，为了从给定TLE的B*值中提取SI单位（m²/kg）的弹道系数，以便正确代入半长轴衰减公式，p₀的正确值和单位是什么？

内容的提问来源于stack exchange，提问作者Aastha Upadhyay