C++ Solar eclipse program(Shoushi integrated) 1644 - 1785 from Ideone( Date: August 14, 2014 ), C++ (gcc)

C++ Solar eclipse program(Shoushi integrated) 1644 - 1785 from Ideone(...

Language: Layout:

// C++ Solar eclipse program(Shoushi integrated) input as start_year and end_year
    // author : Mr Lam/ Yee Din email : yeedlam@yahoo.com
    // Copyright reserved if you want to use the program, please contact the author
    // August 7, 2012
     
    #include <iostream>
    using namespace std;
    #include <cmath> // or math.h
     
    const int start_year = 1644; // enter start_year
    const int end_year = 1785; // enter end_year
     
    const double Tropical_year = 365.2425;
    const double half_year = Tropical_year/2;
    const int Base_year = 1281;
    const double pi = 3.14159265;
    //===========================================================
    // 4 Paramters for adjustments
    const double Qiying_P = 55.06, Runying_P = 20.2050;
    const double Jiaoying_P = 26.0388, Zhuanying_P = 13.0205;
    //===========================================================
    // Constants for Moon to determine solar eclipse purpose
    const double Shuoshi_C = 29.530593, Zhuanzhong_C = 27.5546;
    const double Jiaozhong_C = 27.212224;
    const double Jiao_eclipse_C = 2.318369;
    const double half_Zhuanzhong = Zhuanzhong_C/2;
    //===================================================
    // Constant for Moon positions : Lead or Lag
    const double Moon_Ref_C = 1.975993;
    const double MoonMeanSpeed = 13.36875;
    //===================================================
     
    typedef struct{
    int year; int month; int day;
    double eclipseFrac; double eclipseBegin; double eclipseMax; double eclipseEnd;
    double eclipseBeginS; double eclipseEndS;
    double Winter_solitice;
    int SunStatus; int MoonStatus;
    int SunStatusS; int MoonStatusS;
    double observed_time; double half_day_hours;
    double ReqSun; double ReqMoon;
    double half_dayoff;
    double Gon;
     
    double EclipseFracS;
    double MoonSpeedFactor;
    double noon_offset;
     
    double MDayT_hour;
    double observed_timeS;
    double SunXS; double MoonXS; double SunMoonXS;
    } EclipseTime;
     
    EclipseTime eclipse_shoushi[100];
     
    EclipseTime eclipse_theo[]
    = {
    {1644, 9, 1, 0.303, 10.941, 12.152, 13.354},
    {1648, 6, 21, 0.982, 6.169, 7.213, 8.374},
    {1650, 10, 25, 0.904, 10.267, 11.603, 12.991},
    {1657, 6, 12, 0.611, 4.301, 5.143, 6.049},
    {1658, 6, 11, 0.53, 9.084, 10.501, 12.049},
    {1665, 1, 16, 0.908, 14.822, 16.267, 17.518},
    {1666, 7, 2, 0.988, 15.563, 16.764, 17.835},
    {1669, 4, 30, 0.538, 13.117, 14.386, 15.578},
    {1671, 9, 3, 0.122, 16.426, 16.968, 17.491},
    {1676, 6, 11, 0.092382, 18.3401, 18.779, 19.207}, //eclipseFrac & eclipseBegin from www.eclipse.org.uk
     
    {1681, 9, 12, 0.383, 8.201, 9.225, 10.333},
    {1685, 11, 26, 0.264, 14.663, 15.478, 16.247},
    {1688, 4, 30, 0.96, 8.182, 9.291, 10.495},
    {1690, 9, 3, 0.216, 6.787, 7.476, 8.206},
    {1691, 2, 28, 0.375, 12.289, 13.526, 14.681},
    {1692, 2, 17, 0.532, 11.276, 12.787, 14.223},
    {1695, 12, 6, 0.973, 15.347, 16.518, 17.570},
    {1697, 4, 21, 0.996, 7.795, 8.899, 10.109},
    {1704, 11, 27, 0.577, 12.241, 13.624, 14.938},
    {1706, 5, 12, 0.574, 18.261, 19.074, 19.839},
     
    {1708, 9, 14, 0.541, 16.541, 17.465, 18.328},
    {1709, 9, 4, 0.523, 5.859, 6.681, 7.56},
    {1712, 7, 4, 0.626, 3.784, 4.684, 5.659},
    {1715, 5, 3, 0.604, 18.146, 18.955, 19.715},
    {1719, 2, 19, 0.741, 14.586, 15.969, 17.186},
    {1720, 8, 4, 0.736, 10.691, 12.225, 13.717},
    {1721, 7, 24, 0.39, 17.33, 18.202, 19.012},
    {1730, 7, 15, 0.829, 11.097, 12.85, 14.604},
    {1731, 12, 29, 0.867, 6.3417, 7.409, 8.602},// eclipseBegin from eclipse.org.uk
    {1735, 10, 16, 0.825, 7.696, 8.869, 10.168},
     
    {1742, 6, 3, 0.702, 6.593, 7.512, 8.511},
    {1745, 4, 2, 0.082, 10.852, 11.496, 12.143},
    {1746, 3, 22, 0.765, 9.315, 10.832, 12.43},
    {1747, 8, 16, 0.249, 17.001, 17.7, 18.361},
    {1751, 5, 25, 0.477, 6.714, 7.543, 8.442},
    {1758, 12, 30, 0.9, 15.005, 16.265, 17.3861}, // eclipseEnd from eclipse.org.uk
    {1760, 6, 13, 0.985, 16.312, 17.343, 18.284},
    {1762, 10, 17, 0.756, 16.68, 9999, 9999},// eclipse.org.uk has no data on Beijing
    {1763, 10, 7, 0.734, 6.0403, 6.933, 7.891},// eclipseBegin from eclipse.org.uk
    {1769, 6, 4, 0.333, 16.984, 17.686, 18.348},
     
    {1770, 5, 25, 0.414, 7.52, 8.394, 9.339},
    {1773, 3, 23, 0.415, 13.015, 14.39, 15.644},
    {1774, 9, 16, 0.387, 7.259, 8.251, 9.337},
    {1775, 8, 26, 0.468, 11.405, 12.923, 14.334},
    {1776, 1, 21, 0.145, 9.526, 10.197, 10.889},
    {1784, 8, 16, 0.187, 5.574, 6.27, 7.019},
    {1785, 8, 5, 0.459, 6.634, 7.691, 8.861},
     
    {99999, 0, 0, 0, 0, 0, 0}
    };
     
    enum Status{LagStart, LagEnd, LeadStart, LeadEnd};
    struct Status_type{char *sun; char *moon;};
     
    class Ref_data
    {
    private:
    Status SunStatus, MoonStatus;
     
    protected:
    int calculated_year; int total_months;
     
    double total_days; double Winter_solitice; double LeapDay;
    double RefEcl; double RefSun; double RefMoonDay; double RefMoon;
    double ReqEcl; double ReqSun; double ReqMoonDay; double ReqMoon;
    double sun_days; int intSun_days;
    double sun_add_minus; double sun_sum; double SunX;
    double moon_add_minus; double moon_sum; double moon_difference; double solar_days;
    double SolarEclipse_angle; int Front;
     
    public:
    Ref_data(int m,int n)
    {
    calculated_year = m;total_months = n;
     
    RefEcl = 0; RefSun = 0; RefMoonDay = 0; RefMoon = 0;
    ReqEcl = 0; ReqSun = 0; ReqMoonDay = 0; ReqMoon = 0;
     
    total_days = 0; LeapDay = 0; Winter_solitice = 0; sun_days = 0;
     
    sun_add_minus = 0; sun_sum = 0; SunX = 0;
    }
    ~Ref_data(){}
     
    void get_basic(double& total_days, double& Winter_solitice, double& LeapDay)
    {
    total_days = (calculated_year - Base_year)* Tropical_year;
    Winter_solitice = fmod((total_days + Qiying_P + int(abs(total_days))* 60), 60);
    LeapDay = fmod((total_days + Runying_P + int(abs(total_days)) * Shuoshi_C), Shuoshi_C);
    }
     
    void get_ref(double total_days, double Winter_solitice, double LeapDay,
    double& BaseSun,double& RefEcl,double& RefSun,double& RefMoonDay,double& RefMoon,double& MinSun)
    {
    RefMoonDay = Winter_solitice - LeapDay;
    if (RefMoonDay < 0){RefMoonDay = RefMoonDay + 60;}
     
    BaseSun = Tropical_year/2 - LeapDay;
     
    RefMoon = fmod((total_days + Zhuanying_P - LeapDay + int(abs(total_days))* Zhuanzhong_C),Zhuanzhong_C);
     
    RefEcl = fmod((total_days + Jiaoying_P - LeapDay + int(abs(total_days))*Jiaozhong_C), Jiaozhong_C);
     
    //--------------------------------------------------------------------------------------
    // Shoushi drift problems referring to base year 1628 from Adam Schall 200 eclipse table
    // calculated data from http://ideone.com/BlkMI
    // for MoonDay = -(0.005499/347) => approx -0.0000158473
    // for sun = -(5.99089/347) => -0.0172648
    // for moon = (0.162147/347) => 0.000467281
    // for eclipse = (0.0122486/347) => 0.0000352987
    //--------------------------------------------------------------------------------------
    // slightly modified according to data record in Chongzhen Lishu(page 788-34) as :
     
    RefMoonDay = -0.0081 + RefMoonDay + 0.01041 - 0.01041/347*(calculated_year - 1281);
    if (RefMoonDay < 0){RefMoonDay = RefMoonDay + 60;}
     
    MinSun = 6.572 - 6.1419 + 6.1419/347 *(calculated_year - 1281);
    RefSun = BaseSun - MinSun;
     
    RefMoon = 0.0402 + RefMoon + 0.16309/347 *(calculated_year - 1281);
    RefEcl = RefEcl + 0.01735/347*(calculated_year - 1281);
    }
     
    void get_required(double RefEcl, double RefSun, double RefMoonDay, double RefMoon,
    double& ReqEcl, double& ReqSun, double& ReqMoonDay, double& ReqMoon)
    {
    ReqEcl = fmod(((Jiao_eclipse_C * total_months)+ RefEcl),Jiaozhong_C);
    ReqSun = RefSun + (Shuoshi_C * total_months);
    ReqSun = fmod(ReqSun, Tropical_year);
    ReqMoonDay = fmod((RefMoonDay + Shuoshi_C * total_months),60);
    ReqMoon = fmod(( RefMoon + Moon_Ref_C * total_months),Zhuanzhong_C);
    }
     
    void checkSunStatus(double ReqSun, double& sun_days, int& SunStatus)
    {
    SunStatus = LagStart; sun_days = 0;
     
    if (ReqSun < half_year)
    {
    SunStatus = LagStart; sun_days = ReqSun;
    if (ReqSun > 93.712025){sun_days = half_year - sun_days; SunStatus = LagEnd;}
    }
     
    if (ReqSun > half_year )
    {
    SunStatus = LeadStart; sun_days = ReqSun - half_year ;
    if (sun_days > 88.909225){sun_days = half_year - sun_days; SunStatus = LeadEnd;}
    }
    }
     
    void calculateSunOffset(int SunStatus, double sun_days, int intSun_days,
    double& sun_add_minus, double& sun_sum, double& SunX)
    {
    sun_add_minus = sun_sum = SunX = 0;
     
    if (fmod(SunStatus, 3.0) == 0) // SunStaus = LagStart & LeadEnd
    {
    sun_add_minus = 484.8473 -(4.4362 + 0.0162 *(intSun_days - 1)/2)* intSun_days;
    sun_sum = (487.06 -(2.21 + 0.0027 * intSun_days)* intSun_days )* intSun_days;
    }
     
    if (fmod(SunStatus, 3.0) > 0) // SunStatus = LagEnd & LeadStart
    {
    sun_add_minus = 510.8569 -(4.9386 + 0.0186 *(intSun_days - 1)/2)*intSun_days;
    sun_sum = (513.32 -(2.46 + 0.0031 * intSun_days)*intSun_days)*intSun_days;
    }
    SunX =((sun_days - intSun_days)*sun_add_minus + sun_sum)/10000;
    }
     
    // Calculate Solar eclipse angle
     
    void get_SolarEclipse_angle(int SunStatus,double ReqEcl_angle,double SunX,double& SolarEclipse_angle,int& Front)
    {
    if (SunStatus < LeadStart){SolarEclipse_angle = ReqEcl_angle - abs(SunX);}
    if (SunStatus > LagEnd){SolarEclipse_angle = ReqEcl_angle + abs(SunX);}
     
    if (SolarEclipse_angle < 7){Front = 1;}
    if (SolarEclipse_angle > 342){Front = 1;}
     
    if(SolarEclipse_angle > 202)
    {
    if (SolarEclipse_angle < 342){
    // cout << endl;
    // cout << "==============================" << endl;
    // cout << " no solar eclipse observed" << endl;
    // cout << "==============================" << endl;
    // cout << endl;
    }
    }
     
    if(SolarEclipse_angle < 175)
    {
    if (SolarEclipse_angle > 7){
    // cout << endl;
    // cout << "=================================" << endl;
    // cout << " no solar eclipse observed" << endl;
    // cout << "=================================" << endl;
    // cout << endl;
    }
    }
     
    if (SolarEclipse_angle > 175)
    {
    if ( SolarEclipse_angle < 202){Front = -1;}
    }
    // cout << endl;
    // cout << "=================" << endl;
    // cout << "Front = " << Front << endl;
    // cout << "=================" << endl;
    // cout << endl;
    }
     
    void get_calculatedMoon_data(double ReqMoon, double& CalculatedMoon, double& solar_days,
    double& moon_x12_2, double& moon_minus_solar_days,
    int& MoonStatus, double& MoonLeadLag, double& moon_add_minus)
    {
    double moon_84_check;
    if (ReqMoon > half_Zhuanzhong){CalculatedMoon = ReqMoon - half_Zhuanzhong;MoonStatus = LagStart;}
    if (ReqMoon < half_Zhuanzhong){CalculatedMoon = ReqMoon;MoonStatus = LeadStart;}
     
    moon_x12_2 = CalculatedMoon * 12.2;
    solar_days = int(moon_x12_2) * 0.082008;
     
    moon_84_check = moon_x12_2;
    if (moon_x12_2 > 84)
    {
    if (MoonStatus == LeadStart){MoonStatus = LeadEnd;}
    if (MoonStatus == LagStart){MoonStatus = LagEnd;}
     
    moon_84_check = 167 - moon_84_check;
    }
    MoonLeadLag = moon_84_check;
     
    moon_minus_solar_days = CalculatedMoon - solar_days;
     
    if( moon_minus_solar_days < 0 )
    {
    solar_days = solar_days - 0.082008;
    moon_minus_solar_days = CalculatedMoon - solar_days;
     
    if (fmod(MoonStatus, 2.0) == 0){MoonLeadLag = int(MoonLeadLag);}
    if (fmod(MoonStatus, 2.0) > 0){MoonLeadLag = 167 - int(moon_x12_2);}
    }
    moon_add_minus = 11.081575 -(0.05815 + 0.00195 *(MoonLeadLag-1)/2)* MoonLeadLag;
    }
     
     
    void getMoon_adjusted_data(int MoonStatus, double MoonLeadLag, double moon_add_minus,
    double moon_minus_solar_days,
    double& moon_sum, double& MoonX, double& MoonSpeedFactor)
    {
    moon_sum = 0, MoonX = 0, MoonSpeedFactor = 0;
     
    switch(MoonStatus)
    {
    case LagStart:
    MoonSpeedFactor = 0.9855 +(0.05815 + 0.00195 *(MoonLeadLag-1)/2)* MoonLeadLag/100;
    moon_sum =(0.1111 -(0.000281 + 0.00000325 *MoonLeadLag)*MoonLeadLag)*MoonLeadLag;
    MoonX = moon_sum + moon_minus_solar_days * moon_add_minus/820;
     
    break;
     
    case LagEnd:
    MoonSpeedFactor = 1.2071 -(0.05815 + 0.00195 *(MoonLeadLag-1)/2)*MoonLeadLag/100;
    MoonLeadLag = MoonLeadLag + 1;
    moon_sum =(0.1111 -(0.000281 + 0.00000325 * MoonLeadLag)*MoonLeadLag)*MoonLeadLag;
    MoonX = moon_sum - moon_minus_solar_days * moon_add_minus/820;
     
    break;
     
    case LeadStart:
    MoonSpeedFactor = 1.2071 -(0.05815 + 0.00195 *(MoonLeadLag-1)/2)*MoonLeadLag/100;
    moon_sum =(0.1111 -(0.000281 + 0.00000325 *MoonLeadLag)*MoonLeadLag)*MoonLeadLag;
    MoonX = moon_sum + moon_minus_solar_days * moon_add_minus/820;
     
    break;
     
    case LeadEnd:
    MoonSpeedFactor = 0.9855 +(0.05815 + 0.00195*(MoonLeadLag-1)/2)*MoonLeadLag/100;
    MoonLeadLag = MoonLeadLag + 1;
    moon_sum =(0.1111 -(0.000281 + 0.00000325 * MoonLeadLag)*MoonLeadLag)*MoonLeadLag;
    MoonX = moon_sum - moon_minus_solar_days * moon_add_minus/820;
     
    break;
    default:
    cout << " no MoonStatus";
     
    break;
    }
    }
     
    };
     
    class base2 {
    protected:
    int k;
    public:
    base2(int x) {k=x;}
    ~base2() {}
    };
     
    class derived: public Ref_data, public base2 {int j;
    public:
    derived(int m, int n, int x, int z): Ref_data(m,n), base2(x){j=z;}
    ~derived() {}
     
    void show() {cout << calculated_year << " " << k << " " << j << "\n";}
    };
     
    int main()
    {
    cout << "\tt" << "\tt" << "\tt" << "\tt" << "\tt" << "\tt" << "\tt" << "\tt" << "\tt" << "\tt" << "\tt" << "\tt" << "\tt"
    << "\tt" << "\tt" << "\tt" << endl;
     
    cout << "\tSolar eclipse time is in 24 hours\n"
    << "\tBeijing local noon time is 12 as reference\n"
    << "\tmonth means total months from winter solitice\n"
    << "\tReqMDay=ReqMoonDay\n"
    << "year(month)\n";
     
    int index = 0;int Pindex = 0;int calculated_year, total_months;
     
    for (calculated_year = start_year; calculated_year - 1< end_year; calculated_year ++)
    {
    for (total_months = 0; total_months < 14; total_months ++)
    {
    Ref_data year(calculated_year,total_months);
     
    double total_days, Winter_solitice, LeapDay;
    year.get_basic(total_days, Winter_solitice, LeapDay);
     
    double BaseSun, RefEcl, RefSun, RefMoonDay, RefMoon, MinSun;
    year.get_ref(total_days,Winter_solitice,LeapDay,BaseSun,RefEcl,RefSun,RefMoonDay,RefMoon,MinSun);
     
    double ReqEcl, ReqEcl_angle, ReqSun, ReqMoonDay, ReqMoon;
    year.get_required(RefEcl,RefSun,RefMoonDay,RefMoon,ReqEcl,ReqSun,ReqMoonDay,ReqMoon);
     
    ReqEcl_angle = ReqEcl * MoonMeanSpeed;
     
    double ReqSun_days = ReqSun;
    if (ReqSun > half_year){ReqSun_days = ReqSun - half_year;};
     
    double sun_days;int SunStatus, intSun_days;
    year.checkSunStatus(ReqSun,sun_days, SunStatus);
     
    intSun_days = int(sun_days);
     
    double sun_add_minus, sun_sum, SunX;
    year.calculateSunOffset(SunStatus, sun_days, intSun_days, sun_add_minus, sun_sum, SunX);
     
    //-----------------
    double ReqSun360 = fmod(ReqSun * 360/Tropical_year + 180, 360);
    double e = 0.0167;
    double SunX_Th = 180/pi *(2 *e* sin(ReqSun360 *pi/180)+ 1.25 * e*e * sin(2*ReqSun360 *pi/180));
    double SFactor = MoonMeanSpeed * 0.082 + 0.082; // = 1.1782375
     
    double SunX360 = SunX_Th * SFactor;
    SunX = Tropical_year/360* SunX360;
    //-----------------
    double MinSun360 = MinSun * 360/Tropical_year;
    double RSunYel = ReqSun360 + SunX_Th + MinSun360 + 0.098;
     
    // 求正交平均
    double PosExtra = -570 *SunX_Th/6973;
    double PosX360 = PosExtra * SFactor;
    double PosX = PosX360 * Tropical_year/360;
     
    //^^^^^^^^^^^^^^(Convert Shoushi ReqEcl as SunPos)^^^^^^^^^^
    double SunPos = (ReqEcl - Jiaozhong_C/2)*360/Jiaozhong_C - PosExtra + SunX_Th;
    if (SunPos < 0){SunPos = SunPos + 360;}
    //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
     
    double UsePos = RSunYel - SunPos;
     
    // 求正交實均
    double SunPosReal = 180/pi *atan(56.0/59 *tan(SunPos * pi/180));
    if ((SunPos > 90)&(SunPos < 270)){SunPosReal = SunPosReal + 180;}
    if (SunPos > 270){SunPosReal = SunPosReal + 360;}
     
    // 正交實均
    double PosRealAdd = SunPos - SunPosReal;
     
    //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
    double PosRX360 = PosRealAdd * SFactor;
    double PosRX = PosRX360 * Tropical_year/360;
    //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
     
    // 求正交實行
    double PosReal = UsePos + PosRealAdd;
    double MoonPos = RSunYel - PosReal; // same as MoonPos = SunPos - PosRealAdd
    if (MoonPos < 0){MoonPos = MoonPos + 360;}
    //----------
    double YelWh = 5 + 17.0/60 + 20.0/3600;
    double YMoonPos = 180/pi * atan(cos(YelWh * pi/180)* tan(MoonPos * pi/180));
    if ((MoonPos > 90)&(MoonPos < 270)){YMoonPos = YMoonPos + 180;}
     
    double YelD = YMoonPos - MoonPos;
    if(YelD < -350){YelD = YelD + 360;}
    //------
    double SolarEclipse_angle;int Front = 0;
    year.get_SolarEclipse_angle(SunStatus,ReqEcl_angle,SunX, SolarEclipse_angle, Front);
    if (SolarEclipse_angle < 7){SolarEclipse_angle = SolarEclipse_angle + 363.7934196;}
     
    double SolEA = SolarEclipse_angle - PosX - PosRX;
    // if(NS_f < 0){SolEA = SolarEclipse_angle - NS_Hu/60;}
     
    double SolEA360 = fmod(SolEA + 363.7934196/2, 363.7934196)* 360/363.7934196;
     
    double LatYel = 180/pi *asin(sin(YelWh *pi/180)*sin(MoonPos *pi/180));
    LatYel = 180/pi *asin(sin(YelWh *pi/180)*sin(SolEA360 *pi/180));
     
    //RRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRR
    // Convert Shoushi ReqMoon to consider Ee_o
    //---------------------------------------------------------------
    double OReqMoon = ReqMoon;
    double OReqMoon360 = fmod(OReqMoon * 360/Zhuanzhong_C + 180, 360);
    if (OReqMoon360 < 0){OReqMoon360 = fmod(360 + OReqMoon360, 360);}
    //------------------------------
    // apply Earth Center Offset for Moon adjustment
    double Ee_o = 0;
    if (cos(OReqMoon360 * pi/180)< 0){Ee_o = 0;}
    double y_scale = 1;
    double x_cir = cos(OReqMoon360 * pi/180);
    double y_cir = sin(OReqMoon360 * pi/180);
    double r_elip = 1/sqrt(x_cir*x_cir + y_cir*y_cir/(y_scale*y_scale));
    double x_o = r_elip * x_cir;
    double y_o = r_elip * y_cir;
     
    double EMoon360 = 180/pi *atan(y_o/(x_o - Ee_o));
    if ((OReqMoon360 > 90) & (OReqMoon360 < 270)){EMoon360 = 180 + EMoon360;}
    if (OReqMoon360 > 270){EMoon360 = 360 + EMoon360;}
    if ((OReqMoon360 < 90) & (EMoon360 < 0)){EMoon360 = 180 + EMoon360;}
    //------------------------------
    double MOff = EMoon360 - OReqMoon360;
    double EMoon = fmod((EMoon360 + 180)*Zhuanzhong_C/360, Zhuanzhong_C);
    double MF = 0;
    ReqMoon = EMoon;
    double ReqMoon360 = fmod(ReqMoon * 360/Zhuanzhong_C + 180, 360);
    //RRRRRRRRRRRRRRRRRRRRRRRRRRRRRR
     
    int MoonStatus = 0;double CalculatedMoon = 0, MoonLeadLag, moon_add_minus, solar_days;
    double moon_x12_2, moon_minus_solar_days;
    year.get_calculatedMoon_data(ReqMoon,CalculatedMoon,solar_days,moon_x12_2, moon_minus_solar_days,
    MoonStatus, MoonLeadLag, moon_add_minus);
     
    double moon_sum, MoonX, MoonSpeedFactor;
    year.getMoon_adjusted_data(MoonStatus, MoonLeadLag, moon_add_minus, moon_minus_solar_days,
    moon_sum, MoonX, MoonSpeedFactor);
    //--------------
    double SunXS = -abs(SunX);
    if (SunStatus > LagEnd){SunXS = abs(SunX);}
     
    double MoonXS = -abs(MoonX);
    if (MoonStatus > LagEnd){MoonXS = abs(MoonX);}
     
    double SunMoonX = SunXS - MoonXS;
    double SunMoonT = SunMoonX * 0.082/MoonSpeedFactor;
    //--------------
    double RSunEnd = RSunYel + SunMoonT * 360/Tropical_year;
     
    double MoonT24 = SunMoonT * 24;
     
    //=========================================
    // Use HouBian steps to simulate MidMoon
    //=========================================
     
    const double Qiying_HB = 32.12254;
    double SumDay_HB = (calculated_year - 1723)* 365.24233442;
    double Winter_HB = fmod((SumDay_HB + Qiying_HB +(int)(abs(SumDay_HB))*60),60);
    double UseDays = SumDay_HB +(Qiying_HB -(int)Qiying_HB)-(Winter_HB -(int)Winter_HB);
     
    double DayMoon = fmod(UseDays* 47435.0234086, 1296000)/3600;
    double MoonYearBase = DayMoon + 176 + 27.0/60 + 48.0/3600 + 53.0/21600;
     
    const double Shuoying_HB = 15.12633;
    const double Shuoshi_HB = 29.53059053;
    double days_for_moon = UseDays - Shuoying_HB;
    double Remainder_moon = fmod(days_for_moon, Shuoshi_HB);
    if (Remainder_moon < 0){Remainder_moon = Shuoshi_HB + Remainder_moon;}
     
    double SumMonth_HB = total_months - 1;
    if (SumMonth_HB < 0){SumMonth_HB = 0;}
    if (calculated_year == 1719){SumMonth_HB = SumMonth_HB - 1;}
    if (calculated_year == 1643){SumMonth_HB = SumMonth_HB - 1;}
     
    double Next_moon_days = Shuoshi_HB - Remainder_moon;
    double MoonDay = (Next_moon_days + SumMonth_HB * Shuoshi_HB)*47435.0234086/3600;
    double MidMoon = fmod(MoonYearBase + MoonDay, 360);
     
    //=========================================
    // Use HouBian steps to simulate MaxMoon
    //=========================================
    // 求積日最高平行
    double BaseMaxMoon = fmod(UseDays * 401.070226, 1296000)/3600;
    // 求最高年根
    double MaxYearBase = BaseMaxMoon +(241 + 15.0/60 + 45.0/3600 + 38.0/216000);
    // 求最高日數
    double MaxDay = (Next_moon_days + SumMonth_HB * Shuoshi_HB)* 401.070226/3600;
    // 求最高平行
    double MaxMoon =fmod(MaxYearBase + MaxDay, 360);
    if(MaxMoon < 0){MaxMoon = 360 + MaxMoon;}
    //-----------------------------------------
    double ReqMoon_HB = MidMoon - MaxMoon;
    if (ReqMoon_HB < 0){ReqMoon_HB = 360 + ReqMoon_HB;}
    //-----------------------------------------
     
    double MoonX360 = MoonXS * 360/Tropical_year;
    double MoonX_Th = MoonX360 + SunX_Th - SunX360;
    double UseRMoon = MidMoon + MoonX_Th;
    double RMoonEnd = fmod(UseRMoon + SunMoonT *(MoonSpeedFactor/0.082 + 1)*360/Tropical_year, 360);
    //------------
    double SunMoonXS = SunMoonX;
    int SunStatusS = SunStatus;
    int MoonStatusS = MoonStatus;
    //------------
     
    double MDay0 = ReqMoonDay - (int)ReqMoonDay;
    double MDay0_24 = MDay0 * 24;
     
    double MDayT = MDay0 + SunMoonT;
    if(MDayT < 0){MDayT = 1 + MDayT;}
    if(MDayT > 1){MDayT = MDayT - 1;}
    double MDayT_day = int(MDayT);
    double MDayT_hour = MDayT - int(MDayT);
    double MDayT_24hour = MDayT_hour * 24;
     
    //****************************
    // Calculate half day hours
     
    double sun_rise_time
    = 2068.3 +(((0.122928 -((intSun_days-53)* 0.000608))*(intSun_days-1))/2)*intSun_days + 0.06 * intSun_days;
     
    double sun_set_time = 10000 - sun_rise_time;
     
    double half_day_hours = sun_set_time - 5000;
    if (SunStatus == LagEnd){half_day_hours = 5000 - half_day_hours;}
    if (SunStatus == LeadStart){half_day_hours = 5000 - half_day_hours;}
    //*****************
    //----------------------
    double time_from_noon = (MDayT - int(MDayT)) - 0.5;
    double observed_timeS = (0.5 - abs(time_from_noon))* (time_from_noon)/0.96;
    //----------------------
     
    int SeeTimeTable[][19]
    = {
    //h_deg-120 -105 -90 -75 -60 -45 -30 -15 0 15 30 45 60 75 90 105 120 <--
    //Gon h20 h19 h18 h17 h16 h15 h14 h13 h12 h11 h10 h9 h8 h7 h6 h5 h4 Gon
    //--------------------------------------------------------------------------------------------
    {9, -61, -72, -81, -85, -85, -78, -61, -34, -0, 34, 61, 78, 85, 85, 81, 72, 61, 9},
    {10, -51, -61, -68, -70, -66, -54, -32, -2, 33, 63, 84, 95, 99, 96, 90, 80, 67, 8},
    {11, -40, -49, -55, -55, -50, -37, -14, 15, 48, 75, 93, 103, 106, 103, 96, 86, 74, 7},
    {0, -33, -42, -48, -50, -46, -35, -15, 13, 44, 72, 92, 103, 108, 106, 100, 91, 80, 6},
    {1, -34, -44, -52, -55, -54, -45, -28, -3, 29, 60, 84, 98, 105, 105, 100, 92, 81, 5},
    {2, -45, -56, -64, -69, -69, -62, -46, -20, 13, 46, 73, 90, 97, 98, 94, 86, 75, 4},
    {3, -61, -72, -81, -85, -85, -78, -61, -34, 0, 34, 61, 78, 85, 85, 81, 72, 61, 3},
    //---------------------------------------------------------------------------------------------
    //Gon h4 h5 h6 h7 h8 h9 h10 h11 h12 h13 h14 h15 h16 h17 h18 h19 h20 Gon
    };
    //-->
     
    double degree = fmod(RSunEnd + 270, 360);
    double Gon = degree/30.0;
    double Hour = MDayT_hour * 24;
     
    int LRow = 9;
    int RRow = 9;
    if (fmod(Gon + 3,12.0) < 6){LRow = (int)(fmod((int)Gon + 3, 12.0));}
    if (fmod(Gon + 3,12.0) > 5){RRow = (int)(fmod(9 - (int)Gon, 12.0));}
     
    int TableTime = 0;
    int LHour = 0;
    double t0 = 0;
    double t1 = 0;
    double SeeT = 0;
     
    if (LRow < 6)
    {
    LHour = (int)Hour - 3;
    TableTime = SeeTimeTable[LRow][LHour];
    t0 = TableTime +(SeeTimeTable[LRow][LHour+1] - TableTime)*(Hour -(int)Hour);
     
    TableTime = SeeTimeTable[LRow+1][LHour];
    t1 = TableTime +(SeeTimeTable[LRow+1][LHour+1] - TableTime)*(Hour -(int)Hour);
    }
     
    int RHour = 0;
    if (LRow > 5)
    {
    if (RRow < 7)
    {
    RHour = 21 - (int)Hour;
    TableTime = SeeTimeTable[RRow][RHour];
    t0 = TableTime +(SeeTimeTable[RRow][RHour-1] - TableTime)*(Hour -(int)Hour);
    t0 = -t0;
     
    TableTime = SeeTimeTable[RRow-1][RHour];
    t1 = TableTime +(SeeTimeTable[RRow-1][RHour-1] - TableTime)*(Hour -(int)Hour);
    t1 = -t1;
    }
    }
    SeeT = t0 + (t1 - t0)*(Gon -(int)Gon);
    double SeeHour = SeeT/60;
    double observed_time = SeeHour/24;
    double eclipseMax_hour = MDayT_hour + observed_time;
     
    //-------------------
    double eclipseMax_24hour = eclipseMax_hour * 24;
     
    // cout << "===============================================" << endl;
    // cout << " eclipseMax_24hour = " << eclipseMax_24hour << endl;
    // cout << "===============================================" << endl;
    // cout << endl;
    //********************
     
    double time_difference_from_noon = time_from_noon + observed_time;
    time_difference_from_noon = abs(time_difference_from_noon)* 10000;
     
    double MaxPhase_lead_lag = ReqSun + int(MDayT) + eclipseMax_hour;
     
    double MaxPhaseSun_ref = MaxPhase_lead_lag - ReqMoonDay;
     
    double MaxPhaseSun_days;
    year.checkSunStatus(MaxPhaseSun_ref, MaxPhaseSun_days, SunStatus);
     
    intSun_days = int(MaxPhaseSun_days);
     
    double MaxPhaseSun_add_minus, MaxPhaseSun_sum, MaxPhaseSunX;
    year.calculateSunOffset(SunStatus, MaxPhaseSun_days, intSun_days,MaxPhaseSun_add_minus, MaxPhaseSun_sum, MaxPhaseSunX);
     
    //*************
    double MaxPhaseSun_adjust = 0;
    if (SunStatus < LeadStart){MaxPhaseSun_adjust = MaxPhaseSun_ref - MaxPhaseSunX;}
    if (SunStatus > LagEnd){MaxPhaseSun_adjust = MaxPhaseSun_ref + MaxPhaseSunX;}
     
    //*************
    double MaxPhaseSun_NS;
    year.checkSunStatus(MaxPhaseSun_adjust, MaxPhaseSun_NS, SunStatus);
     
    double NS_difference = 4.46 -((MaxPhaseSun_NS * MaxPhaseSun_NS)/1870);
    double NS_Sh = NS_difference -(NS_difference * time_difference_from_noon)/half_day_hours;
     
    double MaxPhaseSun_EW = MaxPhaseSun_adjust;
    if(MaxPhaseSun_EW > half_year){MaxPhaseSun_EW = MaxPhaseSun_EW - half_year;}
     
    double EW_difference = ((half_year - MaxPhaseSun_EW)* MaxPhaseSun_EW)/1870;
     
    double EW_level = (EW_difference * time_difference_from_noon)/2500;
    double EW_Sh = EW_level;
    if(EW_level > EW_difference ){EW_Sh = 2 * EW_difference - EW_level;}
     
    int AfterNoonSign;
    AfterNoonSign = 1;
    if (time_from_noon < 0){AfterNoonSign = -1;}
     
    switch(SunStatus)
    {
    case LagStart:
    NS_Sh = Front * NS_Sh;
    EW_Sh = -AfterNoonSign * Front * EW_Sh;
    break;
     
    case LagEnd:
    NS_Sh = -Front * NS_Sh;
    EW_Sh = -AfterNoonSign * Front * EW_Sh;
    break;
     
    case LeadStart:
    NS_Sh = -Front * NS_Sh;
    EW_Sh = AfterNoonSign * Front * EW_Sh;
    break;
     
    case LeadEnd:
    NS_Sh = Front * NS_Sh;
    EW_Sh = AfterNoonSign * Front * EW_Sh;
    break;
    default:
    cout << " no SunStatus" << endl;
    }
    //---------------------------
    double eclipse_ref_angle = 0;
    if (Front == -1){eclipse_ref_angle = 188.05 + NS_Sh + EW_Sh;}
    if (Front == 1){eclipse_ref_angle = 357.64 + NS_Sh + EW_Sh;}
     
    double postiveEclipse_angle = SolarEclipse_angle - eclipse_ref_angle;
    double SexEcl = Front * postiveEclipse_angle;
     
    double EclipseFracS = 0;
    double EclipseRef = 6;
    if (SexEcl < 0){EclipseRef = 8;}
     
    EclipseFracS = (EclipseRef - abs(postiveEclipse_angle))/EclipseRef;
    //**************
     
    int Lat_LongTable[][19]
    = {
    //h_deg-120 -105 -90 -75 -60 -45 -30 -15 0 15 30 45 60 75 90 105 120 <--
    //Gon h20 h19 h18 h17 h16 h15 h14 h13 h12 h11 h10 h9 h8 h7 h6 h5 h4 Gon
    //--------------------------------------------------------------------------------------------
    {9, 10, 15, 20, 25, 27, 30, 35, 40, 42, 40, 35, 30, 27, 25, 20, 15, 10, 9},//Lat(min)
    {9, 0, 2, 0, 47, 47, 41, 28, 26, 25, 27, 28, 41, 47, 47, 0, 2, 1, 9},//Lat(sec)
    {9, 30, 35, 40, 41, 42, 39, 31, 17, 0, 17, 31, 39, 42, 41, 40, 35, 30, 9},//Long(min)
    {9, 0, 0, 0, 20, 17, 30, 58, 6, 0, 6, 58, 30, 17, 21, 0, 0, 0, 9},//Long(sec)
     
    {10, 32, 33, 34, 35, 36, 40, 41, 41, 38, 32, 26, 20, 16, 14, 14, 13, 12, 8},
    {10, 0, 0, 40, 33, 38, 26, 29, 24, 25, 36, 40, 51, 48, 39, 0, 0, 0, 8},
    {10, 0, 0, 0, 35, 33, 30, 19, 3, 13, 30, 40, 46, 46, 46, 47, 48, 49, 8},
    {10, 0, 0, 0, 7, 27, 8, 49, 35, 10, 43, 37, 30, 42, 37, 0, 0, 0, 8},
     
    {11, 36, 37, 38, 39, 40, 41, 39, 36, 29, 21, 16, 11, 9, 8, 7, 6, 5, 7},
    {11, 0, 0, 26, 37, 39, 29, 34, 30, 29, 44, 44, 53, 57, 48, 51, 0, 0, 7},
    {11, 32, 31, 30, 30, 29, 23, 12, 1, 21, 35, 43, 47, 49, 49, 48, 47, 46, 7},
    {11, 0, 0, 13, 14, 16, 22, 29, 27, 15, 2, 1, 50, 44, 45, 47, 0, 0, 7},
     
    {0, 39, 40, 41, 41, 40, 38, 35, 28, 21, 13, 9, 7, 6, 5, 5, 5, 4, 6},
    {0, 0, 0, 27, 28, 30, 31, 25, 36, 41, 56, 53, 55, 54, 57, 53, 0, 0, 6},
    {0, 28, 29, 30, 30, 28, 21, 13, 0, 15, 32, 41, 47, 49, 49, 48, 48, 47, 6},
    {0, 0, 0, 3, 5, 4, 11, 13, 36, 32, 7, 57, 51, 44, 47, 47, 0, 0, 6},
     
    {1, 38, 39, 41, 40, 38, 33, 27, 21, 13, 8, 7, 6, 8, 8, 9, 10, 11, 5},
    {1, 0, 47, 26, 28, 21, 30, 36, 40, 55, 46, 52, 56, 48, 52, 27, 40, 0, 5},
    {1, 26, 27, 28, 30, 31, 27, 19, 7, 10, 26, 37, 45, 49, 49, 48, 48, 47, 5},
    {1, 0, 33, 26, 15, 9, 7, 11, 13, 23, 11, 5, 49, 46, 57, 53, 13, 0, 5},
     
    {2, 35, 35, 36, 34, 30, 25, 18, 10, 9, 6, 7, 10, 13, 16, 17, 16, 17, 4},
    {2, 17, 7, 33, 31, 30, 39, 54, 18, 50, 57, 52, 48, 47, 10, 11, 46, 0, 4},
    {2, 27, 28, 29, 31, 34, 34, 27, 13, 3, 18, 34, 43, 46, 47, 46, 45, 45, 4},
    {2, 41, 53, 52, 54, 36, 12, 22, 29, 20, 20, 38, 22, 27, 27, 34, 34, 0, 4},
     
    {3, 26, 27, 27, 26, 21, 17, 13, 8, 6, 8, 13, 17, 21, 26, 27, 27, 26, 3},
    {3, 40, 15, 55, 42, 48, 45, 45, 52, 56, 52, 45, 45, 48, 42, 55, 15, 40, 3},
    {3, 39, 40, 41, 42, 43, 39, 31, 18, 0, 18, 31, 39, 43, 42, 41, 40, 39, 3},
    {3, 13, 20, 27, 29, 35, 29, 37, 21, 0, 38, 37, 29, 25, 39, 27, 20, 13, 3},
     
    //---------------------------------------------------------------------------------------------
    //Gon h4 h5 h6 h7 h8 h9 h10 h11 h12 h13 h14 h15 h16 h17 h18 h19 h20 Gon
    //-->
    };
     
    enum TRow{Lat_min, Lat_sec, Long_min, Long_sec};
     
    LRow = 9;
    RRow = 9;
    if (fmod(Gon + 3,12.0) < 6){LRow = (int)(fmod((int)Gon + 3, 12.0));}
    if (fmod(Gon + 3,12.0) > 5){RRow = (int)(fmod(9 - (int)Gon, 12.0));}
     
    double NS_0 = 0;
    double NS_1 = 0;
    double NS_01 = 0;
    double NS_N0 = 0;
    double NS_N1 = 0;
    double NS_N01 = 0;
     
    double EW_0 = 0;
    double EW_1 = 0;
    double EW_01 = 0;
    double EW_N0 = 0;
    double EW_N1 = 0;
    double EW_N01 = 0;
     
    LHour = 0;
    double NS_Hu = 0;
    double EW_Hu = 0;
     
    if (LRow < 6)
    {
    LHour = (int)Hour - 3;
    NS_0 = Lat_LongTable[LRow *4+Lat_min][LHour] + Lat_LongTable[LRow*4+Lat_sec][LHour]/60.0;
    NS_1 = Lat_LongTable[LRow *4+Lat_min][LHour+1] + Lat_LongTable[LRow*4+Lat_sec][LHour+1]/60.0;
    NS_01 = NS_0 + (NS_1 - NS_0)*(Hour -(int)Hour);
     
    NS_N0 = Lat_LongTable[(LRow+1)*4+Lat_min][LHour] + Lat_LongTable[(LRow+1)*4+Lat_sec][LHour]/60.0;
    NS_N1 = Lat_LongTable[(LRow+1)*4+Lat_min][LHour+1] + Lat_LongTable[(LRow+1)*4+Lat_sec][LHour+1]/60.0;
    NS_N01 = NS_N0 + (NS_N1 - NS_N0)*(Hour -(int)Hour);
    //--------------
    EW_0 = Lat_LongTable[LRow *4+Long_min][LHour] + Lat_LongTable[LRow*4+Long_sec][LHour]/60.0;
    EW_1 = Lat_LongTable[LRow *4+Long_min][LHour+1] + Lat_LongTable[LRow*4+Long_sec][LHour+1]/60.0;
    EW_01 = EW_0 + (EW_1 - EW_0)*(Hour -(int)Hour);
     
    EW_N0 = Lat_LongTable[(LRow+1)*4+Long_min][LHour] + Lat_LongTable[(LRow+1)*4+Long_sec][LHour]/60.0;
    EW_N1 = Lat_LongTable[(LRow+1)*4+Long_min][LHour+1] + Lat_LongTable[(LRow+1)*4+Long_sec][LHour+1]/60.0;
    EW_N01 = EW_N0 + (EW_N1 - EW_N0)*(Hour -(int)Hour);
    }
     
    RHour = 0;
    if (LRow > 5)
    {
    if (RRow < 7)
    {
    RHour = 21 - (int)Hour;
    NS_0 = Lat_LongTable[RRow *4+Lat_min][RHour] + Lat_LongTable[RRow*4+Lat_sec][RHour]/60.0;
    NS_1 = Lat_LongTable[RRow *4+Lat_min][RHour-1] + Lat_LongTable[RRow*4+Lat_sec][RHour-1]/60.0;
    NS_01 = NS_0 + (NS_1 - NS_0)*(Hour -(int)Hour);
     
    NS_N0 = Lat_LongTable[(RRow-1)*4+Lat_min][RHour] + Lat_LongTable[(RRow-1)*4+Lat_sec][RHour]/60.0;
    NS_N1 = Lat_LongTable[(RRow-1)*4+Lat_min][RHour-1] + Lat_LongTable[(RRow-1)*4+Lat_sec][RHour-1]/60.0;
    NS_N01 = NS_N0 + (NS_N1 - NS_N0)*(Hour -(int)Hour);
    //---------------------------------
    EW_0 = Lat_LongTable[RRow *4+Long_min][RHour] + Lat_LongTable[RRow*4+Long_sec][RHour]/60.0;
    EW_1 = Lat_LongTable[RRow *4+Long_min][RHour-1] + Lat_LongTable[RRow*4+Long_sec][RHour-1]/60.0;
    EW_01 = EW_0 + (EW_1 - EW_0)*(Hour -(int)Hour);
     
    EW_N0 = Lat_LongTable[(RRow-1)*4+Long_min][RHour] + Lat_LongTable[(RRow-1)*4+Long_sec][RHour]/60.0;
    EW_N1 = Lat_LongTable[(RRow-1)*4+Long_min][RHour-1] + Lat_LongTable[(RRow-1)*4+Long_sec][RHour+1]/60.0;
    EW_N01 = EW_N0 + (EW_N1 - EW_N0)*(Hour -(int)Hour);
    }
    }
    NS_Hu = NS_01 + (NS_N01 - NS_01)*(Gon -(int)Gon);
    EW_Hu = EW_01 + (EW_N01 - EW_01)*(Gon -(int)Gon);
    //=============
     
    // 求最高平均
    double MaxExtra = 1196 * SunX_Th/6973;
    // 求用最高
    double UseMax = MaxMoon + MaxExtra;
     
    // 求日距月最高
    double MaxSMoon = RSunYel - UseMax;
    if (MaxSMoon < 0){MaxSMoon = MaxSMoon + 360;}
     
    //mmmmmmmmmm
    // 所夾之角
    double UseAngle = 180 - fmod(MaxSMoon *2, 360);
    // 最高實均
    double MaxRealExtra
    = 180/pi *atan((117315 *sin(UseAngle *pi/180))/(550505 - 117315 *cos(UseAngle *pi/180)));
    //mmmmmmmmmm
    //==========
    // From HouBian
    // ------------
    // 求本天心距地數(本時两心差)
    double C_Earth = 117315.0/1000000 *sin(MaxSMoon *2*pi/180)/sin(MaxRealExtra *pi/180);
    //=========
    // From HouBian Step 3 :
    // 推地平高下差及日月視徑第三
     
    // 求太陰距地
    double DistCos = cos(ReqMoon360 * pi/180)* 2 * C_Earth;
    double DistSin = sin(ReqMoon360 * pi/180)* 2 * C_Earth;
    double AddDistCos = 20 - DistCos;
    double DistRef = DistSin * DistSin/AddDistCos;
    double CalRef = (DistRef + AddDistCos)/2;
     
    // 月距地心數(太陰距地)
    double MoonToEarth = 20 - CalRef;
     
    //求日距地心數
    // 太陽實引
    double ReqRealSun = fmod(ReqSun360 + SunX360, 360);
    // 分股
    double SideCos = cos(ReqRealSun * pi/180)* 0.338;
    // 勾
    double SideSin = sin(ReqRealSun * pi/180)* 0.338;
    // 勾弦和
    double AddSideCos = 20 + SideCos;
    // 勾弦較
    double SideRef = SideSin * SideSin/AddSideCos;
    // 弦
    double BaseRef =(SideRef + AddSideCos)/2;
    // 日距地心數
    double SunToEarth = 20 - BaseRef;
     
    //=======
    // From HouBian Step 2 :
     
    // 一小時太陰白道實行
    double HourWhiteMoon = MoonSpeedFactor/0.082 + 1; // simulate HouBian
    double HourRSunYel = 1; // simulate HouBian
     
    // 求斜距交角差
    double XA = 180/pi*atan(HourRSunYel*sin(YelWh*pi/180)/(HourWhiteMoon-HourRSunYel*cos(YelWh*pi/180)));
     
    // 求斜距黃道交角
    double XYelA = YelWh + XA;
    // 求食甚實緯
    double MaxLat = LatYel * cos(XYelA *pi/180);
    //^^^^^^^^^^
    double MaxMoonX =(3450.0 * 10/MoonToEarth)/60;
    double HighLow = MaxMoonX - 10.0/60;
     
    double SunSeeSemiDiameter =(966 * 10/SunToEarth)/60;
    double SunRSemiDiameter = SunSeeSemiDiameter - 15.0/60;
    double MoonSeeSemiDiameter =(940.5 * 10/MoonToEarth)/60;
    double NetDiameter = SunRSemiDiameter + MoonSeeSemiDiameter;
    //^^^^^^^^^
    double MaxSunYel = RSunEnd + observed_time/24;
    double MaxPhaseUseT = eclipseMax_hour;
    //=========
    // Step 4 : 推食甚太陽黃赤經緯宿度及黃赤二經交角第四
    //-----------------
    // 求食甚太陽赤道經度
    double SprYel = MaxSunYel - 90;
    double WinYel = MaxSunYel;
    double MaxSunRed = 180/pi*atan(cos((23+29.0/60)*pi/180)*tan(SprYel*pi/180));
     
    double SprRed = MaxSunRed;
    if ((SprYel > 90)&(SprYel < 270)){SprRed = 180 + MaxSunRed;}
    if (SprYel > 270){SprRed = 360 - MaxSunRed;}
    //-----------------
    // 求食甚太陽赤道緯度
    double ESunLat = 180/pi*asin(sin((23+29.0/60)*pi/180)*sin(SprYel*pi/180));
    // 求太陽距北極
    double SunNPole = 90 - ESunLat;
     
    // 求黃赤二經交角
    double YelRedXA = 90-180/pi*atan(tan((90-(23 + 29.0/60))*pi/180)/cos(SprYel*pi/180));
    if (cos(SprYel*pi/180)< 0){YelRedXA=YelRedXA-180;}
     
    // 求黄白二經交角
    double YelWhXA = XYelA;
    if (cos(MoonPos * pi/180) < 0){YelWhXA = -XYelA;}
     
    // 求赤白二經交角
    double RedWhXA = YelRedXA + YelWhXA;
     
    //=========
    // Step 5 : 又法 page(792-248)推食甚近時第五
     
    // 求用時太陽距午赤道度
    double T_12Red =(MaxPhaseUseT -(int)MaxPhaseUseT - 0.5)* 24 * 15;
    //-----------------
    // 求用時赤經高弧交角
     
    // 北極距天頂 = 90 -(Beijing Lat)
    double NPoleV = 50.08333333;
    // 距極分邊
    double T_12X=180/pi*atan(cos(T_12Red*pi/180)*tan(NPoleV*pi/180));
    // 距日分邊
    double T_12NP = SunNPole - T_12X;
    // 垂弧之正切線
    double T_12Arc = tan(T_12Red*pi/180)*sin(T_12X*pi/180);
    // 用時赤經高弧交角
    double URedXA = 180/pi*atan(T_12Arc/sin(T_12NP*pi/180));
    //---------------
    // 求用時太陽距天頂
    double UseSunHL = 180/pi*asin(sin(NPoleV*pi/180)*sin(T_12Red*pi/180)/sin(URedXA*pi/180));
     
    // 求用時白經高弧交角
    double UWhXA = URedXA + RedWhXA;
    // 求用時高下差
    double UseHL = HighLow*sin(UseSunHL*pi/180);
    // 求用時東西差
    double UseEW = sin(UWhXA*pi/180)*UseHL;
    double EW_f = UseEW/60;
    // 求用時南北差
    double UseNS = -cos(UWhXA*pi/180)*UseHL;
    double NS_f = UseNS/60;
     
    // 求用時視緯
    double UseSeeLat = MaxLat*60 + UseNS;
    //^^^^^^^^^^^^^^^
    double ESeeLat = UseSeeLat/60;
    double SunHalfD = SunRSemiDiameter/60;
    double addD = NetDiameter/60;
     
    double eclipseFrac =(addD - abs(ESeeLat))/(2* SunHalfD);
    //***************
    double RefEA_Sh = 0;
    if (Front == -1){RefEA_Sh = 188.05;}
    if (Front == 1){RefEA_Sh = 357.64;}
    double TestEA = SolarEclipse_angle - RefEA_Sh - PosX - PosRX;
    double NetEA = TestEA -(NS_Sh + EW_Sh);
    double SexEA = NetEA;
    if (Front == -1){SexEA = -1.0 * NetEA;}
     
    double EARef = 6;
    if (SexEA < 0){EARef = 8;}
    EclipseFracS = (EARef - abs(SexEA))/EARef;
    //***************
    // Calculate Eclipse begin time and end time
     
    double AdjMoonSpeedFactor = 0;double AdjMoonLeadLag;int AdjMoonStatus = 0;
     
    double AdjMoon = ReqMoon + SunMoonT;
     
    year.get_calculatedMoon_data(AdjMoon, CalculatedMoon, solar_days, moon_x12_2,
    moon_minus_solar_days, AdjMoonStatus, AdjMoonLeadLag, moon_add_minus);
     
    year.getMoon_adjusted_data(AdjMoonStatus, AdjMoonLeadLag, moon_add_minus,
    moon_minus_solar_days, moon_sum, MoonX, AdjMoonSpeedFactor);
     
    double relativeMoonSun_speed_factor = AdjMoonSpeedFactor - 0.082;
    //---------------------------------
    double relative_time_difference
    = (sqrt((20 - eclipseFrac * 10)* eclipseFrac * 10)* 57.40)/relativeMoonSun_speed_factor;
     
    // Hui Hui method here
    // relative_time_difference
    //= 24*60*sqrt(addD*addD - ESeeLat*ESeeLat)/AdjMoonSpeedFactor;
     
    double eclipseBegin_hour = eclipseMax_hour - relative_time_difference/10000;
    double eclipseBegin_24hour = eclipseBegin_hour * 24;
     
    double eclipseEnd_hour = eclipseMax_hour + relative_time_difference/10000;
    double eclipseEnd_24hour = eclipseEnd_hour * 24;
    //----------------------------------
    double MidEclipse_time
    = (sqrt((20 - EclipseFracS * 10)* EclipseFracS * 10)* 57.40)/relativeMoonSun_speed_factor;
     
    double eclipseBegin_hourS = eclipseMax_hour - MidEclipse_time/10000;
    double eclipseBegin_24hourS = eclipseBegin_hourS * 24;
     
    double eclipseEnd_hourS = eclipseMax_hour + MidEclipse_time/10000;
    double eclipseEnd_24hourS = eclipseEnd_hourS * 24;
    //------------------------------------
    int Print = 0;
    if (eclipseFrac < 0){eclipseFrac = EclipseFracS;}
    if (eclipseFrac > 0)
    {
    if (calculated_year >= eclipse_theo[index].year)
    {
    if (calculated_year > eclipse_theo[index].year)
    {eclipse_shoushi[index].eclipseFrac = -99;index = index + 1;}
     
    if (calculated_year == eclipse_theo[index].year)
    {
    if (total_months >= eclipse_theo[index].month)
    {
    Pindex = index;
    eclipse_shoushi[index].year = calculated_year;
    eclipse_shoushi[index].month = total_months;
     
    eclipse_shoushi[index].eclipseFrac = eclipseFrac;
    eclipse_shoushi[index].MoonSpeedFactor = MoonSpeedFactor;
    eclipse_shoushi[index].EclipseFracS = EclipseFracS;
    eclipse_shoushi[index].SunXS = SunXS;
    eclipse_shoushi[index].MoonXS = MoonXS;
     
    eclipse_shoushi[index].eclipseBegin = eclipseBegin_24hour;
    eclipse_shoushi[index].eclipseBeginS = eclipseBegin_24hourS;
    eclipse_shoushi[index].eclipseMax = eclipseMax_24hour;
    eclipse_shoushi[index].eclipseEnd = eclipseEnd_24hour;
    eclipse_shoushi[index].eclipseEndS = eclipseEnd_24hourS;
     
    eclipse_shoushi[index].Winter_solitice = Winter_solitice;
    eclipse_shoushi[index].SunStatus = SunStatus;
     
    eclipse_shoushi[index].observed_time = observed_time;
    eclipse_shoushi[index].half_day_hours = half_day_hours;
    eclipse_shoushi[index].ReqSun = ReqSun;
    eclipse_shoushi[index].ReqMoon = ReqMoon;
     
    eclipse_shoushi[index].Gon = Gon;
    eclipse_shoushi[index].MoonStatus = MoonStatus;
    eclipse_shoushi[index].SunMoonXS = SunMoonXS;
    eclipse_shoushi[index].MDayT_hour = MDayT_hour;
    eclipse_shoushi[index].observed_timeS = observed_timeS;
     
    eclipse_shoushi[index].SunStatusS = SunStatusS;
    eclipse_shoushi[index].MoonStatusS = MoonStatusS;
     
    index = index + 1;
    Print = 1;
    }
    }
    }
     
    if (Print == 1){//if_Print
    if (calculated_year == eclipse_theo[Pindex].year)
    {
    if (total_months >= eclipse_theo[Pindex].month)
    {
    if (calculated_year <= end_year)
    {
    printf("---------\n");
    printf("%4d(%1d) \t",calculated_year,total_months);
    printf("SprRed\tWinYel\tT_12Red\tYelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA\tHighLow*sin(UseSunHL)=UseHL\n");
    printf("\t\t%3.2f\t%3.2f\t%3.2f\t%3.2f + %3.2f = %3.2f[+]%3.2f=%3.2f\t%3.2f*sin(%3.2f)=%3.2f\n",
    SprRed,WinYel,T_12Red,YelRedXA,YelWhXA,RedWhXA,URedXA,UWhXA,HighLow,UseSunHL,UseHL);
    printf("MidMoon-MaxMoon\tLeapDay ReqMDay\tReqSun()360\tReqMoon()360\tReqEcl");
    printf(" SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd\n");
    printf("%3.2f-%3.2f\t%3.2f\t%3.2f\t%3.2f(%1d)%3.2f",
    MidMoon,MaxMoon,LeapDay,ReqMoonDay, ReqSun,SunStatus,ReqSun360);
    printf("\t%3.2f(%1d)%3.2f\t%3.2f", ReqMoon,MoonStatus,ReqMoon360, ReqEcl);
    printf("\t%3.2f+%3.2f=%3.2f\t%3.2f\t%3.2f/%3.2f\n", SunPos,UsePos,RSunYel,UseRMoon, RSunEnd,RMoonEnd);
    printf("=%3.2f\n",ReqMoon_HB);
    printf("MinSun(360)\tT_12X\tT_12NP\tT_12Arc\tOReqMoon(360)\tSun_Th/SunXS/MoonXS\tEe_o\tEMoon360/MOff/MF\n");
    printf("%3.2f(%3.2f)\t%3.2f\t%3.2f\t%3.2f\t%3.2f(%3.2f)",MinSun,MinSun360,T_12X,T_12NP,T_12Arc,OReqMoon,OReqMoon360);
    printf("\t%3.2f/%3.2f/%3.2f \t%4.3f\t%3.2f/%3.2f/%3.3f\n",SunX_Th,SunXS,MoonXS,Ee_o,EMoon360, MOff,MF);
     
    printf("Eclipse)Frac\tMaxPh(MDayT)\tTestEA\tSexEA\n");
    printf("\t%4.3f\t%4.3f(%4.3f)",eclipseFrac,eclipseMax_24hour, MDayT_24hour);
    if (MoonLeadLag >= 82){cout << "(>)";}
    printf("\t%3.2f\t%3.2f\n",TestEA,SexEA);
    printf("MDay0_24 MoonT24 SeeT(min)\tXA+YelWh=XYelA\tLatYel\tMoonPos/YelD\tSolEA(360)\n");
    printf("%3.3f\t%3.3f\t%3.2f(%3.2f)\t%3.2f+%2.2f=%3.2f\t%3.2f\t%3.2f/%3.3f\t%3.2f(%3.2f)\n",
    MDay0_24,MoonT24,SeeHour,SeeT,XA,YelWh,XYelA,LatYel,MoonPos,YelD,SolEA,SolEA360);
    }
     
    printf("\tNS_Sh/EW_Sh\tNS_Hu/EW_Hu\tESunLat\tEW_f\tMaxLat+NS_f=ESeeLat\taddD\n");
    printf("\t%3.2f/%3.2f\t%3.2f/%3.2f\t%3.2f\t%3.2f\t%3.2f + %3.2f=%3.3f\t%3.3f\n\n",
    NS_Sh,EW_Sh, NS_Hu,EW_Hu,ESunLat, EW_f,MaxLat,NS_f,ESeeLat,addD);
    }
    }
    }//if_Print end
    }// if_1 end
    }
    }
    cout << endl;
    cout << "========================================\n"
    << "\tShoushi calendar method referring to Chonzhen Lishu 200 year eclipse table for drift problem\n"
    << "\t\tdrift rate = offset/347 and 347 = 1628 - 1281(base year for Shoushi)\n"
    << "\t\tobserved_time is integrated with theorectical HuiHui calendar method table\n"
    << "\t\teclipse fraction is integrated from HouBian for comparison with Shoushi\n"
    << endl
     
    << "\tData record in Chongzhen Lishu(page 788-34):\n"
    << "\t\t65 fan late => 0.0065\n"
    << "\t\t6 days late => more than 6\n"
    << "\t\t2 degree late => (2/365.2425)* 27.5546 => more than 0.150883864\n\n"
     
    << "\tShoushi drift from Adam Schall 200 eclipse table for year 1628\n"
    << "\t\tcalculated data from http://ideone.com/BlkMI\n"
    << "\t\tfor moon_days = -(0.005499/347) => approx -0.0000158473\n"
    << "\t\tfor sun = -(5.99089/347) => -0.0172648\n"
    << "\t\tfor moon = (0.162147/347) => 0.000467281\n"
    << "\t\tfor eclipse = (0.0122486/347) => 0.0000352987\n"
    << "-----------------------------\n"
    << endl
    << " slightly modified according to data record in Chongzhen Lishu(page 788-34) as :\n"
    << "RefMoonDay = -0.0081 + RefMoonDay + 0.01041 - 0.01041/347 *(calculated_year - 1281)\n"
    << endl
    << "MinSun = 6.572 - 6.1419 +(6.1419/347)*(calculated_year - 1281)\n"
    << "RefSun = BaseSun - MinSun\n"
    << endl
    << "RefMoon = 0.0402 + RefMoon + 0.16309/347 *(calculated_year - 1281)\n"
    << "RefEcl = RefEcl + 0.01735/347*(calculated_year - 1281)\n"
    << endl
    << "----------------------------\n"
    << "\t\tInverted Kepler elliptic model is implemented as :\n"
    << "\tReqSun360 = fmod(ReqSun * 360/Tropical_year + 180, 360)\n"
    << "\te = 0.0167\n"
    << "\tSunX_Th = 180/pi*(2 *e* sin(ReqSun360 *pi/180)+ 1.25* e*e * sin(2*ReqSun360 *pi/180))\n"
    << "\tSFactor = 13.36875 * 0.082 + 0.082// = 1.1782375\n"
    << "\tSunX360 = SFactor * SunX_Th\n"
    << "\tSunX = Tropical_year/360* SunX360\n"
    << "\t=================================================\n"
    << "\tRSunYel = ReqSun360 + SunX_Th + MinSun360 + 0.098\n"
    << "\t=================================================\n"
    << endl
    << "----------------------\n"
    << "\t\tConvert Shoushi ReqMoon to consider Ee_o\n"
    << "\tOReqMoon = ReqMoon\n"
    << "\tOReqMoon360 = fmod(OReqMoon * 360/Zhuanzhong_C + 180, 360)\n"
    << "\tif (OReqMoon360 < 0){OReqMoon360 = fmod(360 + OReqMoon360, 360);}\n\n"
     
    << "\t\tapply Earth Center Offset for Moon adjustment\n"
    << "\tEe_o = 0\n"
    << "\tif (cos(OReqMoon360 * pi/180)< 0){Ee_o = 0}\n"
    << "\ty_scale = 1\n"
    << "\tx_cir = cos(OReqMoon360 * pi/180)\n"
    << "\ty_cir = sin(OReqMoon360 * pi/180)\n"
    << "\tr_elip = 1/sqrt(x_cir*x_cir + y_cir*y_cir/(y_scale*y_scale))\n"
    << "\tx_o = r_elip * x_cir\n"
    << "\ty_o = r_elip * y_cir\n"
     
    << "\t\tEMoon360 = 180/pi *atan(y_o/(x_o - Ee_o))\n"
    << "\t\tif ((OReqMoon360 > 90)&(OReqMoon360 < 270)){EMoon360 = 180 + EMoon360;}\n"
    << "\t\tif (OReq_moon360 > 270){EMoon360 = 360 + EMoon360;}\n"
    << "----------------------\n"
    << "\tMOff = EMoon360 - OReqMoon360\n"
    << "\tEMoon = fmod((EMoon360 + 180)*Zhuanzhong_C/360, Zhuanzhong_C)\n\n"
    << "\tReqMoon = EMoon\n\n"
     
    << "\tMDay0 = ReqMoonDay -(int)ReqMoonDay\n\n"
     
    << "\t\tHuiHui observed time table is implemented as :\n"
    << "\tdegree = fmod(RSunEnd + 270, 360)\n"
    << "\tGon = degree/30.0\n"
    << "\tHour = MDayT_hour * 24\n"
    << endl
    << "int SeeTimeTable[][19]\n"
    << " = {\n"
    << "h_deg -120 -105 -90 -75 -60 -45 -30 -15 0 15 30 45 60 75 90 105 120 <--\n"
    << "Gon h20 h19 h18 h17 h16 h15 h14 h13 h12 h11 h10 h9 h8 h7 h6 h5 h4 Gon\n"
    << "-----------------------------------------------------------------------------------------------\n"
    << "{9, -61, -72, -81, -85, -85, -78, -61, -34, -0, 34, 61, 78, 85, 85, 81, 72, 61, 9},\n"
    << "{10, -51, -61, -68, -70, -66, -54, -32, -2, 33, 63, 84, 95, 99, 96, 90, 80, 67, 8},\n"
    << "{11, -40, -49, -55, -55, -50, -37, -14, 15, 48, 75, 93, 103, 106, 103, 96, 86, 74, 7},\n"
    << "{0, -33, -42, -48, -50, -46, -35, -15, 13, 44, 72, 92, 103, 108, 106, 100, 91, 80, 6},\n"
    << "{1, -34, -44, -52, -55, -54, -45, -28, -3, 29, 60, 84, 98, 105, 105, 100, 92, 81, 5},\n"
    << "{2, -45, -56, -64, -69, -69, -62, -46, -20, 13, 46, 73, 90, 97, 98, 94, 86, 75, 4},\n"
    << "{3, -61, -72, -81, -85, -85, -78, -61, -34, 0, 34, 61, 78, 85, 85, 81, 72, 61, 3},\n"
    <<"-----------------------------------------------------------------------------------------------\n"
    << "Gon h4 h5 h6 h7 h8 h9 h10 h11 h12 h13 h14 h15 h16 h17 h18 h19 h20 Gon\n"
    << "-->\n" << endl
     
    << "\tobserved_time = SeeT/(24*60)\n"
    << endl
    << "\teclipseMax_hour = MDayT_hour + observed_time\n\n"
    << "\tt\tt\tt\tt\tt\tt\tt\tt\tt\tt\tt\tt\tt\tt" << endl
    << "======================================================================================================================\n"
    << "Shoushi(year)(months)"
    << "\t\teclipseFraction \teclipseBegin\t\tMaxPhase\t\teclipseEnd" << endl;
    cout << "Theorectical(month/day)"
    << "\t\t\t(EclipseFracS)\t(eclipseBeginS)\t(observed_time x 24 x 60)\t(eclipse_endS)" << endl;
    cout << " half_day_hours/ReqSun/ReqMoon"
    << "\t\t\t\t\t\t(observed_timeS x 24 x 60)" << endl
    << " SunXS(SunStatusS)/MoonSpeedFactor" << endl
    << " MoonXS(MoonStatusS)/SunMoonXS" << endl
    << " Gon/MDayT" << endl;
    cout << "___________________________________________________________"
    << "___________________________________________________________" << endl
    << endl;
     
    int record;
    for (record = 0; record < index; record ++)
    {
    printf("%4d(%1d)",eclipse_shoushi[record].year,eclipse_shoushi[record].month);
    printf(" %4.2f",eclipse_shoushi[record].half_day_hours);
    printf("/%2.1f/%2.1f", eclipse_shoushi[record].ReqSun, eclipse_shoushi[record].ReqMoon);
    printf("\t %5.4f(%5.4f)",eclipse_shoushi[record].eclipseFrac, eclipse_shoushi[record].EclipseFracS);
    cout << "\t " << eclipse_shoushi[record].eclipseBegin << "(" << eclipse_shoushi[record].eclipseBeginS << ")"
    << "\t " << eclipse_shoushi[record].eclipseMax << "(" << eclipse_shoushi[record].observed_time *24*60 << ")\t "
    << eclipse_shoushi[record].eclipseEnd << "(" << eclipse_shoushi[record].eclipseEndS << ")\n";
     
    cout << eclipse_theo[record].month << "/" << eclipse_theo[record].day;
    printf("\t%4.3f(%1d)/", eclipse_shoushi[record].SunXS, eclipse_shoushi[record].SunStatusS );
    printf("%4.3f ", eclipse_shoushi[record].MoonSpeedFactor);
     
    printf("\t-)%5.4f", eclipse_theo[record].eclipseFrac);
    printf("\t\t-)%5.4f", eclipse_theo[record].eclipseBegin);
    printf("\t\t-)%5.4f", eclipse_theo[record].eclipseMax);
    printf("(%5.4f)", eclipse_shoushi[record].observed_timeS * 24 * 60);
    printf(" \t-)%6.4f\n", eclipse_theo[record].eclipseEnd);
     
    printf("\t%4.3f(%1d)/", eclipse_shoushi[record].MoonXS, eclipse_shoushi[record].MoonStatusS);
    printf("%4.3f ", eclipse_shoushi[record].SunMoonXS);
    cout << "\t--------------------------------------------------------------------------------------" << endl;
    printf("\t%3.2f/%3.2f", eclipse_shoushi[record].Gon, eclipse_shoushi[record].MDayT_hour * 24);
     
    printf("\t\t %5.4f",eclipse_shoushi[record].eclipseFrac - eclipse_theo[record].eclipseFrac);
    printf("(%5.4f)", eclipse_shoushi[record].EclipseFracS - eclipse_theo[record].eclipseFrac);
    printf("\t %5.4f", eclipse_shoushi[record].eclipseBegin - eclipse_theo[record].eclipseBegin);
    printf("(%5.4f)\t ",eclipse_shoushi[record].eclipseBeginS - eclipse_theo[record].eclipseBegin);
    printf("%5.4f\t\t", eclipse_shoushi[record].eclipseMax - eclipse_theo[record].eclipseMax);
    printf("%5.4f", eclipse_shoushi[record].eclipseEnd - eclipse_theo[record].eclipseEnd);
    printf("(%5.4f)\n\n", eclipse_shoushi[record].eclipseEndS - eclipse_theo[record].eclipseEnd);
    };
     
    cout << endl;
    return 0;
    }
    //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
    // 求正交平均
    // PosExtra = -570 *SunX_Th/6973;
    //======================================
    // Convert Shoushi ReqEcl as SunPos
    //======================================
    // 求日距正交
    // SunPos = (ReqEcl - Jiaozhong_C/2)* 360/Jiaozhong_C - PosExtra + SunX_Th;
    // if (SunPos < 0){SunPos = SunPos + 360;}
     
    // 求用正交
    // UsePos = RSunYel - SunPos;
     
    // MoonPos = RSunYel - PosReal
     
    // LatYel = 180/pi *asin(sin(YelWh *pi/180)*sin(MoonPos *pi/180))
    // Here use LatYel = 180/pi *asin(sin(YelWh *pi/180)*sin(SolEA360 *pi/180));
    //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
    // ESeeLat = UseSeeLat/60;
    // SunHalfD = SunRSemiDiameter/60;
    // addD = NetDiameter/60;
     
    // eclipseFrac =(addD - abs(ESeeLat))/(2* SunHalfD);
    //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
    // PosRX360 = PosRealAdd * SFactor;
    // PosRX = PosRX360 * Tropical_year/360;
     
    // SolEA = SolarEclipse_angle - PosX - PosRX
    // SolEA360 = fmod(SolEA + 363.7934196/2, 363.7934196)*360/363.7934196
    //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
    // EclipseFracS modified to include PosX and PosRX as :
    //--------------------------------------------
    // if (Front == -1){RefEA_Sh = 188.05}
    // if (Front == 1){RefEA_Sh = 357.64}
    // TestEA = SolarEclipse_angle - RefEA_Sh - PosX - PosRX
    // NetEA = TestEA -(NS_Sh + EW_Sh)
    // SexEA = NetEA
    // if (Front == -1){SexEA = -1*NetEA}
     
    // EARef = 6;
    // if (SexEA < 0){EARef = 8}
    // EclipseFracS = (EARef - abs(SexEA))/EARef

Show compiler warnings [ + ] Compiler args [ + ] Show input
Compilation time: 1.23 sec, absolute running time: 0.15 sec, cpu time: 0.01 sec, memory peak: 3 Mb, absolute service time: 1.41 sec

t t t t t t t t t t t t t t t t Solar eclipse time is in 24 hours Beijing local noon time is 12 as reference month means total months from winter solitice ReqMDay=ReqMoonDay year(month) --------- 1644(9) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 160.76 249.17 -1.12 -22.10 + 5.71 = -16.39[+]-1.63=-18.02 55.17*sin(31.79)=29.06 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 251.35-294.98 10.87 52.98 65.43(0)244.49 10.44(3)316.37 13.91 2.20+247.43=249.63 255.43 249.17/249.17 =316.37 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 6.86(6.76) 50.08 31.78 -0.01 10.44(316.37) -1.71/-2.05/3.83 0.000 316.37/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.403 11.925(12.317) -4.41 0.09 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 23.624 -11.306 -0.39(-23.52) 0.42+5.29=5.71 0.16 2.09/-0.009 183.64(1.73) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD -4.33/0.01 18.67/9.76 8.15 -0.15 0.16 + -0.46=-0.302 0.514 --------- 1648(7) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 90.03 180.02 -75.49 -0.01 + 5.68 = 5.67[+]-55.41=-49.74 57.43*sin(64.41)=51.80 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 179.36-89.61 24.84 0.92 357.57(3)172.44 20.65(0)89.76 14.03 5.83+173.77=179.61 173.97 180.03/180.03 =89.76 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 6.93(6.83) 16.67 49.84 -1.11 20.65(89.76) 0.25/0.29/-5.42 0.000 89.76/0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.906 6.967(8.344)(>) -0.56 -0.08 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 22.130 10.214 -1.38(-82.60) 0.39+5.29=5.68 0.51 5.54/-0.023 187.49(5.53) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD -1.26/0.62 20.40/42.07 23.48 -0.66 0.51 + -0.56=-0.051 0.521 --------- 1650(11) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 209.74 301.92 -7.41 -20.24 + -5.67 = -25.91[+]-7.15=-33.06 58.82*sin(52.53)=46.69 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 303.46-185.02 17.06 17.31 118.19(1)296.50 22.84(1)118.44 26.84 173.14+128.58=301.73 299.00 301.92/301.93 =118.44 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 6.96(6.86) 49.85 52.32 -0.10 22.84(118.44) -1.73/-2.07/-4.84 0.000 118.44/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.946 11.506(12.187) -0.71 0.31 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 7.421 4.766 -0.68(-40.84) 0.38+5.29=5.67 0.62 173.49/0.027 356.93(173.21) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD -0.65/-0.36 31.32/17.15 -12.16 -0.42 0.62 + -0.65=-0.031 0.534 --------- 1657(6) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 80.42 171.20 -101.79 3.80 + -5.69 = -1.89[+]-49.05=-50.94 56.80*sin(83.76)=56.46 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 170.21-94.79 4.60 38.82 348.12(3)163.13 19.55(0)75.42 26.46 170.59+0.16=170.75 164.90 171.20/171.20 =75.42 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 7.09(6.98) -13.73 80.54 -1.14 19.55(75.42) 0.54/0.65/-5.30 0.000 75.42/0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.601 5.214(6.518) -2.69 0.45 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 19.623 10.895 -1.30(-78.24) 0.40+5.29=5.69 0.80 171.06/0.037 354.95(171.25) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 0.08/-3.22 29.70/38.74 23.19 -0.73 0.80 + -0.59=0.206 0.519 --------- 1658(6) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 69.03 160.64 -16.30 8.20 + -5.72 = 2.48[+]-34.11=-31.63 54.27*sin(22.57)=20.83 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 159.49-134.27 15.47 33.18 337.23(3)152.39 15.71(0)25.22 27.06 178.99+-18.63=160.36 156.95 160.64/160.64 =25.22 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 7.10(7.00) 48.92 18.99 -0.22 15.71(25.22) 0.87/1.04/-2.42 0.000 25.22/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.544 10.913(11.212) 5.37 3.86 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 4.433 6.778 -0.30(-17.90) 0.43+5.29=5.72 0.07 179.04/0.004 363.01(179.22) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 2.21/-0.70 9.58/12.44 22.08 -0.18 0.07 + -0.30=-0.224 0.508 --------- 1665(1) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL -61.16 26.80 62.74 11.08 + -5.72 = 5.36[+]43.26=48.61 54.00*sin(84.26)=53.73 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 26.25-44.05 3.01 54.69 201.92(1)19.02 12.41(3)342.20 26.69 173.72+-146.85=26.88 27.85 26.80/26.80 =342.20 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 7.23(7.12) 28.70 82.14 0.93 12.41(342.20) 0.64/0.76/1.75 0.000 342.20/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.841 16.183(14.690) 0.31 2.92 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 16.635 -1.945 1.49(89.57) 0.43+5.29=5.72 0.53 174.04/0.025 357.95(174.22) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD -0.63/-1.98 23.67/43.87 -20.84 0.67 0.53 + -0.59=-0.063 0.515 --------- 1666(7) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 101.38 190.46 69.68 -4.51 + 5.68 = 1.17[+]55.95=57.13 57.31*sin(60.23)=49.74 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 190.17-103.27 13.88 46.24 2.96(0)182.92 20.43(0)86.90 13.99 5.00+185.07=190.07 184.83 190.46/190.46 =86.90 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 7.24(7.14) 22.55 44.38 1.04 20.43(86.90) -0.10/-0.11/-5.43 0.000 86.90/0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.962 16.645(15.392) -1.42 -5.66 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 5.851 9.542 1.25(75.17) 0.39+5.29=5.68 0.43 4.75/-0.020 186.63(4.68) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD -1.22/-5.86 22.20/38.71 23.07 0.70 0.43 + -0.45=-0.021 0.521 --------- 1669(5) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 37.83 130.25 34.80 18.34 + 5.68 = 24.03[+]43.71=67.74 57.42*sin(39.30)=36.37 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 128.90-218.42 16.98 59.81 306.00(2)121.61 6.93(3)270.49 13.50 0.32+130.19=130.51 133.96 130.25/130.26 =270.49 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 7.30(7.19) 44.47 30.61 0.49 6.93(270.49) 1.61/1.93/5.43 0.000 270.49/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.669 14.320(13.282)(>) -5.53 5.29 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 19.547 -6.265 1.04(62.27) 0.39+5.29=5.68 0.06 0.30/-0.001 182.52(0.62) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 1.36/-1.60 8.07/27.09 14.92 0.56 0.06 + -0.23=-0.173 0.523 --------- 1671(9) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 162.23 250.73 71.58 -22.30 + -5.70 = -28.00[+]50.30=22.29 56.09*sin(71.06)=53.05 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 253.00-313.83 9.20 16.20 66.62(0)245.67 9.12(3)299.17 26.31 166.15+85.11=251.26 258.04 250.73/250.74 =299.17 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 7.33(7.23) 20.69 61.76 1.06 9.12(299.17) -1.73/-2.07/4.80 0.000 299.17/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.227 16.772(15.931) -7.38 -5.40 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 4.838 -12.907 0.84(50.47) 0.41+5.29=5.70 1.22 166.83/0.054 350.26(166.61) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 0.82/-2.79 38.86/29.75 7.56 0.34 1.22 + -0.82=0.399 0.518 --------- 1676(6) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 80.28 171.07 107.19 3.86 + -5.72 = -1.86[+]47.18=45.31 54.16*sin(87.46)=54.11 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 170.30-147.93 4.52 18.51 347.87(3)162.87 15.49(0)22.36 27.03 178.18+-7.34=170.84 168.06 171.07/171.08 =22.36 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 7.42(7.32) -19.45 86.27 1.08 15.49(22.36) 0.55/0.66/-2.17 0.000 22.36/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.037 19.146(17.712) 4.54 2.45 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 12.154 5.558 1.43(86.02) 0.43+5.29=5.72 0.15 178.27/0.007 362.18(178.40) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD -0.09/2.17 24.39/43.29 23.18 0.64 0.15 + -0.63=-0.487 0.507 --------- 1681(10) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 170.66 259.83 -43.57 -23.15 + 5.68 = -17.47[+]-41.33=-58.81 57.87*sin(53.17)=46.32 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 262.23-1.78 29.36 17.99 75.81(0)254.72 6.16(2)260.45 14.45 9.18+251.21=260.39 267.87 259.84/259.85 =260.45 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 7.51(7.40) 40.90 45.07 -0.62 6.16(260.45) -1.84/-2.19/5.38 0.000 260.45/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.323 9.096(10.431) 2.21 -5.54 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 23.875 -13.445 -1.33(-80.10) 0.39+5.29=5.68 0.76 8.72/-0.037 190.26(8.28) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD -1.49/-1.84 10.48/36.12 4.03 -0.66 0.76 + -0.40=0.357 0.527 --------- 1685(12) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 242.82 334.78 55.31 -10.49 + -5.65 = -16.14[+]39.84=23.70 61.16*sin(79.87)=60.21 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 335.78-172.86 13.80 53.59 150.36(1)328.20 26.25(1)162.93 26.16 164.93+169.82=334.75 334.31 334.78/334.79 =162.93 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 7.58(7.47) 34.23 76.90 0.81 26.25(162.93) -1.03/-1.23/-1.68 0.000 162.93/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.294 15.687(14.796) -8.40 -5.57 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 14.053 0.743 0.89(53.49) 0.36+5.29=5.65 1.31 165.66/0.059 349.24(165.59) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD -0.10/-2.74 39.22/29.56 -21.13 0.40 1.31 + -0.92=0.389 0.547 --------- 1688(5) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 38.06 130.49 -45.93 18.28 + 5.66 = 23.94[+]-48.84=-24.89 60.22*sin(47.05)=44.08 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 128.99-271.56 16.90 39.50 305.74(2)121.35 2.87(2)217.43 14.07 7.90+122.69=130.59 132.05 130.49/130.50 =217.43 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 7.63(7.52) 39.74 35.26 -0.66 2.87(217.43) 1.62/1.93/3.40 0.000 217.43/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.942 8.938(9.616) 1.68 0.84 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 12.077 -2.461 -0.68(-40.69) 0.37+5.29=5.66 0.71 7.51/-0.032 189.73(7.75) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 0.25/2.27 26.94/24.94 15.00 -0.31 0.71 + -0.67=0.043 0.536 --------- 1690(9) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 162.32 250.84 -68.52 -22.31 + -5.67 = -27.98[+]-49.97=-77.96 58.60*sin(68.75)=54.62 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 253.08-6.97 9.12 55.89 66.37(0)245.41 5.06(2)246.12 26.88 173.74+77.61=251.35 258.34 250.84/250.85 =246.12 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 7.67(7.56) 23.64 58.84 -1.02 5.06(246.12) -1.73/-2.06/5.02 0.000 246.12/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.284 7.432(9.079) -0.14 -2.14 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 21.369 -12.290 -1.65(-98.79) 0.38+5.29=5.67 0.57 174.05/0.025 357.50(173.77) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 1.28/0.72 7.39/45.94 7.52 -0.89 0.57 + -0.19=0.381 0.529 --------- 1691(3) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL -18.61 69.84 26.57 22.19 + 5.71 = 27.90[+]25.18=53.08 55.03*sin(53.75)=44.38 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 67.72-26.71 19.99 53.07 243.53(2)60.04 16.92(0)41.02 13.58 1.44+67.95=69.39 63.81 69.84/69.85 =41.02 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 7.69(7.58) 46.91 50.98 0.37 16.92(41.02) 1.68/2.00/-3.66 0.000 41.02/0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.426 13.771(12.682) -4.45 7.20 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 1.773 10.909 1.09(65.34) 0.42+5.29=5.71 0.16 1.37/-0.006 183.60(1.68) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 2.99/-0.24 21.64/29.41 -7.89 0.59 0.15 + -0.44=-0.290 0.517 --------- 1692(2) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL -29.31 58.53 17.23 20.33 + 5.73 = 26.06[+]16.23=42.29 53.80*sin(54.33)=43.71 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 57.00-66.18 1.34 47.44 232.64(2)49.30 13.07(3)350.82 14.19 9.26+49.20=58.46 57.65 58.53/58.54 =350.82 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 7.70(7.59) 48.78 53.22 0.23 13.07(350.82) 1.47/1.76/0.92 0.000 350.82/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.428 13.148(12.248) 2.95 0.49 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 10.583 1.665 0.90(54.01) 0.44+5.29=5.73 0.83 8.80/-0.037 191.00(9.01) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 3.50/-0.06 28.02/24.32 -12.01 0.49 0.82 + -0.54=0.284 0.512 --------- 1695(12) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 253.19 344.51 67.30 -6.62 + 5.66 = -0.95[+]45.06=44.11 59.16*sin(88.46)=59.14 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 345.02-220.81 4.43 55.38 159.55(1)337.26 23.28(1)124.21 14.30 8.37+335.88=344.26 340.66 344.51/344.52 =124.21 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 7.76(7.65) 24.76 87.82 1.00 23.28(124.21) -0.75/-0.90/-4.56 0.000 124.21/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 1.000 16.487(15.339) 1.66 1.74 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 9.090 6.249 1.15(68.86) 0.38+5.29=5.66 0.71 7.95/-0.033 189.71(7.73) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 0.12/3.28 34.58/35.81 -22.58 0.69 0.71 + -0.71=-0.001 0.538 --------- 1697(5) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 29.41 121.57 -49.61 20.31 + -5.65 = 14.66[+]-48.12=-33.46 60.93*sin(51.69)=47.81 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 119.83-276.74 26.18 17.40 296.30(3)112.04 1.77(2)203.10 26.50 172.49+-50.91=121.58 121.72 121.57/121.59 =203.09 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 7.79(7.68) 37.76 40.20 -0.72 1.77(203.10) 1.76/2.10/2.23 0.000 203.10/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.933 8.692(9.357) -0.63 1.46 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 9.571 -0.214 -0.66(-39.85) 0.36+5.29=5.65 0.62 172.87/0.030 357.01(173.29) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 0.35/-2.44 30.97/24.68 12.04 -0.44 0.61 + -0.66=-0.050 0.540 --------- 1704(12) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 243.16 335.10 26.05 -10.37 + -5.67 = -16.04[+]21.68=5.64 58.37*sin(65.76)=53.22 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 335.87-226.00 13.72 33.27 150.11(1)327.95 22.19(1)109.87 26.73 172.50+162.32=334.82 330.97 335.10/335.11 =109.87 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 7.92(7.80) 47.04 64.15 0.36 22.19(109.87) -1.03/-1.24/-5.15 0.000 109.87/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.587 13.737(13.411) -1.18 1.37 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 6.583 6.827 0.33(19.56) 0.38+5.29=5.67 0.67 172.88/0.030 356.46(172.75) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD -1.34/-1.20 40.92/12.39 -21.19 0.09 0.66 + -0.88=-0.219 0.534 --------- 1706(5) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 48.64 141.08 104.00 15.27 + 5.66 = 20.92[+]48.11=69.03 60.38*sin(88.71)=60.37 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 139.79-285.22 5.94 24.82 316.38(3)131.84 2.65(2)214.58 14.04 7.21+133.97=141.18 142.68 141.08/141.09 =214.57 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 7.95(7.84) -16.13 88.07 1.11 2.65(214.58) 1.41/1.68/3.18 0.000 214.58/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.479 18.934(17.286) 0.98 -4.49 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 19.798 -2.512 1.65(98.85) 0.37+5.29=5.66 0.65 6.85/-0.029 189.03(7.06) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 0.07/-3.57 14.25/47.93 18.06 0.94 0.65 + -0.36=0.286 0.536 --------- 1708(10) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 172.21 261.52 84.55 -23.25 + -5.67 = -28.92[+]50.19=21.27 58.76*sin(83.66)=58.40 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 263.88-20.63 27.69 41.21 77.01(0)255.90 4.84(2)243.26 26.84 173.14+88.88=262.03 269.05 261.52/261.53 =243.26 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 7.99(7.87) 6.47 80.16 1.18 4.84(243.26) -1.85/-2.21/4.91 0.000 243.26/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.470 17.637(16.791) -0.73 1.48 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 5.090 -12.299 0.85(50.75) 0.38+5.29=5.67 0.63 173.49/0.027 356.91(173.18) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 0.38/-2.60 40.91/29.94 3.37 0.35 0.62 + -0.91=-0.283 0.531 --------- 1709(9) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 162.66 251.20 -80.33 -22.36 + -5.65 = -28.01[+]-50.66=-78.66 61.30*sin(77.87)=59.93 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 253.16-60.11 9.04 35.58 66.11(0)245.16 1.00(2)193.06 0.24 181.32+70.11=251.43 254.76 251.20/251.22 =193.06 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 8.01(7.89) 11.35 71.27 -1.16 1.00(193.06) -1.72/-2.06/1.31 0.000 193.06/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.500 6.645(8.374) 7.07 4.24 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 13.900 -5.525 -1.73(-103.78) 0.36+5.29=5.65 -0.08 181.26/-0.005 364.71(180.90) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 0.67/2.16 7.79/48.56 7.38 -0.98 -0.08 + -0.20=-0.279 0.541 --------- 1712(7) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 102.92 191.88 -103.81 -5.11 + -5.72 = -10.83[+]-48.36=-59.19 53.97*sin(85.32)=53.79 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 191.90-175.25 12.13 49.15 3.90(0)183.85 15.05(0)16.65 26.96 176.52+15.24=191.76 190.30 191.88/191.90 =16.65 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 8.06(7.94) -15.92 82.97 -1.12 15.05(16.65) -0.13/-0.15/-1.65 0.000 16.65/0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.671 5.080(6.552) 2.82 8.65 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 3.596 2.957 -1.47(-88.36) 0.43+5.29=5.72 0.30 176.70/0.014 360.46(176.70) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 0.06/-5.89 23.04/44.10 22.95 -0.77 0.30 + -0.46=-0.157 0.506 --------- 1715(5) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 39.77 132.22 102.95 17.84 + -5.65 = 12.18[+]48.37=60.56 61.04*sin(89.64)=61.04 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 130.64-290.40 15.23 2.72 306.94(3)122.53 1.55(2)200.24 26.47 171.85+-39.63=132.22 132.30 132.21/132.23 =200.23 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 8.11(8.00) -14.99 89.46 1.13 1.55(200.24) 1.60/1.91/1.98 0.000 200.24/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.702 18.863(17.174) -1.27 -4.46 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 17.291 -0.117 1.69(101.35) 0.36+5.29=5.65 0.68 172.26/0.033 356.37(172.65) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD -0.13/3.31 11.97/48.76 15.53 0.89 0.67 + -0.50=0.172 0.540 --------- 1719(3) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL -27.83 60.08 59.98 20.63 + -5.72 = 14.91[+]43.28=58.20 54.30*sin(75.61)=52.59 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 58.65-85.03 29.20 10.66 233.83(2)50.48 11.76(3)333.63 26.58 173.26+-113.12=60.14 60.86 60.08/60.09 =333.62 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 8.18(8.07) 30.88 70.59 0.89 11.76(333.63) 1.50/1.79/2.51 0.000 333.63/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.792 15.999(14.385) 0.05 -1.21 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 15.798 -1.413 1.61(96.85) 0.43+5.29=5.72 0.56 173.60/0.027 357.69(173.96) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD -0.88/2.14 14.85/44.87 -11.46 0.74 0.55 + -0.46=0.091 0.514 --------- 1720(8) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 134.27 221.80 5.69 -16.15 + 5.69 = -10.46[+]11.15=0.68 56.93*sin(23.17)=22.40 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 222.57-144.25 10.54 2.21 34.88(0)214.38 19.77(0)78.33 13.89 2.56+218.94=221.50 217.49 221.80/221.82 =78.32 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 8.20(8.08) 49.94 22.78 0.08 19.77(78.33) -1.06/-1.27/-5.35 0.000 78.33/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.668 12.379(12.454) -3.93 0.39 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 5.013 7.441 -0.07(-4.49) 0.40+5.29=5.69 0.20 2.43/-0.010 184.12(2.20) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD -3.70/0.16 12.95/8.33 17.28 0.00 0.20 + -0.37=-0.172 0.520 --------- 1721(8) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 123.47 211.23 95.87 -12.70 + 5.72 = -6.98[+]50.45=43.48 54.38*sin(81.67)=53.80 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 211.86-183.73 21.42 56.58 23.99(0)203.64 15.93(0)28.13 14.49 10.94+200.14=211.09 209.36 211.23/211.24 =28.13 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 8.22(8.10) -6.97 77.05 1.18 15.93(28.13) -0.75/-0.90/-2.67 0.000 28.13/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.447 18.391(17.273) 4.11 -8.06 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 13.823 3.450 1.12(67.08) 0.43+5.29=5.72 0.93 10.40/-0.043 192.16(10.15) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD -0.32/-3.64 35.30/31.28 19.92 0.62 0.93 + -0.65=0.276 0.508 --------- 1730(7) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 114.08 202.29 19.63 -9.36 + -5.72 = -15.08[+]37.98=22.90 53.90*sin(24.76)=22.57 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 202.70-188.91 1.17 34.47 14.54(0)194.33 14.83(0)13.79 26.92 175.69+26.53=202.22 201.43 202.29/202.31 =13.79 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 8.38(8.26) 48.39 19.98 0.27 14.83(13.79) -0.46/-0.55/-1.38 0.000 13.79/0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.909 13.309(12.946) 1.97 -5.03 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 11.316 1.629 0.36(21.79) 0.44+5.29=5.72 0.38 175.91/0.017 359.61(175.86) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 3.48/3.51 9.89/13.96 21.64 0.15 0.38 + -0.35=0.033 0.505 --------- 1731(13) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 442.50 366.88 -67.81 2.98 + 5.67 = 8.65[+]-45.25=-36.61 58.84*sin(89.30)=58.84 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 6.62-248.13 12.05 26.02 180.83(1)358.24 22.85(1)118.49 14.23 8.20+358.35=366.55 1.87 366.89/6.90 =118.49 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 8.40(8.27) 24.30 89.01 -1.01 22.85(118.49) -0.06/-0.07/-4.83 0.000 118.49/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.844 7.479(8.725) 1.63 -3.05 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 0.532 8.193 -1.25(-74.73) 0.38+5.29=5.67 0.71 7.78/-0.033 189.68(7.70) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 0.36/-1.78 32.07/38.15 -23.30 -0.58 0.70 + -0.79=-0.083 0.537 --------- 1735(11) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 200.56 292.24 -52.04 -21.91 + 5.67 = -16.23[+]-40.60=-56.83 58.29*sin(68.31)=54.16 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 294.64-42.76 26.02 33.96 107.73(0)286.18 5.50(2)251.88 14.35 7.78+285.00=292.78 300.10 292.24/292.26 =251.88 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 8.47(8.34) 36.33 62.35 -0.76 5.50(251.88) -1.85/-2.21/5.21 0.000 251.88/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.733 8.531(10.058) 0.87 -3.99 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 23.038 -12.979 -1.53(-91.66) 0.38+5.29=5.67 0.64 7.39/-0.031 188.92(6.95) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD -0.66/-2.46 15.24/42.28 -8.67 -0.76 0.64 + -0.49=0.142 0.531 --------- 1742(6) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 70.62 162.12 -71.95 7.60 + 5.65 = 13.25[+]-55.37=-42.12 60.68*sin(62.40)=53.78 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 161.39-312.54 13.56 55.47 337.66(3)152.81 2.21(2)208.86 13.97 5.69+156.54=162.23 163.93 162.12/162.14 =208.86 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 8.59(8.47) 20.33 47.39 -1.07 2.21(208.86) 0.86/1.03/2.73 0.000 208.86/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.726 7.204(8.413) -0.57 3.78 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 11.240 -2.827 -1.21(-72.57) 0.37+5.29=5.65 0.51 5.40/-0.023 187.48(5.52) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD -0.31/3.51 25.12/37.45 22.29 -0.60 0.51 + -0.66=-0.159 0.537 --------- 1745(4) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 11.45 102.45 -7.65 22.99 + 5.72 = 28.71[+]-10.08=18.62 54.60*sin(35.66)=31.83 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 100.13-67.69 16.65 9.04 275.45(2)91.50 16.26(0)32.45 13.47 0.30+101.72=102.02 96.81 102.45/102.46 =32.44 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 8.64(8.52) 49.83 35.24 -0.10 16.26(32.45) 1.91/2.29/-3.02 0.000 32.45/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.134 11.490(11.252) -5.48 6.27 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 0.936 10.316 0.24(14.28) 0.43+5.29=5.72 0.06 0.29/-0.001 182.57(0.66) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 0.41/0.38 23.89/5.88 4.93 0.17 0.06 + -0.50=-0.442 0.513 --------- 1746(4) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 1.22 91.33 -16.72 23.48 + 5.72 = 29.20[+]-19.16=10.04 53.99*sin(42.23)=36.29 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 89.41-107.17 27.53 3.41 264.56(2)80.76 12.42(3)342.25 14.08 8.33+82.96=91.29 90.79 91.33/91.35 =342.24 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 8.66(8.54) 48.86 40.61 -0.23 12.42(342.25) 1.90/2.27/1.74 0.000 342.25/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.672 10.885(10.782) 2.14 -0.58 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 9.746 1.036 0.10(6.21) 0.43+5.29=5.72 0.75 7.91/-0.033 190.19(8.21) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 0.74/0.82 29.77/3.64 0.53 0.11 0.75 + -0.60=0.155 0.511 --------- 1747(8) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 135.86 223.39 87.88 -16.62 + -5.70 = -22.32[+]51.67=29.35 56.08*sin(77.71)=54.80 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 224.22-163.09 8.87 25.43 36.07(0)215.56 18.46(0)61.13 26.28 166.50+56.62=223.11 219.67 223.38/223.40 =61.13 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 8.68(8.55) 2.53 70.64 1.19 18.46(61.13) -1.09/-1.31/-4.82 0.000 61.13/0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.260 17.859(16.814) -6.92 -6.57 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 10.227 6.587 1.04(62.67) 0.41+5.29=5.70 1.18 167.16/0.053 350.72(167.07) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 0.62/-0.97 36.58/31.47 16.83 0.45 1.18 + -0.80=0.380 0.516 --------- 1751(6) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 61.30 153.33 -68.44 11.03 + -5.65 = 5.38[+]-55.01=-49.63 61.24*sin(60.53)=53.32 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 152.24-317.72 22.84 33.36 328.21(3)143.50 1.11(2)194.52 26.40 170.41+-17.07=153.34 153.47 153.34/153.35 =194.52 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 8.75(8.62) 23.72 45.42 -1.02 1.11(194.52) 1.12/1.34/1.45 0.000 194.52/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.596 7.438(8.553) -2.74 -0.24 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 8.733 -0.180 -1.12(-66.91) 0.36+5.29=5.65 0.81 170.89/0.038 354.90(171.20) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 0.66/-3.17 26.90/35.15 20.86 -0.68 0.80 + -0.58=0.229 0.540 --------- 1758(13) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 440.61 368.62 68.37 3.73 + -5.68 = -1.95[+]45.48=43.53 57.95*sin(89.59)=57.95 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 8.27-266.98 10.38 49.24 182.02(1)359.41 21.53(1)101.30 26.62 172.21+196.03=368.24 3.00 368.62/8.64 =101.30 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 8.87(8.75) 23.78 89.42 1.02 21.53(101.30) -0.02/-0.02/-5.36 0.000 101.30/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.936 16.558(15.160) -1.26 1.72 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 5.746 9.414 1.40(83.87) 0.39+5.29=5.68 0.67 172.60/0.031 356.38(172.67) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD -0.25/-2.73 27.34/41.84 -23.20 0.67 0.67 + -0.70=-0.029 0.533 --------- 1760(6) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 81.93 172.59 80.46 3.21 + 5.65 = 8.86[+]54.51=63.37 60.82*sin(68.27)=56.50 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 172.20-326.20 2.60 40.79 348.30(3)163.30 1.99(2)206.00 13.93 4.88+167.83=172.72 174.55 172.59/172.61 =206.00 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 8.91(8.78) 11.20 55.52 1.16 1.99(206.00) 0.54/0.64/2.49 0.000 206.00/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 1.017 17.364(15.909) -1.41 -2.61 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 18.961 -3.052 1.46(87.32) 0.36+5.29=5.65 0.43 4.64/-0.020 186.64(4.70) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD -0.66/-3.36 19.47/43.83 23.28 0.84 0.43 + -0.42=0.008 0.537 --------- 1762(11) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 202.18 293.97 85.29 -21.65 + -5.66 = -27.32[+]49.91=22.59 59.26*sin(87.61)=59.21 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 296.29-61.61 24.35 57.18 108.92(0)287.36 4.19(2)234.69 26.74 171.76+122.68=294.43 301.06 293.96/293.98 =234.68 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 8.94(8.82) 5.60 93.71 1.19 4.19(234.69) -1.84/-2.20/4.51 0.000 234.69/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.690 17.686(16.802) -2.05 1.34 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 4.252 -11.450 0.88(53.04) 0.38+5.29=5.66 0.75 172.17/0.033 355.59(171.88) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD -0.01/-3.38 40.11/30.06 -9.31 0.38 0.74 + -0.91=-0.169 0.535 --------- 1763(10) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 192.38 283.46 -77.36 -22.91 + -5.65 = -28.55[+]-48.83=-77.38 61.43*sin(83.82)=61.08 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 285.57-101.09 5.69 51.54 98.03(0)276.62 0.34(2)184.49 0.14 179.73+103.91=283.65 286.37 283.46/283.48 =184.49 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 8.96(8.83) 14.66 80.67 -1.13 0.34(184.49) -1.91/-2.28/0.46 0.000 184.49/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.718 6.843(8.587) 5.52 2.27 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 13.062 -4.474 -1.74(-104.69) 0.36+5.29=5.65 0.06 179.75/0.001 363.16(179.37) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 0.40/2.85 9.11/48.60 -5.32 -0.99 0.06 + -0.22=-0.165 0.544 --------- 1769(6) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 72.46 163.84 85.04 6.90 + -5.65 = 1.25[+]53.40=54.65 61.31*sin(72.13)=58.35 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 163.04-331.38 11.89 18.69 338.85(3)153.99 0.89(2)191.67 26.36 169.62+-5.78=163.85 164.05 163.83/163.85 =191.66 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 9.07(8.94) 5.90 61.59 1.18 0.89(191.67) 0.82/0.98/1.17 0.000 191.67/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.433 17.669(16.143) -3.55 -7.20 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 16.454 -0.310 1.53(91.54) 0.36+5.29=5.65 0.88 170.14/0.041 354.09(170.40) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 0.63/3.02 17.99/45.08 22.50 0.79 0.88 + -0.56=0.314 0.539 --------- 1770(6) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 61.74 153.76 -54.58 10.87 + -5.66 = 5.22[+]-54.83=-49.61 60.20*sin(49.87)=46.03 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 152.32-10.86 22.76 13.05 327.96(3)143.25 24.61(1)141.47 26.97 178.00+-24.57=153.43 148.70 153.76/153.78 =141.46 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 9.09(8.95) 34.71 34.34 -0.80 24.61(141.47) 1.13/1.35/-3.48 0.000 141.47/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.364 8.362(9.338) 4.48 5.07 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 1.264 8.074 -0.98(-58.56) 0.37+5.29=5.66 0.15 178.10/0.008 362.12(178.34) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 1.50/-2.10 22.50/32.51 20.94 -0.58 0.15 + -0.50=-0.345 0.535 --------- 1773(4) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 2.63 92.87 35.63 23.46 + -5.72 = 17.74[+]35.38=53.12 54.68*sin(50.50)=42.19 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 91.06-126.01 25.86 26.62 265.75(2)81.93 11.10(3)325.06 26.48 172.31+-79.37=92.94 93.88 92.87/92.89 =325.05 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 9.14(9.01) 44.17 44.69 0.50 11.10(325.06) 1.90/2.27/3.21 0.000 325.06/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.582 14.375(13.145) -0.78 -0.79 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 14.960 -1.815 1.23(73.81) 0.43+5.29=5.72 0.63 172.69/0.031 356.86(173.14) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD -1.53/1.55 12.90/32.99 1.14 0.56 0.63 + -0.42=0.206 0.514 --------- 1774(9) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 164.76 253.46 -54.33 -22.61 + 5.70 = -16.92[+]-46.72=-63.63 56.52*sin(58.86)=48.38 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 254.98-185.23 7.20 18.17 66.80(0)245.84 19.12(0)69.76 13.78 0.45+252.78=253.23 250.22 253.46/253.48 =69.75 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 9.16(9.02) 34.88 48.61 -0.80 19.12(69.76) -1.73/-2.07/-5.15 0.000 69.76/0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.322 8.378(9.866) -6.07 2.58 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 4.175 5.691 -1.49(-89.29) 0.41+5.29=5.70 0.01 0.42/-0.002 181.98(0.08) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD -1.43/-2.06 8.58/40.18 6.51 -0.72 0.01 + -0.36=-0.351 0.520 --------- 1775(9) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 154.66 242.70 17.68 -21.11 + 5.72 = -15.39[+]25.10=9.71 54.05*sin(33.31)=29.68 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 244.26-224.71 18.08 12.54 55.90(0)235.10 15.27(0)19.56 14.39 8.69+234.00=242.69 242.65 242.70/242.72 =19.55 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 9.17(9.04) 48.71 30.75 0.24 15.27(19.56) -1.55/-1.85/-1.92 0.000 19.56/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.528 13.179(13.107) 1.80 -4.05 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 12.985 0.122 0.07(4.27) 0.43+5.29=5.72 0.72 8.25/-0.035 189.85(7.87) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD -2.88/0.62 22.93/7.92 10.53 0.08 0.72 + -0.49=0.232 0.508 --------- 1776(2) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL -57.06 30.72 -18.03 12.51 + -5.65 = 6.86[+]-15.56=-8.70 60.58*sin(62.27)=53.62 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 29.80-241.16 28.95 40.19 203.54(1)20.62 25.15(1)148.65 25.98 164.46+-134.00=30.46 26.79 30.72/30.74 =148.64 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 9.19(9.06) 48.66 61.38 -0.24 25.15(148.65) 0.69/0.82/-2.92 0.000 148.65/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.207 10.798(10.883) -8.51 -5.62 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 4.656 6.227 -0.09(-5.11) 0.37+5.29=5.65 1.32 165.22/0.060 349.13(165.49) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD -3.24/0.34 41.30/5.42 -20.03 -0.14 1.32 + -0.88=0.433 0.544 --------- 1784(9) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 145.98 233.64 -80.34 -19.28 + -5.73 = -25.01[+]-51.89=-76.90 53.75*sin(73.94)=51.66 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 235.11-229.89 27.36 50.44 46.46(0)225.79 14.18(0)5.22 26.82 173.34+60.40=233.74 234.82 233.65/233.67 =5.22 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 9.33(9.20) 11.35 64.99 -1.16 14.18(5.22) -1.35/-1.62/-0.54 0.000 5.22/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.198 6.644(8.328) -0.45 -0.96 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 10.479 -2.151 -1.68(-101.01) 0.44+5.29=5.73 0.60 173.67/0.027 357.19(173.47) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 0.90/-0.39 9.17/47.86 13.66 -0.84 0.60 + -0.20=0.403 0.506 --------- 1785(8) SprRed WinYel T_12Red YelRedXA+YelWhXA=RedWhXA[+]URedXA=UWhXA HighLow*sin(UseSunHL)=UseHL 135.35 222.87 -63.11 -16.47 + -5.71 = -22.18[+]-53.16=-75.34 55.23*sin(58.72)=47.20 MidMoon-MaxMoon LeapDay ReqMDay ReqSun()360 ReqMoon()360 ReqEcl SunPos+UsePos=RSunYel UseRMoon RSunEnd/RMoonEnd 224.39-269.37 8.71 44.80 35.56(0)215.05 10.33(3)315.02 0.22 181.68+41.61=223.29 228.45 222.88/222.90 =315.02 MinSun(360) T_12X T_12NP T_12Arc OReqMoon(360) Sun_Th/SunXS/MoonXS Ee_o EMoon360/MOff/MF 9.35(9.22) 28.40 44.62 -0.94 10.33(315.02) -1.08/-1.29/3.92 0.000 315.02/-0.00/0.000 Eclipse)Frac MaxPh(MDayT) TestEA SexEA 0.361 7.793(9.277) 7.53 6.95 MDay0_24 MoonT24 SeeT(min) XA+YelWh=XYelA LatYel MoonPos/YelD SolEA(360) 19.289 -10.011 -1.48(-89.07) 0.42+5.29=5.71 -0.13 181.59/-0.007 365.17(181.36) NS_Sh/EW_Sh NS_Hu/EW_Hu ESunLat EW_f MaxLat+NS_f=ESeeLat addD 1.70/-1.12 8.79/42.08 16.98 -0.76 -0.12 + -0.20=-0.324 0.512 ======================================== Shoushi calendar method referring to Chonzhen Lishu 200 year eclipse table for drift problem drift rate = offset/347 and 347 = 1628 - 1281(base year for Shoushi) observed_time is integrated with theorectical HuiHui calendar method table eclipse fraction is integrated from HouBian for comparison with Shoushi Data record in Chongzhen Lishu(page 788-34): 65 fan late => 0.0065 6 days late => more than 6 2 degree late => (2/365.2425)* 27.5546 => more than 0.150883864 Shoushi drift from Adam Schall 200 eclipse table for year 1628 calculated data from http://ideone.com/BlkMI for moon_days = -(0.005499/347) => approx -0.0000158473 for sun = -(5.99089/347) => -0.0172648 for moon = (0.162147/347) => 0.000467281 for eclipse = (0.0122486/347) => 0.0000352987 ----------------------------- slightly modified according to data record in Chongzhen Lishu(page 788-34) as : RefMoonDay = -0.0081 + RefMoonDay + 0.01041 - 0.01041/347 *(calculated_year - 1281) MinSun = 6.572 - 6.1419 +(6.1419/347)*(calculated_year - 1281) RefSun = BaseSun - MinSun RefMoon = 0.0402 + RefMoon + 0.16309/347 *(calculated_year - 1281) RefEcl = RefEcl + 0.01735/347*(calculated_year - 1281) ---------------------------- Inverted Kepler elliptic model is implemented as : ReqSun360 = fmod(ReqSun * 360/Tropical_year + 180, 360) e = 0.0167 SunX_Th = 180/pi*(2 *e* sin(ReqSun360 *pi/180)+ 1.25* e*e * sin(2*ReqSun360 *pi/180)) SFactor = 13.36875 * 0.082 + 0.082// = 1.1782375 SunX360 = SFactor * SunX_Th SunX = Tropical_year/360* SunX360 ================================================= RSunYel = ReqSun360 + SunX_Th + MinSun360 + 0.098 ================================================= ---------------------- Convert Shoushi ReqMoon to consider Ee_o OReqMoon = ReqMoon OReqMoon360 = fmod(OReqMoon * 360/Zhuanzhong_C + 180, 360) if (OReqMoon360 < 0){OReqMoon360 = fmod(360 + OReqMoon360, 360);} apply Earth Center Offset for Moon adjustment Ee_o = 0 if (cos(OReqMoon360 * pi/180)< 0){Ee_o = 0} y_scale = 1 x_cir = cos(OReqMoon360 * pi/180) y_cir = sin(OReqMoon360 * pi/180) r_elip = 1/sqrt(x_cir*x_cir + y_cir*y_cir/(y_scale*y_scale)) x_o = r_elip * x_cir y_o = r_elip * y_cir EMoon360 = 180/pi *atan(y_o/(x_o - Ee_o)) if ((OReqMoon360 > 90)&(OReqMoon360 < 270)){EMoon360 = 180 + EMoon360;} if (OReq_moon360 > 270){EMoon360 = 360 + EMoon360;} ---------------------- MOff = EMoon360 - OReqMoon360 EMoon = fmod((EMoon360 + 180)*Zhuanzhong_C/360, Zhuanzhong_C) ReqMoon = EMoon MDay0 = ReqMoonDay -(int)ReqMoonDay HuiHui observed time table is implemented as : degree = fmod(RSunEnd + 270, 360) Gon = degree/30.0 Hour = MDayT_hour * 24 int SeeTimeTable[][19] = { h_deg -120 -105 -90 -75 -60 -45 -30 -15 0 15 30 45 60 75 90 105 120 <-- Gon h20 h19 h18 h17 h16 h15 h14 h13 h12 h11 h10 h9 h8 h7 h6 h5 h4 Gon ----------------------------------------------------------------------------------------------- {9, -61, -72, -81, -85, -85, -78, -61, -34, -0, 34, 61, 78, 85, 85, 81, 72, 61, 9}, {10, -51, -61, -68, -70, -66, -54, -32, -2, 33, 63, 84, 95, 99, 96, 90, 80, 67, 8}, {11, -40, -49, -55, -55, -50, -37, -14, 15, 48, 75, 93, 103, 106, 103, 96, 86, 74, 7}, {0, -33, -42, -48, -50, -46, -35, -15, 13, 44, 72, 92, 103, 108, 106, 100, 91, 80, 6}, {1, -34, -44, -52, -55, -54, -45, -28, -3, 29, 60, 84, 98, 105, 105, 100, 92, 81, 5}, {2, -45, -56, -64, -69, -69, -62, -46, -20, 13, 46, 73, 90, 97, 98, 94, 86, 75, 4}, {3, -61, -72, -81, -85, -85, -78, -61, -34, 0, 34, 61, 78, 85, 85, 81, 72, 61, 3}, ----------------------------------------------------------------------------------------------- Gon h4 h5 h6 h7 h8 h9 h10 h11 h12 h13 h14 h15 h16 h17 h18 h19 h20 Gon --> observed_time = SeeT/(24*60) eclipseMax_hour = MDayT_hour + observed_time t t t t t t t t t t t t t t ====================================================================================================================== Shoushi(year)(months) eclipseFraction eclipseBegin MaxPhase eclipseEnd Theorectical(month/day) (EclipseFracS) (eclipseBeginS) (observed_time x 24 x 60) (eclipse_endS) half_day_hours/ReqSun/ReqMoon (observed_timeS x 24 x 60) SunXS(SunStatusS)/MoonSpeedFactor MoonXS(MoonStatusS)/SunMoonXS Gon/MDayT ______________________________________________________________________________________________________________________ 1644(9) 2687.29/65.4/10.4 0.4034(0.9845) 10.7612(10.4748) 11.9254(-23.5225) 13.0897(13.376) 9/1 -2.046(0)/1.024 -)0.3030 -)10.9410 -)12.1520(9.6583) -)13.3540 3.834(3)/-5.880 -------------------------------------------------------------------------------------- 5.31/12.32 0.1004(0.6815) -0.1798(-0.4662) -0.2266 -0.2643(0.0220) 1648(7) 2928.11/357.6/20.6 0.9055(0.9906) 5.62438(5.6184) 6.96737(-82.6023) 8.31036(8.31633) 6/21 0.295(3)/1.102 -)0.9820 -)6.1690 -)7.2130(-79.4410) -)8.3740 -5.424(0)/5.719 -------------------------------------------------------------------------------------- 3.00/8.34 -0.0765(0.0086) -0.5446(-0.5506) -0.2456 -0.0636(-0.0577) 1650(11) 2306.48/118.2/22.8 0.9458(0.9481) 10.2156(10.2155) 11.5063(-40.8428) 12.797(12.7971) 10/25 -2.066(1)/1.144 -)0.9040 -)10.2670 -)11.6030(5.7533) -)12.9910 -4.836(1)/2.770 -------------------------------------------------------------------------------------- 7.06/12.19 0.0418(0.0441) -0.0514(-0.0515) -0.0967 -0.1940(-0.1939) 1657(6) 2910.99/348.1/19.6 0.6013(0.9251) 3.95455(3.84453) 5.21388(-78.2361) 6.47321(6.58323) 6/12 0.651(3)/1.074 -)0.6110 -)4.3010 -)5.1430(-93.0517) -)6.0490 -5.296(0)/5.946 -------------------------------------------------------------------------------------- 2.71/6.52 -0.0097(0.3141) -0.3464(-0.4565) 0.0709 0.4242(0.5342) 1658(6) 2877.81/337.2/15.7 0.5443(0.3561) 9.58911(9.77514) 10.9133(-17.9025) 12.2375(12.0515) 6/11 1.041(3)/1.004 -)0.5300 -)9.0840 -)10.5010(-23.0159) -)12.0490 -2.419(0)/3.459 -------------------------------------------------------------------------------------- 2.35/11.21 0.0143(-0.1739) 0.5051(0.6911) 0.4123 0.1885(0.0025) 1665(1) 2094.00/201.9/12.4 0.8410(0.5132) 14.6974(14.8685) 16.1828(89.5742) 17.6682(17.497) 1/16 0.760(2)/0.997 -)0.9080 -)14.8220 -)16.2670(65.2165) -)17.5180 1.745(3)/-0.985 -------------------------------------------------------------------------------------- 9.89/14.69 -0.0670(-0.3948) -0.1246(0.0465) -0.0842 0.1502(-0.0210) 1666(7) 2931.43/3.0/20.4 0.9623(0.2924) 15.2886(15.6858) 16.6451(75.1688) 18.0016(17.6044) 7/2 -0.114(0)/1.096 -)0.9880 -)15.5630 -)16.7640(76.0410) -)17.8350 -5.428(0)/5.314 -------------------------------------------------------------------------------------- 3.35/15.39 -0.0257(-0.6956) -0.2744(0.1228) -0.1189 0.1666(-0.2306) 1669(5) 2724.07/306.0/6.9 0.6692(0.1178) 13.0378(13.6801) 14.3197(62.2699) 15.6016(14.9593) 4/30 1.927(3)/1.099 -)0.5380 -)13.1170 -)14.3860(35.7791) -)15.5780 5.426(3)/-3.499 -------------------------------------------------------------------------------------- 1.34/13.28 0.1312(-0.4202) -0.0792(0.5631) -0.0663 0.0236(-0.6187) 1671(9) 2681.01/66.6/9.1 0.2266(0.3249) 15.8782(15.7319) 16.7722(50.4681) 17.6662(17.8125) 9/3 -2.066(0)/1.048 -)0.1220 -)16.4260 -)16.9680(82.6027) -)17.4910 4.804(3)/-6.871 -------------------------------------------------------------------------------------- 5.36/15.93 0.1046(0.2029) -0.5478(-0.6941) -0.1958 0.1752(0.3215) 1676(6) 2910.99/347.9/15.5 0.0373(0.5909) 18.7419(17.7832) 19.1458(86.0217) 19.5498(20.5085) 6/11 0.660(3)/1.002 -)0.0924 -)18.3401 -)18.7790(93.5342) -)19.2070 -2.169(0)/2.829 -------------------------------------------------------------------------------------- 2.70/17.71 -0.0551(0.4986) 0.4018(-0.5569) 0.3668 0.3428(1.3015) 1681(10) 2623.19/75.8/6.2 0.3231(0.3072) 8.12116(8.14096) 9.09559(-80.096) 10.07(10.0502) 9/12 -2.195(0)/1.109 -)0.3830 -)8.2010 -)9.2250(-42.6313) -)10.3330 5.383(2)/-7.577 -------------------------------------------------------------------------------------- 5.66/10.43 -0.0599(-0.0758) -0.0798(-0.0600) -0.1294 -0.2630(-0.2828) 1685(12) 2137.53/150.4/26.2 0.2944(0.3039) 14.8116(14.8) 15.6872(53.4917) 16.5629(16.5745) 11/26 -1.227(1)/1.196 -)0.2640 -)14.6630 -)15.4780(67.0117) -)16.2470 -1.678(1)/0.452 -------------------------------------------------------------------------------------- 8.16/14.80 0.0304(0.0399) 0.1486(0.1370) 0.2092 0.3159(0.3275) 1688(5) 2724.07/305.7/2.9 0.9424(0.8608) 7.68186(7.69203) 8.93814(-40.6853) 10.1944(10.1842) 4/30 1.932(3)/1.175 -)0.9600 -)8.1820 -)9.2910(-59.6951) -)10.4950 3.402(2)/-1.470 -------------------------------------------------------------------------------------- 1.35/9.62 -0.0176(-0.0992) -0.5001(-0.4900) -0.3529 -0.3006(-0.3108) 1690(9) 2681.01/66.4/5.1 0.2835(0.7322) 6.52837(6.18394) 7.43214(-98.7911) 8.33592(8.68034) 9/3 -2.062(0)/1.135 -)0.2160 -)6.7870 -)7.4760(-69.0673) -)8.2060 5.024(2)/-7.086 -------------------------------------------------------------------------------------- 5.36/9.08 0.0675(0.5162) -0.2586(-0.6031) -0.0439 0.1299(0.4743) 1691(3) 2281.95/243.5/16.9 0.4265(-0.1999) 12.58(-nan) 13.7713(65.3359) 14.9625(-nan) 2/28 2.003(2)/1.022 -)0.3750 -)12.2890 -)13.5260(20.1111) -)14.6810 -3.661(0)/5.663 -------------------------------------------------------------------------------------- 11.33/12.68 0.0515(-0.5749) 0.2910( -nan) 0.2453 0.2815( -nan) 1692(2) 2224.12/232.6/13.1 0.4277(0.9187) 11.9038(11.6358) 13.1484(54.0062) 14.393(14.661) 2/17 1.758(2)/0.990 -)0.5320 -)11.2760 -)12.7870(7.5985) -)14.2230 0.920(3)/0.838 -------------------------------------------------------------------------------------- 10.95/12.25 -0.1043(0.3867) 0.6278(0.3598) 0.3614 0.1700(0.4380) 1695(12) 2105.40/159.6/23.3 1.0001(0.7094) 15.2058(15.2611) 16.4868(68.8589) 17.7678(17.7125) 12/6 -0.901(1)/1.153 -)0.9730 -)15.3470 -)16.5180(75.3119) -)17.5700 -4.561(1)/3.660 -------------------------------------------------------------------------------------- 8.48/15.34 0.0271(-0.2636) -0.1412(-0.0859) -0.0312 0.1978(0.1425) 1697(5) 2668.37/296.3/1.8 0.9328(0.7565) 7.45225(7.48687) 8.69238(-39.8497) 9.93251(9.89789) 4/21 2.104(3)/1.190 -)0.9960 -)7.7950 -)8.8990(-64.4105) -)10.1090 2.233(2)/-0.130 -------------------------------------------------------------------------------------- 1.05/9.36 -0.0632(-0.2395) -0.3427(-0.3081) -0.2066 -0.1765(-0.2111) 1704(12) 2137.53/150.1/22.2 0.5875(0.7722) 12.5463(12.4643) 13.7369(19.5595) 14.9274(15.0095) 11/27 -1.235(1)/1.130 -)0.5770 -)12.2410 -)13.6240(38.9060) -)14.9380 -5.155(1)/3.920 -------------------------------------------------------------------------------------- 8.17/13.41 0.0105(0.1952) 0.3053(0.2233) 0.1129 -0.0106(0.0715) 1706(5) 2786.73/316.4/2.6 0.4791(0.4389) 17.8625(17.895) 18.9335(98.8525) 20.0045(19.9721) 5/12 1.681(3)/1.178 -)0.5740 -)18.2610 -)19.0740(92.4223) -)19.8390 3.185(2)/-1.504 -------------------------------------------------------------------------------------- 1.70/17.29 -0.0949(-0.1351) -0.3985(-0.3660) -0.1405 0.1655(0.1331) 1708(10) 2610.09/77.0/4.8 0.4704(0.7527) 16.5424(16.3866) 17.6368(50.7542) 18.7312(18.887) 9/14 -2.207(0)/1.139 -)0.5410 -)16.5410 -)17.4650(89.9430) -)18.3280 4.912(2)/-7.120 -------------------------------------------------------------------------------------- 5.72/16.79 -0.0706(0.2117) 0.0014(-0.1544) 0.1718 0.4032(0.5590) 1709(9) 2681.01/66.1/1.0 0.5002(0.2936) 5.57814(5.77306) 6.64457(-103.78) 7.71101(7.51609) 9/4 -2.058(0)/1.199 -)0.5230 -)5.8590 -)6.6810(-79.0702) -)7.5600 1.308(2)/-3.365 -------------------------------------------------------------------------------------- 5.37/8.37 -0.0228(-0.2294) -0.2809(-0.0859) -0.0364 0.1510(-0.0439) 1712(7) 2931.06/3.9/15.1 0.6714(-0.4421) 3.65923(-nan) 5.07954(-88.3583) 6.49984(-nan) 7/4 -0.150(0)/0.997 -)0.6260 -)3.7840 -)4.6840(-92.9560) -)5.6590 -1.648(0)/1.497 -------------------------------------------------------------------------------------- 3.40/6.55 0.0454(-1.0681) -0.1248( -nan) 0.3955 0.8408( -nan) 1715(5) 2730.05/306.9/1.5 0.7023(0.4423) 17.6794(17.834) 18.8633(101.35) 20.0472(19.8926) 5/3 1.907(3)/1.193 -)0.6040 -)18.1460 -)18.9550(91.9737) -)19.7150 1.978(2)/-0.071 -------------------------------------------------------------------------------------- 1.41/17.17 0.0983(-0.1617) -0.4666(-0.3120) -0.0917 0.3322(0.1776) 1719(3) 2229.64/233.8/11.8 0.7924(0.8487) 14.539(14.5237) 15.9989(96.8461) 17.4587(17.474) 2/19 1.788(2)/1.004 -)0.7410 -)14.5860 -)15.9690(59.7137) -)17.1860 2.509(3)/-0.721 -------------------------------------------------------------------------------------- 11.00/14.38 0.0514(0.1077) -0.0470(-0.0623) 0.0299 0.2727(0.2880) 1720(8) 2854.22/34.9/19.8 0.6685(0.9346) 11.0864(11.0118) 12.3794(-4.48657) 13.6724(13.747) 8/4 -1.270(0)/1.079 -)0.7360 -)10.6910 -)12.2250(13.6558) -)13.7170 -5.351(0)/4.082 -------------------------------------------------------------------------------------- 4.39/12.45 -0.0675(0.1986) 0.3954(0.3208) 0.1544 -0.0446(0.0300) 1721(8) 2894.60/24.0/15.9 0.4467(-0.0081) 17.1537(-nan) 18.3913(67.0821) 19.6289(-nan) 7/24 -0.900(0)/1.007 -)0.3900 -)17.3300 -)18.2020(92.3745) -)19.0120 -2.666(0)/1.766 -------------------------------------------------------------------------------------- 4.04/17.27 0.0567(-0.3981) -0.1763( -nan) 0.1893 0.6169( -nan) 1730(7) 2917.52/14.5/14.8 0.9087(0.3714) 11.8067(12.1357) 13.309(21.7938) 14.8113(14.4823) 7/15 -0.555(0)/0.994 -)0.8290 -)11.0970 -)12.8500(27.2257) -)14.6040 -1.378(0)/0.823 -------------------------------------------------------------------------------------- 3.74/12.95 0.0797(-0.4576) 0.7097(1.0387) 0.4590 0.2073(-0.1217) 1731(13) 2068.36/180.8/22.8 0.8439(0.6191) 6.20608(6.28743) 7.47931(-74.7334) 8.75254(8.67119) 12/29 -0.072(1)/1.144 -)0.8670 -)6.3417 -)7.4090(-74.4143) -)8.6020 -4.834(1)/4.762 -------------------------------------------------------------------------------------- 9.23/8.72 -0.0231(-0.2479) -0.1356(-0.0543) 0.0703 0.1505(0.0692) 1735(11) 2370.28/107.7/5.5 0.7328(0.5008) 7.27224(7.39914) 8.53082(-91.6588) 9.7894(9.66249) 10/16 -2.210(1)/1.125 -)0.8250 -)7.6960 -)8.8690(-50.8564) -)10.1680 5.208(2)/-7.418 -------------------------------------------------------------------------------------- 6.74/10.06 -0.0922(-0.3242) -0.4238(-0.2969) -0.3382 -0.3786(-0.5055) 1742(6) 2881.38/337.7/2.2 0.7256(0.3702) 6.0039(6.23448) 7.20363(-72.5713) 8.40336(8.17278) 6/3 1.026(3)/1.185 -)0.7020 -)6.5930 -)7.5120(-78.5852) -)8.5110 2.727(2)/-1.701 -------------------------------------------------------------------------------------- 2.40/8.41 0.0236(-0.3318) -0.5891(-0.3585) -0.3084 -0.1076(-0.3382) 1745(4) 2530.69/275.4/16.3 0.1341(-0.0447) 10.7542(-nan) 11.4903(14.2837) 12.2263(-nan) 4/2 2.286(3)/1.012 -)0.0820 -)10.8520 -)11.4960(-21.9128) -)12.1430 -3.018(0)/5.304 -------------------------------------------------------------------------------------- 0.41/11.25 0.0521(-0.1267) -0.0978( -nan) -0.0057 0.0833( -nan) 1746(4) 2416.29/264.6/12.4 0.6718(0.9272) 9.46217(9.38271) 10.8855(6.21344) 12.3088(12.3883) 3/22 2.265(2)/0.997 -)0.7650 -)9.3150 -)10.8320(-34.2009) -)12.4300 1.741(3)/0.525 -------------------------------------------------------------------------------------- 0.04/10.78 -0.0932(0.1622) 0.1472(0.0677) 0.0535 -0.1212(-0.0417) 1747(8) 2845.58/36.1/18.5 0.2604(0.1790) 16.9067(17.0512) 17.859(62.6714) 18.8113(18.6668) 8/16 -1.308(0)/1.050 -)0.2490 -)17.0010 -)17.7000(90.0898) -)18.3610 -4.822(0)/3.514 -------------------------------------------------------------------------------------- 4.45/16.81 0.0114(-0.0700) -0.0943(0.0502) 0.1590 0.4503(0.3058) 1751(6) 2841.13/328.2/1.1 0.5962(0.9703) 6.30792(6.20331) 7.43759(-66.9095) 8.56725(8.67186) 5/25 1.338(3)/1.198 -)0.4770 -)6.7140 -)7.5430(-76.7799) -)8.4420 1.448(2)/-0.110 -------------------------------------------------------------------------------------- 2.11/8.55 0.1192(0.4933) -0.4061(-0.5107) -0.1054 0.1253(0.2299) 1758(13) 2068.30/182.0/21.5 0.9364(0.7129) 15.2382(15.2911) 16.5577(83.8676) 17.8773(17.8243) 12/30 -0.024(1)/1.115 -)0.9000 -)15.0050 -)16.2650(72.7450) -)17.3861 -5.356(1)/5.332 -------------------------------------------------------------------------------------- 9.29/15.16 0.0364(-0.1871) 0.2332(0.2861) 0.2927 0.4912(0.4382) 1760(6) 2913.29/348.3/2.0 1.0169(0.6736) 16.1199(16.1879) 17.3641(87.3247) 18.6083(18.5403) 6/13 0.644(3)/1.187 -)0.9850 -)16.3120 -)17.3430(82.3605) -)18.2840 2.486(2)/-1.841 -------------------------------------------------------------------------------------- 2.75/15.91 0.0319(-0.3114) -0.1921(-0.1241) 0.0211 0.3243(0.2563) 1762(11) 2363.78/108.9/4.2 0.6901(0.7761) 16.4727(16.4423) 17.6863(53.0448) 18.8999(18.9303) 10/17 -2.197(1)/1.152 -)0.7560 -)16.6800 -)9999.0000(90.0138) -)9999.0000 4.508(2)/-6.705 -------------------------------------------------------------------------------------- 6.80/16.80 -0.0659(0.0201) -0.2073(-0.2377) -9981.3137 -9980.1001(-9980.0697) 1763(10) 2436.16/98.0/0.3 0.7178(0.6212) 5.66678(5.70829) 6.8426(-104.692) 8.01843(7.97692) 10/7 -2.278(1)/1.205 -)0.7340 -)6.0403 -)6.9330(-76.3150) -)7.8910 0.461(2)/-2.739 -------------------------------------------------------------------------------------- 6.45/8.59 -0.0162(-0.1128) -0.3735(-0.3320) -0.0904 0.1274(0.0859) 1769(6) 2884.85/338.9/0.9 0.4334(0.1004) 16.6536(17.131) 17.6691(91.5436) 18.6845(18.2072) 6/4 0.984(3)/1.200 -)0.3330 -)16.9840 -)17.6860(84.7730) -)18.3480 1.174(2)/-0.189 -------------------------------------------------------------------------------------- 2.46/16.14 0.1004(-0.2326) -0.3304(0.1470) -0.0169 0.3365(-0.1408) 1770(6) 2841.13/328.0/24.6 0.3645(0.1543) 7.39311(7.6921) 8.36158(-58.5619) 9.33006(9.03106) 5/25 1.346(3)/1.175 -)0.4140 -)7.5200 -)8.3940(-64.7404) -)9.3390 -3.476(1)/4.822 -------------------------------------------------------------------------------------- 2.13/9.34 -0.0495(-0.2597) -0.1269(0.1721) -0.0324 -0.0089(-0.3079) 1773(4) 2429.53/265.7/11.1 0.5819(0.9009) 13.0331(12.9051) 14.3753(73.8122) 15.7176(15.8456) 3/23 2.272(2)/1.013 -)0.4150 -)13.0150 -)14.3900(32.3707) -)15.6440 3.206(3)/-0.934 -------------------------------------------------------------------------------------- 0.10/13.15 0.1669(0.4859) 0.0181(-0.1099) -0.0147 0.0736(0.2016) 1774(9) 2681.01/66.8/19.1 0.3219(0.5695) 7.35254(7.11874) 8.37792(-89.2856) 9.4033(9.6371) 9/16 -2.069(0)/1.064 -)0.3870 -)7.2590 -)8.2510(-54.8279) -)9.3370 -5.146(0)/3.077 -------------------------------------------------------------------------------------- 5.45/9.87 -0.0651(0.1825) 0.0935(-0.1403) 0.1269 0.0663(0.3001) 1775(9) 2747.66/55.9/15.3 0.5279(0.4933) 11.8546(11.8838) 13.1786(4.27464) 14.5026(14.4734) 8/26 -1.854(0)/0.999 -)0.4680 -)11.4050 -)12.9230(31.4109) -)14.3340 -1.916(0)/0.062 -------------------------------------------------------------------------------------- 5.09/13.11 0.0599(0.0253) 0.4496(0.4788) 0.2556 0.1686(0.1394) 1776(2) 2096.67/203.5/25.2 0.2069(0.2980) 10.038(9.90952) 10.7977(-5.10649) 11.5573(11.6859) 1/21 0.821(2)/1.183 -)0.1450 -)9.5260 -)10.1970(-31.6621) -)10.8890 -2.923(1)/3.744 -------------------------------------------------------------------------------------- 10.02/10.88 0.0619(0.1530) 0.5120(0.3835) 0.6007 0.6683(0.7969) 1784(9) 2797.30/46.5/14.2 0.1977(0.8798) 5.73641(5.13441) 6.64426(-101.012) 7.55211(8.1541) 8/16 -1.616(0)/0.989 -)0.1870 -)5.5740 -)6.2700(-79.6391) -)7.0190 -0.536(0)/-1.080 -------------------------------------------------------------------------------------- 4.79/8.33 0.0107(0.6928) 0.1624(-0.4396) 0.3743 0.5331(1.1351) 1785(8) 2849.95/35.6/10.3 0.3611(-0.1586) 6.67784(-nan) 7.7929(-89.0741) 8.90797(-nan) 8/5 -1.291(0)/1.025 -)0.4590 -)6.6340 -)7.6910(-65.7765) -)8.8610 3.924(3)/-5.215 -------------------------------------------------------------------------------------- 4.43/9.28 -0.0979(-0.6176) 0.0438( -nan) 0.1019 0.0470( -nan)