已知：日出日落時(shí)太陽(yáng)的位置h＝-0.833°,要計(jì)算地的地理位置,經(jīng)度Long,緯度G1at,時(shí)區(qū)zone,UTo為上次計(jì)算的日出日落時(shí)間,第一次計(jì)算時(shí)UTo＝180°.
(1)先計(jì)算出從格林威治時(shí)間公元2000年1月1日到計(jì)算日天數(shù)days;
(2)計(jì)算從格林威治時(shí)間公元2000年1月1日到計(jì)算日的世紀(jì)數(shù)t,則t＝(daysUTo／360)／36525;
(3)計(jì)算太陽(yáng)的平黃徑L=280.46036000.770×t;
(4)計(jì)算太陽(yáng)的平近點(diǎn)角G＝357.52835999.050×t
(5)計(jì)算太陽(yáng)的黃道經(jīng)度λ＝L1.915×sinG0.020xsin(2G);
(6)計(jì)算地球的傾角ε＝23.4393-0.0130×t;
(7)計(jì)算太陽(yáng)的偏差δ＝arcsin(sinε×sinλ);
(8)計(jì)算格林威治時(shí)間的太陽(yáng)時(shí)間角GHA：GHA=UTo-180-1.915×sinG-0.020×sin(2G)2.466×sin(2λ)-0.053×sin(4λ)
(9)計(jì)算修正值e：e=arcos{[sinh-sin(Glat)sin(δ)]/cos(Glat)cos(δ)}
(10)計(jì)算新的日出日落時(shí)間UT＝UTo-(GHALong±e);其中“”表示計(jì)算日出時(shí)間,“-”表示計(jì)算日落時(shí)間;
(11)比較UTo和UT之差的絕對(duì)值,如果大于0.1°即0.007小時(shí),把UT作為新的日出日落時(shí)間值,重新從第(2)步開始進(jìn)行迭代計(jì)算,如果UTo和UT之差的絕對(duì)值小于0.007小時(shí),則UT即為所求的格林威治日出日落時(shí)間;
(12)上面的計(jì)算以度為單位,即180°=12小時(shí),因此需要轉(zhuǎn)化為以小時(shí)表示的時(shí)間,再加上所在的時(shí)區(qū)數(shù)Zone,即要計(jì)算地的日出日落時(shí)間為T=UT/15Zone上面的計(jì)算日出日落時(shí)間方法適用于小于北緯60°和南緯60°之間的區(qū)域,如果計(jì)算位置為西半球時(shí),經(jīng)度Long為負(fù)數(shù).