Input time, output sun longitude and latitude

1 files changed, 20 insertions, 24 deletions

diff --git a/src/parameters.cpp b/src/parameters.cpp
index 57d98c2..2e2758a 100644
--- a/src/parameters.cpp
+++ b/src/parameters.cpp
@@ -1,5 +1,6 @@
-#include "parameters.h"
 #include "calendar.h"
+#include "parameters.h"
+#include <utility>
 using namespace std;
 
 //凡无特殊说明的函数，全部返回弧度制
@@ -11,7 +12,7 @@ double parameter::get_longitude(vector<double> l, double t) {
         L += l[i] * pow(t, i);
     }
     //返回弧度制，太阳地心黄经=地球日心黄经+180°
-    return L + pi;
+    return L + M_PI;
 }
 
 // 计算地球日心黄纬
@@ -36,19 +37,19 @@ double parameter::get_R(vector<double> r, double t) {
 // 转换FK5误差
 double parameter::delta_FK5(double L, double B, double T) {
     // L转角度制
-    L *= 180 / pi;
+    L *= 180 / M_PI;
     //计算L'，角度制
     double L_p = L - 1.397 * T - 0.00031 * T * T;
     //计算L'，弧度制
-    L_p *= pi / 180;
+    L_p *= M_PI / 180;
     while (L_p < 0)
-        L_p += 2 * pi;
-    while (L_p >= 2 * pi)
-        L_p -= 2 * pi;
+        L_p += 2 * M_PI;
+    while (L_p >= 2 * M_PI)
+        L_p -= 2 * M_PI;
     //计算delta_L，单位为角秒
     double delta_L = -0.09033 + 0.03916 * (cos(L_p) + sin(L_p)) * tan(B);
     //转换为弧度制
-    delta_L *= pi / (180 * 3600);
+    delta_L *= M_PI / (180 * 3600);
     return delta_L;
 }
 
@@ -88,14 +89,14 @@ double parameter::nutation(double T) {
                nutation_parameters[i][2] * ang[2] +
                nutation_parameters[i][3] * ang[3] +
                nutation_parameters[i][4] * ang[4];
-        sita *= pi / 180;
+        sita *= M_PI / 180;
         s += (nutation_parameters[i][5] + nutation_parameters[i][6] * T) *
              sin(sita);
         // c+=(nutation_parameters[i][7]+nutation_parameters[i][8]*T)*cos(sita);
     }
     // 计算得到的s和c单位0.0001秒，转换为弧度制
     s *= 0.0001;
-    s *= pi / (180 * 3600);
+    s *= M_PI / (180 * 3600);
     return s;
 }
 
@@ -104,7 +105,7 @@ double parameter::aberration(double R) {
     //公式，单位0.0001角秒
     double delta = -20.4898 / R;
     //转换为弧度制
-    delta *= pi / (180 * 3600);
+    delta *= M_PI / (180 * 3600);
     return delta;
 }
 
@@ -157,7 +158,7 @@ vector<vector<double>> parameter::get_parameters(double t) {
 }
 
 // 计算太阳地心黄经，外部调用，返回角度制
-double parameter::sun_longitude(double t) {
+pair<double, double> parameter::sun_longitude(double t) {
     //从文件中计算地日运行参数
     vector<vector<double>> p = get_parameters(t);
 
@@ -177,13 +178,8 @@ double parameter::sun_longitude(double t) {
     //光行差修正，返回弧度制
     L += aberration(R);
 
-    // L转角度制
-    L *= 180 / pi;
-    while (L < 0)
-        L += 360;
-    while (L >= 360)
-        L -= 360;
-    return L;
+    // 返回太阳地心黄经和黄纬
+    return make_pair(L, -B);
 }
 
 // 计算月球地心黄经，外部调用，返回角度制
@@ -202,7 +198,7 @@ double parameter::moon_longitude(double t) {
                moon_longitude_parameters[i][1] * angles[1] +
                moon_longitude_parameters[i][2] * angles[2] +
                moon_longitude_parameters[i][3] * angles[3];
-        sita *= pi / 180;
+        sita *= M_PI / 180;
         EI += moon_longitude_parameters[i][4] * sin(sita) *
               pow(E, fabs(moon_longitude_parameters[i][1]));
     }
@@ -212,12 +208,12 @@ double parameter::moon_longitude(double t) {
     a1 = 119.75 + 131.849 * T;
     a2 = 53.09 + 479264.290 * T;
 
-    EI += 3958.0 * sin(a1 * pi / 180);
-    EI += 1962.0 * sin((L - angles[3]) * pi / 180);
-    EI += 318.0 * sin(a2 * pi / 180);
+    EI += 3958.0 * sin(a1 * M_PI / 180);
+    EI += 1962.0 * sin((L - angles[3]) * M_PI / 180);
+    EI += 318.0 * sin(a2 * M_PI / 180);
 
     L += EI / 1000000;
-    L += 180 / pi * nutation(T); //地球章动修正
+    L += 180 / M_PI * nutation(T); //地球章动修正
 
     while (L < 0)
         L += 360;