１）時間引数

　時間引数Tは、位置を求めたい日時のユリウス日JDから、次の式で求めます。

　　　T＝(JD−２４５１５４５)/３６５２５　（J２０００.０）

　この時間引数Tを用いて、位置を求めたい日時の軌道要素の値を計算しておきます。

　＜参考＞　ＰＯＶ−Ｒａｙにおける太陽の軌道要素の式
　//****************************************************************************
　　#local SML=mod(280.466449+0.98564736*36525*T+3.036e-4*pow(T,2),360);
　　#local SPL=mod(282.937348+0.470783e-4*36525*T+4.597e-4*pow(T,2),360);
　　#local SEC=0.01670862-4.204e-5*T-1.24e-7*pow(T,2);
　　#local SAX=1.00000102;
　//****************************************************************************
　上から順番に、平均黄経SML、近日点黄経SPL、軌道離心率SEC、軌道長半径SAXです。