日出日落时间计算程序.docxVIP

实用标准文档

//日出日落时间计算C语言程序

#definePI3.1415926#includemath.h#includeiostreamusingnamespacestd;

intdays_of_month_1[]={31,28,31,30,31,30,31,31,30,31,30,31};intdays_of_month_2[]={31,29,31,30,31,30,31,31,30,31,30,31};

longdoubleh=-0.833;

//定义全局变量

voidinput_date(intc[]){inti;

coutEnterthedate(form:20090310):endl;for(i=0;i3;i++){

cinc[i];

}

}

//输入日期

voidinput_glat(intc[]){inti;

coutEnterthedegreeoflatitude(range:0°-60°,form:404040(means40°40′40″)):endl;

for(i=0;i3;i++){

文案大全

实用标准文档

cinc[i];

}

}

//输入纬度

voidinput_glong(intc[]){inti;

coutEnterthedegreeoflongitude(westisnegativ,form:404040(means40°40′40″)):endl;

for(i=0;i3;i++){cinc[i];

}

}

//输入经度

intleap_year(intyear){

if(((year%400==0)||(year%100!=0)(year%4==0)))return1;elsereturn0;

}

//判断是否为闰年:若为闰年,返回1；若非闰年,返回0

intdays(intyear,intmonth,intdate){inti,a=0;

for(i=2000;iyear;i++){

文案大全

实用标准文档

if(leap_year(i))a=a+366;elsea=a+365;

}

if(leap_year(year)){for(i=0;imonth-1;i++){

a=a+days_of_month_2[i];

}

}

else{

for(i=0;imonth-1;i++){a=a+days_of_month_1[i];

}

}

a=a+date;returna;

}

//求从格林威治时间公元2000年1月1日到计算日天数days

longdoublet_century(intdays,longdoubleUTo){return((longdouble)days+UTo/360)/36525;

}

//求格林威治时间公元2000年1月1日到计算日的世纪数t

文案大全

实用标准文档

longdoubleL_sun(longdoublet_century){return(280.460+36000.770*t_century);

}

//求太阳的平黄径

longdoubleG_sun(longdoublet_century){return(357.528+35999.050*t_century);

}

//求太阳的平近点角

longdoubleecliptic_longitude(longdoubleL_sun,longdoubleG_sun){return(L_sun+1.915*sin(G_sun*PI/180)+0.02*sin(2*G_sun*PI/180));

}

//求黄道经度

longdoubleearth_tilt(longdoublet_century){return(23.4393-0.0130*t_century);

}

//求地球倾角

longdoublesun_deviation(longdoubleearth_tilt,longdoubleecliptic_longitude){return

(180/PI*asin(sin(PI/180*earth_tilt)*sin(PI/180*ecliptic_longitude)));

}

文案大全

实用标准文档

//求太阳偏差

longdoubleGHA(longdoubleUTo,longdoubleG_sun,longdoubleecliptic_longitude){return

(UTo-180-1.915