(49g) (50g) Saison begin and Julian Day conversion

[Gil]

Thanks to John Keith and DavidM, here is a shorter version of SEASON (no more IF... THEN DROP END).

In bold the change — at the very end of the program.

\<< \-> y
\<< y -1000 < y 3000 > OR "Not very accurate" 1 IFTE DEG y 1000. /
[[ 485 ]
[ 203 ]
[ 199 ]
[ 182 ]
[ 156 ]
[ 136 ]
[ 77 ]
[ 74 ]
[ 70 ]
[ 58 ]
[ 52 ]
[ 50 ]
[ 45 ]
[ 44 ]
[ 29 ]
[ 18 ]
[ 17 ]
[ 16 ]
[ 14 ]
[ 12 ]
[ 12 ]
[ 12 ]
[ 9 ]
[ 8 ]]
[[ 324.96 ]
[ 337.23 ]
[ 342.08 ]
[ 27.85 ]
[ 73.14 ]
[ 171.52 ]
[ 222.54 ]
[ 296.72 ]
[ 243.58 ]
[ 119.81 ]
[ 297.17 ]
[ 21.02 ]
[ 247.54 ]
[ 325.15 ]
[ 60.93 ]
[ 155.12 ]
[ 288.79 ]
[ 198.04 ]
[ 199.76 ]
[ 95.39 ]
[ 287.11 ]
[ 320.81 ]
[ 227.73 ]
[ 15.45 ]]
[[ 1934.136 ]
[ 32964.467 ]
[ 20.186 ]
[ 445267.112 ]
[ 45036.886 ]
[ 22518.443 ]
[ 65928.934 ]
[ 3034.906 ]
[ 9037.513 ]
[ 33718.147 ]
[ 150.678 ]
[ 2281.226 ]
[ 29929.562 ]
[ 31555.956 ]
[ 4443.417 ]
[ 67555.328 ]
[ 4562.452 ]
[ 62894.029 ]
[ 31436.921 ]
[ 14577.848 ]
[ 31931.756 ]
[ 34777.259 ]
[ 1222.114 ]
[ 16859.074 ]] 0 { } { } { } { } \-> Y A B C JDE0 T W \GD\Gl S
\<<
IF y 1000 \<=
THEN '1721139.29189+365242.1374*Y+.06134*Y^2+.00111*Y^3-.00071*Y^4' EVAL '1721233.25401+365241.72562*Y-.05323*Y^2+.00907*Y^3+.00025*Y^4' EVAL '1721325.70455+365242.49558*Y-.11677*Y^2-.00297*Y^3+.00074*Y^4' EVAL '1721414.39987+365242.88257*Y-.00769*Y^2-.00933*Y^3-.00006*Y^4' EVAL
ELSE Y 2 - 'Y' STO '2451623.80984+365242.37404*Y+.05169*Y^2-.00411*Y^3-.00057*Y^4' EVAL '2451716.56767+365241.62603*Y+.00325*Y^2+.00888*Y^3-.0003*Y^4' EVAL '2451810.21715+365242.01767*Y-.11575*Y^2+.00337*Y^3+.00078*Y^4' EVAL '2451900.05952+365242.74049*Y-.06223*Y^2-.00823*Y^3+.00032*Y^4' EVAL
END 4 \->LIST 'JDE0' STO 1 4
FOR i T JDE0 i GET 2451545 - 36525 / + 'T' STO W T i GET 35999.373 * 2.47 - + 'W' STO \GD\Gl 1 .0334 W i GET COS * + .0007 W i GET 2 * COS * + + '\GD\Gl' STO S 0 1 24
FOR j A j GET B j GET C j GET T i GET * + COS * +
NEXT + 'S' STO JDE0 i GET S i GET .00001 * \GD\Gl i GET / + DUP IP SWAP FP JD\-> 4 ROLL 4 ROLL DROP2 2 \->LIST
NEXT
\>>
\>> 5 ROLL DUP 1 == DROPN
\>>
'SEASON' STO

Just added here the new full Directory EPHEM.V2.DOC, with the SEASON file/program corrected.

Regards,
Gil