HourAngle Fixed, and sun pos right

1 files changed, 14 insertions, 6 deletions

diff --git a/src/calendar.h b/src/calendar.h
index 810fa0c..248026f 100644
--- a/src/calendar.h
+++ b/src/calendar.h
@@ -10,8 +10,8 @@
 using namespace std;
 
 // 哈尔滨南岗地方经纬度
-#define LONGITUDE 126.68 * M_PI / 180
-#define LATITUDE 45.75 * M_PI / 180
+#define LONGITUDE 126.533 * M_PI / 180
+#define LATITUDE 45.800 * M_PI / 180
 
 extern double delta;
 extern char jieqi[25][10];
@@ -45,15 +45,15 @@ struct point {
     }
 };
 
+void radToDMS(const char *message, int lowerBound, int upperBound,
+              double radians);
+
 class Julian {
   private:
     static double d[23][5];
 
     double dt_ext(int y, double jsd);
 
-    //计算力学时与世界时之差,传入年份
-    double delta_t(int y);
-
   public:
     // 计算儒略日
     double getJulianDay(time_t time);
@@ -63,6 +63,9 @@ class Julian {
 
     //儒略千年数转时间戳
     time_t kiloYearToTime(double t, int year);
+
+    //计算力学时与世界时之差,传入年份
+    double delta_t(int y);
 };
 
 class parameter {
@@ -95,8 +98,13 @@ class parameter {
 
     double moon_longitude(double t);
 
-    double getHourAngle(double julianKiloTime, double longitude,
+    double get_epsilon(double T);
+
+    // 计算时角，返回弧度制
+    double getHourAngle(double julianKiloYear, double longitude,
                         double latitude, double alpha);
 };
 
+extern Julian julian;
+extern parameter p;
 #endif
\ No newline at end of file