/* Almanac v3.0 JavaScript programming and graphics, copyright: Adrian R. Ashford, 3rd January 2011. Used with permission by SkyandTelescope.com Parts of sunrise, sunset, moonrise and moonset from moonup.bas and sunup.bas by Roger W. Sinnott, Sky and Telescope magazine. ISS predictions based on MYBIRDSa.BAS by Roger W. Sinnott, 3rd June 2001. 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("images/moon_phases/mp0.gif","images/moon_phases/mp1.gif","images/moon_phases/mp2.gif","images/moon_phases/mp3.gif","images/moon_phases/mp4.gif","images/moon_phases/mp5.gif","images/moon_phases/mp6.gif","images/moon_phases/mp7.gif","images/moon_phases/mp8.gif","images/moon_phases/mp9.gif","images/moon_phases/mp10.gif","images/moon_phases/mp11.gif","images/moon_phases/mp12.gif","images/moon_phases/mp13.gif","images/moon_phases/mp14.gif","images/moon_phases/mp15.gif","images/moon_phases/mp16.gif","images/moon_phases/mp17.gif","images/moon_phases/mp18.gif","images/moon_phases/mp19.gif","images/moon_phases/mp20.gif","images/moon_phases/mp21.gif","images/moon_phases/mp22.gif","images/moon_phases/mp23.gif","images/moon_phases/mp24.gif","images/moon_phases/mp25.gif","images/moon_phases/mp26.gif","images/moon_phases/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-------------------------------------------------------------------------\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-------------------------------------------------------------------------\n"; phenomena += moon_data; ISS = document.planets.iss.checked; if (ISS == true) { phenomena += satellite_calc(); } else { phenomena +="-------------------------------------------------------------------------\n"; } phenomena += "2011 Phenomena\n"; phenomena += "yyyy mm dd hh UT All dates and times are UT\n" + "2011 01 03 19 UT Earth at perihelion, 91,407,282 miles (147,105,760 km)\n" + "2011 01 04 01 UT Meteor shower peak -- Quadrantids\n" + "2011 01 04 Partial solar eclipse (Europe, northern Africa, western Asia)\n" + "2011 01 08 04 UT Mercury at greatest illuminated extent (morning)\n" + "2011 01 08 16 UT Venus at greatest elongation, 47.0\u00B0 west of Sun (morning)\n" + "2011 01 09 15 UT Mercury at greatest elongation, 23.3\u00B0 west of Sun (morning)\n" + "2011 01 10 06 UT Moon at apogee, 251,641 miles (404,977 km)\n" + "2011 01 22 00 UT Moon at perigee, 225,428 miles (362,792 km)\n" + "2011 02 06 23 UT Moon at apogee, 252,229 miles (405,924 km)\n" + "2011 02 19 07 UT Moon at perigee, 222,604 miles (358,247 km)\n" + "2011 03 06 08 UT Moon at apogee, 252,639 miles (406,583 km)\n" + "2011 03 14 10 UT Mercury at greatest illuminated extent (evening)\n" + "2011 03 19 19 UT Moon at perigee, 221,565 miles (356,575 km)\n" + "2011 03 20 23:21 UT Autumn (Southern Hemisphere) begins at the equinox\n" + "2011 03 23 01 UT Mercury at greatest elongation, 18.6\u00B0 east of Sun (evening)\n" + "2011 04 02 09 UT Moon at apogee, 252,684 miles (406,656 km)\n" + "2011 04 04 00 UT Saturn at opposition\n" + "2011 04 17 06 UT Moon at perigee, 222,507 miles (358,090 km)\n" + "2011 04 22 23 UT Meteor shower peak -- Lyrids\n" + "2011 04 29 18 UT Moon at apogee, 252,301 miles (406,039 km)\n" + "2011 05 06 13 UT Meteor shower peak -- Eta Aquarids\n" + "2011 05 07 19 UT Mercury at greatest elongation, 26.6\u00B0 west of Sun (morning)\n" + "2011 05 08 20 UT Mercury at greatest illuminated extent (morning)\n" + "2011 05 15 11 UT Moon at perigee, 225,020 miles (362,135 km)\n" + "2011 05 27 10 UT Moon at apogee, 251,657 miles (405,003 km)\n" + "2011 06 01 Partial solar eclipse (northern North America, Greenland, Iceland, northeastern Asia)\n" + "2011 06 11 02 UT Mercury at greatest illuminated extent (morning)\n" + "2011 06 12 02 UT Moon at perigee, 228,161 miles (367,189 km)\n" + "2011 06 15 Total lunar eclipse, 19:22 UT to 21:03 UT\n" + "2011 06 21 17:16 UT Winter (Southern Hemisphere) begins at the solstice\n" + "2011 06 24 04 UT Moon at apogee, 251,203 miles (404,271 km)\n" + "2011 07 01 Partial solar eclipse (ocean south of Africa)\n" + "2011 07 04 15 UT Earth at aphelion, 94,511,923 miles (152,102,197 km)\n" + "2011 07 07 14 UT Moon at perigee, 229,640 miles (369,570 km)\n" + "2011 07 19 11 UT Mercury at greatest illuminated extent (evening)\n" + "2011 07 20 05 UT Mercury at greatest elongation, 26.8\u00B0 east of Sun (evening)\n" + "2011 07 21 23 UT Moon at apogee, 251,254 miles (404,355 km)\n" + "2011 07 28 14 UT Meteor shower peak -- Delta Aquarids\n" + "2011 08 02 21 UT Moon at perigee, 227,273 miles (365,761 km)\n" + "2011 08 13 06 UT Meteor shower peak -- Perseids\n" + "2011 08 18 16 UT Moon at apogee, 251,755 miles (405,161 km)\n" + "2011 08 22 23 UT Neptune at opposition\n" + "2011 08 30 18 UT Moon at perigee, 224,227 miles (360,858 km)\n" + "2011 09 03 06 UT Mercury at greatest elongation, 18.1\u00B0 west of Sun (morning)\n" + "2011 09 11 07 UT Mercury at greatest illuminated extent (morning)\n" + "2011 09 15 06 UT Moon at apogee, 252,317 miles (406,065 km)\n" + "2011 09 23 09:05 UT Spring (Southern Hemisphere) begins at the equinox\n" + "2011 09 26 00 UT Uranus at opposition\n" + "2011 09 28 01 UT Moon at perigee, 222,176 miles (357,557 km)\n" + "2011 10 12 12 UT Moon at apogee, 252,546 miles (406,434 km)\n" + "2011 10 21 22 UT Meteor shower peak -- Orionids\n" + "2011 10 26 12 UT Moon at perigee, 221,862 miles (357,052 km)\n" + "2011 10 29 02 UT Jupiter at opposition\n" + "2011 11 05 23 UT Meteor shower peak -- Southern Taurids\n" + "2011 11 08 13 UT Moon at apogee, 252,387 miles (406,177 km)\n" + "2011 11 12 22 UT Meteor shower peak -- Northern Taurids\n" + "2011 11 14 08 UT Mercury at greatest elongation, 22.7\u00B0 east of Sun (evening)\n" + "2011 11 15 10 UT Mercury at greatest illuminated extent (evening)\n" + "2011 11 18 04 UT Meteor shower peak -- Leonids\n" + "2011 11 23 23 UT Moon at perigee, 223,502 miles (359,691 km)\n" + "2011 11 25 Partial solar eclipse (Antarctica, southernmost Africa, Tasmania)\n" + "2011 12 06 01 UT Moon at apogee, 251,912 miles (405,414 km)\n" + "2011 12 10 Total lunar eclipse, 14:06 to 14:58 UT\n" + "2011 12 14 17 UT Meteor shower peak -- Geminids\n" + "2011 12 22 03 UT Moon at perigee, 226,676 miles (364,800 km)\n" + "2011 12 22 05:30 UT Summer (Southern Hemisphere) begins at the solstice\n" + "2011 12 22 21 UT Mercury at greatest illuminated extent (morning)\n" + "2011 12 23 04 UT Mercury at greatest elongation, 21.8\u00B0 west of Sun (morning)\n-------------------------------------------------------------------------\n"; 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; for (c0=1; c0<3; c0++) { 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[c0] = a5; dec[c0] = d5; t += 1; } 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 satellite_calc() { with (Math) { PROCinit(); PROCheader(); T = JDS - JDE; do { PROCsatvec(); PROCrangevec(); PROCsunvec(); if (EL >= SatTest && ECL == " vis") { PROCassemble(); } else { if (Pflag == 1) { dur = floor((Tend - Tstart) * 1440 + 0.5); if (dur >= 1) { output_display += " " + dur + " "; } else { output_display += " <1 "; } output_display += P1 + " " + P2 + " " + P3 + "\n"; Pcount += 1; Pflag = 0; } } T += Stepday; } while(Pcount < Plimit && (JDE + T - JDS) < Looktime); if (Pcount == 0) { output_display += "Sorry, no visible ISS passes in the next " + Looktime + " days.\n" } output_display += "-------------------------------------------------------------------------\n"; return output_display; } } function PROCinit() { with (Math) { if (iss_warning_flag == 0) { alert("Please use this option sparingly and remember to uncheck it when not required!\n\nChoosing to display International Space Station predictions will make the Almanac very sluggish on some computers owing to the huge number of additional calculations.\n\nYou may experience an 'unresponsive script' error (or similar); just click 'continue' if this happens."); iss_warning_flag += 1; } output_display = ""; DR = PI / 180; RD = 1 / DR; azimuth = new Array("N ","NNE","NE ","ENE","E ","ESE","SE ","SSE","S ","SSW","SW ","WSW","W ","WNW","NW ","NNW"); Stepday = 1 / 1440; Looktime = 3; Plimit = 5; Agelimit = 10; SunTest = -6; SatTest = 10; iss_elements(); LA = latitude; LO = -longitude; HT = 10; TZ = t_zone; d_s_time = document.planets.d_s_t.checked; if (d_s_time == true) { TZ += 1; } start_date = new Date(); YR = start_date.getUTCFullYear(); MN = start_date.getUTCMonth() + 1; DY = start_date.getUTCDate(); HR = start_date.getUTCHours(); MIN = start_date.getUTCMinutes(); JDS = caljd(YR,MN,DY + HR / 24 + MIN / 1440); LA = DR * LA; LO = DR * LO; HT = HT / 1000; CL = cos(LA); SL = sin(LA); CO = cos(LO); SO = sin(LO); RE = 6378.140; FL = 1 / 298.257; RP = RE * (1 - FL); XX = RE * RE; ZZ = RP * RP; D = sqrt(XX * CL * CL + ZZ * SL * SL); Rx = XX / D + HT; Rz = ZZ / D + HT; Ux = CL * CO; Ex = -SO; Nx = -SL * CO; Uy = CL * SO; Ey = CO; Ny = -SL * SO Uz = SL; Ez = 0; Nz = CL; Ox = Rx * Ux; Oy = Rx * Uy; Oz = Rz * Uz; RA = DR * RA; INC = DR * INC; WP = DR * WP; MA = DR * MA; MM = MM * 2 * PI; M2 = M2 * 2 * PI; JDE = caljd(YE,1,0) + TE; AGE = JDS - JDE; if (AGE > Agelimit) { alert("Warning! ISS orbital elements are getting out of date.\nAge = " + floor(AGE) + " days."); } GM = 3.986E+5; J2 = 1.08263E-3; N0 = MM / 86400; A0 = pow(GM / (N0 * N0),1 / 3); B0 = A0 * sqrt(1 - EC * EC); SI = sin(INC); CI = cos(INC); PC = RE * A0 / (B0 * B0); PC = 1.5 * J2 * PC * PC * MM; QD = -PC * CI; WD = PC * (5 * CI * CI - 1) / 2; DC = -2 * M2 / MM / 3; Pflag = 0; Pcount = 0; } } function PROCheader() { output_display += "-------------------------------------------------------------------------\n"; output_display += "Current nighttime ISS passes visible at your chosen location:\n"; output_display += " Local Duration\n"; output_display += " Date Time (min) Approach Max. elevation Departure\n"; output_display += "- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -\n"; } function PROCsatvec() { with (Math) { DT = DC * T / 2; KD = 1 + 4 * DT; KDP = 1 - 7 * DT; M = MA + MM * T * (1 - 3 * DT); DR0 = floor(M / (2 * PI)); M = M - DR0 * 2 * PI; RN = RV + DR0; EA = M; for (count = 1; count <= 10; count++) { E1 = M + EC * sin(EA); if (abs(E1-EA) < 1E-8) break; EA = E1; } A = A0 * KD; B = B0 * KD; RS = A * (1 - EC * cos(EA)); Sx = A * (cos(EA) - EC); Sy = B * sin(EA); AP = WP + WD * T * KDP; CW = cos(AP); SW = sin(AP); RAAN = RA + QD * T * KDP; CQ = cos(RAAN); SQ = sin(RAAN); CXx = CW * CQ - SW * CI * SQ; CXy = -SW * CQ - CW * CI * SQ; CXz = SI * SQ; CYx = CW * SQ + SW * CI * CQ; CYy = -SW * SQ + CW * CI * CQ; CYz = -SI * CQ; CZx = SW * SI; CZy = CW * SI; CZz = CI; SATx = Sx * CXx + Sy * CXy; SATy = Sx * CYx + Sy * CYy; SATz = Sx * CZx + Sy * CZy; H = gmst(JDE + T); GHAA = (H * 15) * DR; C = cos(-GHAA); S = sin(-GHAA); Sx = SATx * C - SATy * S; Sy = SATx * S + SATy * C; Sz = SATz; C = cos(GHAA); S = sin(GHAA); Ocx = Ox * C - Oy * S; Ocy = Ox * S + Oy * C; Ocz = Oz; } } function PROCrangevec() { with (Math) { Rx = Sx - Ox; Ry = Sy - Oy; Rz = Sz - Oz; R = sqrt(Rx * Rx + Ry * Ry + Rz * Rz); Rx = Rx / R; Ry = Ry / R; Rz = Rz / R; U = Rx * Ux + Ry * Uy + Rz * Uz; East = Rx * Ex + Ry * Ey; North = Rx * Nx + Ry * Ny + Rz * Nz; AZ0 = RD * arc_tan(East,North); EL = RD * atan(U / sqrt(1 - U * U)); } } function SUNxyz(j_date) { with (Math) { var n, ll, g, le, ob; n = j_date - 2451545.0; ll = (280.460 + 0.9856474 * n) / 360; ll = (ll - floor(ll)) * 2 * PI; g = (357.528 + 0.9856003 * n) / 360; g = (g - floor(g)) * 2 * PI; le = ll + 1.915 * DR * sin(g) + 0.020 * DR * sin(2 * g); ob = (23.439 - 0.0000004 * n) * DR; SUNx = cos(le); SUNy = cos(ob) * sin(le); SUNz = sin(ob) * sin(le); } } function PROCsunvec() { with (Math) { SUNxyz(JDE + T); CUA = -(SATx * SUNx + SATy * SUNy + SATz * SUNz) / RS; UMD = RS * sqrt(1 - CUA * CUA) / RE; if (CUA >= 0) {ECL = " + ";} else {ECL = " - ";} if (UMD <= 1 && CUA >= 0) ECL = " ECL"; C = cos( - GHAA); S = sin( - GHAA); Hx = SUNx * C - SUNy * S; Hy = SUNx * S + SUNy * C; Hz = SUNz; cSunZD = Hx * Ux + Hy * Uy + Hz * Uz; SunZD = atan(sqrt(1 - cSunZD * cSunZD) / cSunZD); if (cSunZD < 0) SunZD += PI; SunAlt = 90 - SunZD * RD; if (SunAlt < SunTest && ECL != " ECL") ECL = " vis"; } } function PROCassemble() { with (Math) { jd_temp = JDE + T + TZ / 24 + 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); dy = bb - dd - floor(30.6001 * ee) + ff; mth = ee; if (ee < 13.5) mth = ee - 1; if (ee > 13.5) mth = ee - 13; yr = cc; if (mth > 2.5) yr = cc - 4716; if (mth < 2.5) yr = cc - 4715; HR1 = (dy - floor(dy)) * 24; HR = floor(HR1); MIN = floor(60 * (HR1 - HR)); PM = " pm"; if (HR < 12) PM = " am"; if (HR > 12) HR = HR - 12; if (HR == 0) HR = 12; if (HR < 10) HR = "0" + HR; if (MIN < 10) MIN = "0" + MIN; time = " " + HR + ":" + MIN + PM; dy = floor(dy); if (mth < 10) mth = "0" + mth; if (dy < 10) dy = "0" + dy; date = yr + "/" + mth + "/" + dy + " "; AZ0 += 11.25; if (AZ0 >= 360) AZ0 -= 360; AZI = floor(AZ0 / 22.5); AZ = azimuth[AZI]; P0 = floor(EL + 0.5) + "\u00B0" + " above " + AZ; if (EL < 9.5) P0 = " " + P0; if (Pflag == 1) { if (EL > EL2) { EL2 = EL; P2 = P0; } P3 = P0; Tend = T; } else { output_display += date + time; P1 = P0; Tstart = T; P2 = P0; EL2 = EL; P3 = P0; Tend = T; } Pflag = 1; } } function caljd(year,month,day) { with (Math) { var y0, m0, a, b, b1; y0 = year; m0 = month; if (m0 < 2.5) { y0 = y0 - 1; m0 = m0 + 12; } a = floor(y0 / 100); b = 0; b1 = 2 - a + floor(a / 4); if (y0 > 1582) b = b1; if (y0 == 1582 && month > 10) b = b1; if (y0 == 1582 && month == 10 && day >= 15.0) b = b1; return floor(365.25 * (y0 + 4716)) + floor(30.6001 * (m0 + 1)) + day + b - 1524.5; } } function gmst(jd) { with (Math) { var a, b, c, d, k, g; a = 0.7790572732639; b = 1.002737909350795; c = 8.07745E-16; d = 1.47E-24; k = jd - 2451545.0; g = a + k * (b + k * (c + k * d)); return 24 * (g - floor(g)); } } function arc_tan(y,x) { with (Math) { var a; if (x != 0) {a = atan(y / x);} else {a = PI / 2 * sign_of(y);} if (x < 0) a += PI; if (a < 0) a += 2 * PI; return a; } } 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 = "custom.html?"; 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,'_self'); } }