summaryrefslogtreecommitdiffstats
path: root/org
diff options
context:
space:
mode:
Diffstat (limited to 'org')
-rw-r--r--org/madore/ephem/Test.java10
1 files changed, 9 insertions, 1 deletions
diff --git a/org/madore/ephem/Test.java b/org/madore/ephem/Test.java
index 0551837..cb9dc82 100644
--- a/org/madore/ephem/Test.java
+++ b/org/madore/ephem/Test.java
@@ -5,7 +5,11 @@ import java.util.Date;
public final class Test {
public static void main(String[] args) {
- Date now = new Date();
+ Date now;
+ if ( args.length >= 1 )
+ now = new Date(Long.parseLong(args[0]));
+ else
+ now = new Date();
long nowTime = now.getTime();
int nowUtcMjd = (int)(nowTime/86400000) + 40587;
double nowUtcSeconds = (nowTime%86400000)/1000.;
@@ -32,6 +36,10 @@ public final class Test {
Math.atan(sun.v[2]/Math.sqrt(sun.v[0]*sun.v[0]+sun.v[1]*sun.v[1]))/Comput.degree,
Math.atan2(sun.v[1], sun.v[0])/Comput.degree);
sun = Precession.matrix(Ephem.fromJd(nowTtJd)).apply(sun0);
+ sun0 = Frames.Rotation.rotx(Precession.epsilon.v(Ephem.fromJd(nowTtJd))*Comput.arcsecond).apply(sun);
+ System.out.format("☉ ecliptic(mean) lat=%.6f° long=%.6f°\n",
+ Math.atan(sun0.v[2]/Math.sqrt(sun0.v[0]*sun0.v[0]+sun0.v[1]*sun0.v[1]))/Comput.degree,
+ Math.atan2(sun0.v[1], sun0.v[0])/Comput.degree);
System.out.format("☉ equatorial(mean) lat=%.6f° long=%.6f°\n",
Math.atan(sun.v[2]/Math.sqrt(sun.v[0]*sun.v[0]+sun.v[1]*sun.v[1]))/Comput.degree,
Math.atan2(sun.v[1], sun.v[0])/Comput.degree);