/// <summary>
/// 月亮到中升降时刻计算,传入jd含义与St()函数相同
/// </summary>
/// <param name="jd"></param>
/// <returns></returns>
public static LunarInfoListT <double> Mt(double jd)
{
SZJ.dt = JD.deltatT2(jd);
SZJ.E = ZB.hcjj(jd / 36525);
jd -= LunarHelper.mod2(0.1726222 + 0.966136808032357 * jd - 0.0366 * SZJ.dt - SZJ.L / LunarHelper.pi2, 1); //查找最靠近当日中午的月上中天,mod2的第1参数为本地时角近似值
LunarInfoListT <double> r = new LunarInfoListT <double>(3, 0d);
double sv = LunarHelper.pi2 * 0.966;
r.z__ = r.s__ = r.j__ = r.c__ = r.h__ = jd;
SZJ.Mcoord(jd, 1, r); //月亮坐标
r.s__ += (-r.H0 - r.H) / sv;
r.j__ += (r.H0 - r.H) / sv;
r.z__ += (0 - r.H) / sv;
SZJ.Mcoord(r.s__, 1, r); r.s__ += (-r.H0 - r.H) / sv;
SZJ.Mcoord(r.j__, 1, r); r.j__ += (+r.H0 - r.H) / sv;
SZJ.Mcoord(r.z__, 0, r); r.z__ += (0 - r.H) / sv;
return(r);
}
/// <summary>
/// 太阳到中升降时刻计算,传入jd是当地中午12点时间对应的2000年首起算的格林尼治时间UT
/// </summary>
/// <param name="jd"></param>
/// <returns></returns>
public static LunarInfoListT <double> St(double jd)
{
SZJ.dt = JD.deltatT2(jd);
SZJ.E = ZB.hcjj(jd / 36525);
jd -= LunarHelper.mod2(jd - SZJ.L / LunarHelper.pi2, 1); //查找最靠近当日中午的日上中天,mod2的第1参数为本地时角近似值
LunarInfoListT <double> r = new LunarInfoListT <double>(3, 0d);
double sv = LunarHelper.pi2;
r.z__ = r.s__ = r.j__ = r.c__ = r.h__ = jd;
SZJ.Scoord(jd, 1, 1, r); //太阳坐标
r.s__ += (-r.H0 - r.H) / sv; //升起
r.j__ += (r.H0 - r.H) / sv; //降落
r.c__ += (-r.H1 - r.H) / sv; //民用晨
r.h__ += (r.H1 - r.H) / sv; //民用昏
r.z__ += (0 - r.H) / sv; //中天
SZJ.Scoord(r.s__, 1, 0, r); r.s__ += (-r.H0 - r.H) / sv;
SZJ.Scoord(r.j__, 1, 0, r); r.j__ += (+r.H0 - r.H) / sv;
SZJ.Scoord(r.c__, 0, 1, r); r.c__ += (-r.H1 - r.H) / sv;
SZJ.Scoord(r.h__, 0, 1, r); r.h__ += (+r.H1 - r.H) / sv;
SZJ.Scoord(r.z__, 0, 0, r); r.z__ += (0 - r.H) / sv;
return(r);
}
