日出日落时间计算程序(C语言).doc

/日出日落时间计算C语言程序#define PI 3.1415926 #include #include using namespace std; int days_of_month_1=31,28,31,30,31,30,31,31,30,31,30,31; int days_of_month_2=31,29,31,30,31,30,31,31,30,31,30,31; long double h=-0.833; /定义全局变量 void input_date(int c) int i; coutEnter the date (form: 2009 03 10):endl; for(i=0;ici; /输入日期 void input_glat(int c) int i; coutEnter the degree of latitude(range: 0- 60,form: 40 40 40 (means 404040):endl; for(i=0;ici; /输入纬度 void input_glong(int c) int i; coutEnter the degree of longitude(west is negativ,form: 40 40 40 (means 404040):endl; for(i=0;ici; /输入经度 int leap_year(int year) if(year%400=0) | (year%100!=0) & (year%4=0) return 1; else return 0; /判断是否为闰年:若为闰年,返回1；若非闰年,返回0 int days(int year, int month, int date) int i,a=0; for(i=2000;iyear;i+) if(leap_year(i) a=a+366; else a=a+365; if(leap_year(year) for(i=0;imonth-1;i+) a=a+days_of_month_2i; else for(i=0;i=UTo) d=UT-UTo; else d=UTo-UT; if(d=0.1) UTo=UT; UT=UT_rise(UTo,GHA(UTo,G_sun(t_century(days(year,month,date),UTo),ecliptic_longitude(L_sun(t_century(days(year,month,date),UTo),G_sun(t_century(days(year,month,date),UTo),glong,e(h,glat,sun_deviation(earth_tilt(t_century(days(year,month,date),UTo),ecliptic_longitude(L_sun(t_century(days(year,month,date),UTo),G_sun(t_century(days(year,month,date),UTo); result_rise(UT,UTo,glong,glat,year,month,date); return UT; /判断并返回结果（日出） long double result_set(long double UT, long double UTo, long double glong, long double glat, int year, int month, int date) long double d; if(UT=UTo) d=UT-UTo; else d=UTo-UT; if(d=0.1) UTo=UT; UT=UT_set(UTo,GHA(UTo,G_sun(t_century(days(year,month,date),UTo),ecliptic_longitude(L_sun(t_century(days(year,month,date),UTo),G_sun(t_century(days(year,month,date),UTo),glong,e(h,glat,sun_deviation(earth_tilt(t_century(days(year,month,date),UTo),ecliptic_longitude(L_sun(t_century(days(year,month,date),UTo),G_sun(t_century(days(year,month,date),UTo); result_set(UT,UTo,glong,glat,year,month,date); return UT; /判断并返回结果（日落） int Zone(long double glong) if(glong=0) return (int)(int)(glong/15.0)+1); else return (int)(int)(glong/15.0)-1); /求时区 void output(long double rise, long double set, long double glong) if(int)(60*(rise/15+Zone(glong)-(int)(rise/15+Zone(glong)10) coutThe time at which the sun rises is (int)(rise/15+Zone(glong):0(int)(60*(rise/15+Zone(glong)-(int)(rise/15+Zone(glong) .n; else coutThe time at which the sun rises is (int)(rise/15+Zone(glong):(int)(60*(rise/15+Zone(glong)-(int)(rise/15+Zone(glong) .n; if(int)(60*(set/15+Zone(glong)-(int)(set/15+Zone(glong)10) coutThe time at which the sun sets is (int)(set/15+Zone(glong): (int)(60*(set/15+Zone(glong)-(int)(set/15+Zone(glong) .n; else coutThe time at which the sun sets is (int)(set/15+Zone(glong):(int)(60*(set/15+Zone(glong)-(int)(set/15+Zone(glong) .n; /打印结果 int main() long double UTo=180.0; int year,month,date; long double glat,glong; int c3; input_date(c); year=c0; month=c1; date=c2; input_glat(c); glat=c0+c1/60+c2/3600; input_glong(c); glong=c0+c1/60+c2/3600; long double rise,set; rise=result_rise(UT_rise(UTo,GHA(UTo,G_sun(t_century(days(year,month,date),UTo),ecliptic_longitude(L_sun(t_century(days(year,month,date),UTo),G_sun(t_century(days(year,month,date),UTo),glong,e(h,glat,sun_deviation(earth_tilt(t_century(days(year,month,date),UTo),ecliptic_longitude(L_sun(t_century(days(year,month,date),UTo),G_sun(t_century(days(year,month,date),UTo),UTo,glong,glat,year,month,date); set=result_set(UT_set(UTo,GHA(UTo,G_sun(t_century(days(year,month,date),UTo),ecliptic_longitude(L_sun(t_century(days(year,month,date),UTo),G_sun(t_century(days(year,month,date),UTo),glong,e(h,glat,sun_deviation(earth_tilt(t_century(days(year,month,date)