0.31,1.9,14390.93,-3,0,0,0.31,4.6,8852.22,4,51,0,
0.31,0.6,6551.97,-2,0,0,0.30,4.7,-7595.61,-4,-51,0,
0.30,1.9,7143.05,0,0,0,0.29,3.2,-1428.95,3,0,0,
0.27,4.9,-2476.01,-1,0,0,0.26,3.2,41748.23,5,100,0,
0.25,3.4,-1009.85,-6,-70,0,0.24,1.9,32514.46,-4,0,0,
0.24,3.3,32933.21,0,0,0,0.21,3.6,22233.29,-5,0,0,
0.21,4.4,47847.67,0,0,0,0.21,3.9,23418.91,-3,0,0,
0.17,5.8,14951.65,-2,0,0,0.16,2.0,38890.67,-2,0,0),
new Array(//MB1
0.074,4.10,6480.986,0,7,0,0.030,0.9,7737.590,0,7,0,
0.022,5.0,15019.227,-3,-19,0,0.020,1.1,23347.918,-1,6,0,
0.019,1.7,-1847.705,-2,-18,0,0.017,5.6,16133.856,1,24,0,
0.016,0,9061.768,-1,0,0,0.014,3.9,733.077,-2,-26,0,
0.013,2.4,17390.460,1,25,0,0.013,5.6,8399.685,0,0,0,
0.013,0.9,-523.527,-2,-26,0,0.012,3.2,7805.164,-1,0,0,
0.011,3.7,8435.767,-2,0,0,0.011,5.9,8431.165,1,0,0)
),

MR:new Array( //精度3千米
new Array(//MR0
385000.5,0,0,0,0,0,
20905.4,5.497147,8328.6914246,1.52292,25.07,-0.1236,3699.1,4.89979,7214.062865,-2.1848,-18.9,0.083,
2956.0,0.97216,15542.754290,-0.6618,6.2,-0.041,569.9,1.5695,16657.382849,3.046,50,-0.25,
246.2,5.6858,-1114.62856,-3.708,-44,0.21,204.6,1.0153,14914.45233,-0.635,6,-0.04,
170.7,3.3277,23871.44571,0.86,31,-0.2,152.1,4.943,6585.76091,-2.16,-19,0.1,
129.6,0.743,-7700.38947,-1.55,-25,0.1,108.7,5.198,7771.37714,-0.33,3,0,
104.8,2.312,8956.99338,1.50,25,-0.1,79.7,5.383,-8538.24089,2.80,26,-0.1,
48.9,6.240,628.3020,-0.03,0,0,34.8,2.730,22756.8172,-2.85,-13,0,
30.8,4.071,16171.0562,-0.69,6,0,24.2,1.715,7842.3648,-2.21,-19,0.1,
23.2,3.925,24986.0743,4.57,75,-0.4,21.6,0.375,14428.1257,-4.37,-38,0.2,
16.7,2.014,8399.6791,-0.4,3,0,14.4,3.33,-9443.3200,-5.2,-69,0,
12.8,3.37,23243.1438,0.9,31,0,11.6,5.09,31085.5086,-1.3,12,0,
10.4,5.68,32200.1371,2.4,56,0,10.3,0.86,-1324.1780,0.6,7,0,
10.1,5.73,-1742.9305,-3.7,-44,0,9.9,1.06,14286.1504,-0.6,6,0,
8.8,4.79,-9652.8694,-0.9,-18,0,8.4,5.98,-557.3143,-1.9,-22,0,
7.0,4.67,-16029.0809,-3.1,-50,0,6.3,1.27,16100.0686,1.2,28,0,
5.8,4.67,17285.6848,3.0,50,0,5.0,4.99,5957.459,-2.1,-19,0,
4.4,4.60,-209.549,4.3,51,0,4.1,3.21,7004.513,2.1,32,0,
4.0,2.77,22128.515,-2.8,-13,0,3.3,0.67,14985.440,-2.5,-16,0,
3.1,0.11,16866.932,-1.3,0,0,2.6,0.14,24499.748,0.8,31,0,
2.4,1.67,8470.667,-2.2,-19,0,2.1,0.70,-7072.088,-1.6,-25,0,
1.9,0.42,13799.824,-4.3,-38,0,1.7,3.63,-8886.006,-3,-47,0,
1.6,5.1,30457.207,-1,12,0,1.4,1.2,39414.200,0,37,0,
1.4,6.2,23314.131,-1,9,0,1.2,2.3,9585.295,1,25,0,
1.1,6.3,33314.766,6,100,0,1.1,6.2,1256.604,0,0,0,
1.1,4.1,8364.740,-2,-19,0,0.9,4.4,16728.371,1,28,0,
0.9,4.6,6656.749,-4,-41,0,0.9,2.8,70.988,-2,-22,0,
0.8,5.7,31571.835,2,56,0,0.8,5.1,-9095.555,1,0,0,
0.8,1.0,-17772.011,-7,-94,0,0.8,2.7,15752.304,-5,-45,0,
0.7,0.3,8326.390,3,50,0,0.7,1.3,8330.993,0,0,0,
0.7,1.8,40528.829,4,81,0,0.7,3.4,22614.842,1,31,0,
0.7,0.9,-1952.480,1,7,0,0.6,6.0,8393.126,-2,-19,0,
0.6,5.0,24080.995,-3,-20,0,0.6,5.8,23385.119,-3,0,0,
0.5,4.3,6099.43,-6,-63,0,0.5,1.8,14218.58,0,0,0,
0.5,5.2,7143.08,0,0,0,0.5,3.4,-10071.62,-5,-69,0,
0.5,2.4,-17981.56,-2,0,0,0.5,4.9,-8294.93,-2,0,0,
0.5,0.2,8362.45,1,0,0,0.4,4.5,29970.88,-5,0,0,
0.4,2.3,-24357.77,-5,-75,0,0.4,1.1,13657.85,-1,0,0,
0.4,0.5,8311.77,-2,0,0,0.4,3.6,24428.76,3,53,0,
0.4,0.7,25614.38,5,75,0,0.3,5.8,-2371.23,-4,0,0,
0.3,0.9,9166.54,-3,0,0,0.3,0.4,-8257.70,-3,0,0,
0.3,4.8,-10281.17,-1,0,0,0.3,5.8,5889.88,-2,0,0,
0.3,0.6,38299.57,-4,0,0,0.3,5.6,15333.20,4,57,0,
0.3,2.8,21500.21,-3,0,0,0.3,0.7,14357.14,-2,0,0),
new Array(//MR1
0.514,4.16,14914.4523,-0.6,6,0,0.382,1.80,6585.7609,-2.2,-19,0,
0.327,2.40,7700.3895,1.5,25,0,0.264,5.45,8956.9934,1.5,25,0,
0.123,3.10,628.302,0,0,0,0.078,0.93,16171.056,-0.7,6,0,
0.061,4.86,7842.365,-2.2,-19,0,0.050,4.2,14286.150,-1,6,0,
0.042,5.2,8399.679,0,0,0,0.032,0.2,23243.144,1,31,0,
0.025,2.6,-1742.931,-4,-44,0,0.025,1.8,5957.459,-2,-19,0,
0.018,4.8,16029.081,3,50,0,0.014,1.5,17285.68,3,50,0,
0.014,1.0,15542.75,-1,0,0,0.013,5.0,8326.39,3,50,0,
0.012,4.8,8470.67,-2,0,0,0.012,2.8,8330.99,0,0,0,
0.011,2.4,7072.09,2,0,0,0.010,5.9,22128.52,-3,0,0),
new Array(//MR2
0.0015,4.2,14914.452,-1,6,0,0.0011,1.8,6585.761,-2,-19,0,
0.0009,2.4,7700.389,2,25,0,0.0008,5.5,8956.993,1,25,0)
),

//=====================

//星歷函數(日月球面坐標計算)
Enn:function(ob,t,n){ //計算E_L0或E_L1或E_L2等
  var i,j,F,N,v=0,tn=1,c;
  if(ob==this.EL){
    var t2=t*t, t3=t2*t, t4=t3*t, t5=t4*t; //千年數的各次方
    v += 1753469512 + 6283319653318*t + 529674*t2 + 432*t3 - 1124*t4 - 9*t5 + 630 * Math.cos(6+3*t); //地球平黃經(已擬合DE406)
  }
  n*=3; if(n<0) n = ob[0].length;
  for(i=0;i<ob.length;i++,tn*=t){
   F = ob[i];
   N = int2( n*F.length/ob[0].length+0.5 );  if(i) N+=3;  if(N >= F.length) N=F.length;
   for(j=0,c=0;j<N;j+=3) c+=F[j]*Math.cos(F[j+1] +t*F[j+2]);
   v += c*tn;
  }
  return v/1000000000;
},
E_coord:function(t,re,n1,n2,n3){ //返回地球坐標,t為世紀數
  t/=10;
  re[0]= this.Enn( this.EL, t, n1);
  re[1]= this.Enn( this.EB, t, n2);
  re[2]= this.Enn( this.ER, t, n3);
},

E_Lon:function(t,n){ return this.Enn(this.EL,t/10,n); }, //地球經度計算,返回Date分點黃經,傳入世紀數、取項數

Mnn:function(ob,t,n){ //計算ML0或ML1或ML2
  var i,j,F,N,v=0,tn=1,c;
  var t2=t*t,t3=t2*t,t4=t3*t,t5=t4*t,tx=t-10;
  if(ob==this.ML){
    v += (3.81034409 + 8399.684730072*t -3.319e-05*t2 + 3.11e-08*t3 - 2.033e-10*t4)*rad; //月球平黃經(弧度)
    v += 5028.792262*t + 1.1124406*t2 + 0.00007699*t3 - 0.000023479*t4 -0.0000000178*t5;  //歲差(角秒)
    if(tx>0) v += -0.866 +1.43*tx +0.054*tx*tx; //對公元3000年至公元5000年的擬合,最大誤差小于10角秒
  }
  t2/=1e4,t3/=1e8,t4/=1e8;
  n*=6; if(n<0) n = ob[0].length;
  for(i=0;i<ob.length;i++,tn*=t){
   F=ob[i];
   N = int2( n*F.length/ob[0].length+0.5 );  if(i) N+=6;  if(N >= F.length) N=F.length;
   for(j=0,c=0;j<N;j+=6) c+=F[j]*Math.cos(F[j+1] +t*F[j+2] +t2*F[j+3] +t3*F[j+4] +t4*F[j+5]);
   v += c*tn;
  }
  if(ob!=this.MR) v/=rad;
  return v;
},
M_coord:function(t,re,n1,n2,n3){ //返回月球坐標,n1,n2,n3為各坐標所取的項數
  re[0] = this.Mnn( this.ML, t, n1 );
  re[1] = this.Mnn( this.MB, t, n2 );
  re[2] = this.Mnn( this.MR, t, n3 );
},
M_Lon:function(t,n){ return this.Mnn(this.ML,t,n); }, //月球經度計算,返回Date分點黃經,傳入世紀數,n是項數比例

//=========================
E_v:function(t){ //地球速度,t是世紀數,誤差小于萬分3
 var f=628.307585*t;
 return 628.332 +21 *Math.sin(1.527+f) +0.44 *Math.sin(1.48+f*2) 
         +0.129*Math.sin(5.82+f)*t +0.00055*Math.sin(4.21+f)*t*t;
},
M_v:function(t){ //月球速度計算,傳入世經數
  var v = 8399.71 - 914*Math.sin( 0.7848 + 8328.691425*t + 0.0001523*t*t ); //誤差小于5%
  v-=179*Math.sin( 2.543  +15542.7543*t )  //誤差小于0.3%
    +160*Math.sin( 0.1874 + 7214.0629*t )
    +62 *Math.sin( 3.14   +16657.3828*t )
    +34 *Math.sin( 4.827  +16866.9323*t )
    +22 *Math.sin( 4.9    +23871.4457*t )
    +12 *Math.sin( 2.59   +14914.4523*t )
    +7  *Math.sin( 0.23   + 6585.7609*t )
    +5  *Math.sin( 0.9    +25195.624 *t )
    +5  *Math.sin( 2.32   - 7700.3895*t )
    +5  *Math.sin( 3.88   + 8956.9934*t )
    +5  *Math.sin( 0.49   + 7771.3771*t );
   return v;
},

//=========================
MS_aLon:function(t,Mn,Sn){ //月日視黃經的差值
  return this.M_Lon(t,Mn) + ZB.gxc_moonLon(t) - ( this.E_Lon(t,Sn) + ZB.gxc_sunLon(t) + Math.PI );
},
S_aLon:function(t,n){  //太陽視黃經
 return this.E_Lon(t,n) + ZB.nutationLon(t) + ZB.gxc_sunLon(t) + Math.PI; //注意，這里的章動計算很耗時
},

//=========================

E_Lon_t:function(W){ //已知地球真黃經求時間
  var t,v= 628.3319653318;
  t  = ( W - 1.75347          )/v; v=this.E_v(t);   //v的精度0.03%，詳見原文
  t += ( W - this.E_Lon(t,10) )/v; v=this.E_v(t);   //再算一次v有助于提高精度,不算也可以
  t += ( W - this.E_Lon(t,-1) )/v;
  return t;
},
M_Lon_t:function(W){ //已知真月球黃經求時間
  var t,v= 8399.70911033384;
  t  = ( W - 3.81034          )/v;
  t += ( W - this.M_Lon(t,3 ) )/v; v=this.M_v(t);  //v的精度0.5%，詳見原文
  t += ( W - this.M_Lon(t,20) )/v;
  t += ( W - this.M_Lon(t,-1) )/v;
  return t;
},
MS_aLon_t:function(W){ //已知月日視黃經差求時間
  var t,v= 7771.37714500204;
  t  = ( W + 1.08472               )/v;
  t += ( W - this.MS_aLon(t, 3, 3) )/v; v=this.M_v(t)-this.E_v(t);  //v的精度0.5%，詳見原文
  t += ( W - this.MS_aLon(t,20,10) )/v;
  t += ( W - this.MS_aLon(t,-1,60) )/v;
  return t;
},
S_aLon_t:function(W){ //已知太陽視黃經反求時間
  var t,v= 628.3319653318;
  t  = ( W - 1.75347-Math.PI   )/v; v=this.E_v(t); //v的精度0.03%，詳見原文
  t += ( W - this.S_aLon(t,10) )/v; v=this.E_v(t); //再算一次v有助于提高精度,不算也可以
  t += ( W - this.S_aLon(t,-1) )/v;
  return t;
},
/****
MS_aLon_t1:function(W){ //已知月日視黃經差求時間,高速低精度,誤差不超過40秒
  var t,v = 7771.37714500204;
  t  = ( W + 1.08472               )/v;
  t += ( W - this.MS_aLon(t, 3, 3) )/v;  v=this.M_v(t)-this.E_v(t);  //v的精度0.5%，詳見原文
  t += ( W - this.MS_aLon(t,50,20) )/v;
  return t;
},
S_aLon_t1:function(W){ //已知太陽視黃經反求時間,高速低精度,最大誤差不超過50秒,平均誤差15秒
  var t,v= 628.3319653318;
  t  = ( W - 1.75347-Math.PI   )/v; v = 628.332 + 21*Math.sin( 1.527+628.307585*t );
  t += ( W - this.S_aLon(t,3) )/v;
  t += ( W - this.S_aLon(t,40))/v;
  return t;
},
****/
MS_aLon_t2:function(W){ //已知月日視黃經差求時間,高速低精度,誤差不超過600秒
  var t,v = 7771.37714500204;
  t  = ( W + 1.08472 )/v;
  var L,t2 = t*t;
  t -= ( -0.00003309*t2 + 0.10976*Math.cos( 0.784758 + 8328.6914246*t + 0.000152292*t2 ) + 0.02224 *Math.cos( 0.18740  + 7214.0628654*t - 0.00021848 *t2 ) - 0.03342 *Math.cos( 4.669257 + 628.307585*t ) )/v;
  L = this.M_Lon(t,20) - (4.8950632+ 628.3319653318*t + 0.000005297*t*t + 0.0334166*Math.cos(4.669257+628.307585*t) + 0.0002061*Math.cos(2.67823+628.307585*t)*t + 0.000349*Math.cos(4.6261+1256.61517*t) - 20.5/rad);
  v = 7771.38 - 914*Math.sin( 0.7848 + 8328.691425*t + 0.0001523*t*t ) - 179*Math.sin( 2.543 + 15542.7543*t ) - 160*Math.sin( 0.1874 + 7214.0629*t );
  t += ( W - L )/v;
  return t;
},
S_aLon_t2:function(W){ //已知太陽視黃經反求時間,高速低精度,最大誤差不超過600秒
  var t,L,v= 628.3319653318;
  t =  ( W - 1.75347-Math.PI )/v;
  t -= ( 0.000005297*t*t + 0.0334166 * Math.cos( 4.669257 + 628.307585*t) + 0.0002061 * Math.cos( 2.67823  + 628.307585*t)*t )/v;
  t += ( W - this.E_Lon(t,8) - Math.PI + (20.5+17.2*Math.sin(2.1824-33.75705*t))/rad )/v;
  return t;
},
moonIll:function(t){ //月亮被照亮部分的比例
 var t2 = t*t, t3 =t2*t, t4 = t3*t;
 var D,M,m,a, dm=Math.PI/180;
 D = 297.8502042 + 445267.1115168*t - 0.0016300*t2 + t3/545868 - t4/113065000; //日月距角
 M = 357.5291092 +  35999.0502909*t - 0.0001536*t2 + t3/24490000;              //太陽平近點
 m = 134.9634114 + 477198.8676313*t + 0.0089970*t2 + t3/69699 - t4/14712000;   //月亮平近點
 D*=dm, M*=dm, m*=dm;
 a = Math.PI - D + (-6.289*Math.sin(m) +2.100*Math.sin(M) -1.274*Math.sin(D*2-m) -0.658*Math.sin(D*2) -0.214*Math.sin(m*2) -0.110*Math.sin(D))*dm;
 return (1+Math.cos(a))/2;
},
moonRad:function(r,h){ //轉入地平緯度及地月質心距離,返回站心視半徑(角秒)
 return 358473400/r*(1+Math.sin(h)*cs_rEar/r);
},
shiCha:function(t){ //時差計算(高精度)
  var t2=t*t,t3=t2*t,t4=t3*t,t5=t4*t;
  var L = ( 1753469512 + 628331965331.8*t + 5296.74*t2 + 0.432*t3 - 0.1124*t4 - 0.00009*t5 + 630 * Math.cos(6+0.3*t) )/1000000000 + Math.PI -20.5/rad;

  var E,dE,dL,f,z=new Array();
  dL= -17.2*Math.sin(2.1824-33.75705*t)/rad; //黃經章
  dE=   9.2*Math.cos(2.1824-33.75705*t)/rad; //交角章
  E = ZB.hcjj(t)+dE; //真黃赤交角

  //地球坐標
  z[0]=this.E_Lon(t,50)+Math.PI+ZB.gxc_sunLon(t)+dL;
  z[1]=-( 2796*Math.cos(3.1987+8433.46616*t) +1016*Math.cos(5.4225+550.75532*t) +804*Math.cos(3.88  +522.3694*t) )/1000000000;

  ZB.llrConv( z, E ); //z太陽地心赤道坐標
  z[0]-=dL*Math.cos(E);

  L = rad2mrad(L-z[0]);
  if(L>Math.PI) L-=pi2;
  return L/pi2; //單位是周(天)
},
shiCha2:function(t){ //時差計算(低精度),誤差約在1秒以內
  var L = ( 1753469512 + 628331965331.8*t + 5296.74*t*t )/1000000000 + Math.PI;
  var z=new Array();
  var E= (84381.4088 -46.836051*t)/rad;
  z[0]=this.E_Lon(t,5)+Math.PI, z[1]=0; //地球坐標
  ZB.llrConv( z, E ); //z太陽地心赤道坐標
  L = rad2mrad(L-z[0]);
  if(L>Math.PI) L-=pi2;
  return L/pi2; //單位是周(天)
}

};