Sunrise/Sunset AlgorithmSource: Alm

Sunrise/Sunset Algorithm

Source:
Almanac for Computers, 1990
published by Nautical Almanac Office
United States Naval Observatory
Washington, DC 20392

Inputs:
day, month, year: date of sunrise/sunset
latitude, longitude: location for sunrise/sunset
zenith: Sun's zenith for sunrise/sunset
offical = 90 degrees 50'
civil = 96 degrees
nautical = 102 degrees
astronomical = 108 degrees

NOTE: longitude is positive for East and negative for West
NOTE: the algorithm assumes the use of a calculator with the
trig functions in "degree" (rather than "radian") mode. Most
programming languages assume radian arguments, requiring back
and forth convertions. The factor is 180/pi. So, for instance,
the equation RA = atan(0.91764 * tan(L)) would be coded as RA
= (180/pi)*atan(0.91764 * tan((pi/180)*L)) to give a degree
answer with a degree input for L.


1. first calculate the day of the year

N1 = floor(275 * month / 9)
N2 = floor((month + 9) / 12)
N3 = (1 + floor((year - 4 * floor(year / 4) + 2) / 3))
N = N1 - (N2 * N3) + day - 30

2. convert the longitude to hour value and calculate an approximate time

lngHour = longitude / 15

if rising time is desired:
t = N + ((6 - lngHour) / 24)
if setting time is desired:
t = N + ((18 - lngHour) / 24)

3. calculate the Sun's mean anomaly

M = (0.9856 * t) - 3.289

4. calculate the Sun's true longitude

L = M + (1.916 * sin(M)) + (0.020 * sin(2 * M)) + 282.634
NOTE: L potentially needs to be adjusted into the range [0,360) by adding/subtracting 360

5a. calculate the Sun's right ascension

RA = atan(0.91764 * tan(L))
NOTE: RA potentially needs to be adjusted into the range [0,360) by adding/subtracting 360

5b. right ascension value needs to be in the same quadrant as L

Lquadrant = (floor( L/90)) * 90
RAquadrant = (floor(RA/90)) * 90
RA = RA + (Lquadrant - RAquadrant)

5c. right ascension value needs to be converted into hours

RA = RA / 15

6. calculate the Sun's declination

sinDec = 0.39782 * sin(L)
cosDec = cos(asin(sinDec))

7a. calculate the Sun's local hour angle

cosH = (cos(zenith) - (sinDec * sin(latitude))) / (cosDec * cos(latitude))

if (cosH > 1)
the sun never rises on this location (on the specified date)
if (cosH < -1)
the sun never sets on this location (on the specified date)

7b. finish calculating H and convert into hours

if if rising time is desired:
H = 360 - acos(cosH)
if setting time is desired:
H = acos(cosH)

H = H / 15

8. calculate local mean time of rising/setting

T = H + RA - (0.06571 * t) - 6.622

9. adjust back to UTC

UT = T - lngHour
NOTE: UT potentially needs to be adjusted into the range [0,24) by adding/subtracting 24

10. convert UT value to local time zone of latitude/longitude

localT = UT + localOffset

ผลลัพธ์ (ไทย) 1:

อัลกอริทึมซันไรส์/ซันเซทที่มา: Almanac สำหรับคอมพิวเตอร์ 1990 เผยแพร่ โดยสำนักงาน Almanac ทะเล หอดูดาวกองทัพเรือสหรัฐอเมริกา วอชิงตัน DC 20392อินพุต: วัน เดือน ปี: วันของซันไรส์/ซันเซท ละติจูด ลองจิจูด: ตำแหน่งที่ตั้งสำหรับซันไรส์/ซันเซท สุดยอด: สุดยอดของดวงอาทิตย์สำหรับซันไรส์/ซันเซท อย่างเป็นทาง = 90 องศา 50' แพ่ง = 96 องศา ทะเล = 102 องศา ดาราศาสตร์ = 108 องศา หมายเหตุ: ลองจิจูดเป็นบวกสำหรับตะวันออก และลบในตะวันตก หมายเหตุ: อัลกอริทึมถือว่าการใช้เครื่องคิดเลขพร้อม trig ฟังก์ชันในโหมด "องศา" (แทนที่ "เรเดียน") มากที่สุด เรเดียนอาร์กิวเมนต์ ต้องการกลับไปถือว่าภาษาเขียนโปรแกรม และมา convertions ตัวเป็น 180/pi ดังนั้น เช่น สมการ RA = atan (0.91764 * tan(L)) จะมีโค้ดเป็น RA = (180/pi) * atan (0.91764 * tan((pi/180)*L)) ให้เป็น ตอบ ด้วยระดับสัญญาณขาเข้าแบบสำหรับ L.1. คำนวณวันที่ของปี N1 =ชั้น (275 * เดือน/9) N2 =ชั้น ((month + 9) / 12) N3 = (1 + ชั้น ((ปี - 4 * floor(year / 4) + 2 / 3)) N = N1 - (N2 * N3) + วัน - 302. แปลงลองจิจูดค่าชั่วโมง และคำนวณเป็นเวลาโดยประมาณ lngHour =ลองจิจูด / 15 ถ้าต้องการเวลาเพิ่มขึ้น: t = N + ((6-lngHour) / 24) ถ้าตั้งค่าเวลา ที่ถูกต้อง: t = N + ((18-lngHour) / 24)3. คำนวณอนอมัลลีเฉลี่ยของดวงอาทิตย์ M = (0.9856 * t) - 3.2894. คำนวณลองจิจูดจริงของดวงอาทิตย์ L = M + (1.916 * sin(M)) + (0.020 * sin(2 * M)) + 282.634 หมายเหตุ: L อาจต้องปรับปรุงในช่วง [0,360) โดยการเพิ่ม/ลบ 3605 a. คำนวณของดวงอาทิตย์ไรต์แอสเซนชัน RA = atan(0.91764 * tan(L)) หมายเหตุ: RA อาจต้องปรับปรุงในช่วง [0,360) โดยการเพิ่ม/ลบ 3605b ค่าไรต์แอสเซนชันต้องควอดร้อนท์เดียวเป็น L Lquadrant (ชั้น (L 90)) = * 90 RAquadrant = (floor(RA/90)) * 90 RA = RA + (Lquadrant - RAquadrant)ค่าไรต์แอสเซนชันค. 5 ต้องถูกแปลงเป็นชั่วโมง RA = RA / 156. คำนวณปฏิเสธของดวงอาทิตย์ sinDec = 0.39782 * sin(L) cosDec = cos(asin(sinDec))7a การคำนวณมุมชั่วโมงภายในของดวงอาทิตย์ cosH = (cos(zenith) - (sinDec * sin(latitude))) / (cosDec * cos(latitude)) ถ้า (cosH > 1) พระอาทิตย์ไม่เคยขึ้นที่ตำแหน่งนี้ (ในวันที่ระบุ) ถ้า (cosH < -1) พระอาทิตย์ไม่เคยตกที่ตำแหน่งนี้ (ในวันที่ระบุ)7b เสร็จคำนวณ H และแปลงเป็นชั่วโมง ถ้าหากต้องการเวลาเพิ่มขึ้น: H = 360 - acos(cosH) ถ้าตั้งค่าเวลา ที่ถูกต้อง: H = acos(cosH) H = H / 158. คำนวณเวลาท้องถิ่นของไรซิ่ง / T = H + RA - (0.06571 * t) - 6.6229. ปรับไป UTC UT = T - lngHour หมายเหตุ: UT อาจต้องปรับปรุงในช่วง [0,24) โดยการเพิ่ม/ลบ 2410. ค่า UT แปลงโซนเวลาท้องถิ่นของละติจูด/ลองจิจูด localT = UT + localOffset

ผลลัพธ์ (ไทย) 2:

พระอาทิตย์ขึ้น / ซันเซ็ทอัลกอริทึมที่มา: ปูมสำหรับคอมพิวเตอร์ 1990 ที่ตีพิมพ์โดยปูมการเดินเรือสำนักงานสหรัฐอเมริกานาวีหอคอยกรุงวอชิงตันดีซี 20,392 ปัจจัยการผลิต: วันเดือนปี: วันที่พระอาทิตย์ขึ้น / ละติจูดลองจิจูด: สถานที่สำหรับพระอาทิตย์ขึ้น / สุดยอด: สุดยอดของดวงอาทิตย์ สำหรับพระอาทิตย์ขึ้น / ทางการ = 90 องศา 50 ' ประชา = 96 องศา ทะเล = 102 องศา ดาราศาสตร์ = 108 องศาหมายเหตุ: เส้นแวงเป็นบวกสำหรับตะวันออกและเชิงลบสำหรับเวสต์หมายเหตุ: ขั้นตอนวิธีการจะถือว่าการใช้งานของเครื่องคิดเลขที่มีฟังก์ชั่นในการหนุน "การศึกษาระดับปริญญา "(มากกว่า" เรเดียน ") โหมด ส่วนใหญ่การเขียนโปรแกรมภาษาถือว่าข้อโต้แย้งเรเดียนต้องกลับมา convertions ปัจจัยที่ 180 / ปี่ ดังนั้นสำหรับตัวอย่างเช่นสม RA = Atan (0.91764 * สีน้ำตาล (L)) จะถูกกำหนดเป็น RA = (180 / ปี่) * Atan (0.91764 * สีน้ำตาล ((ปี่ / 180) * L)) เพื่อให้การศึกษาระดับปริญญาคำตอบ ด้วยการป้อนข้อมูลสำหรับการศึกษาระดับปริญญาลิตร1 ครั้งแรกของการคำนวณวันของปีN1 = ชั้น (275 * เดือน / 9) N2 = ชั้น ((เดือน + 9) / 12) N3 = (1 + ชั้น ((ปี - 4 * ชั้น (ปี / 4) + 2) / 3)) N = N1 - (N2 N3 *) + วัน - 30 2. แปลงเส้นแวงเป็นค่าชั่วโมงและคำนวณเวลาโดยประมาณlngHour = แวง / 15 ถ้ามีเวลาเพิ่มขึ้นเป็นที่ต้องการ: t = N + ((6 - lngHour ) / 24) ถ้าการตั้งค่าเวลาที่ต้องการ: t = N + ((18 - lngHour) / 24) 3 คำนวณความผิดปกติของดวงอาทิตย์เฉลี่ยM = (0.9856 * t) - 3.289 4 คำนวณเส้นแวงที่แท้จริงของดวงอาทิตย์L = M + (1.916 * บาป (M)) + (0.020 * บาป (2 * M)) + 282.634 หมายเหตุ: L อาจจะต้องปรับลงในช่วง [0360) โดยการเพิ่ม / ลบ 360 5a . คำนวณขวาขึ้นของดวงอาทิตย์RA = Atan (0.91764 * สีน้ำตาล (L)) หมายเหตุ: RA อาจจะต้องปรับลงในช่วง [0360) โดยการเพิ่ม / ลบ 360 5b ค่าขวาขึ้นจะต้องอยู่ในด้านเดียวกับ L = Lquadrant (ชั้น (L / 90)) * 90 = RAquadrant (ชั้น (RA / 90)) * 90 = RA RA + (Lquadrant - RAquadrant) 5c ค่าขวาขึ้นจะต้องมีการแปลงเป็นชั่วโมงRA = RA / 15 6. การคำนวณของดวงอาทิตย์ปฏิเสธsinDec = 0.39782 * บาป (L) cosDec = cos (สิน (sinDec)) 7a คำนวณมุมชั่วโมงในท้องถิ่นของดวงอาทิตย์กระบอง = (cos (สุดยอด) - (sinDec * บาป (ละติจูด))) / (cosDec * cos (ละติจูด)) ถ้า (กระบอง> 1) ดวงอาทิตย์ไม่เคยขึ้นบนสถานที่นี้ (ในวันที่ระบุ ) ถ้า (กระบอง <-1) ดวงอาทิตย์ไม่เคยตั้งในสถานที่นี้ (ในวันที่ระบุ) 7b เสร็จสิ้นการคำนวณ H และแปลงเป็นชั่วโมงหากถ้ามีเวลาเพิ่มขึ้นเป็นที่ต้องการ: H = 360 - acos (กระบอง) ถ้าการตั้งค่าเวลาที่ต้องการ: H = acos (กระบอง) H = H / 15 8. คำนวณเวลาเฉลี่ยในท้องถิ่นของการเพิ่มขึ้น / การตั้งค่าT = H + RA - (0.06571 * t) - 6.622 9 กลับไปปรับ UTC UT = T - lngHour หมายเหตุ: UT อาจจะต้องปรับลงในช่วง [0,24) โดยการเพิ่ม / ลบ 24 10. แปลงค่า UT โซนเวลาท้องถิ่นของละติจูด / ลองจิจูดlocalT = UT + localOffset

ผลลัพธ์ (ไทย) 3:

เวลาพระอาทิตย์ขึ้น / พระอาทิตย์ตกขั้นตอนวิธี

ที่มา :
ปูมสำหรับคอมพิวเตอร์ 1990 เผยแพร่โดยสำนักงานปูมการเดินเรือ

สหรัฐอเมริกานาวีหอคอย
Washington , DC 20392

กระผม :
วันที่ ปี เดือน วันพระอาทิตย์ขึ้น / พระอาทิตย์ตก
ละติจูด , ลองจิจูด : สถานที่สำหรับพระอาทิตย์ขึ้น / พระอาทิตย์ตก
สุดยอด : ซันสุดยอดสำหรับพระอาทิตย์ขึ้นพระอาทิตย์ตก /
อย่างเป็นทางการ = 90 องศา 50 '
แพ่ง = 96 องศา

= 102 องศาลมดาราศาสตร์ = 108 องศา

หมายเหตุ : เส้นแวงเป็นบวก และตะวันออกเป็นตะวันตก
หมายเหตุ : ขั้นตอนวิธีจะใช้เครื่องคิดเลขกับ
หนุนการทำงานใน " ปริญญา " ( ไม่ใช่ " เรเดียน " ) โหมด ภาษาโปรแกรมส่วนใหญ่
ถือว่าเรเดียนอาร์กิวเมนต์ต้องกลับมา
มา convertions . ปัจจัยคือ 180 / PI ดังนั้นตัวอย่างเช่น
สมการ รา = ATAN ( 091764 * ตาล ( L ) จะถูกเข้ารหัสเป็นรา
= ( 180 / PI ) * ATAN ( 0.91764 * แทน ( ( pi / 180 ) * L ) ) เพื่อให้ระดับ
ตอบกับปริญญาที่ใส่ L


1 แรกคำนวณวัน

( / 9 ชั้น 1 = 275 * เดือน )
2 = ชั้น ( เดือน 9 ) / 12 )
3 = 1 ชั้น ( ( ปี 4 * ชั้น ( ปี 4 ) 2 ) 3 )
n = 1 - ( 2 * 3 ) วันนี้ - 30

2 แปลงเส้นลองจิจูดค่าชั่วโมงคำนวณ

เวลาโดยประมาณlnghour = ลองจิจูด / 15

ถ้าเพิ่มขึ้นเวลาที่ต้องการ :
t = N ( 6 - lnghour ) / 24 )
ถ้าตั้งค่าเวลาที่ต้องการ :
t = N ( 18 - lnghour ) / 24 )

3 คำนวณดวงอาทิตย์หมายถึงมิติ

M = ( 0.9856 * t ) - 3.289

4 . คำนวณดวงอาทิตย์จริงลองจิจูด

L = M ( 1.916 * sin ( M ) ( 0.020 * sin ( 2 * m ) 282.634
หมายเหตุ : ผมอาจต้องมีการปรับเปลี่ยนในช่วง [ 0 ,360 ) โดยการเพิ่ม / ลบ 360

5a . คำนวณพระอาทิตย์ขึ้นเยื้องไปทางขวา

รา = ATAN ( 0.91764 * ตาล ( L )
หมายเหตุ : ราอาจต้องมีการปรับเปลี่ยนในช่วง [ 0360 ) โดยการเพิ่ม / ลบ 360

5B ไรต์แอสเซนชันค่าต้องอยู่ในเครือข่ายเดียวกันฉัน

lquadrant = ( ชั้น ( L / 90 ) ) * 90 = (
raquadrant พื้น ( Ra / 90 ) ) * 90 =
รา รา ( lquadrant - raquadrant )

5 .ไรต์แอสเซนชันมีค่าต้องแปลงเป็นชั่วโมง

รา = รา / 15

6 คำนวณดวงอาทิตย์ปฏิเสธ

sindec = 0.39782 * sin ( L )
cosdec = cos ( เกลือ ( sindec ) )

7A . คำนวณดวงอาทิตย์เป็นมุมชั่วโมงท้องถิ่น

Cos = ( cos ( Zenith ) - ( sindec * sin ( ละติจูด ) ) ) ) ) ) ) / ( cosdec * cos ( ละติจูด )

ถ้า ( Cos > 1 )
พระอาทิตย์ไม่เคยขึ้นในตำแหน่งนี้ ( วันที่กำหนด )
ถ้า ( Cos < - 1 )
ดวงอาทิตย์ไม่เคยตกดินในตำแหน่งนี้ ( ในวันที่ที่ระบุไว้ )

7b จบคำนวณ H และแปลงเป็นชั่วโมง

ถ้าเพิ่มขึ้นเวลาที่ต้องการ :
H = 360 - ACOS ( Cos )
ถ้าตั้งค่าเวลาที่ต้องการ :
H = ACOS ( Cos )

H
= H / 15
8 คำนวณเวลาที่เพิ่มขึ้น / การตั้งค่าท้องถิ่นหมายถึง
t = H รา - ( 0.06571 * t ) - 6.622

9 ปรับกลับไปอยู่

แต่ = t - lnghour
หมายเหตุ :แต่อาจต้องมีการปรับเปลี่ยนในช่วง [ 0,24 ) โดยการเพิ่ม / ลบ 24

10 แปลงโซนเวลาท้องถิ่นแต่ค่าของละติจูด / ลองจิจูด

localt = UT localoffset