90){write_db(places[ndx][NAME]+": invalid latitude ("+latitude+")");
errors++;
}
longitude=places[ndx][LONG];
if(Math.abs(lon_val=longitude.DMSto())>180){write_db(places[ndx][NAME]+": invalid longitude ("+longitude+")");
errors++;
}
timezone=places[ndx][TZ];
if(get_tz_ndx(timezone)<0){write_db(places[ndx][NAME]+": invalid time zone ("+timezone+")");
errors++;
}
else if(lon_val*timezone<0){write_db(places[ndx][NAME]+": warning: longitude-time zone mismatch ("
+longitude+" "+"<->"+" "+timezone+")");
warnings++;
}
dst_flag=places[ndx][DST];
if(dst_flag!==0&&dst_flag!=1){write_db(places[ndx][NAME]+": invalid DST flag ("+dst_flag+")");
errors++;
}
}
if(from_database===true)
form.loc_type[0].checked=true;
else
form.loc_type[1].checked=true;
form.place.selectedIndex=loc_ndx;
write_db(form.place.options.length+' locations');
if(warnings==1)
write_db(warnings+' location warning');
else if(warnings>1)
write_db(warnings+' location warnings');
if(errors==1)
write_db(errors+' location error');
else
write_db(errors+' location errors');
}
function checkDSTRules(year)
{var
errors=0,
loc,
ndx;
for(ndx=0; places[ndx]; ndx++){if(places[ndx][DST]==1){
loc=new Place(places[ndx]);
if(loc.checkDSTLimits(year)===false)
errors++;
}
}
if(errors==1)
write_db(errors+' DST rule error');
else
write_db(errors+' DST rule errors');
}
Place.prototype.checkDSTLimits=function(year)
{var errors=0;
if(this.uses_dst===0||flags.show_dst===false)
return true;
if(this.getDSTRule(year)===false){write_db('No DST rule for '+this.name);
return false;
}
if(this.checkDSTChangeDate(this.DSTStartStr, "start")===false)
errors++;
if(this.checkDSTChangeDate(this.DSTEndStr, "end")===false)
errors++;
if(this.getDSTLimits(year)===false){write_db(this.name+": cannot get DST start/end dates");
errors++;
}
if(this.dst_start!==null&&this.dst_end!==null){var dst_diff=this.dst_end-this.dst_start;
if(this.latitude/dst_diff<0){write_db(this.name+": DST start/end reversed?");
errors++;
}
}
if(errors>0)
return false;
else
return true;
};
Place.prototype.checkDSTChangeDate=function(changeStr, tag)
{var
errors=0,
change=changeStr.split(":");
var
change_month=(change[0]-0),
change_day=(change[1]-0),
change_week=(change[2]-0),
change_DoW=(change[3]-0),
change_time=(change[4]-0),
change_type=change[5];
if(change_month===0)
return true;
if(isNaN(change[0])||change_month<0||change_month>12){write_db(this.name+": invalid "+tag+" month ("+change[0]+")");
errors++;
}
if(isNaN(change[1])||change_day<1||
change_day>monthdays[1][change_month-1]){write_db(this.name+": invalid "+tag+" day ("+change[1]+" "
+months[change_month-1]+")");
errors++;
}
if(isNaN(change[2])||change_week<0||change_week>5){write_db(this.name+": invalid "+tag+" week ("+change[2]+")");
errors++;
}
if(isNaN(change[3])||change_DoW<0||change_DoW>6){write_db(this.name+": invalid "+tag+" DoW ("+change[3]+")");
errors++;
}
if(isNaN(change[4])||change_time<0||change_time>24){write_db(this.name+": invalid "+tag+" time ("+change[4]+")");
errors++;
}
if(change_type.search(/[suw]/)<0){write_db(this.name+": invalid "+tag+" time type ("+change_type+")");
errors++;
}
if(errors>0)
return false;
else
return true;
};
function searchKeyFilter(event)
{var
BACKSPACE=8,
RETURN=13,
e=window.event||event,
code=e.keyCode||e.charCode;
switch (e.type){case "keyup":
if(code==BACKSPACE)
search_places("backspace");
else if(code>=32)
search_places("new");
break;
case "keypress":
if(code==RETURN){search_places("next");
return false;
}
break;
default:
break;
}
return true;
}
// incoporates many ideas from http://www.geonames.org/export/ajax-postalcode-autocomplete.html
GeoNames={};
GeoNames.isProvince=function(str)
{var result;
result=str.match(/\b(AB|BC|MB|NB|NL|NS|NT|NU|NV|ON|PE|PQ|QC|SK|YT)\b/i);
if(result!==null)
return result[1];
else
return null;
};
GeoNames.provinces={"AB": "Alberta",
"BC": "British Columbia",
"MB": "Manitoba",
"NB": "New Brunswick",
"NL": "Newfoundland and Labrador",
"NS": "Nova Scotia",
"NT": "Northwest Territories",
"NU": "Nunavut",
"NV": "Nunavut",
"ON": "Ontario",
"PE": "Prince Edward Island",
"QC": "Quebec",
"PQ": "Quebec",
"SK": "Saskatchewan",
"YT": "Yukon"
};
GeoNames.closeSelectList=function()
{document.getElementById('selectListElement').innerHTML='';
document.getElementById('selectListElement').style.visibility='hidden';
};
GeoNames.selectListMouseOut=function(ndx)
{document.getElementById('locId'+ndx).className='selectListTR';
};
GeoNames.selectListMouseOver=function(ndx)
{document.getElementById('locId'+ndx).className='selectListMouseOver';
};
var geoData;
GeoNames.selectListMouseDown=function(ndx)
{GeoNames.closeSelectList();
GeoNames.setLocation(form, GeoNames.geoData[ndx]);
};
GeoNames.enableDisplayPrint=function(timer)
{clearTimeout(timer);
HTTP.Display.disabled=false;
HTTP.Print.disabled=false;
};
GeoNames.getLocation=function(data)
{var i, g;
if((GeoNames.geoData=data.geonames)===undefined){alert("Query failed");
return;
}
clearTimeout(HTTP.queryTimer);
HTTP.lookUpGoButton.disabled=false;
var
TH='| ',
TD=' | ',
TD_L=' | ',
TD_C=' | ';
TD_R=' | ';
var selectList;
var selectListHTML='';
switch (GeoNames.geoData.length){case 0:
alert("No matches");
break;
case 1:
GeoNames.setLocation(form, GeoNames.geoData[0]);
break;
default:
selectListHTML+="";
selectListHTML+=TR+
TH+"Name"+THE+
TH+"Country"+THE+
TH+"State/Province"+THE+
TH+"County/Other"+THE+
TH+"Feat. Type"+THE+
TH+"Lat, °"+THE+
TH+"Long, °"+THE;
if(flags.elevation_in_feet===true)
selectListHTML+=TH+"Elev, ft"+THE;
else
selectListHTML+=TH+"Elev, m"+THE;
selectListHTML+=TRE;
for(i=0; i"+
TD_L+g.name+TDE+
TD+g.countryName+TDE+
TD_L+g.adminName1+TDE+
TD_L+g.adminName2+TDE+
TD_L+g.fcodeName+TDE+
TD_R+g.lat.toDMS(0).replace(/^-/, MINUS)+TDE+
TD_R+g.lng.toDMS(0).replace(/^-/, MINUS)+TDE;
if(g.elevation===undefined)
selectListHTML+=TD_C+SP+TDE;
else if(flags.elevation_in_feet===true)
selectListHTML+=TD_R+g.elevation.MtoFT().toFixed(0)+TDE;
else
selectListHTML+=TD_R+g.elevation+TDE;
selectListHTML+=TRE;
}
selectListHTML+=" ";
selectList=document.getElementById('selectListElement');
selectList.innerHTML=selectListHTML;
selectList.style.visibility='visible';
break;
}
};
GeoNames.setLocation=function(form, gdata)
{var admin_name;
var latitude, longitude, timezone;
var query, type;
if(gdata.elevation===undefined){query="lat="+gdata.lat+"&lng="+gdata.lng;
if(gdata.lat>=-56&&gdata.lat<=60)
type="srtm3";
else if(gdata.lat>=-65&&gdata.lat<=83)
type="astergdem";
else
type="gtopo30";
GeoNames.query(type, query);
}
else if(flags.elevation_in_feet===true)
form.elevation.value=addCommas(gdata.elevation.MtoFT().toFixed(0));
else
form.elevation.value=addCommas(gdata.elevation);
if(gdata.timezone.gmtOffset===undefined||
gdata.timezone.dstOffset===undefined){
query="lat="+gdata.lat+"&lng="+gdata.lng;
GeoNames.query("timezone", query);
}
else{
if(gdata.lat<0)
timezone=gdata.timezone.dstOffset;
else
timezone=gdata.timezone.gmtOffset;
if((tz_ndx=get_tz_ndx(timezone))<0){tz_ndx=-1;
alert("Invalid time zone: "+timezone);
}
form.timezone.selectedIndex=tz_ndx+1;
if(gdata.timezone.dstOffset!=gdata.timezone.gmtOffset)
form.uses_dst.selectedIndex=1;
else
form.uses_dst.selectedIndex=0;
}
switch (gdata.countryCode){case "US":
admin_name=gdata.adminCode1;
break;
case "CA":
for(var code in GeoNames.provinces){if(GeoNames.provinces[code]==gdata.adminName1){
admin_name=code;
break;
}
}
if(admin_name===undefined)
admin_name=gdata.adminName1;
break;
default:
admin_name=gdata.countryName;
break;
}
form.place_name.value=gdata.name+", "+admin_name;
if((latitude=gdata.lat)<0){form.latitude.value=(-latitude).toDMS(-1);
form.lat_sign.selectedIndex=1;
}
else{form.latitude.value=latitude.toDMS(-1);
form.lat_sign.selectedIndex=0;
}
if((longitude=gdata.lng)<0){form.longitude.value=(-longitude).toDMS(-1);
form.long_sign.selectedIndex=1;
}
else{form.longitude.value=longitude.toDMS(-1);
form.long_sign.selectedIndex=0;
}
SaveState.write(form);
return true;
};
GeoNames.setElevation=function(data, obj)
{if(data[obj]===undefined){
alert("Could not get elevation");
GeoNames.enableDisplayPrint(HTTP.elevationTimer);
return false;
}
else if(obj=="srtm3"&&data[obj]==-32768)
form.elevation.value=0;
else if(obj=="astergdem"&&data[obj]==-9999)
form.elevation.value=0;
else if(obj=="gtopo30"&&data[obj]==-9999)
form.elevation.value=0;
else if(flags.elevation_in_feet===true)
form.elevation.value=addCommas(data[obj].MtoFT().toFixed(0));
else
form.elevation.value=addCommas(data[obj]);
GeoNames.enableDisplayPrint(HTTP.elevationTimer);
SaveState.write(form);
return true;
};
GeoNames.setTimeZone=function(data)
{var timezone, tz_ndx;
if(data.gmtOffset===undefined||
data.dstOffset===undefined){
alert("Could not get time zone");
GeoNames.enableDisplayPrint(HTTP.timezoneTimer);
return false;
}
else{
if(data.lat<0)
timezone=data.dstOffset;
else
timezone=data.gmtOffset;
}
if((tz_ndx=get_tz_ndx(timezone))<0){tz_ndx=-1;
alert("Invalid time zone: "+timezone);
}
form.timezone.selectedIndex=tz_ndx+1;
if(data.dstOffset!=data.gmtOffset)
form.uses_dst.selectedIndex=1;
else
form.uses_dst.selectedIndex=0;
GeoNames.enableDisplayPrint(HTTP.timezoneTimer);
SaveState.write(form);
return true;
};
GeoNames.query=function(type, query)
{var url;
var callback, obj;
var name, name_param, rest, province;
var timer=null;
var buttons=[];
switch (type){case "search":
url = "http://ws.geonames.org/searchJSON?";
url+="style=FULL";
callback=GeoNames.getLocation;
obj=null;
if(query.indexOf("~")>=0){name_param="&name=";
query=query.replace(/\~/g, "");
}
else
name_param="&name_equals=";
if(query.indexOf(",") >= 0) { name=query.replace(/,.*/, "");
rest=query.replace(/.*, */, "");
if((province=GeoNames.isProvince(rest))!==null)
rest=rest.replace(province, GeoNames.provinces[province.toUpperCase()]);
query=name_param+encodeURIComponent(name)+
"&q="+encodeURIComponent(rest);
}
else
query=name_param+encodeURIComponent(query);
buttons=[ "lookUpGoButton" ];
timer="queryTimer";
break;
case "srtm3":
url = "http://ws.geonames.org/srtm3JSON?";
callback=GeoNames.setElevation;
obj="srtm3";
buttons=[ "Display", "Print" ];
timer="elevationTimer";
break;
case "astergdem":
url = "http://ws.geonames.org/astergdemJSON?";
callback=GeoNames.setElevation;
obj="astergdem";
buttons=[ "Display", "Print" ];
timer="elevationTimer";
break;
case "gtopo30":
url = "http://ws.geonames.org/gtopo30JSON?";
callback=GeoNames.setElevation;
obj="gtopo30";
buttons=[ "Display", "Print" ];
timer="elevationTimer";
break;
case "timezone":
url = "http://ws.geonames.org/timezoneJSON?";
callback=GeoNames.setTimeZone;
obj=null;
buttons=[ "Display", "Print" ];
timer="timezoneTimer";
break;
default:
break;
}
HTTP.getJSON(url, query, callback, obj, timer, buttons);
};
var HTTP={};
HTTP.initialize_buttons=function()
{
HTTP.lookUpGoButton=document.getElementById("cstmLocLookupGo");
HTTP.Display=document.getElementsByName("Display")[0];
HTTP.Print=document.getElementsByName("Print")[0];
};
HTTP.getJSON=function(url, query, callback, obj, timer, buttons)
{var i;
var timeout=16000;
var script=document.createElement("script");
document.body.appendChild(script);
var funcname="func"+HTTP.getJSON.counter++;
HTTP.getJSON[funcname]=function(text)
{
callback(text, obj);
document.body.removeChild(script);
delete HTTP.getJSON[funcname];
};
HTTP[timer]=setTimeout(function()
{
document.body.removeChild(script);
delete HTTP.getJSON[funcname];
alert("Request timed out");
for(i=0; i=0){if(calc_date>this.dst_start_date&&calc_datethis.dst_end_date&&calc_date0){this.snorm[n][m]=this.snorm[n][m-1]*Math.sqrt(((n-m+1)*j)/(n+m));
j=1;
this.c_main[m-1][n]=this.snorm[n][m]*this.c_main[m-1][n];
this.c_sec[m-1][n]=this.snorm[n][m]*this.c_sec[m-1][n];
}
this.c_main[n][m]=this.snorm[n][m]*this.c_main[n][m];
this.c_sec[n][m]=this.snorm[n][m]*this.c_sec[n][m];
}
this.fn[n]=n+1;
this.fm[n]=n;
}
for(n=0; n<=this.maxord; n++){for(m=0; m<=n; m++)
this.p[n][m]=this.snorm[n][m];
}
this.k[1][1]=0;
this.start_date=Date.UTC(this.epoch, 0, 0).MsecToDay();
this.end_date=Date.UTC(this.epoch+this.t_valid, 0, 0).MsecToDay();
return true;
};
MagModel.prototype.calcMagDec=function (loc, date)
{
var m, n;
var
dt,
aor,
ar,
br,
bt,
bp,
bpp,
par,
parp,
c2,
c4,
dyear,
rlon,
rlat,
sin_rlon,
sin_rlat,
cos_rlon,
cos_rlat,
sin_rlon2,
sin_rlat2,
cos_rlon2,
cos_rlat2,
temp1,
temp2,
t_long;
dyear=getDecimalYear(date);
var
glat=loc.latitude,
glon=loc.longitude,
elev=0;
dt=dyear-this.epoch;
rlon=glon.DtoR();
rlat=glat.DtoR();
sin_rlon=Math.sin(rlon);
sin_rlat=Math.sin(rlat);
cos_rlon=Math.cos(rlon);
cos_rlat=Math.cos(rlat);
sin_rlon2=Math.square(sin_rlon);
sin_rlat2=Math.square(sin_rlat);
cos_rlon2=Math.square(cos_rlon);
cos_rlat2=Math.square(cos_rlat);
this.sin_p[1]=sin_rlon;
this.cos_p[1]=cos_rlon;
c2=this.a2-this.b2;
c4=this.a4-this.b4;
if(elev!=this.old_elev||glat!=this.old_lat){this.q=Math.sqrt(this.a2-c2*sin_rlat2);
this.q1=elev*this.q;
this.q2=Math.square((this.q1+this.a2)/(this.q1+this.b2));
this.cos_t=sin_rlat/Math.sqrt(this.q2*cos_rlat2+sin_rlat2);
this.sin_t=Math.sqrt(1.0-Math.square(this.cos_t));
this.r2=Math.square(elev)+2.0*this.q1+(this.a4-c4*sin_rlat2)/Math.square(this.q);
this.r=Math.sqrt(this.r2);
this.d=Math.sqrt(this.a2*cos_rlat2+this.b2*sin_rlat2);
this.cos_a=(elev+this.d)/this.r;
this.sin_a=c2*cos_rlat*sin_rlat/(this.r*this.d);
}
if(this.reset===true||glon!=this.old_lon){for(m=2; m<=this.maxord; m++){
this.sin_p[m]=this.sin_p[1]*this.cos_p[m-1]+this.cos_p[1]*this.sin_p[m-1];
this.cos_p[m]=this.cos_p[1]*this.cos_p[m-1]-this.sin_p[1]*this.sin_p[m-1];
}
}
aor=MagModel.MEAN_RAD/this.r;
ar=Math.square(aor);
br=0;
bt=0;
bp=0;
bpp=0;
for(n=1; n<=this.maxord; n++){ar*=aor;
for(m=0; m<=n; m++){
if(this.reset===true||elev!=this.old_elev||glat!=this.old_lat){if(n==m){
this.p[n][m]=this.sin_t*this.p[n-1][m-1];
this.dp[n][m]=this.sin_t*this.dp[n-1][m-1]+this.cos_t*this.p[n-1][m-1];
}
else if(n==1&&m===0){this.p[n][m]=this.cos_t*this.p[n-1][m];
this.dp[n][m]=this.cos_t*this.dp[n-1][m]-this.sin_t*this.p[n-1][m];
}
else if(n>1&&n!=m){if(m>n-2){
this.p[n-2][m]=0.0;
this.dp[n-2][m]=0.0;
}
this.p[n][m]=this.cos_t*this.p[n-1][m]-this.k[n][m]*this.p[n-2][m];
this.dp[n][m]=this.cos_t*this.dp[n-1][m]-this.sin_t*this.p[n-1][m]
-this.k[n][m]*this.dp[n-2][m];
}
}
if(this.reset===true||dyear!=this.old_dyear){this.c_time[n][m]=this.c_main[n][m]+dt*this.c_sec[n][m];
if(m!==0)
this.c_time[m-1][n]=this.c_main[m-1][n]+dt*this.c_sec[m-1][n];
}
par=ar*this.p[n][m];
if(m===0){temp1=this.c_time[n][m]*this.cos_p[m];
temp2=this.c_time[n][m]*this.sin_p[m];
}
else{temp1=this.c_time[n][m]*this.cos_p[m]+this.c_time[m-1][n]*this.sin_p[m];
temp2=this.c_time[n][m]*this.sin_p[m]-this.c_time[m-1][n]*this.cos_p[m];
}
bt-=ar*temp1*this.dp[n][m];
bp+=this.fm[m]*temp2*par;
br+=this.fn[n]*temp1*par;
if(Math.abs(this.sin_t)=30){loc.mag_dec=mag_model.calcMagDec(loc, calc_date);
mag_dec_date=calc_date;
}
if(sun.az_rise!==null)
sun.az_rise=(sun.az_rise-loc.mag_dec).setRange(360);
if(sun.az_set!==null)
sun.az_set=(sun.az_set-loc.mag_dec).setRange(360);
if(moon.az_rise!==null)
moon.az_rise=(moon.az_rise-loc.mag_dec).setRange(360);
if(moon.az_set!==null)
moon.az_set=(moon.az_set-loc.mag_dec).setRange(360);
}
if(rs_opts.select(sun, moon, loc)===true||calc_params.calc_type==POSITIONS){found_one=true;
this.times[ndx++]=format_times(CalcDate, calc_params.calc_type,
sun, moon, rs_opts, t_offset, loc);
}
}
rs_opts.matches=ndx;
return found_one;
}
function calculate_positions(calc_params, sun, moon, loc)
{var
altitude,
calc_time,
ndx=0,
path_angle,
places,
spos=new SunPos(),
mpos=new MoonPos(),
CalcDate=new Date(),
jday,
start_date=calc_params.start_date;
var
ncalcs,
start_time=calc_params.start_time,
end_time=calc_params.end_time,
time_step=calc_params.time_step;
end_time+=1e-4;
var line;
CalcDate.setTime(start_date.DayToMsec());
jday=CalcDate.getJulian();
calcObliquity(jday, spos, mpos);
calcNutation(jday, spos, mpos);
calcDeltaT(jday, spos, mpos);
this.positions.length=0;
ncalcs=0;
places=2;
for(calc_time=start_time; calc_time<=end_time; calc_time+=time_step){if(++ncalcs>2880)
break;
calc_time=Math.round(calc_time*60)/60;
spos.calcPos(jday, calc_time);
spos.geoToTopo(loc);
spos.calcAltAz(loc);
spos.addRefraction(loc);
spos.calcDecl(loc);
mpos.calcPos(jday, calc_time);
mpos.calcPhase(spos.dist, loc);
mpos.geoToTopo(loc);
mpos.calcAltAz(loc);
mpos.addRefraction(loc);
mpos.calcDecl(loc);
if(flags.show_mag_north===true){spos.azimuth=(spos.azimuth-loc.mag_dec).setRange(360);
mpos.azimuth=(mpos.azimuth-loc.mag_dec).setRange(360);
}
CalcDate.setTime((start_date+calc_time/24).DayToMsec());
line=TR+LF;
line+=TD;
if(time_step<0.5&&Math.rnd(calc_time+loc.ut_offset, 3) % 1===0)
line+=calc_time.showPosTime(CalcDate, loc).bold();
else
line+=calc_time.showPosTime(CalcDate, loc);
line+=TDE+LF;
line+=TD+spos.azimuth.toFixed(places)+SP2+TDE+LF;
if((altitude=spos.altitude)<0)
line+=TD+MINUS+(-altitude).toFixed(places);
else
line+=TD+altitude.toFixed(places);
line+=SP2+TDE+LF;
if((path_angle=spos.getPathAngle(loc))>90)
line+=TD+SP+(180-path_angle).toFixed(1)+MINUS+SP+TDE+LF;
else
line+=TD+SP+path_angle.toFixed(1)+SP2+TDE+LF;
if(altitude>0){var shadow_len=(1/Math.tan(altitude.DtoR() ));
if(shadow_len>100)
line+=TD_C+">"+" 100"+TDE+LF;
else
line+=TD+shadow_len.toFixed(2)+SP+TDE+LF;
}
else
line+=TD_C+MDASH+TDE+LF;
line+=TD+mpos.azimuth.toFixed(places)+SP2+TDE+LF;
if((altitude=mpos.altitude)<0)
line+=TD+MINUS+(-altitude).toFixed(places);
else
line+=TD+altitude.toFixed(places);
line+=SP2+TDE+LF;
if((path_angle=mpos.getPathAngle(loc))>90)
line+=TD+SP+(180-path_angle).toFixed(1)+MINUS+SP+TDE+LF;
else
line+=TD+SP+path_angle.toFixed(1)+SP2+TDE+LF;
line+=TD+mpos.phase.toFixed(3)+SP+TDE+LF;
line+=TD+SP+mpos.sd.toFixed(3)+SP+TDE+LF;
line+=TRE+LF;
this.positions[ndx++]=line;
}
return true;
}
Place.prototype.calcObserverCoord=function()
{var
earth_rad=EARTH_RADIUS*1e3,
elevation=this.elevation,
latitude=this.latitude.DtoR(),
ratio=0.99664719,
u;
u=Math.atan(ratio*Math.tan(latitude));
this.rhoSinPhi=ratio*Math.sin(u)+elevation/earth_rad*Math.sin(latitude);
this.rhoCosPhi=Math.cos(u)+elevation/earth_rad*Math.cos(latitude);
this.rho=Math.sqrt(Math.square(this.rhoSinPhi)+Math.square(this.rhoCosPhi));
};
Place.prototype.calcDensities=function()
{var
elevation=this.elevation,
height=this.height,
lapse_rate=6.5e-3,
ndx=4.256,
std_temp=273.15+10;
this.density=calcDensity(elevation);
this.density_gnd=calcDensity(elevation-height);
};
calcDensity=function(elevation)
{var
lapse_rate=6.5e-3,
ndx=4.256,
temp0=273.15+10;
return Math.pow(1-lapse_rate/temp0*elevation, ndx);
};
function calc_k(elevation)
{var
lapse_rate=6.5e-3,
ndx=3.256,
temp0=273.15+10;
return K_REFR*Math.pow(1.0-lapse_rate/temp0*elevation, ndx);
}
Place.prototype.calcRSAlt=function()
{var
refr_0h_gnd,
horizon_true;
refr_0h_gnd=this.refr_0h_elev=refraction_app(0, this);
this.refr_0h=this.refr_0h_elev/this.density;
if(this.height>0)
refr_0h_gnd*=this.density_gnd/this.density;
this.rs_alt=-refr_0h_gnd;
this.dip=0;
if(this.height>0){var
density,
k,
r_E;
density=this.density;
k=calc_k(this.elevation-2*this.height/3);
r_E=EARTH_RADIUS*1e3*this.rho;
this.dip=Math.sqrt(2*(1-k)*this.height/r_E).RtoD();
horizon_true=this.dip/(1-k);
this.rs_alt-=horizon_true;
}
};
function calcObliquity(jday, sun, moon)
{var
T, T2, T3,
obliquity;
T=(jday-J2000)/36525;
T2=T*T;
T3=T2*T;
obliquity=-46.8150*T
-0.00059*T2
+0.001813*T3;
obliquity/=3600.0;
obliquity+=23.43929111;
obliquity=obliquity.DtoR();
sun.obliquity=moon.obliquity=obliquity;
}
function calcNutation(jday, sun, moon)
{var
T,
T2,
L_s,
G_s,
L_m,
N_m,
nut_lon,
nut_obl,
eq_equinoxes;
T=(jday-J2000)/36525;
T2=T*T;
L_s=0.77907196+100.00213759*T+0.00000084*T2;
G_s=0.99312619+99.99735956*T-0.00000044*T2;
L_m=0.60643382+1336.85522267*T-0.00000313*T2;
N_m=0.34734264-5.37260713*T+0.00000579*T2;
L_s=L_s.setRange(1)*Math.PI2;
G_s=G_s.setRange(1)*Math.PI2;
L_m=L_m.setRange(1)*Math.PI2;
N_m=N_m.setRange(1)*Math.PI2;
nut_lon=-17.200*Math.sin(N_m)
-1.319*Math.sin(2*L_s)
-0.227*Math.sin(2*L_m)
+0.206*Math.sin(2*N_m)
+0.143*Math.sin(G_s);
nut_obl=9.203*Math.cos(N_m)
+0.574*Math.cos(2*L_s)
+0.098*Math.cos(2*L_m)
-0.090*Math.cos(2*N_m);
nut_lon=(nut_lon/3600).DtoR();
nut_obl=(nut_obl/3600).DtoR();
eq_equinoxes=nut_lon*Math.cos(sun.obliquity+nut_obl);
sun.nut_lon=moon.nut_lon=nut_lon;
sun.nut_obl=moon.nut_obl=nut_obl;
sun.eq_equinoxes=moon.eq_equinoxes=eq_equinoxes;
}
function calcDeltaT(jday, sun, moon)
{
var
deltaT,
T,
t,
u,
year;
T=(jday-J2000)/36525;
year=2000+100*T;
if(year>=2005){if(year<2050){
t=year-2000;
deltaT=62.92+0.32217*t+0.005589*Math.square(t);
}
else if(year<2150){t=(year-1820)/100;
deltaT=-20+32*Math.square(t)-0.5628*(2150-year);
}
else{u=(year-1820)/100;
deltaT=-20+32*Math.square(u);
}
}
else if(year>=1986){t=year-2000;
deltaT=63.86+0.3345*t-0.060374*Math.square(t)
+0.0017275*Math.pow(t, 3)+0.000651814*Math.pow(t, 4)
+0.00002373599*Math.pow(t, 5);
}
else if(year>=1961){t=year-1975;
deltaT=45.45+1.067*t-Math.square(t)/260-Math.pow(t, 3)/718;
}
else if(year>=1941){t=year-1950;
deltaT=29.07+0.407*t-Math.square(t)/233+Math.pow(t, 3)/2547;
}
else if(year>=1920){t=year-1920;
deltaT=21.20+0.84493*t-0.076100*Math.square(t)
+0.0020936*Math.pow(t, 3);
}
else if(year>=1900){t=year-1900;
deltaT=-2.79+1.494119*t-0.0598939*Math.square(t)
+0.0061966*Math.pow(t, 3)-0.000197*Math.pow(t, 4);
}
else if(year>=1860){t=year-1860;
deltaT=7.62+0.5737*t-0.251754*Math.square(t)
+0.01680668*Math.pow(t, 3)-0.0004473624*Math.pow(t, 4)
+Math.pow(t, 5)/233174;
}
else if(year>=1800){t=year-1800;
deltaT=13.72-0.332447*t+0.0068612*Math.square(t)
+0.0041116*Math.pow(t, 3)-0.00037436*Math.pow(t, 4)
+0.0000121272*Math.pow(t, 5)-0.0000001699*Math.pow(t, 6)
+0.000000000875*Math.pow(t, 7);
}
else if(year>=1700){t=year-1700;
deltaT=8.83+0.1603*t-0.0059285*Math.square(t)
+0.00013336*Math.pow(t, 3)-Math.pow(t, 4)/1174000;
}
else if(year>=1600){t=year-1600;
deltaT=120-0.9808*t-0.01532*Math.square(t)+Math.pow(t, 3)/7129;
}
else{u=(year-1000)/100;
deltaT=1574.2-556.01*u+71.23472*Math.square(u)
+0.319781*Math.pow(u, 3)-0.8503463*Math.pow(u, 4)
-0.005050998*Math.pow(u, 5)+0.0083572073*Math.pow(u, 6);
}
sun.deltaT=moon.deltaT=deltaT/86400;
}
function GSTofUT(jday, ut)
{var
T,
T2, T3,
gst0,
gst;
T=(jday-J2000)/36525.0;
T2=T*T;
T3=T2*T;
gst0=6.697374558+2400.05133691*T+0.000025862*T2-1.72e-9*T3;
return (gst0+ut.UTtoST()).setRange(24);
}
function refraction_app(alt, loc)
{var refraction;
if(alt<-4.3425||alt>89.9225)
return 0;
refraction=1/Math.tan((alt+7.31/(alt+4.4)).DtoR());
refraction*=loc.density;
return refraction/60;
}
function refraction_true(alt, loc)
{var refr0, refr;
if(alt<-5.0016||alt>89.891)
refr=0;
else if(Math.abs(loc.elevation)>1e-6&&alt<=5){refr0=-1;
refr=0;
while(Math.abs(refr-refr0)>0.01){refr0=refr;
if(alt<-4.3425||alt>89.9225)
return 0;
refr=1/Math.tan((alt+refr+7.31/(alt+refr+4.4)).DtoR());
refr*=loc.density;
refr/=60;
}
}
else{refr=1.02/Math.tan((alt+10.3/(alt+5.11)).DtoR());
refr*=loc.density;
refr/=60;
}
return refr;
}
function parallax_geo(alt, sin_hp, loc)
{var num, denom;
num=loc.rho*sin_hp*Math.cos(alt.DtoR());
denom=1.0-loc.rho*sin_hp*Math.sin(alt.DtoR());
return Math.atan2(num, denom).RtoD();
}
function parallax_topo(alt, sin_hp, loc)
{return Math.asin(loc.rho*sin_hp*Math.cos(alt.DtoR())).RtoD();
}
function getTopoSD()
{return this.sd0*(1+Math.sin(this.altitude.DtoR())*this.sin_hp);
}
function getGeoAlt(alt, loc)
{var
rs_alt_app,
sin_hp=this.sin_hp;
alt=appAltToTrueAlt(alt, loc);
if(sin_hp>0)
alt+=parallax_topo(alt, sin_hp, loc);
return alt;
}
function getTopoAlt(loc)
{var
alt=this.altitude,
sin_hp=this.sin_hp;
var alt0=alt;
if(sin_hp>0)
alt-=parallax_geo(alt, sin_hp, loc);
return alt;
}
function getAppAlt(loc)
{var
alt=this.altitude,
refraction,
sin_hp=this.sin_hp;
if(sin_hp>0)
alt-=parallax_geo(alt, sin_hp, loc);
refraction=refraction_true(alt, loc);
if(alt+refraction>=0)
alt+=refraction;
return alt;
}
function getRefrAlt(loc)
{var
alt=this.altitude,
refraction;
refraction=refraction_true(alt, loc);
if(alt+refraction>=0)
alt+=refraction;
return alt;
}
addRefraction=function(loc)
{this.altitude=trueAltToAppAlt(this.altitude, loc);
};
function trueAltToAppAlt(altitude, loc)
{if(loc.height>0){
if(altitude>=loc.rs_alt&&altitude<=-loc.refr_0h_elev){var
alt,
density_vh,
elev_vh,
elevation,
k,
r_E,
y_diff;
elevation=loc.elevation;
r_E=(EARTH_RADIUS*1e3*loc.rho);
alt=altitude+loc.refr_0h*loc.density_gnd;
var alt0=-99;
elev_vh=loc.elevation-loc.height;
y_diff=loc.height;
while(Math.abs(alt-alt0)>0.01){alt0=alt;
k=calc_k(elev_vh+y_diff/3);
if(alt<0)
y_diff=(1-k)*r_E/2*Math.square(alt.DtoR());
else
y_diff=0;
elev_vh=elevation-y_diff;
density_vh=calcDensity(elev_vh);
alt=altitude+loc.refr_0h*density_vh;
}
altitude=alt;
if(altitude<0){
k=calc_k(elev_vh+y_diff/3);
altitude*=(1-k);
}
}
else if(altitude>-loc.refr_0h_elev)
altitude+=refraction_true(altitude, loc);
}
else{var refraction=refraction_true(altitude, loc);
if(altitude+refraction>=0)
altitude+=refraction;
}
return altitude;
}
function appAltToTrueAlt(alt_app, loc)
{var alt_true;
if(Math.abs(alt_app)<1e-9)
alt_true=loc.rs_alt;
else if(loc.height>0&&alt_app>=-loc.dip&&alt_app<0){
var
density_vh,
elev_vh,
elevation,
k,
r_E,
y_diff;
elevation=loc.elevation;
elev_vh=loc.elevation-loc.height;
r_E=EARTH_RADIUS*1e3*loc.rho;
y_diff=loc.height;
k=calc_k(elev_vh+y_diff/3);
y_diff=Math.square(alt_app.DtoR())*r_E/(2*(1-k));
alt_true=alt_app/(1-k);
elev_vh=elevation-y_diff;
density_vh=calcDensity(elev_vh);
alt_true-=loc.refr_0h*density_vh;
}
else if(alt_app>=0)
alt_true=alt_app-refraction_app(alt_app, loc);
else
alt_true=alt_app;
return alt_true;
}
function getPathAngle(loc)
{var
alt=this.altitude.DtoR(),
argument,
declination=this.declination;
argument=(loc.sin_lat-Math.sin(declination)*Math.sin(alt))
/(Math.cos(declination)*Math.cos(alt));
if(argument>1.0)
return 0;
else if(argument<-1)
return 180;
else
return Math.acos(argument).RtoD();
}
function calcDeclination(altitude, azimuth, loc)
{var alt, az, arg;
alt=altitude.DtoR();
az=azimuth.DtoR();
arg=Math.sin(alt)*loc.sin_lat+Math.cos(alt)*loc.cos_lat*Math.cos(az);
return Math.asin(arg);
}
function calcDecl(loc)
{this.declination=calcDeclination(this.altitude, this.azimuth, loc);
}
function eclipticToTopo(loc)
{var
num,
denom,
cos_obl,
sin_obl,
rhoCosPhi=loc.rhoCosPhi,
rhoSinPhi=loc.rhoSinPhi,
lambda_g=this.lambda,
beta_g=this.beta,
lambda_t,
beta_t,
cos_beta,
lst,
cos_lst,
sin_lst,
sin_hp=this.sin_hp;
lst=(this.gst+loc.longitude/15).setRange(24);
cos_obl=Math.cos(this.obliquity+this.nut_obl);
sin_obl=Math.sin(this.obliquity+this.nut_obl);
cos_beta=Math.cos(beta_g);
cos_lst=Math.cos(lst.HtoR());
sin_lst=Math.sin(lst.HtoR());
denom=Math.cos(lambda_g)*cos_beta-rhoCosPhi*sin_hp*cos_lst;
num=Math.sin(lambda_g)*cos_beta
-sin_hp*(rhoSinPhi*sin_obl+rhoCosPhi*cos_obl*sin_lst);
lambda_t=Math.atan2(num, denom);
if(lambda_t<0.0)
lambda_t+=Math.PI2;
num=Math.cos(lambda_t)*(Math.sin(beta_g)
-sin_hp*(rhoSinPhi*cos_obl-rhoCosPhi*sin_obl*sin_lst));
beta_t=Math.atan(num/denom);
this.lambda_t=lambda_t;
this.beta_t=beta_t;
}
function eclipticToEquat(obliquity)
{var
sin_long,
cos_obl,
sin_obl;
sin_long=Math.sin(this.lambda);
sin_obl=Math.sin(obliquity);
cos_obl=Math.cos(obliquity);
this.ra=Math.atan2(sin_long*cos_obl
-Math.tan(this.beta)*sin_obl, Math.cos(this.lambda));
if(this.ra<0.0)
this.ra+=Math.PI2;
this.declination=Math.asin(Math.sin(this.beta)*cos_obl
+Math.cos(this.beta)*sin_obl*sin_long);
}
function geoToTopo(loc)
{var
lha=(this.gha+loc.longitude).DtoR(),
cos_dec=Math.cos(this.declination),
a=cos_dec*Math.sin(lha),
b=cos_dec*Math.cos(lha)-loc.rhoCosPhi*this.sin_hp,
c=Math.sin(this.declination)-loc.rhoSinPhi*this.sin_hp,
q=Math.sqrt(Math.square(a)+Math.square(b)+Math.square(c));
this.lha=Math.atan2(a, b);
this.ra=(this.gst.HtoR()+loc.longitude.DtoR())-this.lha;
this.declination=Math.asin(c/q);
this.sd=this.sd0/q;
}
function calcAltAz(loc)
{var
arg, num, denom,
lha=this.gst.HtoR()+loc.longitude.DtoR()-this.ra,
decl=this.declination,
sin_lat=loc.sin_lat,
cos_lat=loc.cos_lat;
num=Math.sin(lha);
denom=Math.cos(lha)*sin_lat-Math.tan(decl)*cos_lat;
this.azimuth=Math.atan2(num, denom).RtoD()+180;
arg=sin_lat*Math.sin(decl)+cos_lat*Math.cos(decl)*Math.cos(lha);
this.altitude=Math.asin(arg).RtoD();
this.zenith=90-this.altitude;
}
function getAzRS(loc, type)
{var
alt_rs=this.altitude.DtoR(),
azimuth,
cos_az,
declination=this.declination;
cos_az=
(Math.sin(declination)-loc.sin_lat*Math.sin(alt_rs))
/(loc.cos_lat*Math.cos(alt_rs));
if(cos_az>1.0)
azimuth=null;
else if(cos_az<-1.0)
azimuth=null;
else
azimuth=Math.acos(cos_az).RtoD();
if(azimuth===null)
return null;
else
return (type=="R" ? azimuth : 360-azimuth).setRange(360);
}
function getMaxAlt(loc)
{var
altitude;
altitude=90.0-Math.abs(this.declination.RtoD()-loc.latitude);
this.altitude=altitude;
return this.getAppAlt(loc);
}
function getEventTime(jday, loc, type, alt_rs, alt_rs_type)
{var
cos_H,
longitude=loc.longitude,
ntries,
sign,
sin_h,
lha0,
lha,
sd,
ut,
ut0=-1;
var
sin_lat=loc.sin_lat,
cos_lat=loc.cos_lat;
switch (type){case "M":
case "R":
sign=1;
break;
case "T":
sign=0;
break;
case "S":
case "E":
sign=-1;
break;
default:
break;
}
switch (type){case "M":
case "E":
sin_h=Math.sin(alt_rs.DtoR());
break;
case "T":
break;
default:
if(this.body=="Sun"){
this.altitude=this.getGeoAlt(alt_rs, loc);
if(alt_rs_type=="top")
this.altitude-=this.getTopoSD();
sin_h=Math.sin(this.altitude.DtoR());
}
break;
}
if(loc.timezone>0&&longitude<0)
longitude+=360;
else if(loc.timezone<0&&longitude>0)
longitude-=360;
ut=12-loc.longitude/15;
lha=0;
ut0=ut-12;
ntries=0;
while(Math.abs(ut-ut0)>0.005){if(++ntries>10)
break;
ut0=ut;
if(flags.interpolate_rs_positions===true)
this.interpPos(jday, ut, loc);
else
this.calcPos(jday, ut);
lha0=this.gha+longitude;
if(this.body=="Moon"){
this.altitude=this.getGeoAlt(alt_rs, loc);
if(alt_rs_type=="top")
this.altitude-=this.getTopoSD();
sin_h=Math.sin(this.altitude.DtoR());
}
if(type=="T")
cos_H=1;
else{cos_H=(sin_h-sin_lat*Math.sin(this.declination))/
(cos_lat*Math.cos(this.declination));
}
if(cos_H>1)
lha=0;
else if(cos_H<-1)
lha=180;
else
lha=Math.acos(cos_H).RtoD();
ut=ut0-(lha0+sign*lha)/15;
while(utloc.ut_max)
ut-=24;
}
if(cos_H>1)
ut=Number.NEGATIVE_INFINITY;
else if(cos_H<-1)
ut=Number.POSITIVE_INFINITY;
else if(ntries>10)
ut=F_SET_NO_EVENT;
return (ut);
}
function Array2D2(rows, cols)
{var a=new Array(rows);
for(var i=0; i0){if(this.pos[i][0]<=this.pos[i-1][0])
this.pos[i][0]+=Math.PI2;
if(this.pos[i][2]<=this.pos[i-1][2])
this.pos[i][2]+=Math.PI2;
if(this.pos[i][4]<=this.pos[i-1][4])
this.pos[i][4]+=Math.PI2;
}
}
}
function interpPos(jday, ut, loc)
{var n;
n=(ut+loc.ut_offset_start)/24;
if(this.body=="Moon"){
this.lambda=interpolate(this.pos[0][0], this.pos[1][0], this.pos[2][0], n);
this.beta=interpolate(this.pos[0][1], this.pos[1][1], this.pos[2][1], n);
this.lambda_s=interpolate(this.pos[0][2], this.pos[1][2], this.pos[2][2], n);
this.sin_hp=interpolate(this.pos[0][3], this.pos[1][3], this.pos[2][3], n);
this.sd0=interpolate(this.pos[0][6], this.pos[1][6], this.pos[2][6], n);
}
else{
this.sin_hp=this.pos[1][3];
this.sd0=this.pos[1][6];
}
this.ra=interpolate(this.pos[0][4], this.pos[1][4], this.pos[2][4], n);
this.declination=interpolate(this.pos[0][5], this.pos[1][5], this.pos[2][5], n);
this.gst=GSTofUT(jday, ut)+this.eq_equinoxes.RtoH();
this.gha=this.gst-this.ra.RtoH();
this.gha=this.gha.setRange(24)*15;
}
function interpolate(y1, y2, y3, p)
{var a, c;
a=y2-y1;
c=y3-y2-a;
return y1+p*(2*a+c*(2*p-1));
}
function SunPos()
{this.body="Sun";
this.HP0=(8.794/3600).DtoR();
this.SIN_HP0=Math.sin(this.HP0);
this.SD0=959.63/3600;
}
SunPos.prototype.eclipticToEquat=eclipticToEquat;
SunPos.prototype.geoToTopo=geoToTopo;
SunPos.prototype.getTopoSD=getTopoSD;
SunPos.prototype.calcAltAz=calcAltAz;
SunPos.prototype.getGeoAlt=getGeoAlt;
SunPos.prototype.getAppAlt=getAppAlt;
SunPos.prototype.addRefraction=addRefraction;
SunPos.prototype.getRefrAlt=getRefrAlt;
SunPos.prototype.calcDecl=calcDecl;
SunPos.prototype.getAzRS=getAzRS;
SunPos.prototype.getMaxAlt=getMaxAlt;
SunPos.prototype.getPathAngle=getPathAngle;
SunPos.prototype.getEventTime=getEventTime;
SunPos.prototype.initPosArray=initPosArray;
SunPos.prototype.interpPos=interpPos;
SunPos.prototype.calcPos=function(jday, ut)
{var
T,
T1900,
L_s,
G_s,
G_2,
G_4,
G_5,
L_m,
N_m,
lambda,
dist;
T=(jday+ut/24.0+this.deltaT-J2000)/36525.0;
T2=T*T;
T1900=T+1;
L_m=0.60643382+1336.85522267*T-0.00000313*T2;
N_m=0.34734264-5.37260713*T+0.00000579*T2;
L_s=0.77907196+100.00213759*T+0.00000084*T2;
G_s=0.99312619+99.99735956*T-0.00000044*T2;
G_2=0.140023+162.54946219*T;
G_4=0.053856+53.16627469*T;
G_5=0.056531+8.43029617*T;
L_m=L_m.setRange(1)*Math.PI2;
N_m=N_m.setRange(1)*Math.PI2;
L_s=L_s.setRange(1)*Math.PI2;
G_s=G_s.setRange(1)*Math.PI2;
G_2=G_2.setRange(1)*Math.PI2;
G_4=G_4.setRange(1)*Math.PI2;
G_5=G_5.setRange(1)*Math.PI2;
lambda=6910*Math.sin(G_s)
+72*Math.sin(2*G_s)
-17*T1900*Math.sin(G_s)
-7*Math.cos(G_s-G_5)
+6*Math.sin(L_m-L_s)
+5*Math.sin(4*G_s-8*G_4+3*G_5)
-5*Math.cos(2*G_s-2*G_2)
-4*Math.sin(G_s-G_2)
+4*Math.cos(4*G_s-8*G_4+3*G_5)
+3*Math.sin(2*G_s-2*G_2)
-3*Math.sin(G_5)
-3*Math.sin(2*G_s-2*G_5);
lambda-=20;
dist=1.00014
-0.01675*Math.cos(G_s)
-0.00014*Math.cos(2*G_s);
lambda*=(1).DtoR()/3600.0;
this.lambda=(L_s+lambda+this.nut_lon).setRange(Math.PI2);
this.beta=0;
this.dist=dist;
this.sin_hp=this.SIN_HP0/dist;
this.sd0=this.SD0/dist;
this.eclipticToEquat(this.obliquity+this.nut_obl);
this.gst=GSTofUT(jday, ut)+this.eq_equinoxes.RtoH();
this.gha=this.gst-this.ra.RtoH();
this.gha=this.gha.setRange(24)*15;
};
SunPos.prototype.getDayLen=function(rise, set)
{var
day,
rise_event=rise.testEvent(),
set_event=set.testEvent();
if(rise_event===true&&set_event===true){if(rise=F_CIRCUMPOLAR||set>=F_CIRCUMPOLAR){day=24;
}
else
day=0;
}
else if(rise_event===true&&set_event===false)
day=24-rise;
else if(rise_event===false&&set_event===true)
day=set;
return day;
};
function MoonPos()
{this.body="Moon";
this.RADIUS=1738;
}
MoonPos.prototype.eclipticToEquat=eclipticToEquat;
MoonPos.prototype.geoToTopo=geoToTopo;
MoonPos.prototype.getTopoSD=getTopoSD;
MoonPos.prototype.calcAltAz=calcAltAz;
MoonPos.prototype.getGeoAlt=getGeoAlt;
MoonPos.prototype.getTopoAlt=getTopoAlt;
MoonPos.prototype.getAppAlt=getAppAlt;
MoonPos.prototype.addRefraction=addRefraction;
MoonPos.prototype.getRefrAlt=getRefrAlt;
MoonPos.prototype.getAzRS=getAzRS;
MoonPos.prototype.getMaxAlt=getMaxAlt;
MoonPos.prototype.calcAltAz=calcAltAz;
MoonPos.prototype.calcDecl=calcDecl;
MoonPos.prototype.eclipticToTopo=eclipticToTopo;
MoonPos.prototype.getPathAngle=getPathAngle;
MoonPos.prototype.getEventTime=getEventTime;
MoonPos.prototype.initPosArray=initPosArray;
MoonPos.prototype.interpPos=interpPos;
MoonPos.prototype.calcPos=function(jday, ut)
{var
T,
T2,
T1900,
L_m,
G_m,
F_m,
N_m,
D_ms,
G_s,
L_s,
G_2,
d_lambda,
lambda,
lambda_s,
beta,
dist;
T=(jday+ut/24+this.deltaT-J2000)/36525.0;
T2=T*T;
T1900=T+1;
L_m=0.60643382+1336.85522467*T-0.00000313*T2;
G_m=0.37489701+1325.55240982*T+0.00002565*T2;
F_m=0.25909118+1342.22782980*T-0.00000892*T2;
D_ms=0.82736186+1236.85308708*T-0.00000397*T2;
N_m=0.34734264-5.37260513*T+0.00000579*T2;
L_s=0.77907196+100.00213759*T+0.00000084*T2;
G_s=0.99312619+99.99735956*T-0.00000044*T2;
L_2=0.50549945+162.55336953*T+8.628e-07*T2;
L_m=L_m.setRange(1)*Math.PI2;
G_m=G_m.setRange(1)*Math.PI2;
F_m=F_m.setRange(1)*Math.PI2;
D_ms=D_ms.setRange(1)*Math.PI2;
N_m=N_m.setRange(1)*Math.PI2;
L_s=L_s.setRange(1)*Math.PI2;
G_s=G_s.setRange(1)*Math.PI2;
L_2=L_2.setRange(1)*Math.PI2;
d_lambda=6910*Math.sin(G_s)
+72*Math.sin(2*G_s)
-17*T1900*Math.sin(G_s);
d_lambda-=20;
d_lambda*=(1).DtoR()/3600.0;
lambda_s=(L_s+d_lambda+this.nut_lon).setRange(Math.PI2);
d_lambda=22640*Math.sin(G_m)
-4586*Math.sin(G_m-2*D_ms)
+2370*Math.sin(2*D_ms)
+769*Math.sin(2*G_m)
-668*Math.sin(G_s)
-412*Math.sin(2*F_m)
-212*Math.sin(2*G_m-2*D_ms)
-206*Math.sin(G_m-2*D_ms+G_s)
+192*Math.sin(G_m+2*D_ms)
+165*Math.sin(2*D_ms-G_s)
+148*Math.sin(G_m-G_s)
-125*Math.sin(D_ms)
-110*Math.sin(G_m+G_s)
-55*Math.sin(2*F_m-2*D_ms)
-45*Math.sin(G_m+2*F_m)
+40*Math.sin(G_m-2*F_m)
-38*Math.sin(G_m-4*D_ms)
+36*Math.sin(3*G_m)
-31*Math.sin(2*G_m-4*D_ms)
+28*Math.sin(G_m-2*D_ms-G_s)
-24*Math.sin(2*D_ms+G_s)
+19*Math.sin(G_m-D_ms)
+18*Math.sin(D_ms+G_s)
+15*Math.sin(G_m+2*D_ms-G_s)
+14*Math.sin(2*G_m+2*D_ms)
+14*Math.sin(4*D_ms)
-13*Math.sin(3*G_m-2*D_ms)
-11*Math.sin(G_m+16*L_s-18*L_2)
+10*Math.sin(2*G_m-G_s)
+9*Math.sin(G_m-2*F_m-2*D_ms)
+9*Math.cos(G_m+16*L_s-18*L_2)
-9*Math.sin(2*G_m-2*D_ms+G_s)
-8*Math.sin(G_m+D_ms)
+8*Math.sin(2*D_ms-2*G_s)
-8*Math.sin(2*G_m+G_s)
-7*Math.sin(2*G_s)
-7*Math.sin(G_m-2*D_ms+2*G_s)
+7*Math.sin(N_m)
-6*Math.sin(G_m-2*F_m+2*D_ms)
-6*Math.sin(2*F_m+2*D_ms)
-4*Math.sin(G_m-4*D_ms+G_s)
+4*T1900*Math.cos(G_m+16*L_s-18*L_2)
-4*Math.sin(2*G_m+2*F_m)
+4*T1900*Math.sin(G_m+16*L_s-18*L_2)
+3*Math.sin(G_m-3*D_ms)
-3*Math.sin(G_m+2*D_ms+G_s)
-3*Math.sin(2*G_m-4*D_ms+G_s)
+3*Math.sin(G_m-2*G_s)
+3*Math.sin(G_m-2*D_ms-2*G_s)
-2*Math.sin(2*G_m-2*D_ms-G_s)
-2*Math.sin(2*F_m-2*D_ms+G_s)
+2*Math.sin(G_m+4*D_ms)
+2*Math.sin(4*G_m)
+2*Math.sin(4*D_ms-G_s)
+2*Math.sin(2*G_m-D_ms);
beta=18461*Math.sin(F_m)
+1010*Math.sin(G_m+F_m)
+1000*Math.sin(G_m-F_m)
-624*Math.sin(F_m-2*D_ms)
-199*Math.sin(G_m-F_m-2*D_ms)
-167*Math.sin(G_m+F_m-2*D_ms)
+117*Math.sin(F_m+2*D_ms)
+62*Math.sin(2*G_m+F_m)
+33*Math.sin(G_m-F_m+2*D_ms)
+32*Math.sin(2*G_m-F_m)
-30*Math.sin(F_m-2*D_ms+G_s)
-16*Math.sin(2*G_m+F_m-2*D_ms)
+15*Math.sin(G_m+F_m+2*D_ms)
+12*Math.sin(F_m-2*D_ms-G_s)
-9*Math.sin(G_m-F_m-2*D_ms+G_s)
-8*Math.sin(F_m+N_m)
+8*Math.sin(F_m+2*D_ms-G_s)
-7*Math.sin(G_m+F_m-2*D_ms+G_s)
+7*Math.sin(G_m+F_m-G_s)
-7*Math.sin(G_m+F_m-4*D_ms)
-6*Math.sin(F_m+G_s)
-6*Math.sin(3*F_m)
+6*Math.sin(G_m-F_m-G_s)
-5*Math.sin(F_m+D_ms)
-5*Math.sin(G_m+F_m+G_s)
-5*Math.sin(G_m-F_m+G_s)
+5*Math.sin(F_m-G_s)
+5*Math.sin(F_m-D_ms)
+4*Math.sin(3*G_m+F_m)
-4*Math.sin(F_m-4*D_ms)
-3*Math.sin(G_m-F_m-4*D_ms)
+3*Math.sin(G_m-3*F_m)
-2*Math.sin(2*G_m-F_m-4*D_ms)
-2*Math.sin(3*F_m-2*D_ms)
+2*Math.sin(2*G_m-F_m+2*D_ms)
+2*Math.sin(G_m-F_m+2*D_ms-G_s)
+2*Math.sin(2*G_m-F_m-2*D_ms)
+2*Math.sin(3*G_m-F_m);
dist=60.36298
-3.27746*Math.cos(G_m)
-0.57994*Math.cos(G_m-2*D_ms)
-0.46357*Math.cos(2*D_ms)
-0.08904*Math.cos(2*G_m)
+0.03865*Math.cos(2*G_m-2*D_ms)
-0.03237*Math.cos(2*D_ms-G_s)
-0.02688*Math.cos(G_m+2*D_ms)
-0.02358*Math.cos(G_m-2*D_ms+G_s)
-0.02030*Math.cos(G_m-G_s)
+0.01719*Math.cos(D_ms)
+0.01671*Math.cos(G_m+G_s)
+0.01247*Math.cos(G_m-2*F_m)
+0.00704*Math.cos(G_s)
+0.00529*Math.cos(2*D_ms+G_s)
-0.00524*Math.cos(G_m-4*D_ms)
+0.00398*Math.cos(G_m-2*D_ms-G_s)
-0.00366*Math.cos(3*G_m)
-0.00295*Math.cos(2*G_m-4*D_ms)
-0.00263*Math.cos(D_ms+G_s)
+0.00249*Math.cos(3*G_m-2*D_ms)
-0.00221*Math.cos(G_m+2*D_ms-G_s)
+0.00185*Math.cos(2*F_m-2*D_ms)
-0.00161*Math.cos(2*D_ms-2*G_s)
+0.00147*Math.cos(G_m+2*F_m-2*D_ms)
-0.00142*Math.cos(4*D_ms)
+0.00139*Math.cos(2*G_m-2*D_ms+G_s)
-0.00118*Math.cos(G_m-4*D_ms+G_s)
-0.00116*Math.cos(2*G_m+2*D_ms)
-0.00110*Math.cos(2*G_m-G_s);
this.dist=dist;
this.sin_hp=1/dist;
d_lambda*=(1).DtoR()/3600;
lambda=(L_m+d_lambda+this.nut_lon).setRange(Math.PI2);
beta*=(1).DtoR()/3600;
this.sd0=Math.asin(this.RADIUS/EARTH_RADIUS*this.sin_hp).RtoD();
this.lambda=lambda;
this.beta=beta;
this.lambda_s=lambda_s;
this.eclipticToEquat(this.obliquity+this.nut_obl);
this.gst=GSTofUT(jday, ut)+this.eq_equinoxes.RtoH();
this.gha=this.gst-this.ra.RtoH();
this.gha=this.gha.setRange(24)*15;
};
MoonPos.prototype.calcPhase=function(sun_dist, loc)
{var
lambda,
beta,
cos_elong,
elong,
phase_angle,
moon_dist;
sun_dist*=Math.AU;
moon_dist=this.dist*EARTH_RADIUS;
this.eclipticToTopo(loc);
lambda=this.lambda_t;
beta=this.beta_t;
cos_elong=Math.cos(beta)*Math.cos(lambda-this.lambda_s);
elong=Math.acos(cos_elong);
phase_angle=Math.atan2(sun_dist*Math.sin(elong),
moon_dist-sun_dist*cos_elong);
this.phase=(1+Math.cos(phase_angle))/2;
this.phase_angle=Math.abs(phase_angle.RtoD());
if(elong.setRange(Math.PI2)'+LF
+''+LF;
}
function add_help(name, href)
{if(!href)
return name;
return ""+name+"";
}
function add_place_link(name, ndx)
{var
tag=''
+name+'';
return tag;
}
function add_date_link(tipStr, dateStr, date)
{var
tag=''
+dateStr+'';
return tag;
}
function begin_document(title)
{var
d,
handle;
var
wid=0.98*screen.availWidth,
ht=0.86*screen.availHeight,
xpos=5,
ypos=0.005*screen.availHeight,
win_size,
win_attribs,
win_features;
if(flags.use_popups===true){win_size=",height="+ht+",width="+wid+","+win_attribs;
win_attribs="resizable=1,menubar=1,toolbar=0,status=1,scrollbars=1";
win_pos="left="+xpos+","+"top="+ypos;
win_features=win_size+","+win_pos+","+win_attribs;
handle=window.open("",title,win_features);
}
else
handle=window.open("",title);
handle.document.close();
handle.document.open("text/html");
return handle;
}
function add_buttons()
{var form;
if(flags.print!==true){form=""+LF;
this.buttons=form;
}
else
this.buttons=null;
}
function build_positions_header(hdr)
{var header;
if(arguments.length>0)
header=hdr;
else
header=BR+LF;
header+=POS_TABLE;
header+=POS_TABLE_COLSPEC;
header+=THD+LF;
header+=TR+LF;
header+=TH21+"Time";
if(flags.show_ut===true)
header+=BR+" (UT)";
header+=THE+LF;
header+=TH14+"Sun"+THE+LF;
header+=TH15+"Moon"+THE+LF;
header+=TRE+LF;
header+=TR+LF;
if(flags.show_mag_north===true)
header+=TH11+add_help("Azimuth*", "posAzimuth")+THE+LF;
else
header+=TH11+add_help("Azimuth", "posAzimuth")+THE+LF;
header+=TH11+add_help("Altitude", "posAltitude")+THE+LF;
header+=TH11+add_help("Path"+BR+"Angle", "posPathAngle")+THE+LF;
header+=TH11+add_help("Shadow"+BR+"Length", "posShadowLength")+THE+LF;
if(flags.show_mag_north===true)
header+=TH11+add_help("Azimuth*", "posAzimuth")+THE+LF;
else
header+=TH11+add_help("Azimuth", "posAzimuth")+THE+LF;
header+=TH11+add_help("Altitude", "posAltitude")+THE+LF;
header+=TH11+add_help("Path"+BR+"Angle", "posPathAngle")+THE+LF;
header+=TH11+add_help("Phase", "posPhase")+THE+LF;
header+=TH11+add_help("SD", "posSD")+THE+LF;
header+=TRE+LF;
header+=THDE+LF;
header+=TBD+LF;
this.positions_header=header;
}
function build_results_doc_head(title)
{var head;
head = DOCTYPE + LF;
head+=""+LF;
head+=""+LF;
head+=CONTENT+LF;
head+="Sun/Moon "+title+""+LF;
head+=SCRIPT;
head+=STYLE;
if(flags.print===true&&flags.use_font_size===true)
head+=PRINT_STYLE;
head+=""+LF;
this.doc_head=head;
}
function output_results(type)
{var i,
handle=begin_document(this.name),
d=handle.document;
d.write(this.doc_head);
d.write(this.page_top);
if(this.rs_criteria)
d.write(this.rs_criteria);
d.write(this.times_header);
if(this.times.length>0){for(i=0; i0){for(i=0; i', LF, '', LF);
d.close();
handle.focus();
if(handle.print&&flags.print===true)
handle.print();
}
function build_page_top(calc_params, rs_opts, loc)
{var header, title, time_zone, elev_units;
if(flags.elevation_in_feet===true)
elev_units="ft";
else
elev_units="m";
header=""+LF;
title="Location: " + loc.name.replace(/'/g, "’");
header+=SP+title.big().bold()+BR+LF;
if(loc.latitude<0)
title="Latitude: "+(-loc.latitude).toDMS(2)+" S";
else
title="Latitude: "+loc.latitude.toDMS(2)+" N";
header+=SP+title.bold()+",";
if(loc.longitude<0)
title="Longitude: "+(-loc.longitude).toDMS(2)+" W";
else
title="Longitude: "+loc.longitude.toDMS(2)+" E";
header+=SP+title.bold();
if(flags.elevation_in_feet===true)
loc.elevation=loc.elevation.MtoFT();
if(loc.elevation<0)
title="Elevation: "+MINUS+addCommas(Math.round(-loc.elevation));
else if(loc.elevation>0)
title="Elevation: "+addCommas(Math.round(loc.elevation));
title+=" "+elev_units;
if(loc.elevation!==0)
header+=","+SP+title.bold();
header+=BR+LF;
title="Time Zone: ";
if(Math.floor(loc.timezone)!=loc.timezone){if(loc.timezone<0)
title+="UTC"+MINUS+(-loc.timezone).toDMS();
else
title+="UTC+"+loc.timezone.toDMS();
}
else{if(loc.timezone<0)
title+="UTC"+MINUS+(-loc.timezone);
else
title+="UTC+"+loc.timezone;
}
if((time_zone=loc.getTimezone()))
title+=" ("+time_zone+")";
header+=SP+title.bold()+BR+LF;
if(loc.mag_dec_date_valid!=INVALID||
flags.show_mag_north===true||flags.show_az_adjust===true){title=add_help("Magnetic Declination: ", "magdec");
if(loc.mag_dec_mean<0)
title+=(-loc.mag_dec_mean).toFixed(1)+DEG+" W";
else
title+=loc.mag_dec_mean.toFixed(1)+DEG+" E";
title=title.bold();
if(calc_params.end_date-calc_params.start_date>180)
title+=" on "+loc.mag_date_str;
if(loc.mag_dec_date_valid!=VALID)
title+=SP2+"(Caution: outside magnetic model range of "
+mag_model.valid_range_str+")";
header+=SP+title+BR+LF;
}
title=(flags.show_mag_north===true ? "Magnetic" : "True")+" North";
title="Azimuths relative to "+title.bold();
if(flags.show_az_adjust===true&&Math.abs(loc.mag_dec_mean)>=0.05){var az_adjustment;
if(loc.mag_dec_mean<0){az_adjustment=(-loc.mag_dec_mean).toFixed(1);
if(flags.show_mag_north===true)
title+=" (subtract "+az_adjustment+DEG+" for True North)";
else
title+=" (add "+az_adjustment+DEG+" for Magnetic North)";
}
else{az_adjustment=loc.mag_dec_mean.toFixed(1);
if(flags.show_mag_north===true)
title+=" (add "+az_adjustment+DEG+" for True North)";
else
title+=" (subtract "+az_adjustment+DEG+" for Magnetic North)";
}
}
if(calc_params.calc_type==TIMES&&flags.show_mag_north===true&&
loc.mag_dec_valid===false){title+=SP2+"(Caution: ";
if(calc_params.start_date>mag_model.end_date+MagModel.GRACE||
calc_params.end_date0){if(flags.elevation_in_feet===true)
loc.height=loc.height.MtoFT();
title=add_help("Height above Horizon:", "HeightAboveHorizon")
+SP+addCommas(loc.height.toFixed(0))+" "+elev_units
+SP+"(elevation at horizon:"
+SP+addCommas((loc.elevation-loc.height).toFixed(0))+" "+elev_units+")"
+BR+LF;
title+=SP+add_help("Dip of the horizon: ", "dip")+loc.dip.toFixed(2)+"°";
header+=SP+title+BR+LF;
}
this.page_top=header;
}
function add_rs_criteria(calc_params, rs_opts)
{if(rs_opts.selections===false){
this.rs_criteria=null;
if(rs_opts.rs_alts===false)
return;
}
var
header="",
date_step_units=[ "days", "weeks", "months (30 day)", "months (1/12 year)", "years" ];
if(rs_opts.selections===true){var
StartDate=new Date(calc_params.start_date.DayToMsec()),
EndDate=new Date(calc_params.end_date.DayToMsec());
header=SP+'Dates: '+StartDate.toLocDateStr(0)
+' to '+EndDate.toLocDateStr(0)+BR+LF;
if(calc_params.date_step!=1)
header+=SP+'Date Interval: '
+(calc_params.date_step/calc_params.date_step_units_value)
+' '+date_step_units[calc_params.date_step_units_ndx]+BR+LF;
}
if(rs_opts.alt_sr_range===true){header+=SP+"Sunrise Altitude "
+(rs_opts.alt_sr_type=="center" ? "(center): " : "(upper limb): ");
if(rs_opts.alt_sr_min<0)
header+=MINUS+-Math.rnd(rs_opts.alt_sr_min, 2);
else
header+=Math.rnd(rs_opts.alt_sr_min, 2);
header+=DEG+" to ";
if(rs_opts.alt_sr_max<0)
header+=MINUS+-Math.rnd(rs_opts.alt_sr_max, 2);
else
header+=Math.rnd(rs_opts.alt_sr_max, 2);
header+=DEG+BR+LF;
}
else if(rs_opts.alt_sr!==0||rs_opts.alt_sr_type=="center"){header+=SP+"Sunrise Altitude: "
+(rs_opts.alt_sr_type=="center" ? "center at " : "upper limb at ");
if(rs_opts.alt_sr<0)
header+=MINUS+-Math.rnd(rs_opts.alt_sr, 2);
else
header+=Math.rnd(rs_opts.alt_sr, 2);
header+=DEG+BR+LF;
}
if(rs_opts.az_sr_check===true){header+=SP+"Sunrise Azimuth: "+Math.rnd(rs_opts.az_sr_min, 2)+DEG;
if(rs_opts.az_sr_max!=rs_opts.az_sr_min)
header+=" to "+Math.rnd(rs_opts.az_sr_max, 2)+DEG;
header+=BR+LF;
}
if(rs_opts.alt_ss_range===true){header+=SP+"Sunset Altitude "
+(rs_opts.alt_ss_type=="center" ? "(center): " : "(upper limb): ");
if(rs_opts.alt_ss_min<0)
header+=MINUS+-Math.rnd(rs_opts.alt_ss_min, 2);
else
header+=Math.rnd(rs_opts.alt_ss_min, 2);
header+=DEG+" to ";
if(rs_opts.alt_ss_max<0)
header+=MINUS+-Math.rnd(rs_opts.alt_ss_max, 2);
else
header+=Math.rnd(rs_opts.alt_ss_max, 2);
header+=DEG+BR+LF;
}
else if(rs_opts.alt_ss!==0||rs_opts.alt_ss_type=="center"){header+=SP+"Sunset Altitude: "
+(rs_opts.alt_ss_type=="center" ? "center at " : "upper limb at ");
if(rs_opts.alt_ss<0)
header+=MINUS+-Math.rnd(rs_opts.alt_ss, 2);
else
header+=Math.rnd(rs_opts.alt_ss, 2);
header+=DEG+BR+LF;
}
if(rs_opts.az_ss_check===true){header+=SP+"Sunset Azimuth: "+Math.rnd(rs_opts.az_ss_min, 2)+DEG;
if(rs_opts.az_ss_max!=rs_opts.az_ss_min)
header+=" to "+Math.rnd(rs_opts.az_ss_max, 2)+DEG;
header+=BR+LF;
}
if(rs_opts.alt_mr_range===true){header+=SP+"Moonrise Altitude "
+(rs_opts.alt_mr_type=="center" ? "(center): " : "(upper limb): ")
+Math.rnd(rs_opts.alt_mr_min, 2)+DEG
+" to "+Math.rnd(rs_opts.alt_mr_max, 2)+DEG+BR+LF;
}
else if(rs_opts.alt_mr!==0||rs_opts.alt_mr_type=="center"){header+=SP+"Moonrise Altitude:"+SP
+(rs_opts.alt_mr_type=="center" ? "center at " : "upper limb at ")
+Math.rnd(rs_opts.alt_mr, 2)+DEG+BR+LF;
}
if(rs_opts.az_mr_check===true){header+=SP+"Moonrise Azimuth: "+Math.rnd(rs_opts.az_mr_min, 2)+DEG;
if(rs_opts.az_mr_max!=rs_opts.az_mr_min)
header+=" to "+Math.rnd(rs_opts.az_mr_max, 2)+DEG;
header+=BR+LF;
}
if(rs_opts.alt_ms_range===true){header+=SP+"Moonset Altitude "
+(rs_opts.alt_ms_type=="center" ? "(center): " : "(upper limb): ")
+Math.rnd(rs_opts.alt_ms_min, 2)+DEG
+" to "+Math.rnd(rs_opts.alt_ms_max, 2)+DEG+BR+LF;
}
else if(rs_opts.alt_ms!==0||rs_opts.alt_ms_type=="center"){header+=SP+"Moonset Altitude: "
+(rs_opts.alt_ms_type=="center" ? "center at " : "upper limb at ")
+Math.rnd(rs_opts.alt_ms, 2)+DEG+BR+LF;
}
if(rs_opts.az_ms_check===true){header+=SP+"Moonset Azimuth: "+Math.rnd(rs_opts.az_ms_min, 2)+DEG;
if(rs_opts.az_ms_max!=rs_opts.az_ms_min)
header+=" to "+Math.rnd(rs_opts.az_ms_max, 2)+DEG;
header+=BR+LF;
}
if(rs_opts.phase_mr_check===true&&
(rs_opts.phase_mr_min>0||rs_opts.phase_mr_max<1)){header+=SP+"Moonrise Phase: "+Math.rnd(rs_opts.phase_mr_min, 3)
+" to "+Math.rnd(rs_opts.phase_mr_max, 3)+BR+LF;
}
if(rs_opts.phase_ms_check===true&&
(rs_opts.phase_ms_min>0||rs_opts.phase_ms_max<1)){header+=SP+"Moonset Phase: "+Math.rnd(rs_opts.phase_ms_min, 3)
+" to "+Math.rnd(rs_opts.phase_ms_max, 3)+BR+LF;
}
if(rs_opts.tdiff_mr_check===true){header+=SP+"Moonrise ";
if(rs_opts.tdiff_mr_min<0)
header+=-rs_opts.tdiff_mr_min+" min before "+rs_opts.tdiff_mr_type;
else
header+=rs_opts.tdiff_mr_min+" min after "+rs_opts.tdiff_mr_type;
header+=" to ";
if(rs_opts.tdiff_mr_max<0)
header+=-rs_opts.tdiff_mr_max+" min before "+rs_opts.tdiff_mr_type;
else
header+=rs_opts.tdiff_mr_max+" min after "+rs_opts.tdiff_mr_type;
header+=BR+LF;
}
if(rs_opts.tdiff_ms_check===true){header+=SP+"Moonset ";
if(rs_opts.tdiff_ms_min<0)
header+=-rs_opts.tdiff_ms_min+" min before "+rs_opts.tdiff_ms_type;
else
header+=rs_opts.tdiff_ms_min+" min after "+rs_opts.tdiff_ms_type;
header+=" to ";
if(rs_opts.tdiff_ms_max<0)
header+=-rs_opts.tdiff_ms_max+" min before "+rs_opts.tdiff_ms_type;
else
header+=rs_opts.tdiff_ms_max+" min after "+rs_opts.tdiff_ms_type;
header+=BR+LF;
}
if(flags.show_ut===true)
header+=SP+"Times in "+"UT".bold()+BR+LF;
if(rs_opts.selections===true){if(rs_opts.matches==1)
header+=BR+SP+rs_opts.matches+' date matches criteria:';
else if(rs_opts.matches>0)
header+=BR+SP+rs_opts.matches+' dates match criteria:';
header+=BR+LF;
}
this.rs_criteria=header;
}
function build_times_header(rs_opts)
{var header;
if(this.times.length===0){this.times_header=BR+SP2+"****No dates meet criteria****".big().bold()+BR;
return;
}
header=BR+LF;
if(rs_opts.show_tdiff_mr===true||rs_opts.show_tdiff_ms===true){header+=TIMES_TABLE2;
header+=TIMES_TABLE_COLSPEC2;
}
else{header+=TIMES_TABLE1;
header+=TIMES_TABLE_COLSPEC1;
}
header+=THD+LF;
header+=TR+LF;
header+=TH31+"Day"+THE+LF;
header+=TH31+"Date"+THE+LF;
header+=TH19+"Sun"+THE+LF;
if(rs_opts.show_tdiff_mr===true||rs_opts.show_tdiff_ms===true)
header+=TH19+"Moon"+THE+LF;
else
header+=TH18+"Moon"+THE+LF;
header+=TRE+LF;
header+=TR+LF;
header+=TH21+add_help("Dawn", "timeDawn")+THE+LF;
header+=TH21+add_help("Rise", "timeRise");
if(rs_opts.alt_sr_type=="center")
header+=RS_CENTER;
if(rs_opts.alt_sr!==0){header+=BR;
if(rs_opts.alt_sr<0)
header+=MINUS+-Math.rnd(rs_opts.alt_sr, 1)+DEG;
else
header+=Math.rnd(rs_opts.alt_sr, 1)+DEG;
}
header+=THE+LF;
header+=TH21+add_help("Transit", "timeTransit")+THE+LF;
header+=TH21+add_help("Set", "timeSet");
if(rs_opts.alt_ss_type=="center")
header+=RS_CENTER;
if(rs_opts.alt_ss!==0){header+=BR;
if(rs_opts.alt_ss<0)
header+=MINUS+-Math.rnd(rs_opts.alt_ss, 1)+DEG;
else
header+=Math.rnd(rs_opts.alt_ss, 1)+DEG;
}
header+=THE+LF;
header+=TH21+add_help("Dusk", "timeDusk")+THE+LF;
header+=TH21+add_help("Day"+BR+"Length", "dayLength")+THE+LF;
if(flags.show_mag_north===true)
header+=TH12+add_help("Azimuth*", "timeRiseAzimuth")+THE+LF;
else
header+=TH12+add_help("Azimuth", "timeRiseAzimuth")+THE+LF;
header+=TH21+add_help("Max"+BR+"Alt", "timeMaxAlt")+THE+LF;
header+=TH21+add_help("Rise", "timeRise");
if(rs_opts.alt_mr_type=="center")
header+=RS_CENTER;
if(rs_opts.alt_mr!==0){header+=BR;
header+=Math.rnd(rs_opts.alt_mr, 1)+DEG;
}
header+=THE+LF;
header+=TH21+add_help("Transit", "timeTransit")+THE+LF;
header+=TH21+add_help("Set", "timeSet");
if(rs_opts.alt_ms_type=="center")
header+=RS_CENTER;
if(rs_opts.alt_ms!==0){header+=BR;
header+=Math.rnd(rs_opts.alt_ms, 1)+DEG;
}
header+=THE+LF;
if(rs_opts.show_tdiff_mr===true||rs_opts.show_tdiff_ms===true)
header+=TH21+add_help("Time"+BR+"Diff", "timeTimeDiff")+THE+LF;
if(flags.show_mag_north===true)
header+=TH12+add_help("Azimuth*", "timeRiseAzimuth")+THE+LF;
else
header+=TH12+add_help("Azimuth", "timeRiseAzimuth")+THE+LF;
header+=TH21+add_help("Max"+BR+"Alt", "timeMaxAlt")+THE+LF;
switch(rs_opts.show_phase){case "rise":
header+=TH12+add_help("At Rise", "timePhase")+THE+LF;
break;
case "transit":
header+=TH12+add_help("At Transit", "timePhase")+THE+LF;
break;
case "set":
header+=TH12+add_help("At Set", "timePhase")+THE+LF;
break;
default:
break;
}
header+=TRE+LF;
header+=TR+LF;
header+=TH11+add_help("Rise", "timeRiseAzimuth")+THE+LF;
header+=TH11+add_help("Set", "timeSetAzimuth")+THE+LF;
header+=TH11+add_help("Rise", "timeRiseAzimuth")+THE+LF;
header+=TH11+add_help("Set", "timeSetAzimuth")+THE+LF;
header+=TH11+add_help("Phase", "timePhase")+THE+LF;
header+=TH11+add_help("SD", "timeSD")+THE+LF;
header+=TRE+LF;
header+=THDE+LF;
header+=TBD+LF;
this.times_header=header;
}
function format_times(cdate, ctype, s, m, rs_opts, t_offset, loc, d)
{var
cdate_time,
dateStr,
t_diff,
rpad_d,
rpad_t,
NO_EVENT=MDASH,
line="";
if(flags.show_am_pm_time===true&&flags.show_ut===false)
rpad_d=0;
else if(flags.show_ut===true){rpad_d=0;
rpad_t=0;
}
else{rpad_d=1;
rpad_t=1;
}
line+=TR_R+LF;
line+=TD_L+cdate.getLocDayName(loc)+TDE+LF;
if(ctype==TIMES)
line+=TD_L
+add_date_link(cdate.toLocDateStr(0), cdate.toLocDateStr(1),
cdate.getTime().MsecToDay());
else
line+=TD_L+cdate.toLocDateStr(1);
if(flags.show_ut===true)
line+=" (UT)";
else if(cdate.getDayDST(loc)===true)
line+="*";
line+=TDE+LF;
line+=s.dawn.showEventTime(cdate, loc, "T", rpad_d)+LF;
line+=s.rise.showEventTime(cdate, loc, "", 0)+LF;
line+=s.transit.showEventTime(cdate, loc, "", rpad_t)+LF;
line+=s.set.showEventTime(cdate, loc, "", 0)+LF;
line+=s.dusk.showEventTime(cdate, loc, "T", 0)+LF;
if(s.day.testEvent()===true)
line+=TD+(s.day).toDMS()+SP2+TDE+LF;
else
line+=TD_C+NO_EVENT+TDE+LF;
if(s.az_rise===null)
line+=TD_C+NO_EVENT+TDE+LF;
else
line+=TD+Math.round(s.az_rise)+SP+TDE+LF;
if(s.az_set===null)
line+=TD_C+NO_EVENT+TDE+LF;
else
line+=TD+Math.round(s.az_set)+TDE+LF;
if(s.transit.testEvent()===false)
line+=TD_C+NO_EVENT+TDE+LF;
else
line+=s.maxalt.showMaxAlt()+LF;
line+=m.rise.showEventTime(cdate, loc, "", 0)+LF;
line+=m.transit.showEventTime(cdate, loc, "", rpad_t)+LF;
line+=m.set.showEventTime(cdate, loc, "", 0)+LF;
if(rs_opts.show_tdiff_mr===true){if(m.az_rise===null)
line+=TD_C+NO_EVENT+TDE+LF;
else{if(rs_opts.tdiff_mr_type=="sunset")
t_diff=m.rise-s.set;
else
t_diff=m.rise-s.rise;
if(t_diff<0)
line+=TD+MINUS+(-t_diff).toDMS()+TDE+LF;
else
line+=TD+t_diff.toDMS()+TDE+LF;
}
}
else if(rs_opts.show_tdiff_ms===true){if(m.az_set===null)
line+=TD_C+NO_EVENT+TDE+LF;
else{if(rs_opts.tdiff_ms_type=="sunset")
t_diff=m.set-s.set;
else
t_diff=m.set-s.rise;
if(t_diff<0)
line+=TD+MINUS+(-t_diff).toDMS()+TDE+LF;
else
line+=TD+t_diff.toDMS()+TDE+LF;
}
}
if(m.az_rise===null)
line+=TD_C+NO_EVENT+TDE+LF;
else
line+=TD+Math.round(m.az_rise)+SP+TDE+LF;
if(m.az_set===null)
line+=TD_C+NO_EVENT+TDE+LF;
else
line+=TD+Math.round(m.az_set)+TDE+LF;
if(m.transit.testEvent()===false)
line+=TD_C+NO_EVENT+TDE+LF;
else
line+=m.maxalt.showMaxAlt()+LF;
switch(rs_opts.show_phase){case "rise":
if(m.az_rise===null){line+=TD_C+NO_EVENT+TDE+LF;
line+=TD_C+NO_EVENT+TDE+LF;
}
else{
line+=TD_L+m.phase_r.toFixed(3);
if(rs_opts.phase_check===true)
line+=(m.waxing_r===true ? PLUS : MINUS);
line+=TDE+LF;
line+=TD+m.sd_r.toFixed(3)+TDE+LF;
}
break;
case "transit":
if(m.transit.testEvent()===false){line+=TD_C+NO_EVENT+TDE+LF;
line+=TD_C+NO_EVENT+TDE+LF;
}
else{
line+=TD_L+m.phase_t.toFixed(3);
if(rs_opts.phase_check===true)
line+=(m.waxing_t===true ? PLUS : MINUS);
line+=TDE+LF;
line+=TD+m.sd_t.toFixed(3)+TDE+LF;
}
break;
case "set":
if(m.az_set===null){line+=TD_C+NO_EVENT+TDE+LF;
line+=TD_C+NO_EVENT+TDE+LF;
}
else{
line+=TD_L+m.phase_s.toFixed(3);
if(rs_opts.phase_check===true)
line+=(m.waxing_s===true ? PLUS : MINUS);
line+=TDE+LF;
line+=TD+m.sd_s.toFixed(3)+TDE+LF;
}
break;
default:
break;
}
line+=TRE+LF;
return line;
}
function build_places_doc_head()
{
var PLACES_STYLE=''+LF;
var header;
header = DOCTYPE + LF;
header+=""+LF;
header+=""+LF;
header+=CONTENT+LF;
header+="Sun/Moon Calculator Locations"+LF;
header+=SCRIPT;
header+=STYLE;
header+=PLACES_STYLE;
if(flags.print===true&&flags.use_font_size===true)
header+=PRINT_STYLE;
header+=""+LF;
this.doc_head=header;
}
function build_places_header(year, date)
{var elev_units, header;
if(flags.elevation_in_feet===true)
elev_units="ft";
else
elev_units="m";
header=""+LF;
header+="Sun/Moon Calculator Locations".big().bold()+VSPACE+LF;
if(flags.show_mag_dec===true){var mag_date=new Date();
mag_date.setTime(date.DayToMsec());
header+=SP+"Magnetic declinations on "+mag_date.toLocDateStr(0);
if(datemag_model.end_date+MagModel.GRACE)
header+=SP2+"(Caution: outside magnetic model range of "
+mag_model.valid_range_str+")";
header+=VSPACE+LF;
}
this.page_top=header;
header=PLACES_TABLE;
header+=PLACES_TABLE_COLSPEC1;
if(flags.show_dst_rules===true)
header+=PLACES_TABLE_COLSPEC2;
if(flags.show_dst_dates===true)
header+=PLACES_TABLE_COLSPEC2;
header+=THD+LF;
header+=TR+LF;
if(flags.show_dst_rules===true||flags.show_dst_dates===true){header+=TH21+"Name"+THE+LF;
header+=TH21+"Latitude"+THE+LF;
header+=TH21+"Longitude"+THE+LF;
header+=TH21+"Elevation,"+BR+elev_units+THE+LF;
if(flags.show_mag_dec===true)
header+=TH21+"Magnetic"+BR+"Declination"+THE+LF;
}
else{header+=TH11+"Name"+THE+LF;
header+=TH11+"Latitude"+THE+LF;
header+=TH11+"Longitude"+THE+LF;
header+=TH11+"Elevation,"+BR+elev_units+THE+LF;
if(flags.show_mag_dec===true)
header+=TH11+"Magnetic"+BR+"Declination"+THE+LF;
}
if(flags.show_mag_dec===true&&flags.show_dst_rules===true&&flags.show_dst_dates===true)
header+=TH21+"Time"+BR+"Zone"+THE+LF;
else if(flags.show_dst_rules===true||flags.show_dst_dates===true)
header+=TH21+"Time Zone"+THE+LF;
else
header+=TH11+"Time Zone"+THE+LF;
if(flags.show_dst_rules===true)
header+=TH12+"DST Rule for "+year+THE+LF;
if(flags.show_dst_dates===true)
header+=TH12+"DST Dates for "+year+THE+LF;
if(flags.show_dst_rules===true||flags.show_dst_dates===true){header+=TRE+LF;
header+=TR+LF;
}
if(flags.show_dst_rules===true){header+=TH11+"Start"+THE+LF;
header+=TH11+"End"+THE+LF;
}
if(flags.show_dst_dates===true){header+=TH11+"Start"+THE+LF;
header+=TH11+"End"+THE+LF;
}
header+=TRE+LF;
header+=THDE+LF;
this.places_header=header;
}
function build_place_list(year, date)
{var
ERRMSG="***Missing DST rule***",
NO_DST_MSG="DST not observed",
TD_ERRMSG,
TD_L_ERR=PAD8+'',
TD_MSG,
TD_NAME;
var
errors=0,
line,
ndx,
i=0,
msg_wid=0;
var
loc=new Place();
if(flags.show_dst_rules===true)
msg_wid++;
if(flags.show_dst_dates===true)
msg_wid++;
switch (msg_wid){case 1:
TD_MSG=TD12;
TD_ERRMSG=PAD8+' | ';
break;
case 2:
TD_MSG=TD14;
TD_ERRMSG=PAD8+' | ';
break;
default:
break;
}
place_list.length=0;
for(ndx=0; places[ndx]; ndx++){loc.Init(places[ndx]);
TD_NAME=TD_L;
name=add_place_link(loc.name.replace(/'/g, "’"), ndx);
line=TD_L+SP;
if(loc.latitude<0)
line+=(-loc.latitude).toDMSLeftPadded(2)+SP+'S';
else
line+=(loc.latitude).toDMSLeftPadded(2)+SP+'N';
line+=TDE+LF;
line+=TD_L+SP;
if(loc.longitude<0)
line+=(-loc.longitude).toDMSLeftPadded(3)+SP+'W';
else
line+=(loc.longitude).toDMSLeftPadded(3)+SP+'E';
line+=TDE+LF;
if(flags.elevation_in_feet===true)
loc.elevation=loc.elevation.MtoFT();
if(loc.elevation<0)
line+=TD_R+SP+MINUS+(-loc.elevation).toFixed(0)+SP;
else if(loc.elevation>0)
line+=TD_R+SP+(loc.elevation).toFixed(0)+SP;
else
line+=TD_R+MDASH+SP2;
line+=TDE+LF;
if(flags.show_mag_dec===true){line+=TD_L+SP2;
loc.mag_dec=mag_model.calcMagDec(loc, date);
if(loc.mag_dec<0)
line+=(-loc.mag_dec).toFixedLeftPadded(1, 3)+' W';
else
line+=loc.mag_dec.toFixedLeftPadded(1, 3)+' E';
line+=TDE+LF;
}
line+=TD_L+SP;
if(loc.timezone<0)
line+='UTC'+MINUS+(-loc.timezone);
else
line+='UTC'+PLUS+loc.timezone;
if(! (flags.show_dst_rules===true&&flags.show_dst_dates===true)){var time_zone;
if((time_zone=loc.getTimezone()))
line+=" ("+time_zone+")";
}
line+=TDE+LF;
if(flags.show_dst_rules===true||flags.show_dst_dates===true){if(loc.uses_dst==1){
var
start_month,
start_day,
start_week,
start_DoW,
end_month,
end_day,
end_week,
end_DoW;
if(loc.getDSTLimits(year)===false){name=name.bold();
TD_NAME=TD_L_ERR;
line+=TD_ERRMSG+ERRMSG.bold()+TDE+LF;
}
else if(loc.dst_start===null||loc.dst_end===null)
line+=TD_MSG+NO_DST_MSG+TDE+LF;
else{if(flags.show_dst_rules===true){
line+=TD_L+SP+showDSTRule(loc.DSTStartStr)+TDE+LF;
line+=TD_L+SP+showDSTRule(loc.DSTEndStr)+TDE+LF;
}
if(flags.show_dst_dates===true){line+=TD_L+SP+showDSTDate(loc.DSTStartStr, year)+TDE+LF;
line+=TD_L+SP+showDSTDate(loc.DSTEndStr, year)+TDE+LF;
}
}
}
else
line+=TD_MSG+NO_DST_MSG+TDE+LF;
}
name=TD_NAME+name+TDE+LF;
line=TR+LF+name+line+TRE+LF;
this.places[i++]=line;
}
}
function output_places()
{var
handle=begin_document("Places"),
d=handle.document;
d.writeln(this.doc_head);
d.writeln(this.page_top);
d.writeln(this.places_header);
d.writeln(this.places.join(LF), TABLE_END);
d.writeln(BR, SP2, this.copyright);
if(this.prog_time)
d.writeln(this.prog_time, BR);
if(this.buttons)
d.writeln(this.buttons);
d.writeln('', LF, '', LF);
d.close();
handle.focus();
if(handle.print&&flags.print===true)
handle.print();
}
function reformat_results(type)
{var title;
if(type=="positions")
title="Positions";
else
title="Times";
flags.print=true;
if(flags.use_font_size===true)
set_font_size();
if(type=="positions")
results=positions_results;
else
results=times_results;
if(flags.use_sep_windows===true&&type=="positions")
results.name="Positions";
else
results.name="Results";
results.build_doc_head(title);
results.add_buttons();
results.output();
flags.print=false;
}
function Result(type)
{this.type=type;
this.doc_head=null;
this.page_top=null;
this.times_header=null;
this.times=new Array();
this.positions_header=null;
this.positions=new Array();
}
Result.prototype.calculate_times=calculate_times;
Result.prototype.calculate_positions=calculate_positions;
Result.prototype.build_doc_head=build_results_doc_head;
Result.prototype.build_page_top=build_page_top;
Result.prototype.add_rs_criteria=add_rs_criteria;
Result.prototype.build_times_header=build_times_header;
Result.prototype.build_positions_header=build_positions_header;
Result.prototype.add_buttons=add_buttons;
Result.prototype.output=output_results;
function PlaceList()
{this.type="places";
this.doc_head=null;
this.page_top=null;
this.places_header=null;
this.places=new Array();
}
PlaceList.prototype.build_doc_head=build_places_doc_head;
PlaceList.prototype.build_places_header=build_places_header;
PlaceList.prototype.build_place_list=build_place_list;
PlaceList.prototype.add_buttons=add_buttons;
PlaceList.prototype.output=output_places;
Number.prototype.showEventTime=function(cdate, loc, type, r_pad)
{var
NO_EVENT=MDASH,
ABOVE="****",
BELOW="-----",
LIGHT = "/////",
DARK="===";
var column=null;
if(type=="T"){if(this<=-F_CIRCUMPOLAR)
column=DARK;
else if(this>=F_CIRCUMPOLAR)
column=LIGHT;
else if(this>F_TEST_NO_EVENT)
column=NO_EVENT;
}
else{if(this<=-F_CIRCUMPOLAR)
column=BELOW;
else if(this>=F_CIRCUMPOLAR)
column=ABOVE;
else if(this>F_TEST_NO_EVENT)
column=NO_EVENT;
}
if(column)
return TD_C+column+TDE;
var
date=new Date(),
time=this,
ut_offset;
date.setTime(cdate.getTime()+time.HtoMsec());
ut_offset=loc.getTimeOffset(date);
if(flags.show_ut===false){time+=ut_offset;
if(time<0||time>24)
column=NO_EVENT;
}
if(column)
return TD_C+column+TDE;
if(time<0)
column=(time+24).formatTime("T")+MINUS;
else if(time>24)
column=(time-24).formatTime("T")+PLUS;
else{column=time.formatTime("T");
if(flags.show_ut===true)
column+=SP2;
}
while(r_pad-->0)
column+=SP;
return TD+column+TDE;
};
Number.prototype.showPosTime=function(date, loc)
{var
time=this,
timeStr,
ut_offset=loc.getTimeOffset(date);
if(flags.show_ut===false)
time=this+ut_offset;
if(time<0)
timeStr=(time+24).formatTime("P")+MINUS;
else if(time>=24){while(time>=24)
time-=24;
timeStr=(time).formatTime("P")+PLUS;
}
else{timeStr=time.formatTime("P");
if(flags.show_ut===true||flags.next_day===true)
timeStr+=SP2;
}
return timeStr;
};
Number.prototype.showMaxAlt=function()
{var column;
if(this<-9.5)
column=MINUS+-Math.round(this);
else if(this<0)
column=MINUS+(-this).toFixed(1);
else if(this<9.5)
column=this.toFixed(1);
else
column=Math.round(this);
return TD+column+SP+TDE;
};
| |