/* Almanac v2.2 Commented out 2009-specific event list Jan 27, 2010 -- DAF JavaScript programming and graphics, copyright: Adrian R. Ashford, December 29th, 2008. Used with permission by Sky Publishing Corporation. http://www.skyandtelescope.com/templates/almanac_v2.js Parts of sunrise, sunset, moonrise and moonset from moonup.bas and sunup.bas by Roger W. Sinnott, Sky and Telescope magazine. Free date & time validation scripts: Sandeep Tamhankar stamhankar@hotmail.com http://javascript.internet.com Free dynamic form event script: http://www.hscripts.com Free variable passing in URL script: http://javascript.internet.com Public domain cookie scripts: Dustin Diaz http://www.dustindiaz.com/top-ten-javascript/ */ myImage = new Array("http://media.skyandtelescope.com/designimages/mp0.gif","http://media.skyandtelescope.com/designimages/mp1.gif","http://media.skyandtelescope.com/designimages/mp2.gif","http://media.skyandtelescope.com/designimages/mp3.gif","http://media.skyandtelescope.com/designimages/mp4.gif","http://media.skyandtelescope.com/designimages/mp5.gif","http://media.skyandtelescope.com/designimages/mp6.gif","http://media.skyandtelescope.com/designimages/mp7.gif","http://media.skyandtelescope.com/designimages/mp8.gif","http://media.skyandtelescope.com/designimages/mp9.gif","http://media.skyandtelescope.com/designimages/mp10.gif","http://media.skyandtelescope.com/designimages/mp11.gif","http://media.skyandtelescope.com/designimages/mp12.gif","http://media.skyandtelescope.com/designimages/mp13.gif","http://media.skyandtelescope.com/designimages/mp14.gif","http://media.skyandtelescope.com/designimages/mp15.gif","http://media.skyandtelescope.com/designimages/mp16.gif","http://media.skyandtelescope.com/designimages/mp17.gif","http://media.skyandtelescope.com/designimages/mp18.gif","http://media.skyandtelescope.com/designimages/mp19.gif","http://media.skyandtelescope.com/designimages/mp20.gif","http://media.skyandtelescope.com/designimages/mp21.gif","http://media.skyandtelescope.com/designimages/mp22.gif","http://media.skyandtelescope.com/designimages/mp23.gif","http://media.skyandtelescope.com/designimages/mp24.gif","http://media.skyandtelescope.com/designimages/mp25.gif","http://media.skyandtelescope.com/designimages/mp26.gif","http://media.skyandtelescope.com/designimages/mp27.gif"); thisImage = 0; images = myImage.length - 1; function getParams() { var idx = document.URL.indexOf('?'); var params = new Array(); if (idx != -1) { var pairs = document.URL.substring(idx+1, document.URL.length).split('&'); for (var i=0; i 59) utc_mins = 59; if (utc_hours < 10) utc_hours = "0" + utc_hours; if (utc_month < 10) utc_month = "0" + utc_month; if (utc_day < 10) utc_day = "0" + utc_day; if (UTdate == "now") { document.planets.date_txt.value = utc_month + "/" + utc_day + "/" + utc_year; } else { document.planets.date_txt.value = UTdate; } if (UTtime == "now") { document.planets.ut_h_m.value = utc_hours + ":" + utc_mins; } else { document.planets.ut_h_m.value = UTtime; } planets(); } function back_to_now() { var nowdate = new Date(); var utc_day = nowdate.getUTCDate(); var utc_month = nowdate.getUTCMonth() + 1; var utc_year = nowdate.getUTCFullYear(); zone_comp = nowdate.getTimezoneOffset() / 1440; var utc_hours = nowdate.getUTCHours(); var utc_mins = nowdate.getUTCMinutes(); var utc_secs = nowdate.getUTCSeconds(); utc_mins += utc_secs / 60.0; utc_mins = Math.floor((utc_mins + 0.5)); if (utc_mins < 10) utc_mins = "0" + utc_mins; if (utc_mins > 59) utc_mins = 59; if (utc_hours < 10) utc_hours = "0" + utc_hours; if (utc_month < 10) utc_month = "0" + utc_month; if (utc_day < 10) utc_day = "0" + utc_day; document.planets.date_txt.value = utc_month + "/" + utc_day + "/" + utc_year; document.planets.ut_h_m.value = utc_hours + ":" + utc_mins; /* if (UTdate == "now") { document.planets.date_txt.value = utc_month + "/" + utc_day + "/" + utc_year; } else { document.planets.date_txt.value = UTdate; } if (UTtime == "now") { document.planets.ut_h_m.value = utc_hours + ":" + utc_mins; } else { document.planets.ut_h_m.value = UTtime; } */ planets(); } function resizeDocumentTo(setw,seth) { return window.resizeTo(setw,seth),window.resizeTo(setw*2-((typeof window.innerWidth == 'undefined')?document.body.clientWidth:window.innerWidth),seth*2-((typeof window.innerHeight == 'undefined')?document.body.clientHeight:window.innerHeight)); } function sign_of(x) { return x == 0 ? 0 : (x > 0 ? 1 : -1); } function proper_ang(big) { with (Math) { var tmp = 0; if (big > 0) { tmp = big / 360.0; tmp = (tmp - floor(tmp)) * 360.0; } else { tmp = ceil(abs(big / 360.0)); tmp = big + tmp * 360.0; } } return tmp; } function proper_ang_rad(big) { with (Math) { var tmp = 0; if (big > 0) { tmp = big / 2 / PI; tmp = (tmp - floor(tmp)) * 2 * PI; } else { tmp = ceil(abs(big / 2 / PI)); tmp = big + tmp * 2 * PI; } } return tmp; } function range_1(num) { var temp = num; if (num > 1) temp = num - 1; if (num < 0) temp = num + 1; return temp; } function range_24(big) { with (Math) { var tmp = 0; if (big > 0) { tmp = big / 24.0; tmp = (tmp - floor(tmp)) * 24.0; } else { tmp = ceil(abs(big / 24.0)); tmp = big + tmp * 24.0; } } return tmp; } function round_10(num) { return Math.floor((num + 0.05) * 10) / 10; } function round_100(num) { return Math.floor((num + 0.005) * 100) / 100; } function round_1000(num) { return Math.floor((num + 0.0005) * 1000) / 1000; } function hundred_miles(num) { return Math.floor((num / 100) + 0.5) * 100; } function isValidDate(dateStr) { if (IsValidTime(document.planets.ut_h_m.value) == true) { var datePat = /^(\d{1,2})(\/|-)(\d{1,2})\2(\d{4})$/; var matchArray = dateStr.match(datePat); if (matchArray == null) { alert("Date is not in a valid format.") return false; } var month = matchArray[1]; var day = matchArray[3]; var year = matchArray[4]; if (month < 1 || month > 12) { alert("Month must be between 1 and 12."); return false; } if (year < 1600 || year > 2400) { alert("Year must be between 1600 and 2400."); return false; } if (day < 1 || day > 31) { alert("Day must be between 1 and 31."); return false; } if ((month == 4 || month == 6 || month == 9 || month == 11) && day == 31) { alert("Month " + month + " doesn't have 31 days!") return false } if (month == 2) { var isleap = (year % 4 == 0 && (year % 100 != 0 || year % 400 == 0)); if (day > 29 || (day == 29 && !isleap)) { alert("February " + year + " doesn't have " + day + " days!"); return false; } } if (month < 10 && month.length == 1) month = "0" + month; if (day < 10 && day.length == 1) day = "0" + day; document.planets.date_txt.value = month + "/" + day + "/" + year; var dt_str = document.planets.date_txt.value; if ((dt_str.substring(2,3) != "/") || (dt_str.substring(5,6) != "/")) { alert ("Date is not in a valid format."); return false; } planets(); return true; } else { return false; } } function IsValidTime(timeStr) { var timePat = /^(\d{1,2}):(\d{2})(:(\d{2}))?(\s?(AM|am|PM|pm))?$/; var matchArray = timeStr.match(timePat); if (matchArray == null) { alert("Time is not in a valid format."); return false; } var hour = matchArray[1]; var minute = matchArray[2]; if (hour < 0 || hour > 23) { alert("Hour must be between 0 and 23."); return false; } if (minute < 0 || minute > 59) { alert ("Minute must be between 0 and 59."); return false; } if (hour < 10 && hour.length == 1) hour = "0" + hour; if (minute < 10 && minute.length == 1) minute = "0" + minute; document.planets.ut_h_m.value = hour + ":" + minute; var tm_str = document.planets.ut_h_m.value; if ((tm_str.substring(2,3) != ":") && (dt_str.length != 5)) { alert ("Time is not in a valid format."); return false; } return true; } function time_change(tmp) { if (isValidDate(document.planets.date_txt.value) == true) { var jd_temp, zz, ff, alpha, aa, bb, cc, dd, ee; var calendar_day, calendar_month, calendar_year; var int_day, hours, minutes; var tm_as_str, ut_hrs, ut_mns, part_day; var jd = julian_date(); tm_as_str = document.planets.ut_h_m.value; ut_hrs = eval(tm_as_str.substring(0,2)); ut_mns = eval(tm_as_str.substring(3,5)); part_day = ut_hrs / 24.0 + ut_mns / 1440.0; with (Math) { jd_temp = jd + part_day + tmp / 24.0 + 0.5; zz = floor(jd_temp); ff = jd_temp - zz; alpha = floor((zz - 1867216.25) / 36524.25); aa = zz + 1 + alpha - floor(alpha / 4); bb = aa + 1524; cc = floor((bb - 122.1) / 365.25); dd = floor(365.25 * cc); ee = floor((bb - dd) / 30.6001); calendar_day = bb - dd - floor(30.6001 * ee) + ff; calendar_month = ee; if (ee < 13.5) calendar_month = ee - 1; if (ee > 13.5) calendar_month = ee - 13; calendar_year = cc; if (calendar_month > 2.5) calendar_year = cc - 4716; if (calendar_month < 2.5) calendar_year = cc - 4715; int_day = floor(calendar_day); hours = (calendar_day - int_day) * 24; minutes = floor((hours - floor(hours)) * 60 + 0.5); hours = floor(hours); if (minutes > 59) {minutes = 0; hours = hours + 1;} if (calendar_month < 10) calendar_month = "0" + calendar_month; if (int_day < 10) int_day = "0" + int_day; if (hours < 10) hours = "0" + floor(hours); if (minutes < 10) minutes = "0" + minutes; } document.planets.date_txt.value = calendar_month + "/" + int_day + "/" + calendar_year; document.planets.ut_h_m.value = hours + ":" + minutes; planets(); return true; } else { return false; } } function zone_date_time() { var jd_temp, zz, ff, alpha, aa, bb, cc, dd, ee; var calendar_day, calendar_month, calendar_year; var int_day, hours, minutes, daylight_saving; var d_of_w, dow_str, x, y, lat_long = "", dst = 0; var tm_as_str, ut_hrs, ut_mns, part_day; var jd = julian_date(); tm_as_str = document.planets.ut_h_m.value; ut_hrs = eval(tm_as_str.substring(0,2)); ut_mns = eval(tm_as_str.substring(3,5)); part_day = ut_hrs / 24.0 + ut_mns / 1440.0; with (Math) { x = floor(abs(latitude)); y = floor((abs(latitude) - x) * 60 + 0.5); if (x < 10) lat_long += "0"; lat_long = lat_long + x + "\u00B0" + " "; if (y < 10) lat_long += "0"; lat_long = lat_long + y + "\'"; if (latitude >= 0) { lat_long += " N "; } else { lat_long += " S "; } x = floor(abs(longitude)); y = floor((abs(longitude) - x) * 60 + 0.5); if (x < 100) lat_long += "0"; if (x < 10) lat_long += "0"; lat_long = lat_long + x + "\u00B0" + " "; if (y < 10) lat_long += "0"; lat_long = lat_long + y + "\'"; if (longitude >= 0) { lat_long += " W"; } else { lat_long += " E"; } document.planets.lat_long.value = lat_long; daylight_saving = document.planets.d_s_t.checked; if (daylight_saving == true) { dst = 1; zone = t_zone + dst; } else { zone = t_zone; } document.planets.timezone.value = zone; jd_temp = jd + part_day + zone / 24 + 0.5; d_of_w = floor(jd_temp + 1.0) % 7; if (d_of_w == 0) dow_str = "Sun."; if (d_of_w == 1) dow_str = "Mon."; if (d_of_w == 2) dow_str = "Tue."; if (d_of_w == 3) dow_str = "Wed."; if (d_of_w == 4) dow_str = "Thu."; if (d_of_w == 5) dow_str = "Fri."; if (d_of_w == 6) dow_str = "Sat."; zz = floor(jd_temp); ff = jd_temp - zz; alpha = floor((zz - 1867216.25) / 36524.25); aa = zz + 1 + alpha - floor(alpha / 4); bb = aa + 1524; cc = floor((bb - 122.1) / 365.25); dd = floor(365.25 * cc); ee = floor((bb - dd) / 30.6001); calendar_day = bb - dd - floor(30.6001 * ee) + ff; calendar_month = ee; if (ee < 13.5) calendar_month = ee - 1; if (ee > 13.5) calendar_month = ee - 13; calendar_year = cc; if (calendar_month > 2.5) calendar_year = cc - 4716; if (calendar_month < 2.5) calendar_year = cc - 4715; int_day = floor(calendar_day); hours = (calendar_day - int_day) * 24; minutes = floor((hours - floor(hours)) * 60 + 0.5); hours = floor(hours); if (minutes > 59) {minutes = 0; hours = hours + 1;} if (calendar_month < 10) calendar_month = "0" + calendar_month; if (int_day < 10) int_day = "0" + int_day; } document.planets.date_zone.value = dow_str + " " + calendar_month + "/" + int_day + "/" + calendar_year + " " + am_pm(hours + minutes / 60.0); } function am_pm(h) { with (Math) { var hours, minutes, ampm, result=""; hours = floor(h); minutes = floor((h - hours) * 60 + 0.5); if (minutes == 60) { hours += 1; minutes = 0; } if (hours > 11) { hours -= 12; ampm = " pm"; } else {ampm = " am";} if (hours < 1) hours = 12; if (hours < 10) {result = result + "0" + floor(hours);} else {result = result + floor(hours);} result = result + ":"; if (minutes < 10) {result = result + "0" + minutes;} else {result = result + minutes;} } result = result + ampm return result; } function julian_date() { var dt_as_str, mm, dd, yy; var yyy, mmm, a, b; dt_as_str = document.planets.date_txt.value; mm = eval(dt_as_str.substring(0,2)); dd = eval(dt_as_str.substring(3,5)); yy = eval(dt_as_str.substring(6,10)); with (Math) { var yyy=yy; var mmm=mm; if (mm < 3) { yyy = yy - 1; mmm = mm + 12; } a = floor(yyy/100); b = 2 - a + floor(a/4); return floor(365.25*yyy) + floor(30.6001*(mmm+1)) + dd + 1720994.5 + b; } } function julian_date_zoned() { var dt_as_str, mm, dd, yy; var yyy, mmm, a, b; dt_as_str = document.planets.date_zone.value; mm = eval(dt_as_str.substring(5,7)); dd = eval(dt_as_str.substring(8,10)); yy = eval(dt_as_str.substring(11,15)); with (Math) { var yyy=yy; var mmm=mm; if (mm < 3) { yyy = yy - 1; mmm = mm + 12; } a = floor(yyy/100); b = 2 - a + floor(a/4); return floor(365.25*yyy) + floor(30.6001*(mmm+1)) + dd + 1720994.5 + b; } } function moon_facts() { with (Math) { var dt_as_str, mm, dd, yy, jd, i, tmp, k, t; var M, M1, F, NM_FM, FQ_LQ, RAD = 180 / PI; var jd_temp, d_of_w, dow_str, zz, ff, alpha; var aa, bb, cc, dd, ee, calendar_day, month; var calendar_month, int_day, hours, minutes; moon_data = ""; phase[0] = "New Moon "; phase[1] = "Moon at first quarter "; phase[2] = "Full Moon "; phase[3] = "Moon at last quarter "; dt_as_str = document.planets.date_txt.value; yy = eval(dt_as_str.substring(6,10)); mm = eval(dt_as_str.substring(0,2)); dd = eval(dt_as_str.substring(3,5)); tmp = floor(((yy + (mm - 1) / 12 + dd / 365.25) - 1900.0) * 12.3685); for (i=0; i<4; i++) { k = tmp + i * 0.25 t = k / 1236.85; M = proper_ang(359.2242 + 29.10535608 * k - 0.0000333 * pow(t,2) - 0.00000347 * pow(t,3)) / RAD; M1 = proper_ang(306.0253 + 385.81691806 * k + 0.0107306 * pow(t,2) + 0.00001236 * pow(t,3)) / RAD; F = proper_ang(21.2964 + 390.67050646 * k - 0.0016528 * pow(t,2) - 0.00000239 * pow(t,3)) / RAD; NM_FM = (0.1734 - 0.000393 * t) * sin(M) + 0.0021 * sin(2 * M) - 0.4068 * sin(M1) + 0.0161 * sin(2 * M1) - 0.0004 * sin(3 * M1) + 0.0104 * sin(2 * F) - 0.0051 * sin(M + M1) - 0.0074 * sin(M - M1) + 0.0004 * sin(2 * F + M) - 0.0004 * sin(2 * F - M) - 0.0006 * sin(2 * F + M1) + 0.0010 * sin(2 * F - M1) + 0.0005 * sin(M + 2 * M1); FQ_LQ = (0.1721 - 0.0004 * t) * sin(M) + 0.0021 * sin(2 * M) - 0.628 * sin(M1) + 0.0089 * sin(2 * M1) - 0.0004 * sin(3 * M1) + 0.0079 * sin(2 * F) - 0.0119 * sin(M + M1) - 0.0047 * sin(M - M1) + 0.0003 * sin(2 * F + M) - 0.0004 * sin(2 * F - M) - 0.0006 * sin(2 * F + M1) + 0.0021 * sin(2 * F - M1) + 0.0003 * sin(M + 2 * M1) + 0.0004 * sin(M - 2 * M1) - 0.0003 * sin(2 * M + M1); phases[i] = 2415020.75933 + 29.53058868 * k + 0.0001178 * pow(t,2) - 0.000000155 * pow(t,3) + 0.00033 * sin(2.907 + 2.319 * t - 1.6e-4 * pow(t,2)); if (i == 0 || i == 2) phases[i] += NM_FM; if (i == 1) phases[i] = phases[i] + FQ_LQ + 0.0028 - 0.0004 * cos(M) + 0.0003 * cos(M1); if (i == 3) phases[i] = phases[i] + FQ_LQ - 0.0028 + 0.0004 * cos(M) - 0.0003 * cos(M1); } new_moon = phases[0]; for (i=0; i<4; i++) { jd_temp = phases[i] + 0.5; d_of_w = floor(jd_temp + 1.0) % 7; if (d_of_w == 0) dow_str = "Sun."; if (d_of_w == 1) dow_str = "Mon."; if (d_of_w == 2) dow_str = "Tue."; if (d_of_w == 3) dow_str = "Wed."; if (d_of_w == 4) dow_str = "Thu."; if (d_of_w == 5) dow_str = "Fri."; if (d_of_w == 6) dow_str = "Sat."; zz = floor(jd_temp); ff = jd_temp - zz; alpha = floor((zz - 1867216.25) / 36524.25); aa = zz + 1 + alpha - floor(alpha / 4); bb = aa + 1524; cc = floor((bb - 122.1) / 365.25); dd = floor(365.25 * cc); ee = floor((bb - dd) / 30.6001); calendar_day = bb - dd - floor(30.6001 * ee) + ff; calendar_month = ee; if (ee < 13.5) calendar_month = ee - 1; if (ee > 13.5) calendar_month = ee - 13; calendar_year = cc; if (calendar_month > 2.5) calendar_year = cc - 4716; if (calendar_month < 2.5) calendar_year = cc - 4715; int_day = floor(calendar_day); hours = (calendar_day - int_day) * 24; minutes = floor((hours - floor(hours)) * 60 + 0.5); hours = floor(hours); if (minutes > 59) {minutes = 0; hours = hours + 1;} if (hours < 10) hours = "0" + hours; if (minutes < 10) minutes = "0" + minutes; if (int_day < 10) int_day = "0" + int_day; if (calendar_month == 1) month = "Jan."; if (calendar_month == 2) month = "Feb."; if (calendar_month == 3) month = "Mar."; if (calendar_month == 4) month = "Apr."; if (calendar_month == 5) month = "May "; if (calendar_month == 6) month = "Jun."; if (calendar_month == 7) month = "Jul."; if (calendar_month == 8) month = "Aug."; if (calendar_month == 9) month = "Sep."; if (calendar_month == 10) month = "Oct."; if (calendar_month == 11) month = "Nov."; if (calendar_month == 12) month = "Dec."; moon_data += dow_str + " " + month + " " + int_day + ", " + calendar_year + " " + phase[i] + "(" + hours + ":" + minutes + " UT)"; if (i == 1 || i == 3) { moon_data += "\n" } else { moon_data += " "; } } } } function planets() { with (Math) { var RAD = 180 / PI; document.planets.size6.value = ""; document.planets.size7.value = ""; document.planets.phase7.value = ""; document.planets.phase6.value = ""; document.planets.dawn.value = ""; document.planets.dusk.value = ""; document.planets.tr0.value = "..."; document.planets.tr1.value = "..."; document.planets.tr2.value = "..."; document.planets.tr3.value = "..."; document.planets.tr4.value = "..."; document.planets.tr5.value = "..."; document.planets.tr6.value = "..."; var planet_dec = new Array(9); var planet_ra = new Array(9); var planet_dist = new Array(9); var planet_phase = new Array(9); var planet_size = new Array(9); var planet_alt = new Array(9); var planet_az = new Array(9); var ang_at_1au = new Array(1,6.68,16.82,9.36,196.94,166.6); var planet_a = new Array(5); var planet_l = new Array(5); var planet_i = new Array(5); var planet_e = new Array(5); var planet_n = new Array(5); var planet_m = new Array(5); var planet_r = new Array(5); var jd, t, l, m, e, c, sl, r, R, Rm, u, p, q, s, v, z, a, b, ee, oe, t0, gt; var g0, x, y, i, indx, count, hl, ha, eg, et, dr, en, de, ra, la, lo, height; var az, al, size, il, dist, rho_sin_lat, rho_cos_lat, h_p, rise_str, set_str; var dawn, dusk, gt_dawn, gt_dusk, dawn_planets, dusk_planets, morn, eve; var sunrise_set, moonrise_set, twilight_begin_end, twilight_start, twilight_end; var sat_lat, phase_angle, dt_as_str, tm_as_str, ut_hrs, ut_mns, part_day, dt; var moon_days, Days, age, N, M_sun, Ec, lambdasun, P0, N0, rise, set, tr, h0; var i_m, e_m, l_m, M_m, N_m, Ev, Ae, A3, A4, lambda_moon, beta_moon, phase; var transit, lambda_transit, beta_transit, de_transit, ra_transit; var de_np, ra_np, d_beta, d_lambda, lambda_moon_D, beta_moon_D, yy, mm, dd; la = latitude / RAD; lo = longitude / RAD; height = 10.0; tm_as_str = document.planets.ut_h_m.value; ut_hrs = eval(tm_as_str.substring(0,2)); ut_mns = eval(tm_as_str.substring(3,5)); dt_as_str = document.planets.date_txt.value; yy = eval(dt_as_str.substring(6,10)); mm = eval(dt_as_str.substring(0,2)); dd = eval(dt_as_str.substring(3,5)); part_day = ut_hrs / 24.0 + ut_mns / 1440.0; jd = julian_date() + part_day; document.planets.phase8.value = round_1000(jd); y = floor(yy / 50 + 0.5) * 50; if (y == 1600) dt = 2; if (y == 1650) dt = 0.8; if (y == 1700) dt = 0.1; if (y == 1750) dt = 0.2; if (y == 1800) dt = 0.2; if (y == 1850) dt = 0.1; if (y == 1900) dt = 0; if (y == 1950) dt = 0.5; if (y == 2000) dt = 1.1; if (y == 2050) dt = 3; if (y == 2100) dt = 5; if (y == 2150) dt = 6; if (y == 2200) dt = 8; if (y == 2250) dt = 11; if (y == 2300) dt = 13; if (y == 2350) dt = 16; if (y == 2400) dt = 19; dt = dt / 1440; jd += dt; zone_date_time(); moon_facts(); t=(jd-2415020.0)/36525; oe=.409319747-2.271109689e-4*t-2.8623e-8*pow(t,2)+8.779e-9*pow(t,3); t0=(julian_date()-2415020.0)/36525; g0=.276919398+100.0021359*t0+1.075e-6*pow(t0,2); g0=(g0-floor(g0))*2*PI; gt=proper_ang_rad(g0+part_day*6.30038808); l=proper_ang_rad(4.88162797+628.331951*t+5.27962099e-6*pow(t,2)); m=proper_ang_rad(6.2565835+628.3019457*t-2.617993878e-6*pow(t,2)-5.759586531e-8*pow(t,3)); e=.01675104-4.18e-5*t-1.26e-7*pow(t,2); c=proper_ang_rad(3.3500896e-2-8.358382e-5*t-2.44e-7*pow(t,2))*sin(m)+(3.50706e-4-1.7e-6*t)*sin(2*m)+5.1138e-6*sin(3*m); sl=proper_ang_rad(l+c); R=1.0000002*(1-pow(e,2))/(1+e*cos(m+c)); de = asin(sin(oe) * sin(sl)); y = sin(sl) * cos(oe); x = cos(sl); ra = proper_ang_rad(atan2(y,x)); dist = R; size = 31.9877 / dist; ha=proper_ang_rad(gt-lo-ra); al=asin(sin(la)*sin(de)+cos(la)*cos(de)*cos(ha)); az=acos((sin(de)-sin(la)*sin(al))/(cos(la)*cos(al))); if (sin(ha)>=0) az=2*PI-az; planet_decl[0] = de; planet_r_a[0] = ra; planet_mag[0] = -26.7; planet_dec[8] = ""; planet_ra[8] = ""; planet_dist[8] = ""; planet_size[8] = ""; planet_alt[8] = ""; planet_az[8] = ""; y = sin(-0.8333/RAD) - sin(la) * sin(de); x = cos(la) * cos(de); if (y/x >= 1) { rise_str = "No rise"; set_str = "No rise"; h0 = 0; } if (y/x <= -1) { rise_str = "No set"; set_str = "No set"; h0 = 0; } if (abs(y/x) < 1) { h0 = acos(y/x) * RAD; tr = range_1((ra * RAD + lo * RAD - g0 * RAD) / 360); document.planets.tr0.value = am_pm(range_1(tr + zone / 24 - dt) * 24); rise_str = am_pm(range_1(tr - h0 / 360 + zone / 24 - dt) * 24); set_str = am_pm(range_1(tr + h0 / 360 + zone / 24 - dt) * 24); } y = sin(-18/RAD) - sin(la) * sin(de); x = cos(la) * cos(de); if (y/x >= 1) { twilight_start = "No twilight"; twilight_end = "No twilight"; h0 = 0; } if (y/x <= -1) { twilight_start = "All night"; twilight_end = "All night"; h0 = 0; } if (abs(y/x) < 1) { h0 = acos(y/x) * RAD; tr = range_1((ra * RAD + lo * RAD - g0 * RAD) / 360); twilight_start = am_pm(range_1(tr - h0 / 360 + zone / 24 - dt) * 24); twilight_end = am_pm(range_1(tr + h0 / 360 + zone / 24 - dt) * 24); } al = al * RAD; az = az * RAD; if (al >= 0) { al += 1.2 / (al + 2); } else { al += 1.2 / (abs(al) + 2); } planet_altitude[0] = al; planet_azimuth[0] = az; de = de * RAD; ra = ra * RAD / 15; x = round_10(abs(de)); if (de < 0) { planet_dec[8] += "-"; } else { planet_dec[8] += "+"; } if (x < 10) planet_dec[8] += "0"; planet_dec[8] += x; if (x == floor(x)) planet_dec[8] += ".0"; planet_dec[8] += "\u00B0"; x = floor(ra); y = floor((ra - x) * 60 + 0.5); if (y == 60) { x += 1; y = 0; } if (x < 10) planet_ra[8] += "0"; planet_ra[8] += x + "h "; if (y < 10) planet_ra[8] += "0"; planet_ra[8] += y + "m"; dist = round_100(dist); if (dist < 10) planet_dist[8] += "0"; planet_dist[8] += dist; size = round_10(size); planet_size[8] += size; if (size == floor(size)) planet_size[8] += ".0"; planet_size[8] += "\'"; sunrise_set = sun_rise_set(la,lo,zone); rise_str = sunrise_set.substring(0,8); set_str = sunrise_set.substring(11,19); planet_alt[8] = rise_str; planet_az[8] = set_str; /* planet_alt[8] = rise_str; planet_az[8] = set_str; sunrise_set = rise_str + " / " + set_str; */ moon_days = 29.53058868; Days = jd - 2444238.5; N = proper_ang(Days / 1.01456167); M_sun = proper_ang(N - 3.762863); Ec = 1.91574168 * sin(M_sun / RAD); lambdasun = proper_ang(N + Ec + 278.83354); l0 = 64.975464; P0 = 349.383063; N0 = 151.950429; i_m = 5.145396; e_m = 0.0549; l_m = proper_ang(13.1763966 * Days + l0); M_m = proper_ang(l_m - 0.111404 * Days - P0); N_m = proper_ang(N0 - 0.0529539 * Days); Ev = 1.2739 * sin((2 * (l_m - lambdasun) - M_m) / RAD); Ae = 0.1858 * sin(M_sun / RAD); A3 = 0.37 * sin(M_sun / RAD); M_m += Ev - Ae - A3; Ec = 6.2886 * sin(M_m / RAD); dist = hundred_miles(238855.7 * (1 - pow(e_m,2)) / (1 + e_m * cos((M_m + Ec) / RAD))); A4 = 0.214 * sin(2 * M_m / RAD); l_m += Ev + Ec - Ae + A4; l_m += 0.6583 * sin(2 * (l_m - lambdasun) / RAD); N_m -= 0.16 * sin(M_sun / RAD); y = sin((l_m - N_m) / RAD) * cos(i_m / RAD); x = cos((l_m - N_m) / RAD); lambda_moon = proper_ang_rad(atan2(y,x) + N_m / RAD); beta_moon = asin(sin((l_m - N_m) / RAD) * sin(i_m / RAD)); d_beta = 8.727e-4 * cos((l_m - N_m) / RAD); d_lambda = 9.599e-3 + 1.047e-3 * cos(M_m / RAD); de = asin(sin(beta_moon) * cos(oe) + cos(beta_moon) * sin(oe) * sin(lambda_moon)); y = sin(lambda_moon) * cos(oe) - tan(beta_moon) * sin(oe); x = cos(lambda_moon); ra = proper_ang_rad(atan2(y,x)); de_np = de; ra_np = ra; size = 2160 / dist * RAD * 60; h_p = 3963.2 / dist * RAD; x = proper_ang(l_m - lambdasun); phase = 50.0 * (1.0 - cos(x / RAD)); age = jd - new_moon; if (age > moon_days) age -= moon_days; if (age < 0) age += moon_days; /* age = (x / 360 * moon_days); */ Rm = R - dist * cos(x / RAD) / 92955800; phase_angle = acos((pow(Rm,2) + pow((dist/92955800),2) - pow(R,2)) / (2 * Rm * dist/92955800)) * RAD / 100; planet_mag[6] = round_10(0.21 + 5 * log(Rm * dist / 92955800)/log(10) + 3.05 * (phase_angle) - 1.02 * pow(phase_angle,2) + 1.05 * pow(phase_angle,3)); if (planet_mag[6] == floor(planet_mag[6])) planet_mag[6] += ".0"; /* thisImage = floor(age / moon_days * 28); */ thisImage = floor(age / moon_days * 27 + 0.5); document.myPicture.src = myImage[thisImage]; age = round_10(age); x = atan(0.996647 * tan(la)); rho_sin_lat = 0.996647 * sin(x) + height * sin(la) / 6378140; rho_cos_lat = cos(x) + height * cos(la) / 6378140; r = dist / 3963.2; ha = proper_ang_rad(gt - lo - ra); x = atan(rho_cos_lat * sin(ha) / (r * cos(de) - rho_cos_lat * cos(ha))); ra -= x; y = ha; ha += x; de = atan(cos(ha) * (r * sin(de) - rho_sin_lat) / (r * cos(de) * cos(y) - rho_cos_lat)); ha = proper_ang_rad(gt - lo - ra); al = asin(sin(de) * sin(la) + cos(de) * cos(la) * cos(ha)); az = acos((sin(de) - sin(la) * sin(al)) / (cos(la) * cos(al))); if (sin(ha) >= 0) az = 2 * PI - az; planet_decl[6] = de; planet_r_a[6] = ra; planet_dec[9] = ""; planet_ra[9] = ""; planet_dist[9] = ""; planet_dist[6] = ""; planet_phase[9] = ""; planet_size[9] = ""; planet_alt[9] = ""; planet_az[9] = ""; lambda_moon_D = lambda_moon - d_lambda * part_day * 24; beta_moon_D = beta_moon - d_beta * part_day * 24; tr = range_1((ra_np * RAD + lo * RAD - g0 * RAD) / 360); transit = tr * 24; lambda_transit = lambda_moon_D + transit * d_lambda; beta_transit = beta_moon_D + transit * d_beta; de_transit = asin(sin(beta_transit) * cos(oe) + cos(beta_transit) * sin(oe) * sin(lambda_transit)); y = sin(lambda_transit) * cos(oe) - tan(beta_transit) * sin(oe); x = cos(lambda_transit); ra_transit = proper_ang_rad(atan2(y,x)); tr = range_1((ra_transit * RAD + lo * RAD - g0 * RAD) / 360); document.planets.tr6.value = am_pm(range_1(tr + zone / 24 - dt) * 24); al = al * RAD; az = az * RAD; if (al >= 0) { al += 1.2 / (al + 2); } else { al += 1.2 / (abs(al) + 2); } planet_altitude[6] = al; planet_azimuth[6] = az; de = de * RAD; ra = ra * RAD / 15; x = round_10(abs(de)); if (de < 0) { planet_dec[9] += "-"; } else { planet_dec[9] += "+"; } if (x < 10) planet_dec[9] += "0"; planet_dec[9] += x; if (x == floor(x)) planet_dec[9] += ".0"; planet_dec[9] += "\u00B0"; x = floor(ra); y = floor((ra - x) * 60 + 0.5); if (y == 60) { x += 1; y = 0; } if (x < 10) planet_ra[9] += "0"; planet_ra[9] += x + "h "; if (y < 10) planet_ra[9] += "0"; planet_ra[9] += y + "m"; if (age < 10) planet_dist[6] += "0"; planet_dist[6] += age; if (age == floor(age)) planet_dist[6] += ".0"; planet_dist[6] += " days"; planet_dist[9] += dist; size = round_10(size); planet_size[9] += size; if (size == floor(size)) planet_size[9] += ".0"; planet_size[9] += "\'"; phase = floor(phase + 0.5); planet_phase[9] += phase + "%"; moonrise_set = moon_rise_set(la,lo,zone); rise_str = moonrise_set.substring(0,8); set_str = moonrise_set.substring(11,19); planet_alt[9] = rise_str; planet_az[9] = set_str; if (twilight_start == "All night") { twilight_begin_end = "twilight all night"; } else if (twilight_start == "No twilight") { twilight_begin_end = "Sky dark for 24h"; } else { twilight_begin_end = twilight_start + " / " + twilight_end; } u=t/5+.1; p=proper_ang_rad(4.14473024+52.96910393*t); q=proper_ang_rad(4.64111846+21.32991139*t); s=proper_ang_rad(4.2501771+7.478172103*t); v=5*q-2*p; z=proper_ang_rad(q-p); if (v >= 0) { v=proper_ang_rad(v); } else { v=proper_ang_rad(abs(v))+2*PI; } planet_a[1]=0.3870986; planet_l[1]=proper_ang_rad(3.109811569+2608.814681*t+5.2552e-6*pow(t,2)); planet_e[1]=.20561421+2.046e-5*t-3e-8*pow(t,2); planet_i[1]=.12222333+3.247708672e-5*t-3.194e-7*pow(t,2); planet_n[1]=.822851951+.020685787*t+3.035127569e-6*pow(t,2); planet_m[1]=proper_ang_rad(1.785111938+2608.787533*t+1.22e-7*pow(t,2)); planet_a[2]=0.7233316; planet_l[2]=proper_ang_rad(5.982413642+1021.352924*t+5.4053e-6*pow(t,2)); planet_e[2]=.00682069-4.774e-5*t+9.1e-8*pow(t,2); planet_i[2]=.059230034+1.75545e-5*t-1.7e-8*pow(t,2); planet_n[2]=1.322604346+.0157053*t+7.15585e-6*pow(t,2); planet_m[2]=proper_ang_rad(3.710626189+1021.328349*t+2.2445e-5*pow(t,2)); planet_a[3]=1.5236883; planet_l[3]=proper_ang_rad(5.126683614+334.0856111*t+5.422738e-6*pow(t,2)); planet_e[3]=.0933129+9.2064e-5*t-7.7e-8*pow(t,2); planet_i[3]=.032294403-1.178097e-5*t+2.199e-7*pow(t,2); planet_n[3]=.851484043+.013456343*t-2.44e-8*pow(t,2)-9.3026e-8*pow(t,3); planet_m[3]=proper_ang_rad(5.576660841+334.0534837*t+3.159e-6*pow(t,2)); planet_l[4]=proper_ang_rad(4.154743317+52.99346674*t+5.841617e-6*pow(t,2)-2.8798e-8*pow(t,3)); a=.0058*sin(v)-1.1e-3*u*cos(v)+3.2e-4*sin(2*z)-5.8e-4*cos(z)*sin(q)-6.2e-4*sin(z)*cos(q)+2.4e-4*sin(z); b=-3.5e-4*cos(v)+5.9e-4*cos(z)*sin(q)+6.6e-4*sin(z)*cos(q); ee=3.6e-4*sin(v)-6.8e-4*sin(z)*sin(q); planet_e[4]=.04833475+1.6418e-4*t-4.676e-7*pow(t,2)-1.7e-9*pow(t,3); planet_a[4]=5.202561+7e-4*cos(2*z); planet_i[4]=.022841752-9.941569952e-5*t+6.806784083e-8*pow(t,2); planet_n[4]=1.735614994+.017637075*t+6.147398691e-6*pow(t,2)-1.485275193e-7*pow(t,3); planet_m[4]=proper_ang_rad(3.932721257+52.96536753*t-1.26012772e-5*pow(t,2)); planet_m[4]=proper_ang_rad(planet_m[4]+a-b/planet_e[4]); planet_l[4]=proper_ang_rad(planet_l[4]+a); planet_e[4]=planet_e[4]+ee; planet_a[5]=9.554747; planet_l[5]=proper_ang_rad(4.652426047+21.35427591*t+5.663593422e-6*pow(t,2)-1.01229e-7*pow(t,3)); a=-.014*sin(v)+2.8e-3*u*cos(v)-2.6e-3*sin(z)-7.1e-4*sin(2*z)+1.4e-3*cos(z)*sin(q); a=a+1.5e-3*sin(z)*cos(q)+4.4e-4*cos(z)*cos(q); a=a-2.7e-4*sin(3*z)-2.9e-4*sin(2*z)*sin(q)+2.6e-4*cos(2*z)*sin(q); a=a+2.5e-4*cos(2*z)*cos(q); b=1.4e-3*sin(v)+7e-4*cos(v)-1.3e-3*sin(z)*sin(q); b=b-4.3e-4*sin(2*z)*sin(q)-1.3e-3*cos(q)-2.6e-3*cos(z)*cos(q)+4.7e-4*cos(2*z)*cos(q); b=b+1.7e-4*cos(3*z)*cos(q)-2.4e-4*sin(z)*sin(2*q); b=b+2.3e-4*cos(2*z)*sin(2*q)-2.3e-4*sin(z)*cos(2*q)+2.1e-4*sin(2*z)*cos(2*q); b=b+2.6e-4*cos(z)*cos(2*q)-2.4e-4*cos(2*z)*cos(2*q); planet_l[5]=proper_ang_rad(planet_l[5]+a); planet_e[5]=.05589232-3.455e-4*t-7.28e-7*pow(t,2)+7.4e-10*pow(t,3); planet_i[5]=.043502663-6.839770806e-5*t-2.703515011e-7*pow(t,2)+6.98e-10*pow(t,3); planet_n[5]=1.968564089+.015240129*t-2.656042056e-6*pow(t,2)-9.2677e-8*pow(t,3); planet_m[5]=proper_ang_rad(3.062463265+21.32009513*t-8.761552845e-6*pow(t,2)); planet_m[5]=proper_ang_rad(planet_m[5]+a-b/planet_e[5]); planet_a[5]=planet_a[5]+.0336*cos(z)-.00308*cos(2*z)+.00293*cos(v); planet_e[5]=planet_e[5]+1.4e-3*cos(v)+1.2e-3*sin(q)+2.7e-3*cos(z)*sin(q)-1.3e-3*sin(z)*cos(q); for (i=1; i<6; i++) { e = planet_m[i]; for (count = 1; count <= 50; count++) { de = e - planet_e[i] * sin(e) - planet_m[i]; if (abs(de) <= 5e-8) break; e = e - de / (1 - planet_e[i] * cos(e)); } v=2*atan(sqrt((1+planet_e[i])/(1-planet_e[i]))*tan(e/2)); planet_r[i]=planet_a[i]*(1-planet_e[i]*cos(e)); u=proper_ang_rad(planet_l[i]+v-planet_m[i]-planet_n[i]); x=cos(planet_i[i])*sin(u); y=cos(u); y=acos(y/sqrt(pow(x,2)+pow(y,2))); if (x<0) y=2*PI-y; hl=proper_ang_rad(y+planet_n[i]); ha=asin(sin(u)*sin(planet_i[i])); x=planet_r[i]*cos(ha)*sin(proper_ang_rad(hl-sl)); y=planet_r[i]*cos(ha)*cos(proper_ang_rad(hl-sl))+R; y=acos(y/sqrt(pow(x,2)+pow(y,2))); if (x<0) y=2*PI-y; eg=proper_ang_rad(y+sl); dr=sqrt(pow(R,2)+pow(planet_r[i],2)+2*planet_r[i]*R*cos(ha)*cos(proper_ang_rad(hl-sl))); size=ang_at_1au[i]/dr; il=floor((pow((planet_r[i]+dr),2)-pow(R,2))/(4*planet_r[i]*dr)*100+.5); phase_angle=acos((pow(planet_r[i],2)+pow(dr,2)-pow(R,2))/(2*planet_r[i]*dr))*RAD; et=asin(planet_r[i]/dr*sin(ha)); en=acos(cos(et)*cos(proper_ang_rad(eg-sl))); de=asin(sin(et)*cos(oe)+cos(et)*sin(oe)*sin(eg)); y=cos(eg); x=sin(eg)*cos(oe)-tan(et)*sin(oe); y=acos(y/sqrt(pow(x,2)+pow(y,2))); if (x<0) y=2*PI-y; ra=y; ha=proper_ang_rad(gt-lo-ra); al=asin(sin(la)*sin(de)+cos(la)*cos(de)*cos(ha)); az=acos((sin(de)-sin(la)*sin(al))/(cos(la)*cos(al))); if (sin(ha)>=0) az=2*PI-az; planet_decl[i] = de; planet_r_a[i] = ra; planet_mag[i]=5*log(planet_r[i]*dr)/log(10); if (i == 1) planet_mag[i] = -0.42 + planet_mag[i] + 0.038*phase_angle - 0.000273*pow(phase_angle,2) + 0.000002*pow(phase_angle,3); if (i == 2) planet_mag[i] = -4.4 + planet_mag[i] + 0.0009*phase_angle + 0.000239*pow(phase_angle,2) - 0.00000065*pow(phase_angle,3); if (i == 3) planet_mag[i] = -1.52 + planet_mag[i] + 0.016*phase_angle; if (i == 4) planet_mag[i] = -9.4 + planet_mag[i] + 0.005*phase_angle; if (i == 5) { sat_lat = asin(0.47063*cos(et)*sin(eg-2.9585)-0.88233*sin(et)); planet_mag[i] = -8.88 + planet_mag[i] + 0.044*phase_angle - 2.6*sin(abs(sat_lat)) + 1.25*pow(sin(sat_lat),2); } planet_mag[i] = round_10(planet_mag[i]); if (planet_mag[i] == floor(planet_mag[i])) planet_mag[i] += ".0"; if (planet_mag[i] >= 0) planet_mag[i] = "+" + planet_mag[i]; planet_dec[i] = ""; planet_ra[i] = ""; planet_dist[i] = ""; planet_size[i] = ""; planet_phase[i] = ""; planet_alt[i] = ""; planet_az[i] = ""; y = sin(-0.5667/RAD) - sin(la) * sin(de); x = cos(la) * cos(de); if (y/x >= 1) { rise_str = "No rise"; set_str = "No rise"; h0 = 0; } if (y/x <= -1) { rise_str = "No set"; set_str = "No set"; h0 = 0; } if (abs(y/x) < 1) { h0 = acos(y/x) * RAD; tr = range_1((ra * RAD + lo * RAD - g0 * RAD) / 360); rise_str = am_pm(range_1(tr - h0 / 360 + zone / 24 - dt) * 24); set_str = am_pm(range_1(tr + h0 / 360 + zone / 24 - dt) * 24); tr = am_pm(range_1(tr + zone / 24 - dt) * 24); if (i == 1) document.planets.tr1.value = tr; if (i == 2) document.planets.tr2.value = tr; if (i == 3) document.planets.tr3.value = tr; if (i == 4) document.planets.tr4.value = tr; if (i == 5) document.planets.tr5.value = tr; } al = al * RAD; az = az * RAD; planet_altitude[i] = al; planet_azimuth[i] = az; de = de * RAD; ra = ra * RAD / 15; x = round_10(abs(de)); if (de < 0) { planet_dec[i] += "-"; } else { planet_dec[i] += "+"; } if (x < 10) planet_dec[i] += "0"; planet_dec[i] += x; if (x == floor(x)) planet_dec[i] += ".0"; planet_dec[i] += "\u00B0"; x = floor(ra); y = floor((ra - x) * 60 + 0.5); if (y == 60) { x += 1; y = 0; } if (x < 10) planet_ra[i] += "0"; planet_ra[i] += x + "h "; if (y < 10) planet_ra[i] += "0"; planet_ra[i] += y + "m"; dr = round_100(dr); if (dr < 10) planet_dist[i] += "0"; planet_dist[i] += dr; size = round_10(size); if (size < 10) planet_size[i] += "0"; planet_size[i] += size; if (size == floor(size)) planet_size[i] += ".0"; planet_size[i] += "\""; planet_phase[i] += il + "%"; planet_alt[i] = rise_str; planet_az[i] = set_str; } dawn_planets = ""; dusk_planets = ""; phenomena = ""; y = sin(-9/RAD) - sin(la) * sin(planet_decl[0]); x = cos(la) * cos(planet_decl[0]); if (y/x >= 1) { rise_str = "No rise"; set_str = "No rise"; dawn_planets = "-----"; dusk_planets = "-----"; } if (y/x <= -1) { rise_str = "No set"; set_str = "No set"; dawn_planets = "-----"; dusk_planets = "-----"; } if (abs(y/x) < 1) { h0 = acos(y/x) * RAD; tr = range_1((planet_r_a[0] * RAD + lo * RAD - g0 * RAD) / 360); dawn = range_1(tr - h0 / 360 - dt); dusk = range_1(tr + h0 / 360 - dt); } gt_dawn = proper_ang_rad(g0 + dawn * 6.30038808); gt_dusk = proper_ang_rad(g0 + dusk * 6.30038808); indx = 0; morn = 0; eve = 0; for (i=0; i<17; i++) { if (i < 6) { ha = proper_ang_rad(gt_dawn - lo - planet_r_a[i]); al = asin(sin(la) * sin(planet_decl[i]) + cos(la) * cos(planet_decl[i]) * cos(ha)); if (al >= 0) { if (morn > 0) dawn_planets += ", "; dawn_planets += planet_name[i]; morn ++; } ha = proper_ang_rad(gt_dusk - lo - planet_r_a[i]); al = asin(sin(la) * sin(planet_decl[i]) + cos(la) * cos(planet_decl[i]) * cos(ha)); if (al >= 0) { if (eve > 0) dusk_planets += ", "; dusk_planets += planet_name[i]; eve ++; } } x = round_10(acos(sin(planet_decl[6]) * sin(planet_decl[i]) + cos(planet_decl[6]) * cos(planet_decl[i]) * cos(planet_r_a[6] - planet_r_a[i])) * RAD); if (x <= 10 && i != 6) { if (indx > 0) phenomena += " and "; if (x < 1) { phenomena += "less than 1"; } else { phenomena += floor(x + 0.5); } phenomena += "\u00B0" + " from " + planet_name[i]; if (x < 0.5) { if (i == 0) { phenomena += " (eclipse possible)"; } else { phenomena += " (occultation possible)"; } } indx ++ } } if (dawn_planets == "") dawn_planets = "-----"; if (dusk_planets == "") dusk_planets = "-----"; if (indx > 0) { phenomena = "The Moon is " + phenomena + ". "; } indx = 0; for (i=1; i<16; i++) { for (count=i+1; count<17; count++) { if (i != 6 && count != 6) { x = round_10(acos(sin(planet_decl[count]) * sin(planet_decl[i]) + cos(planet_decl[count]) * cos(planet_decl[i]) * cos(planet_r_a[count] - planet_r_a[i])) * RAD); if (x <= 5) { if (indx > 0) phenomena += " and "; phenomena += planet_name[count] + " is "; if (x < 1) { phenomena += "less than 1"; } else { phenomena += floor(x + 0.5); } phenomena += "\u00B0" + " from " + planet_name[i]; indx ++ } x = 0; } } } if (indx > 0) phenomena += "." if (jd > 2454048.3007 && jd < 2454048.5078) { phenomena += "\nMercury is transiting the face of the Sun. "; x = round_10((2454048.50704 - jd) * 24); if (x >= 0.1) { phenomena += "The spectacle continues for a further " + x + " hours." } } if (phenomena != "") { phenomena += "\n"; } if (planet_altitude[0] < -18) phenomena += "The sky is dark. "; if (planet_altitude[0] >= -18 && planet_altitude[0] < -12) phenomena += "The Sun is below the horizon (altitude " + floor(planet_altitude[0] + 0.5) + "\u00B0) and the sky is in deep twilight. "; if (planet_altitude[0] >= -12 && planet_altitude[0] < -6) phenomena += "The Sun is below the horizon (altitude " + floor(planet_altitude[0] + 0.5) + "\u00B0) and the sky is in twilight. "; if (planet_altitude[0] >= -6 && planet_altitude[0] < -0.3) phenomena += "The Sun is below the horizon (altitude " + floor(planet_altitude[0] + 0.5) + "\u00B0) and the sky is in bright twilight. "; if (planet_altitude[0] >= -0.3) phenomena += "The Sun is above the horizon (altitude " + floor(planet_altitude[0] + 0.5) + "\u00B0). "; if (planet_altitude[6] >= -0.3) phenomena += "The Moon is above the horizon (altitude " + floor(planet_altitude[6] + 0.5) + "\u00B0). "; if (planet_altitude[6] < -0.3) phenomena += "The Moon is below the horizon. "; phenomena += "\n"; phenomena += moon_data; /* Commented out 2009-specific events 2010/01/27 - DAF phenomena += "-----\n2009 Phenomena\n"; phenomena += "yyyy mm dd hh UT All dates and times are UT\n" + "2009 01 02 04 UT Mercury at greatest illuminated extent (evening)\n" + "2009 01 03 13 UT Meteor shower peak -- Quadrantids\n" + "2009 01 04 14 UT Mercury at greatest elongation, 19.3\u00B0 east of Sun (evening)\n" + "2009 01 04 16 UT Earth at perihelion 91,400,939 miles (147,095,552 km)\n" + "2009 01 10 11 UT Moon at perigee, 222,137 miles (357,495 km)\n" + "2009 01 14 21 UT Venus at greatest elongation, 47.1\u00B0 east of Sun (evening)\n" + "2009 01 23 00 UT Moon at apogee, 252,350 miles (406,118 km)\n" + "2009 01 26 Annular eclipse of the Sun (South Atlantic, Indonesia)\n" + "2009 02 07 20 UT Moon at perigee, 224,619 miles (361,489 km)\n" + "2009 02 09 14:38 UT Penumbral lunar eclipse (midtime)\n" + "2009 02 11 01 UT Mercury at greatest illuminated extent (morning)\n" + "2009 02 13 21 UT Mercury at greatest elongation, 26.1\u00B0 west of Sun (morning)\n" + "2009 02 19 16 UT Venus at greatest illuminated extent (evening)\n" + "2009 02 19 17 UT Moon at apogee, 251,736 miles (405,129 km)\n" + "2009 03 07 15 UT Moon at perigee, 228,053 miles (367,016 km)\n" + "2009 03 08 20 UT Saturn at opposition\n" + "2009 03 19 13 UT Moon at apogee, 251,220 miles (404,299 km)\n" + "2009 03 20 11:44 UT Spring (Northern Hemisphere) begins at the equinox\n" + "2009 03 27 19 UT Venus at inferior conjunction with Sun\n" + "2009 04 02 02 UT Moon at perigee, 229,916 miles (370,013 km)\n" + "2009 04 13 05 UT Mercury at greatest illuminated extent (evening)\n" + "2009 04 16 09 UT Moon at apogee, 251,178 miles (404,232 km)\n" + "2009 04 22 11 UT Meteor shower peak -- Lyrids\n" + "2009 04 26 08 UT Mercury at greatest elongation, 20.4\u00B0 east of Sun (evening)\n" + "2009 04 28 06 UT Moon at perigee, 227,446 miles (366,039 km)\n" + "2009 05 02 14 UT Venus at greatest illuminated extent (morning)\n" + "2009 05 06 00 UT Meteor shower peak -- Eta Aquarids\n" + "2009 05 14 03 UT Moon at apogee, 251,603 miles (404,915 km)\n" + "2009 05 26 04 UT Moon at perigee, 224,410 miles (361,153 km)\n" + "2009 06 05 21 UT Venus at greatest elongation, 45.8\u00B0 west of Sun (morning)\n" + "2009 06 10 16 UT Moon at apogee, 252,144 miles (405,787 km)\n" + "2009 06 13 12 UT Mercury at greatest elongation, 23.5\u00B0 west of Sun (morning)\n" + "2009 06 21 05:46 UT Summer (Northern Hemisphere) begins at the solstice\n" + "2009 06 23 11 UT Moon at perigee, 222,459 miles (358,013 km)\n" + "2009 07 03 21 UT Mercury at greatest illuminated extent (morning)\n" + "2009 07 04 02 UT Earth at aphelion 94,505,048 miles (152,091,132 km)\n" + "2009 07 07 09:39 UT Penumbral lunar eclipse (too slight to see)\n" + "2009 07 07 22 UT Moon at apogee, 252,421 miles (406,232 km)\n" + "2009 07 21 20 UT Moon at perigee, 222,118 miles (357,464 km)\n" + "2009 07 22 Total solar eclipse (India, China, and a few Pacific islands)\n" + "2009 07 28 02 UT Meteor shower peak -- Delta Aquarids\n" + "2009 08 04 01 UT Moon at apogee, 252,294 miles (406,028 km)\n" + "2009 08 06 00:39 UT Weak penumbral lunar eclipse (midtime)\n" + "2009 08 12 17 UT Meteor shower peak -- Perseids\n" + "2009 08 14 18 UT Jupiter at opposition\n" + "2009 08 17 21 UT Neptune at opposition\n" + "2009 08 19 05 UT Moon at perigee, 223,469 miles (359,638 km)\n" + "2009 08 24 16 UT Mercury at greatest elongation, 27.4\u00B0 east of Sun (evening)\n" + "2009 08 26 09 UT Mercury at greatest illuminated extent (evening)\n" + "2009 08 31 11 UT Moon at apogee, 251,823 miles (405,269 km)\n" + "2009 09 16 08 UT Moon at perigee, 226,211 miles (364,052 km)\n" + "2009 09 17 09 UT Uranus at opposition\n" + "2009 09 22 21:19 UT Fall (Northern Hemisphere) begins at the equinox\n" + "2009 09 28 04 UT Moon at apogee, 251,302 miles (404,432 km)\n" + "2009 10 06 02 UT Mercury at greatest elongation, 17.9\u00B0 west of Sun (morning)\n" + "2009 10 11 07 UT Mercury at greatest illuminated extent (morning)\n" + "2009 10 13 12 UT Moon at perigee, 229,327 miles (369,066 km)\n" + "2009 10 21 10 UT Meteor shower peak -- Orionids\n" + "2009 10 25 23 UT Moon at apogee, 251,137 miles (404,166 km)\n" + "2009 11 05 10 UT Meteor shower peak -- Southern Taurids\n" + "2009 11 07 07 UT Moon at perigee, 229,225 miles (368,902 km)\n" + "2009 11 12 10 UT Meteor shower peak -- Northern Taurids\n" + "2009 11 17 15 UT Meteor shower peak -- Leonids\n" + "2009 11 22 20 UT Moon at apogee, 251,489 miles (404,733 km)\n" + "2009 12 04 14 UT Moon at perigee, 225,856 miles (363,481 km)\n" + "2009 12 14 05 UT Meteor shower peak -- Geminids\n" + "2009 12 17 13 UT Mercury at greatest illuminated extent (evening)\n" + "2009 12 18 17 UT Mercury at greatest elongation, 20.3\u00B0 east of Sun (evening)\n" + "2009 12 20 15 UT Moon at apogee, 252,109 miles (405,731 km)\n" + "2009 12 21 17:47 UT Winter (Northern Hemisphere) begins at the solstice\n" + "2009 12 31 Partial lunar eclipse, 18:52 UT to 19:54 UT\n-----\n"; END Comment DAF */ document.planets.ra1.value = planet_ra[1]; document.planets.dec1.value = planet_dec[1]; document.planets.dist1.value = planet_mag[1]; document.planets.size1.value = planet_size[1]; document.planets.phase1.value = planet_phase[1]; document.planets.alt1.value = planet_alt[1]; document.planets.az1.value = planet_az[1]; document.planets.ra2.value = planet_ra[2]; document.planets.dec2.value = planet_dec[2]; document.planets.dist2.value = planet_mag[2]; document.planets.size2.value = planet_size[2]; document.planets.phase2.value = planet_phase[2]; document.planets.alt2.value = planet_alt[2]; document.planets.az2.value = planet_az[2]; document.planets.ra3.value = planet_ra[3]; document.planets.dec3.value = planet_dec[3]; document.planets.dist3.value = planet_mag[3]; document.planets.size3.value = planet_size[3]; document.planets.phase3.value = planet_phase[3]; document.planets.alt3.value = planet_alt[3]; document.planets.az3.value = planet_az[3]; document.planets.ra4.value = planet_ra[4]; document.planets.dec4.value = planet_dec[4]; document.planets.dist4.value = planet_mag[4]; document.planets.size4.value = planet_size[4]; document.planets.phase4.value = planet_phase[4]; document.planets.alt4.value = planet_alt[4]; document.planets.az4.value = planet_az[4]; document.planets.ra5.value = planet_ra[5]; document.planets.dec5.value = planet_dec[5]; document.planets.dist5.value = planet_mag[5]; document.planets.size5.value = planet_size[5]; document.planets.phase5.value = planet_phase[5]; document.planets.alt5.value = planet_alt[5]; document.planets.az5.value = planet_az[5]; document.planets.dist6.value = planet_dist[6]; document.planets.size6.value = sunrise_set; document.planets.phase6.value = phenomena; document.planets.dawn.value = dawn_planets; document.planets.dusk.value = dusk_planets; document.planets.size7.value = moonrise_set; document.planets.phase7.value = twilight_begin_end; document.planets.ra8.value = planet_ra[8]; document.planets.dec8.value = planet_dec[8]; document.planets.dist8.value = planet_mag[0]; document.planets.size8.value = planet_size[8]; document.planets.alt8.value = planet_alt[8]; document.planets.az8.value = planet_az[8]; document.planets.ra9.value = planet_ra[9]; document.planets.dec9.value = planet_dec[9]; document.planets.dist9.value = planet_mag[6]; document.planets.size9.value = planet_size[9]; document.planets.phase9.value = planet_phase[9]; document.planets.alt9.value = planet_alt[9]; document.planets.az9.value = planet_az[9]; } } function moon_rise_set(latitude,longitude,tzone) { with (Math) { var r_a = new Array(3); var dec = new Array(3); var dist = new Array(3); var i, c0, x, y, p1, p2, r1, RAD, k1, b5, l5, h, z0, t, t0; var s, l, m, f, d, n, g, v, u, w, a5, d5, r5, z1, c, z, m8, w8; var a0, d0, p, f0, f1, f2, a, b, a2, d2, l0, l2, h0, h2, d1; var v0=0, v1=0, v2=0, e, t3, rise_str = "", set_str = "", result = ""; p1 = PI; p2 = 2 * p1; r1 = PI / 180; RAD = 1 / r1; k1 = 15 * r1 * 1.002738; b5 = latitude * RAD; l5 = longitude * RAD; h = tzone; l5 = -l5 / 360; z0 = -h / 24; t = (julian_date_zoned() - 2451545.0); t0 = t / 36525; s = (24110.5 + 8640184.813 * t0 + 86636.6 * z0 + 86400 * l5) / 86400; s = s - floor(s); t0 = s * 360 * r1; t = t + z0; for (i=1; i<4; i++) { l = 0.606434 + 0.036601 * t; m = 0.374897 + 0.036292 * t; f = 0.259091 + 0.036748 * t; d = 0.827362 + 0.033863 * t; n = 0.347343 - 1.470939e-04 * t; g = 0.993126 + 2.737778e-03 * t; l = (l - floor(l)) * p2; m = (m - floor(m)) * p2; f = (f - floor(f)) * p2; d = (d - floor(d)) * p2; n = (n - floor(n)) * p2; g = (g - floor(g)) * p2; v = 0.39558*sin(f+n) + 0.082*sin(f) + 0.03257*sin(m-f-n) + 0.01092*sin(m+f+n) + 6.66e-03*sin(m-f) - 6.44e-03*sin(m+f-2*d+n) - 3.31e-03*sin(f-2*d+n) - 3.04e-03*sin(f-2*d) - 2.4e-03*sin(m-f-2*d-n) + 2.26e-03*sin(m+f) - 1.08e-03*sin(m+f-2*d) - 7.9e-04*sin(f-n) + 7.8e-04*sin(f+2*d+n); u = 1 - 0.10828*cos(m) - 0.0188*cos(m-2*d) - 0.01479*cos(2*d) + 1.81e-03*cos(2*m-2*d) - 1.47e-03*cos(2*m) - 1.05e-03*cos(2*d-g) - 7.5e-04*cos(m-2*d+g); w = 0.10478*sin(m) - 0.04105*sin(2*f+2*n) - 0.0213*sin(m-2*d) - 0.01779*sin(2*f+n) + 0.01774*sin(n) + 9.87e-03*sin(2*d) - 3.38e-03*sin(m-2*f-2*n) - 3.09e-03*sin(g) - 1.9e-03*sin(2*f) - 1.44e-03*sin(m+n) - 1.44e-03*sin(m-2*f-n) - 1.13e-03*sin(m+2*f+2*n) - 9.4e-04*sin(m-2*d+g) - 9.2e-04*sin(2*m-2*d); s = w / sqrt(u-v*v); a5 = l + atan(s / sqrt(1-s*s)); s = v / sqrt(u); d5 = atan(s / sqrt(1-s*s)); r5 = 60.40974 * sqrt(u); r_a[i] = a5; dec[i] = d5; dist[i] = r5; t += 0.5; } if (r_a[2] < r_a[1]) r_a[2] += p2; if (r_a[3] < r_a[2]) r_a[3] += p2; z1 = r1 * (90.567 - 41.685 / dist[2]); s = sin(b5*r1); c = cos(b5*r1); z = cos(z1); m8 = 0; w8 = 0; a0 = r_a[1]; d0 = dec[1]; for (c0=0; c0<24; c0++) { p = (c0 + 1) / 24; f0 = r_a[1]; f1 = r_a[2]; f2 = r_a[3]; a = f1 - f0; b = f2 - f1 - a; f = f0 + p * (2 * a + b * (2 * p - 1)); a2 = f; f0 = dec[1]; f1 = dec[2]; f2 = dec[3]; a = f1 - f0; b = f2 - f1 - a; f = f0 + p * (2 * a + b * (2 * p - 1)); d2 = f; l0 = t0 + c0 * k1; l2 = l0 + k1; if (a2 < a0) a2 += p2; h0 = l0 - a0; h2 = l2 - a2; h1 = (h2 + h0) / 2; d1 = (d2 + d0) / 2; if (c0 < 1) v0 = s * sin(d0) + c * cos(d0) * cos(h0) - z; v2 = s * sin(d2) + c * cos(d2) * cos(h2) - z; if (sign_of(v0) != sign_of(v2)) { v1 = s * sin(d1) + c * cos(d1) * cos(h1) - z; a = 2 * v2 - 4 * v1 + 2 * v0; b = 4 * v1 - 3 * v0 - v2; d = b * b - 4 * a * v0; if (d >= 0) { d = sqrt(d); e = (-b + d) / (2 * a); if (e > 1 || e < 0) e = (-b - d) / (2 * a); t3 = am_pm(c0 + e); if (v0 < 0 && v2 > 0) { m8 = 1; rise_str = t3; } if (v0 > 0 && v2 < 0) { w8 = 1; set_str = t3; } } } a0 = a2; d0 = d2; v0 = v2; } if (m8 == 0 && w8 == 0) { if (v2 < 0) { rise_str = "down 24h"; set_str = "down 24h"; } if (v2 > 0) { rise_str = " up 24h "; set_str = " up 24h "; } } else { if (m8 == 0) rise_str = " ------ "; if (w8 == 0) set_str = " ------ "; } result = rise_str + " / " + set_str; return result; } } function sun_rise_set(latitude,longitude,tzone) { with (Math) { var r_a = new Array(2); var dec = new Array(2); var p1 = PI; var p2 = 2 * p1; var dr = p1 / 180; var RAD = 1 / dr; var k1 = 15 * dr * 1.0027379; var b5, l5, h, z0, t, tt, t0, s, l, g, v, u, w; var a5, d5, r5, z1, c, z, m8, w8, a0, d0, da, dd; var c0, p, a2, d2, l0, l2, h0, h1, h2, d1; var v0=0, v1=0, v2=0, a, b, d, e, t3; var rise_str = "", set_str = "", result = ""; b5 = latitude * RAD; l5 = longitude * RAD; h = tzone; l5 = -l5 / 360; z0 = -h / 24; t = (julian_date_zoned() - 2451545.0); tt = t / 36525 + 1; t0 = t / 36525; s = (24110.5 + 8640184.813 * t0 + 86636.6 * z0 + 86400 * l5) / 86400; s = s - floor(s); t0 = s * 360 * dr; t = t + z0; l = .779072+.00273790931*t; g = .993126+.0027377785*t; l = (l-floor(l))*p2; g = (g-floor(g))*p2; v = .39785*sin(l)-.01000*sin(l-g)+.00333*sin(l+g)-.00021*tt*sin(l); u = 1-.03349*cos(g)-.00014*cos(2*l)+.00008*cos(l); w = -.00010-.04129*sin(2*l)+.03211*sin(g)+.00104*sin(2*l-g)-.00035*sin(2*l+g)-.00008*tt*sin(g); s = w/sqrt(u-v*v); a5 = l+atan(s/sqrt(1-s*s)); s = v/sqrt(u); d5 = atan(s/sqrt(1-s*s)); r5 = 1.00021*sqrt(u); r_a[1] = a5; dec[1] = d5; t = t + 1; l = .779072+.00273790931*t; g = .993126+.0027377785*t; l = (l-floor(l))*p2; g = (g-floor(g))*p2; v = .39785*sin(l)-.01000*sin(l-g)+.00333*sin(l+g)-.00021*tt*sin(l); u = 1-.03349*cos(g)-.00014*cos(2*l)+.00008*cos(l); w = -.00010-.04129*sin(2*l)+.03211*sin(g)+.00104*sin(2*l-g)-.00035*sin(2*l+g)-.00008*tt*sin(g); s = w/sqrt(u-v*v); a5 = l+atan(s/sqrt(1-s*s)); s = v/sqrt(u); d5 = atan(s/sqrt(1-s*s)); r5 = 1.00021*sqrt(u); r_a[2] = a5; dec[2] = d5; if (r_a[2] < r_a[1]) r_a[2] += p2; z1 = dr * 90.833; s = sin(b5 * dr); c = cos(b5 * dr); z = cos(z1); m8 = 0; w8 = 0; a0 = r_a[1]; d0 = dec[1]; da = r_a[2] - r_a[1]; dd = dec[2] - dec[1]; for (c0=0; c0<24; c0++) { p = (c0 + 1) / 24; a2 = r_a[1] + p * da; d2 = dec[1] + p * dd; l0 = t0 + c0 * k1; l2 = l0 + k1; h0 = l0 - a0; h2 = l2 - a2; h1 = (h2 + h0) / 2; d1 = (d2 + d0) / 2; if (c0 < 1) v0 = s * sin(d0) + c * cos(d0) * cos(h0) - z; v2 = s * sin(d2) + c * cos(d2) * cos(h2) - z; if (sign_of(v0) != sign_of(v2)) { v1 = s * sin(d1) + c * cos(d1) * cos(h1) - z; a = 2 * v2 - 4 * v1 + 2 * v0; b = 4 * v1 - 3 * v0 - v2; d = b * b - 4 * a * v0; if (d >= 0) { d = sqrt(d); e = (-b + d) / (2 * a); if (e > 1 || e < 0) e = (-b - d) / (2 * a); t3 = am_pm(c0 + e); if (v0 < 0 && v2 > 0) { m8 = 1; rise_str = t3; } if (v0 > 0 && v2 < 0) { w8 = 1; set_str = t3; } } } a0 = a2; d0 = d2; v0 = v2; } if (m8 == 0 && w8 == 0) { if (v2 < 0) { rise_str = "down 24h"; set_str = "down 24h"; } if (v2 > 0) { rise_str = " up 24h "; set_str = " up 24h "; } } else { if (m8 == 0) rise_str = " ------ "; if (w8 == 0) set_str = " ------ "; } result = rise_str + " / " + set_str; return result; } } function getCookie( name ) { var start = document.cookie.indexOf( name + "=" ); var len = start + name.length + 1; if ( ( !start ) && ( name != document.cookie.substring( 0, name.length ) ) ) { return null; } if ( start == -1 ) return null; var end = document.cookie.indexOf( ";", len ); if ( end == -1 ) end = document.cookie.length; return unescape( document.cookie.substring( len, end ) ); } function setCookie( name, value, expires, path, domain, secure ) { var today = new Date(); today.setTime( today.getTime() ); if ( expires ) { expires = expires * 1000 * 60 * 60 * 24; } var expires_date = new Date( today.getTime() + (expires) ); document.cookie = name+"="+escape( value ) + ( ( expires ) ? ";expires="+expires_date.toGMTString() : "" ) + //expires.toGMTString() ( ( path ) ? ";path=" + path : "" ) + ( ( domain ) ? ";domain=" + domain : "" ) + ( ( secure ) ? ";secure" : "" ); } function deleteCookie( name, path, domain ) { if ( getCookie( name ) ) document.cookie = name + "=" + ( ( path ) ? ";path=" + path : "") + ( ( domain ) ? ";domain=" + domain : "" ) + ";expires=Thu, 01-Jan-1970 00:00:01 GMT"; } function save_loc() { deleteCookie('almanac2latdeg','',''); deleteCookie('almanac2longdeg','',''); deleteCookie('almanac2tzonehours','',''); setCookie('almanac2latdeg',round_1000(latitude),365,'','',''); setCookie('almanac2longdeg',round_1000(longitude),365,'','',''); setCookie('almanac2tzonehours',t_zone,365,'','',''); alert("The current location has been saved."); } function load_loc() { var full_custom_URL = "http://www.skyandtelescope.com/observing/almanac/almanacCustom?"; var display_saved_almanac; var lat_deg = getCookie('almanac2latdeg'); var lng_deg = getCookie('almanac2longdeg'); var time_zone_hours = getCookie('almanac2tzonehours'); if (lat_deg == null && lng_deg == null) { alert ("Sorry, but you haven't saved a previous latitude and longitude. If you want the Almanac to remember the currently displayed location, please click the 'Save this?' link. Alternatively, enter a new latitude and longitude by clicking the 'New?' link.") } else { full_custom_URL += "latitude=" + lat_deg + "&"; full_custom_URL += "longitude=" + lng_deg + "&"; full_custom_URL += "tzone=" + time_zone_hours + "&"; full_custom_URL += "UTdate=" + document.planets.date_txt.value + "&"; full_custom_URL += "UTtime=" + document.planets.ut_h_m.value; display_saved_almanac = window.open(full_custom_URL,'chooser'); } }