// Javascript Astrotools // FUNCTIONS FOR CREATING CHARTS // Copyright Ole Nielsen 2002-2004 // Please read copyright notice in astrotools2.html source SUNLIGHT = 0; PLANETLIGHT = 1; MOONLIGHT = 5; // indices 2-4 no longer used // first array stores colours used for visibility diagrams // as RGB triplets of value RRGGBB, eg. #f60 means red=FF, green=66, blue=00 (orange colour) var light = new Array(); light[SUNLIGHT] = ["#ccd","#88f","#44c","#00a","#000"]; light[PLANETLIGHT] = ["#eef","#8b8","#0b0","#0d0","#0f0"]; light[MOONLIGHT]= ["#eef","#bb9","#cc4","#ee0","#ff0"]; // light[MOONLIGHT]= ["#eef","#cc8","#bb0","#dd0","#ff0"]; var psyms = new Array("psym1.png","psym2.png","psym3.png","psym4.png","psym5.png","psym6.png","psym7.png","psym8.png","","psymsol.png","psymmoon.png","","","","","psym.png","","","","","psym.png"); var head1="\n"; var head3="" + "

close window

"; function box(x,y,w,h,bgr) { // draw box with background parameter (colour and/or url as CSS) var str="
"; // the 'font-size' parameter needed due to IE6 bug for small boxes return str; } function writexlbl(leftoff,topoff,lbl,lblpos,xlbl,xpos) { // write x-axis labels var str = "

"+xlbl+"

"; for (var i=0; i"+lbl[i]+"

"; } return str; } function txt_lft(x,y,s) { // left adjusted text return "

" + s + "

"; } function txt_rgt(x,y,w,rpad,s) { // right adjusted text, w is necessary width, rpad space from right var str = "
"; str += "

" + s + "

\n"; return str; } function combineEvents(obj, jday, obs, transit) { // Combines events of Sun and chosen object into a single array // Comprises times as fractions of a day from jday and indices into 'lights' arrays // last two entries store transit times in first index (-1 if not valid) var sunevents = findEvents(SUN,jday,obs); // first entry is state at t=0; var objevents = findEvents(obj,jday,obs); // same here var events = new Array(); var k = 0; var i=1; var n = 1; // counters for events and objevents arrays var tr = [-1,-1]; var trcnt=0; // remember up to two transit times var twl = sunevents[0][1]; // kind of twilight var up = (objevents[0][1]==0); // true if object is up var otype = (up && (obj!=SUN)) ? (obj==MOON?MOONLIGHT:PLANETLIGHT) : SUNLIGHT; events[k++] = new Array(sunevents[i][0],otype,twl); while (true) { var to = objevents[n][0]; var ts = sunevents[i][0] while ((ts = sunevents[i][0]) < to) { if (sunevents[i][1] == 0) { // skip Sun transit i++; continue; } twl = sunevents[i][1]; twl = (twl<0 ? -twl-1 : twl); events[k++] = new Array(ts,otype,twl); i++; } if (to >= 1.0) break; if (objevents[n][1]==0) { // check for object transit tr[trcnt++] = to; n++; continue; } up = !up; // object event must be rise or set otype = (up && (obj!=SUN)) ? (obj==MOON?MOONLIGHT:PLANETLIGHT) : SUNLIGHT; events[k++] = new Array(to,otype,twl); n++; } events[k] = new Array(1.0, -1, -1); events[k+1] = new Array(tr[0],0,0); events[k+2] = new Array(tr[1],0,0); return events; } // end combineEvents() function makeBar(obs, obj, jday, transit,x,y) { var str=""; var events = combineEvents(obj, jday, obs, transit); var el = events.length; var len = 480; // = 3 min per pixel var h = 13; var l1 = 0; var l; for (var i = 0; i < el-3; i++) { l = Math.round( events[i+1][0]*len ); if (l-l1>0) str += box(x+l1,y+1,(l-l1),h,light[events[i][1]][events[i][2]]); l1 = l; } for (var j=0;j<2;j++) { // plot red markers for transits if (transit && events[el-2+j][0] >= 0) str += box(x+events[el-2+j][0]*len,y+1,1,h,"#f00"); } return str; } // end makeBar() // function nextDate(obs1,dstep,origday) *** moved to makelist.js function doVisibility(obs,obj,dspan,dstep,transit) { // Create visibility diagrams for one object // obs is a reference variable, make a copy var obscopy=new Object(); var obsmax=new Object(); for (var i in obs) {obscopy[i] = obs[i]; obsmax[i] = obs[i]; } if (dstep<1.0) dstep=1.0; obscopy.hours=12; obscopy.minutes=0; // graphics only allows for 12:00 var pwin=window.open("","moonlight","menubar,scrollbars,resizable"); var doc=pwin.document; var str = head1 + "Javascript: Visibility" + head3 + "\n

Visibility

of "+bodies[obj].name+"

"; str += "

Observer location: "+sitename() + " (UT "+hmstring(-obs.tz/60.0,true)+")

\n"; nextDate(obsmax, dspan, obs.day); // 'abuse' nextdate to calculate end time jdmax = jd(obsmax); var jday=jd(obscopy); var tmax = Math.floor((jdmax-jday)/dstep); // how many lines? if (dstep<0) tmax = Math.floor(dspan/dstep); var leftoff=130; var topoff=30; str += "
"; str += box(leftoff-1,topoff,482,15*tmax+2,"#aaa url(hourline.png)"); var lbl = ["12","14","16","18","20","22","00","02","04","06","08","10","12"]; var lblpos = [-5,35,75,115,155,195,235,275,315,355,395,435,475]; str += writexlbl(leftoff,topoff,lbl,lblpos,"hours",225); str += txt_lft(leftoff+488,topoff-30,"Illum."); for (var t=0; t  = transit.

"; str += "

close window

"; str += "
\n"; doc.write(str); doc.close(); pwin.focus(); } // end doVisibility() function doPlanetVisibility(obs,transit) { // Create the diagram showing all planets for one date var obscopy=new Object(); for (var i in obs) obscopy[i] = obs[i]; obscopy.hours=12; obscopy.minutes=0; // graphics only allows for 12:00 var pwin=window.open("","planetvis","menubar,scrollbars,resizable"); var doc=pwin.document; var str = head1 + "Javascript: Visibility" + head3 + "

Visibility

of Sun, Moon and Planets

"; str += "

Observer location: " +sitename()+ " (UT "+hmstring(-obs.tz/60.0,true)+")

"; str += "

Date: "+datestring(obscopy)+"

\n"; var objects=[SUN,MOON,0,1,3,4,5,6,7]; var leftoff=120; var topoff=30; str += "
"; str += box(leftoff-1,topoff,482,15*9+2,"#aaa url(hourline.png)"); var lbl = ["12","14","16","18","20","22","00","02","04","06","08","10","12"]; var lblpos = [-5,35,75,115,155,195,235,275,315,355,395,435,475]; str += writexlbl(leftoff,topoff,lbl,lblpos,"hours",225); for (var t=0; t  = transit.

"; str += "

close window

"; str += "\n"; doc.write(str); doc.close(); pwin.focus(); } // end doPlanetVisibility() function doStars(obs,deepsky, minalt, rasort, transit) { // Show visible stars or deep sky objects on selected date var obscopy=new Object(); for (var i in obs) obscopy[i] = obs[i]; obscopy.hours=12; obscopy.minutes=0; // graphics only allows for start at 12:00 // stepsize in Julian day var stepsize=1/96.0; var ord = new Array(); // for storing records of siderial time and index to catalogue var sid = local_sidereal(obscopy) + 12*1.002737; // sidereal time for following midnight var objcnt = 0; // calculate number of objects to show for (var t=0; t<(deepsky?dso.length:stars.length); t++) { var de = parsecol(deepsky?dso[t].de:stars[t].de); if (de >= obs.latitude-90+minalt && de <= obs.latitude+90-minalt) { var h = parsecol(deepsky?dso[t].ra:stars[t].ra) - sid; // negative hour angle if (h<0) h+=24; if (h>12) h-=24; ord[objcnt++] = new Array(h,t); // } } if (rasort) isort(ord); // sort according to transit time starting from north var pwin=window.open("","starvis","menubar,scrollbars,resizable"); var doc=pwin.document; var ostr = (deepsky?"Deep Sky Objects":"Stars"); var str = head1 + "Javascript: Visibility" + head3 + "

Visibility

" + "

of "+ostr+"

"; str += "

Observer location: "+sitename(); str += " (UT "+hmstring(-obs.tz/60.0,true)+")

\n"; str += "

Altitude limit: " + minalt + "°

\n"; str += "

Observation date: " + datestring(obs) + "

"; var leftoff=130; var topoff=30; str += "
"; str += box(leftoff-1,topoff,482,15*objcnt+2,"#aaa url(hourline.png)"); var lbl = ["12","14","16","18","20","22","00","02","04","06","08","10","12"]; var lblpos = [-5,35,75,115,155,195,235,275,315,355,395,435,475]; str += writexlbl(leftoff,topoff,lbl,lblpos,"hours",225); var jday = jd(obscopy); for (var t=0; t  = transit.

"; str += "

close window

"; str += "\n"; doc.write(str); doc.close(); pwin.focus(); } // end doStars() function doDataGrph(obs,obj,dspan,dstep) { var obscopy=new Object(); var obsmax=new Object(); for (var i in obs) {obscopy[i] = obs[i]; obsmax[i] = obs[i]; } var pwin=window.open("","illumdiam","menubar,scrollbars,resizable"); var doc=pwin.document; var str = head1 + "Javascript: Object Data"+ head3; str += "

Object Data

(Angular diameter, magnitude and illumination)

of "+bodies[obj].name+"

\n"; str += "

Observer location: "+sitename() + " (UT "+hmstring(-obs.tz/60.0,true)+")

\n"; nextDate(obsmax, dspan, obs.day); // 'abuse' nextdate to calculate end time jdmax = jd(obsmax); var jday=jd(obscopy); var tmax = Math.floor((jdmax-jday)/dstep); if (dstep<0) tmax = Math.floor(dspan/dstep); var leftoff=130; var topoff=60; str += "
"; str += box(leftoff,topoff,502,15*(tmax-1)+2,"#aaa url(grid.png)"); lblillum = ["0%","10%","20%","30%","40%","50%","60%","70%","80%","90%","100%"]; lblipos= [-3,42,92,142,192,242,292,342,392,442,492]; if (obj==SUN || obj==MOON) { lbldiam = ["1600\"","1700\"","1800\"","1900\"","2000\"","2100\""]; lbldpos = [-12,88,188,288,388,488]; } else { lbldiam = ["0\"","10\"","20\"","30\"","40\"","50\""]; lbldpos = [-3,95,195,295,395,495]; } lblmag = ["14.0","12.0","10.0","8.0","6.0","4.0","2.0","0.0","-2.0","-4.0","-6.0"]; lblmagpos= [-8,42,92,144,194,244,294,344,392,442,492]; str += writexlbl(leftoff,topoff-12,lblmag,lblmagpos,"",180); str += writexlbl(leftoff,topoff,lblillum,lblipos,"",180); str += writexlbl(leftoff,topoff-24,lbldiam,lbldpos,"illumination (%) / magnitude / diameter (arcsec)",140); // str += txt_lft(leftoff+612,topoff-30,ylbl); for (var t=0; t  = Diameter, "; str += "  = Illumination, "; str += "  = Magnitude

"; str += "

close window

\n"; doc.write(str); doc.close(); pwin.focus(); } // end doDataGrph() function doTwilightVisibility(obs,obj,dspan,dstep,sunalt) { // Altitude of object(s) when Sun 6 degrees below horizon var obscopy=new Object(); var obsmax=new Object(); for (var i in obs) {obscopy[i] = obs[i]; obsmax[i] = obs[i]; } if (dstep<1.0) dstep=1.0; var objects = (obj==100 ? [7,6,5,4,3,1,0,10] : [obj]); /* order of planets */ var values=new Array(44); // stepsize in degrees var stepsize=1.0; // Now make the diagram var pwin=window.open("","mercvenus","menubar,scrollbars,resizable"); var doc=pwin.document; if (obj>20) var ostr="Moon and Planets"; else var ostr=bodies[obj].name; var str = head1 + "Javascript: Twilight Altitude" + head3 + "

Twilight Altitude

"; str += "

of "+ostr+"

\n"; str += "

at the time when the Sun is " + (-sunalt) + "° below the horizon

"; str += "

Observer location: " + sitename() + " (UT "+hmstring(-obs.tz/60.0,true)+")

\n"; nextDate(obsmax, dspan, obs.day); // 'abuse' nextdate to calculate end time jdmax = jd(obsmax); var jday=jd(obscopy); var tmax = Math.floor((jdmax-jday)/dstep); // needed for calculating box sizes if (dstep<0) tmax = Math.floor(dspan/dstep); var leftoff=80; var off2 = 320; var topoff=30; str += "
"; str += box(leftoff,topoff,221,15*(tmax-1)+2,"#aaa url(grid.png)"); str += box(leftoff+off2,topoff,221,15*(tmax-1)+2,"#aaa url(grid.png)"); var lbl = new Array("0°","10°","20°","30°","40°"); var lblpos = new Array(-3,45,95,145,195); var lblpos2 = new Array(-3+off2,45+off2,95+off2,145+off2,195+off2); str += writexlbl(leftoff,topoff,lbl,lblpos,"Altitude (morning)",70); str += writexlbl(leftoff,topoff,lbl,lblpos2,"Altitude (evening)",70+off2); // str += txt_lft(leftoff+612,topoff-30,"az"); for (var t=0; t=0 && h<45) { str += "\n"; } if (obj<100) str += txt_lft(leftoff+i*off2+226,topoff-6+15*t,(bodies[p].az<180) ? "rising" : "setting"); } } nextDate(obscopy,dstep, obs.day); } str += "
\n"; if (obj==100) { str += "

M=Moon, 1=Mercury, 2=Venus, 4=Mars, 5=Jupiter, 6=Saturn, 7=Uranus, 8=Neptune

\n"; } str += "

close window

"; doc.write(str); doc.close(); pwin.focus(); } // end doTwilightVisibility() function doAltitude(obs,obj,mstep) { // altitude of one or more objects during one day, if obj=100 plot all planets and Sun/Moon var obscopy=new Object(); // make working copy for (var i in obs) obscopy[i] = obs[i]; obscopy.minutes = 0; // start at full hour for nice display // order of planets, later ones plot on top of earlier ones var objects = (obj==100 ? [7,6,5,4,3,1,0,10,9] : [obj]); // dstep in julian days var dstep=mstep/1440; if (obj>20) var ostr="Sun, Moon and Planets"; else var ostr=bodies[obj].name; var pwin=window.open("","pl_altitude","menubar,scrollbars,resizable"); var doc=pwin.document; var str = head1 + "Javascript: Altitude" + head3 + "

Altitude

of " + ostr + "

"; str += "

Observer location: " + sitename(); str += " (UT " + hmstring(-obs.tz/60.0,true) + ")

\n"; doc.write(str); var leftoff=120; var topoff=30; var htot = 15*1440/mstep; str = "
"; var jday=jd(obscopy); // shade according to sun up/down var ev = findEvents(SUN,jday,obs); var t0=0; var vis = (ev[0][1]==0); var i=1; // vis == true if sun up while (ev[i][0]<1.0) { if ((vis && (ev[i][1]==1)) || (!vis && (ev[i][1]==(-1)))) { str += box(leftoff+100,Math.round(topoff+t0*htot),452,(ev[i][0]-t0)*htot+1,(vis) ? "#ccd" : "#58f"); t0 = ev[i][0]; vis = !vis; } i++; } str += box(leftoff+100,Math.round(topoff+t0*htot),452,(1.0-t0)*htot+1,(vis) ? "#ccd" : "#58f"); str += box(leftoff,topoff,100,htot+2,"#888"); // paint box below horizon // str += box(leftoff+100,topoff,452,(15*1440/mstep+2),"#ddd"); // above horizon str += box(leftoff,topoff,552,htot+2,"url(grid.png)"); // transparent grid var lbl = new Array("-20°","-10°","0°","10°","20°","30°","40°","50°","60°","70°","80°","90°"); var lblpos = new Array(-10,40,97,144,194,244,294,344,394,444,494,544); str += writexlbl(leftoff,topoff,lbl,lblpos,"Altitude",250); str += txt_lft(leftoff+562,topoff-30,"az"); // for each mstep min do for (var t=0;t<1.0/dstep+0.001;t++) { str += txt_rgt(0,(topoff-5+15*t),114,6,datestring(obscopy)+"  "+hmstring2(obscopy.hours,obscopy.minutes,0)); bodies[9].update(jday,obs); for (var n in objects) { var p=objects[n]; bodies[p].update(jday,obs); var h=bodies[p].alt; if (h>=-20.0) { str += "\n"; } } if (obj<21) { // write azimuth str += txt_lft(leftoff+562,topoff-5+15*t,Math.round(bodies[obj].az)+"°"); } jday+=dstep; nextDate(obscopy,dstep, obs.day); } str += "
"; // outer block if (obj==100) { str += "

S=Sun, M=Moon, 1=Mercury, 2=Venus, 4=Mars, 5=Jupiter, 6=Saturn, 7=Uranus, 8=Neptune

"; } str += "

close window

\n"; doc.write(str); doc.close(); pwin.focus(); } // end doAltitude() function doAngles(obs,obj,dspan,dstep,type) { // Common function for declination (0), longitude (1) and elongation (2) var obscopy=new Object(); var obsmax=new Object(); for (var i in obs) {obscopy[i] = obs[i]; obsmax[i] = obs[i]; } var objects=(obj==100 ? [7,6,5,4,3,1,0,10,9] : [obj]); // objects to plot and order, later ones plot on top of earlier ones if (type==0) { var tstr="Declination"; var lbl = new Array("-30°","-25°","-20°","-15°","-10°","-5°","0°","5°","10°","15°","20°","25°","30°"); var lblpos = new Array(-10,40,90,140,190,244,297,347,395,445,495,545,595); var ylbl = "R.A"; } else if (type==1) { var tstr="Ecliptic Longitude"; var lbl = new Array("0°","30°","60°","90°","120°","150°","180°","210°","240°","270°","300°","330°","360°"); var lblpos = new Array(-10,40,97,144,194,244,294,344,394,444,494,544,594); var ylbl = "Lat."; } else { var tstr="Elongation"; var lbl = new Array("0°","15°","30°","45°","60°","75°","90°","105°","120°","135°","150°","165°","180°"); var lblpos = new Array(-3,45,95,145,195,245,295,342,392,442,492,542,592); var ylbl = "P.A."; } var ostr=(obj==100 ? "Sun, Moon and Planets" : bodies[obj].name); if (type==2 && obj==100) ostr = "Moon and Planets"; // Now make the diagram var pwin=window.open("","position","menubar,scrollbars,resizable"); var doc=pwin.document; var str = head1 + "Javascript: " + tstr + head3 + "

" + tstr + "

of "+ostr+"

\n"; str += "

at noon local time (UT " + hmstring(-obs.tz/60.0,true) + ")

\n"; nextDate(obsmax, dspan, obs.day); // 'abuse' nextdate to calculate end time jdmax = jd(obsmax); var jday=jd(obscopy); var tmax = Math.floor((jdmax-jday)/dstep); if (dstep<0) tmax = Math.floor(dspan/dstep); var leftoff=80; var topoff=30; str += "
"; str += box(leftoff,topoff,602,15*(tmax-1)+2,"#aaa url(grid.png)"); str += writexlbl(leftoff,topoff,lbl,lblpos,tstr,270); str += txt_lft(leftoff+612,topoff-30,ylbl); for (var t=0; t\n"; } if (objclose window

\n"; doc.write(str); doc.close(); pwin.focus(); } // end doAngles