《寿星天文历》软件我已用了好几年,的确实用.

感谢您的推广。有空请常来历法知识版块做客，大家一块探讨历法相关知识。

许老师，太佩服你了

今天查日食，发现一个小BUG，
比如2032-11-3有一次日偏食。
数据显示：
*****偏食始： 2032-11-03 03:23:35 + 56°12',+ 52°46'
*****中心始： 2056-11-01 11:20:04 +102°21',- 11°21'
*****中心终： 2056-11-01 11:20:04 +102°21',- 11°21'
*****偏食终： 2032-11-03 07:44:57 +141°14',+ 29°28'
*****//因为月影中心不在地球上，所以软件计算时，就出了差错，显示为2056年了。

这是显示问题，可以把

   re.gk1 = this.qrd(t2+jd,vx,vy,0); //中心始
   re.gk2 = this.qrd(t3+jd,vx,vy,0); //中心终

改为，这样，地外中心线统一显增为2000年1月1日

  if(re.d<1){
   re.gk1 = this.qrd(t2+jd,vx,vy,0); //中心始
   re.gk2 = this.qrd(t3+jd,vx,vy,0); //中心终
  }else{
   re.gk1 = [0,0,0];
   re.gk2 = [0,0,0];
  }

更彻底解决问题，要追朔到显示输出部分。

多谢许兄。其它的日食也有类似“显示问题”。_________比如“食甚: 2020-12-15 00:07:44 UTC，杭州，日偏食。这肯定是观测不到的。”_____________这样的例子很多。因为这些日食，都发生的地球的对面国家或地区。似乎是解方程时，许兄没舍弃其中一个解。直线与球面相交二点。只能有其中一交点是发生日食。另一交点，在地球的背面，不可能同时发生日食的。

这段代码

  if(rmin<=G.sr-G.mr){ //环食计算
   this.LX='环';
   this.sT[3] = this.lineT(G,v,u, G.sr-G.mr, 0); //食既
   this.secXY(this.sT[3],L,fa,high,g);
   this.sT[3] = this.lineT(g,v,u, g.sr-g.mr, 0); //食既再算1次

   this.sT[4] = this.lineT(G,v,u, G.sr-G.mr, 1); //生光
   this.secXY(this.sT[4],L,fa,high,g);
   this.sT[4] = this.lineT(g,v,u, g.sr-g.mr, 1); //生光再算1次
   this.dur = this.sT[4]-this.sT[3];
  }

 

之后补三行代码试试：

 

  if(rmin<=G.sr-G.mr){ //环食计算
   this.LX='环';
   this.sT[3] = this.lineT(G,v,u, G.sr-G.mr, 0); //食既
   this.secXY(this.sT[3],L,fa,high,g);
   this.sT[3] = this.lineT(g,v,u, g.sr-g.mr, 0); //食既再算1次

   this.sT[4] = this.lineT(G,v,u, G.sr-G.mr, 1); //生光
   this.secXY(this.sT[4],L,fa,high,g);
   this.sT[4] = this.lineT(g,v,u, g.sr-g.mr, 1); //生光再算1次
   this.dur = this.sT[4]-this.sT[3];
  }
  if(this.sT[0]<this.sun_s || this.sT[0]>this.sun_j) this.sT[0]=0;
  if(this.sT[1]<this.sun_s || this.sT[1]>this.sun_j) this.sT[1]=0;
  if(this.sT[2]<this.sun_s || this.sT[2]>this.sun_j) this.sT[2]=0;

昨天下载了vsop2013，有打算将寿星更换为vsop2013

这是一个大工程。如果更新版本，可以顺带增加一些常用功能。