1二十四節(jié)氣如何計(jì)算?立春日期的計(jì)算：計(jì)算公式：Y*D+C-L公式解讀：年數(shù)的后2位乘0.2422加3.87取整數(shù)減閏年數(shù)。21世紀(jì)C值=3.87，22世紀(jì)C值=4.15。立春日期的計(jì)算計(jì)算公式：Y*D+C-L公式解讀：年數(shù)的后2位乘0.2422加3.87取整數(shù)減閏年數(shù)。21世紀(jì)C值=3.87，22世紀(jì)C值=4.15。舉例說明：2058年立春日期的計(jì)算步驟58X.0.2422+3.87-(58-1)/4=17-14=3，則2月3日立春。雨水日期的計(jì)算Y*D+C-L公式解讀：年數(shù)的后2位乘0.2422加18.74取整數(shù)減閏年數(shù)。21世紀(jì)雨水的C值18.73。舉例說明：2008年雨水日期=8X.0.2422+18.73-(8-1)/4=20-1=19,2月19日雨水。例外：2026年計(jì)算得出的雨水日期應(yīng)調(diào)減一天為18日。驚蟄日期的計(jì)算Y*D+C-L公式解讀：年數(shù)的后2位乘0.2422加5.63取整數(shù)減閏年數(shù)。21世紀(jì)驚蟄的C值=5.63。舉例說明：2088年驚蟄日期=88X.0.2422+5.63-88/4=26-22=4,3月4日是驚蟄。例外：無。春分日期的計(jì)算Y*D+C-L公式解讀：年數(shù)的后2位乘0.2422加20.646取整數(shù)減閏年數(shù)。21世紀(jì)春分的C值=20.646。舉例說明：2092年春分日期=92X.0.2422+20.646-92/4=42-23=19,3月19日是春分。例外：2084年的計(jì)算結(jié)果加1日。清明節(jié)日期的計(jì)算Y*D+C-L公式解讀：丫丫= =年數(shù)后2位，D=0.2422,L=閏年數(shù)，21世紀(jì)C=4.81,20世紀(jì)=5.59。2舉例說明：2088年清明日期=88X.0.2422+4.81-88/4=26-22=4,4月4日是清明。例外：無。谷雨節(jié)日期的計(jì)算Y*D+C-L公式解讀：丫丫= =年數(shù)后2位，D=0.2422，L=閏年數(shù)，21世紀(jì)C=20.1，20世紀(jì)=20.888。 舉例說明：2088年谷雨日期=88X.0.2422+20.1-88/4=41-22=19,4月19日是谷例外：無。立夏日期的計(jì)算Y*D+C-L公式解讀：丫丫= =年數(shù)后2位，D=0.2422,L=閏年數(shù)，21世紀(jì)C=5.52,20世紀(jì)=6.318。舉例說明：2088年立夏日期=88X.0.2422+5.52-88/4=26-22=4,5月4日是立夏。例外：1911年的計(jì)算結(jié)果加1日。小滿日期的計(jì)算Y*D+C-L公式解讀：丫丫= =年數(shù)后2位，D=0.2422,L=閏年數(shù)，21世紀(jì)C=21.04,20世紀(jì)=21.86。 舉例說明：2088年小滿日期=88X.0.2422+21.04-88/4=42-22=20,5月20日小滿。例外：2008年的計(jì)算結(jié)果加1日。芒種日期的計(jì)算Y*D+C-L公式解讀：丫丫= =年數(shù)后2位，D=0.2422,L=閏年數(shù)，21世紀(jì)C=5.678,20世紀(jì)=6.5。 舉例說明：2088年芒種日期=88X.0.2422+5.678-88/4=26-22=4,6月4日芒種。例外：1902年的計(jì)算結(jié)果加1日。夏至日期的計(jì)算Y*D+C-L公式解讀：丫丫= =年數(shù)后2位，D=0.2422,L=閏年數(shù)，21世紀(jì)C=21.37,20世紀(jì)=22.20。 舉例說明：2088年夏至日期=88X.0.2422+21.37-88/4=42-22=20,6月20日夏至。例外：1928年的計(jì)算結(jié)果加1日。小暑日期的計(jì)算Y*D+C-L公式解讀：丫丫= =年數(shù)后2位，D=0.2422,L=閏年數(shù)，21世紀(jì)C=7.108,20世紀(jì)=7.928。 舉例說明：2088年小暑日期=88X0.2422+7.108-88/4=28-22=6,7月6日是小暑。例外：1925年和2016年的計(jì)算結(jié)果加1日。大暑日期的計(jì)算Y*D+C-L公式解讀：丫丫= =年數(shù)后2位，D=0.2422,L=閏年數(shù)，21世紀(jì)C=22.83,20世紀(jì)=23.65。 舉例說明：2088年大暑日期=88X0.2422+22.83-88/4=44-22=22,7月22日大暑。例外：1922年的計(jì)算結(jié)果加1日。立秋日期的計(jì)算Y*D+C-L公式解讀：丫丫= =年數(shù)后2位，D=0.2422,L=閏年數(shù),21世紀(jì)C=7.5,20世紀(jì)=8.35。 舉例說明：2088年立秋日期=88X0.2422+7.5-88/4=28-22=6,8月6日是立秋。例外：2002年的計(jì)算結(jié)果加1日。雨。3處暑日期的計(jì)算Y*D+C-L公式解讀：丫丫= =年數(shù)后2位，D=0.2422,L=閏年數(shù)，21世紀(jì)C=23.13,20世紀(jì)=23.95。舉例說明：2088年處暑日期=88X0.2422+23.13-88/4=44-22=22,8月22日處暑。例外：無。白露日期的計(jì)算Y*D+C-L公式解讀：丫丫= =年數(shù)后2位，D=0.2422,L=閏年數(shù)，21世紀(jì)C=7.646,20世紀(jì)=8.44。 舉例說明：2088年白露日期=88X0.2422+7.646-88/4=28-22=6,9月6日是白露。例外：1927年的計(jì)算結(jié)果加1日。秋分日期的計(jì)算Y*D+C-L公式解讀：丫丫= =年數(shù)后2位，D=0.2422,L=閏年數(shù)，21世紀(jì)C=23.042,20世紀(jì)=23.822。舉例說明：2088年秋分日期=8X.0.2422+23.042-88/4=44-22=22,9月22日是秋分。例外：1942年的計(jì)算結(jié)果加1日。寒露日期的計(jì)算Y*D+C-L公式解讀：丫丫= =年數(shù)后2位，D=0.2422,L=閏年數(shù)，21世紀(jì)C=8.318,20世紀(jì)=9.098。舉例說明：2088年寒露日期=88X0.2422+8.318-88/4=29-22=7,10月7日是寒露。例外：無。霜降日期的計(jì)算Y*D+C-L公式解讀：丫丫= =年數(shù)后2位，D=0.2422,L=閏年數(shù)，21世紀(jì)C=23.438,20世紀(jì)=24.218。 舉例說明：2088年霜降日期=88X0.2422+23.438-88/4=44-22=22,10月22日霜降。例外：2089年的計(jì)算結(jié)果加1日。立冬日期的計(jì)算Y*D+C-L公式解讀：丫丫= =年數(shù)后2位，D=0.2422,L=閏年數(shù)，21世紀(jì)C=7.438,20世紀(jì)=8.218。舉例說明：2088年立冬日期=88X0.2422+7.438-88/4=28-22=6,11月6日是立冬。例外：2089年的計(jì)算結(jié)果加1日。小雪日期的計(jì)算Y*D+C-L公式解讀：丫丫= =年數(shù)后2位，D=0.2422,L=閏年數(shù)，21世紀(jì)C=22.36,20世紀(jì)=23.08。舉例說明：2088年小雪日期=88X0.2422+22.36-88/4=43-22=21,11月21日小雪。例外：1978年的計(jì)算結(jié)果加1日。大雪日期的計(jì)算Y*D+C-L公式解讀：丫丫= =年數(shù)后2位，D=0.2422,L=閏年數(shù),21世紀(jì)C=7.18,20世紀(jì)=7.9。 舉例說明：2088年大雪日期=88X0.2422+7.18-88/4=28-22=6,12月6日大雪。例外：1954年的計(jì)算結(jié)果加1日。4冬至日期的計(jì)算Y*D+C-L公式解讀：丫丫= =年數(shù)后2位，D=0.2422,L=閏年數(shù)，21世紀(jì)C=21.94,20世紀(jì)=22.60。舉例說明：2088年冬至日期=88X0.2422+21.94-88/4=43-22=21,12月21日冬至。例外：1918年和2021年的計(jì)算結(jié)果減1日。小寒日期的計(jì)算Y*D+C-L公式解讀：丫丫= =年數(shù)后2位，D=0.2422，L=閏年數(shù)，21世紀(jì)C=5.4055，20世紀(jì)=6.11。舉例說明：1988年小寒日期=88X.0.2422+6.11-(88-1)/4=27-21=6,1月6日小寒。例外：1982年計(jì)算結(jié)果加1日,2019年減1日。大寒日期的計(jì)算Y*D+C-L公式解讀：丫丫= =年數(shù)的后2位，D=0.2422,L=閏年數(shù)，21世紀(jì)C=20.12,20世紀(jì)C=20.84。 舉例說明：2089年大寒日期=89X0.2422+20.12-(89-1)/4=41-22=19,1月19日大寒。例外：2082年的計(jì)算結(jié)果加1日,20世紀(jì)無。農(nóng)歷 24 節(jié)氣算法/許劍偉/天星歷算(2012-06-17 18:41:46)轉(zhuǎn)載標(biāo)簽: 八字命理八字算命天文歷法星座農(nóng)歷 24 節(jié)氣算法許劍偉摘要古老而又現(xiàn)代的中國農(nóng)歷，是一種天文學(xué)性質(zhì)的陰陽歷。本文采用 VSOP87 星歷算 法并結(jié)合運(yùn)動(dòng)學(xué)方法以及牛頓求根法得到 24 節(jié)氣的精確時(shí)間，揭開農(nóng)歷計(jì)算的神秘面紗。關(guān)鍵字農(nóng)歷算法、星歷、節(jié)氣正文計(jì)算中國農(nóng)歷，首先要計(jì)算岀二十四節(jié)氣時(shí)刻。在計(jì)算機(jī)問世之前，二十四節(jié)氣的許 算是非常復(fù)雜的。隨著計(jì)算機(jī)及互聯(lián)網(wǎng)的普及，美國航空航天局、法國巴黎天文臺(tái)各自在網(wǎng) 絡(luò)上發(fā)布了精密星歷表的計(jì)算方法，這使得民間計(jì)算農(nóng)歷成為可能。本文以法國巴黎天文臺(tái) 的 VSOP87 算法為基礎(chǔ)，給出中國農(nóng)歷的二十四節(jié)氣算法。在農(nóng)歷中，太陽黃經(jīng)為 0 度時(shí)，對應(yīng)春風(fēng)節(jié)氣。相鄰節(jié)氣對應(yīng)的太陽黃經(jīng)相差 15 度。一周 年內(nèi)，太陽黃經(jīng)從 0 度變化到360 度，共有 24 個(gè)節(jié)氣。一、時(shí)間標(biāo)尺一一儒略日數(shù)計(jì)算計(jì)算星歷之前首先要解決時(shí)間尺問題。公歷規(guī)定平年 365 日，閏年 366 日。1582 年 10 月 4 日以前，公歷規(guī)定每 4 年設(shè)置5一個(gè)閏年，平均年長度 365.25 天，這期間的公歷稱為儒略 歷。在 1582 年 10 月 15 日之后實(shí)行格里高利歷，規(guī)定每 400年 97 閏，平均年長度為 365.2425 天。由于儒略歷存在嚴(yán)重的 多閏”問題，到了 1582 年，公歷跑快了 10 天左右，當(dāng)時(shí)就人為調(diào)整 了 10 天，并從此實(shí)行格里歷。因此務(wù)必注意 1582 年 10 月 4 日(儒略歷)的下一日為 15 82 年 10 月 15 日(格里歷)。就是說 1582 年 10 月份少了 10天。在儒略歷中，能被 4 整除的年份為閏年，這一年有 366 天，其它年份為平年(365 天)。如 900 年和 1236 年為 閏年，而 750 年和 1429 年為平年。格里高利歷法也采用這一規(guī)則， 但下列年份除外： 不能被 100 整除的年份為平年， 如 1700 年， 1800 年， 1900 年和 2100 年。其余能被 400 整除的年份則為閏年，如 1600 年，2000 年和 2400 年。儒略日數(shù)(簡稱儒略日)：儒略日數(shù)是指從公元-4712 年開始連續(xù)計(jì)算日數(shù)得出的天數(shù)及不滿一日的小數(shù)，通常記為JD (*)。傳統(tǒng)上儒略日的計(jì)數(shù)是從格林尼治平午，即世界時(shí) 12 點(diǎn)開始的。若以力學(xué)時(shí)(或 歷書時(shí))為標(biāo)尺，這種計(jì)數(shù)通常表達(dá)為儒略歷書日,即 JDE (*)，其中 E 只是一種表征， 即按每天 86400 個(gè)標(biāo)準(zhǔn)秒長嚴(yán)格地計(jì)日。例如：1977 年 4 月 26.4日UT =JD 2443259.91977 年 4 月 26.4日TD =JDE 2443259.9儒略日的計(jì)算：設(shè) Y 為給定年份，M為月份 ，D 為該月 日期(可以帶小數(shù))。若 M 2，Y 和 M 不變，若 M =1 或 2，以 Y -代 Y,以 M+12 代 M ,換句話說，如果日 期在 1 月或 2 月，則被看作是在前一年的 13 月或 14 月。對格里高利歷有：A = INT ( Y/100 ) B = 2 - A + INT(A/4)對儒略歷，取 B = 0儒略日即為：JD = INT(365.25(Y+4716)+INT(30.6001(M+1)+D+B-1524.5使用數(shù)值 30.6 取代 30.6001 才是正確的，但我們?nèi)允褂?30.6001，以確保總能取得恰當(dāng)?shù)恼麛?shù)。事實(shí)上可用 30.601 甚至 30.61 來取代 30.6001。例如，5 乘 30.6 精確等于 153，然而大多數(shù)計(jì)算機(jī)不能精確表示出 30.6，這導(dǎo)致得出一個(gè) 152.9999998 的結(jié)果，它的整數(shù)部分為 152，如此算出的 JD 就不正確了。由儒略日推算歷日：將 JD 加上 0.5，令 Z 為其整數(shù)部分，F(xiàn) 為尾數(shù)(小數(shù))部分。若 Z 2299161 ，取 A = Z若 Z 大于等于 2299161，計(jì)算a=INT(Z -1867216.25)/36524.25)A=Z+1+a-INT(a/4)然后計(jì)算B = A+1524C = INT(B-122.1)/365.25)D = INT(365.25C)E = INT(B-D)/30.6001)該月 日期(帶小數(shù)部分)則為：d = B - D - INT(30.6001E)+ F月份 m 為：IF E 2 THEN y = C - 4716IF m =1 or m=2 THEN y = C - 4715這個(gè)公式里求 E 時(shí)用的數(shù) 30.6001 不能代之以 30.6 ,哪怕計(jì)算機(jī)沒有先前所說的問題。否 則，你得到的結(jié)果會(huì)是 2 月 0日而不是 1 月 31 日，或者 4 月 0 日而不是 3 月 31 日。值得記住的一個(gè)常數(shù)是：2000 年 1 月 1 日 12:00:00 的儒略日數(shù)是 J2000= 2451545二、力學(xué)時(shí)與世界時(shí)的差值（deltat T）計(jì)算一般的，可以把手表時(shí)（UTC ）近似看作世界時(shí)（UT），二者的主要差別在于時(shí)區(qū)。如北京 手表時(shí) 8 點(diǎn)對應(yīng)世界時(shí) 0 點(diǎn)。世界時(shí)與地球自轉(zhuǎn)嚴(yán)格同步，但有趣的是，我們的手表時(shí)實(shí) 際上稱為協(xié)調(diào)世界時(shí)，它的秒長是原子鐘的秒長，由于地球自轉(zhuǎn)速度不均勻，時(shí)快時(shí)慢，這 就注定手表時(shí)與地球自轉(zhuǎn)不完全同步?，F(xiàn)在，地球自轉(zhuǎn)速度正在變慢，我們不得不在某些年份的年末把手表撥慢 1 秒，使得手表時(shí)更好的與地球自轉(zhuǎn)同步，并美言為 跳秒。力學(xué)時(shí)是 根據(jù)太陽系的動(dòng)力學(xué)原理導(dǎo)岀的，是一種均勻的時(shí)間系統(tǒng)，其秒長與原子鐘的秒長相同。因 此，協(xié)調(diào)世界時(shí)（UTC ）與世界時(shí)（記為 UT ）其本同步，但力學(xué)時(shí)（記作 TD ）與世界時(shí)不 太同步，二者的差值記作 deltat T 或記作 T。利用直接的天文觀測可以得知每年的T， 利用古代的日月食觀測資料可以反推古代的厶 T。所有年份的 T 計(jì)算出來后，可以擬合出 以下多項(xiàng)式表達(dá)，使得T 的計(jì)算更快捷，計(jì)算結(jié)果的單位是秒。我們利用下表可以嚴(yán)格計(jì)算 T （即 T =TD - UT ）年份 a b c d-4000,108371.7,-13036.80,392.000,0.0000-500,17201.0, -627.82,16.170,-0.3413-150,12200.6, -346.41, 5.403,-0.1593150, 9113.8, -328.13,-1.647,0.0377500, 5707.5, -391.41, 0.915,0.3145900, 2203.4, -283.45,13.034,-0.17781300,490.1, -57.35, 2.085,-0.00721600,120.0,-9.81,-1.532,0.14031700,10.2, -0.91, 0.510,-0.03701800,13.4, -0.72, 0.202,-0.01931830,7.8,-1.81, 0.416,-0.02471860,8.3,-0.13,-0.406,0.02921880,-5.4,0.32,-0.183,0.01731900,-2.3,2.06, 0.169,-0.01351920,21.2,1.69,-0.304,0.01671940,24.2,1.22,-0.064,0.00311960,33.2,0.51, 0.231,-0.01091980,51.0,1.29,-0.026,0.00322000,63.87,0.1,0, 0,2005表中每一行適用一定的年代范圍，如第 1 行適用于公元-4000 年到-500 年，第 2 行適用于 公元-500 到-1500 年，其它類推。每行的起始年份記作 Y1，終止年份記作 Y2，如果年份 y 在 Y1 到 Y2 之 間，那么該年的 deltat T 表達(dá)為： T = a + b*t1 + c*t2 + d*t3，單位是秒其中 t1 =( y-Y1 ) / ( Y2-Y1 ) *10 , t2 = t1*t1 , t3 = t1*t1*t1對于 2005 年以后的 deltat T 是未知的，要做外推計(jì)算：2005 至 2014 年建議使用 1995 到 2005 年期間 T 的平均增速計(jì)算，即：7 T = F(y) = 64.7+ (y-2005)* b， 其中速度 b = 0.42114 年以后可以使用二次曲線外推 T = f(y) = -20+ a * (y-1820)/100A2，其中加速度 a = 312114 年到 2014 年之間的外推， 可以在上面兩個(gè)外推算式的基礎(chǔ)上做一次的曲線連接， 使 之連續(xù)即可。 比如可以這么計(jì)算： T = f(y) + (y-2114)* f(2014)- F(2014)/100以下數(shù)值可供程序驗(yàn)證參考2008 年 T = 66.0 秒1950 年 T = 29 秒500 年T = 5710 秒三、太陽視黃經(jīng)(真分點(diǎn)視坐標(biāo))算法基于 VSOP87 半解析法。力學(xué)時(shí) t 為 J2000 起算的儒略世紀(jì)數(shù)，t2 = t*t，t3 = t2*t，t4 = t3*tA、 低精度算法L0(t)= .66+ 53.318*t 弧度B、 中精度算法L1(t)= .66+ 53.318*t+ 53*t*t+ 334116*cos(4.67+628.307585*t)+ 2061*cos(2.678+628.3076*t)*t / 弧度C、 高精度算法L2(t)= .66+ 53.318*t+ 52.9674*t2+ 0.00432*t3- 0.001124*t4+334166* cos( 4.669257+ 628.307585*t)+3489* cos(4.6261+1256.61517*t)+350*cos(2.744 + 575.3385*t)+342*cos(2.829 +0.3523*t)+314*cos(3.628 +7771.3771*t)+268*cos(4.418 + 786.0419*t)+234*cos(6.135 + 393.021*t)+132*cos(0.742 +1150.677*t)+127*cos(2.037 + 52.9691*t)+120*cos(1.11+157.7344*t)+ 99* cos(5.23+ 588.493*t)+ 90*cos(2.05+2.63*t )+ 86*cos(3.51+ 39.815*t)+ 78*cos(1.18+ 522.369*t)+ 75*cos(2.53+ 550.755*t)+ 51*cos(4.58+1884.923*t)+ 49*cos(4.21+ 77.552*t)+ 36*cos(2.92+0.07*t )+ 32*cos(5.85+1179.063*t)+ 28* cos( 1.9+ 79.63*t )+ 27* cos( 0.31+1097.71*t )+2060.6* cos( 2.67823+ 628.307585*t) * t+43.0* cos( 2.635 +1256.6152*t )* t+8.72* cos( 1.072 + 628.3076*t )* t28-994- 834* sin(2.1824-33.75705*t)-64 * sin(3.5069+1256.66393*t) / 弧度最后兩行分別為光行差和章動(dòng)四、太陽黃經(jīng)速度平速度：v0 = 628.18即時(shí) 速度：v1 = 628.332+21* sin(1.527+628.307585*t)速度的單位是 弧度/儒略世紀(jì)”即 弧度/36525 天”注意，平速度比即時(shí)速度的精度要高得多，務(wù)必保留足夠的有效數(shù)字，否則將帶來嚴(yán)重的計(jì) 算誤差。五、節(jié)氣時(shí)刻計(jì)算以上天體黃經(jīng)時(shí)間的函數(shù)，即 L = f(t)，所謂的求節(jié)氣時(shí)刻就是已知 L 求 t，顯然這是在求 解一個(gè)關(guān)于 t 的方程。偉大的英國天文學(xué)家物理學(xué)家牛頓給岀了一種非常有效的迭代算法：牛頓求根法。用這種方法，求 t 所花費(fèi)的時(shí)間僅是求 f(t)花費(fèi)時(shí)間的 1.21.3 倍。設(shè)某個(gè) 節(jié)氣對應(yīng)的黃經(jīng)為 W，那么算法如下。牛頓迭代算法設(shè)計(jì)：第 1 步迭代：t= 0第2 步迭代：t= t+( W- L0(t)/ v0第3 步迭代：t= t+( W- L1(t)/ v1(t)第4 步迭代：t= t+( W- L2(t)/ v1(t)誤差：算法誤差 2 分鐘以內(nèi)，實(shí)際找到的誤差一般在 30 秒以內(nèi)，平均 15 秒注意：W 指的是太陽黃經(jīng)。1999 年春分對應(yīng) W=0，以后每 W 每增加 15 度對應(yīng)下一個(gè)節(jié)氣。 迭代的的結(jié)果是力學(xué)時(shí)，單位是儒略世紀(jì)數(shù)。最后結(jié)果還應(yīng)轉(zhuǎn)換為北京時(shí)間，即：JD = J 2000+ t*36525- T/86400+ 8/24最后使用儒略日數(shù)轉(zhuǎn)公歷”所述方法得到節(jié)氣的日期和時(shí)間。六、計(jì)算結(jié)果比較為了進(jìn)行誤差比較，下文列岀 2007 年的 24 節(jié)氣，并與壽星天文歷比對。壽星天文歷 是筆者制作的一款精度優(yōu)于 1 秒的農(nóng)歷工具，已發(fā)布于互聯(lián)網(wǎng)上，其算法與本文類似。節(jié)氣本文算法壽星天年歷春分 2007-03-2108:06:59 08:07:26清明 2007-04-0512:04:21 12:04:39谷雨 2007-04-2019:06:40 19:07:04立夏 2007-05-0605:20:10 05:20:23小滿 2007-05-2118:11:45 18:11:56芒種 2007-06-0609:27:02 09:27:04夏至 2007-06-2202:06:22 02:06:25小暑 2007-07-0719:41:48 19:41:42大暑 2007-07-2313:00:13 13:00:10立秋 2007-08-0805:31:31 05:31:14節(jié)氣本文算法壽星天年歷秋分 2007-09-2317:51:31 17:51:14寒露 2007-10-0900:12:00 00:11:31霜降 2007-10-2403:15:43 03:15:24立冬 2007-11-0803:24:24 03:24:00小雪 2007-11-2300:50:03 00:49:52大雪 2007-12-0720:14:11 20:14:04冬至 2007-12-2214:07:48 14:07:47小寒 2008-01-0607:24:43 07:24:49大寒 2008-01-2100:43:28 00:43:30立春 2008-02-0419:00:09 19:00:229處暑 2007-08-2320:08:07 20:07:58白露 2007-09-0808:29:56 08:29:29雨水 2008-02-1914:49:27 14:49:32驚蟄 2008-03-0512:58:26 12:58:47
本站僅提供存儲(chǔ)服務(wù)，所有內(nèi)容均由用戶發(fā)布，如發(fā)現(xiàn)有害或侵權(quán)內(nèi)容，請
點(diǎn)擊舉報(bào)。
