//<!-- 
//     ------------------ Das Script beginnt hier ---------------
//     KOMMENTAR ZUM SCRIPT 
//     Wir freuen uns selbstverstaendlich wenn Sie sich
//     fuer die Details unseres kleinen Skripts interessieren. 
//     Es bringt jedoch nichts, wenn Sie dieses Script auf Ihre Seite kopieren. 
//     Ein einfacher Link beweist genau so gut, dass Sie das Skript gefunden haben.
//     Kopieren Sie deshalb diese Skript nicht auf Ihre 'private' Hompage.
//     Arnold Barmettler, astro!nfo  
//         
//         Source code by Arnold Barmettler, www.astronomie.info / www.CalSKY.com
//         based on algorithms by Peter Duffett-Smith's great and easy book
//         'Practical Astronomy with your Calculator'.
//-->


pi = Math.PI;
DEG = pi/180.0;
RAD = 180./pi;


function sqr(x)
{
	return x*x;
}


// return integer value, closer to 0
function Int(x)
{
	if (x<0) 
	{ 
		return(Math.ceil(x)); 
	} else 
		return(Math.floor(x));
}

function frac(x) 
{ 
	return(x-Math.floor(x)); 
}

function Mod(a, b) 
{ 
	return(a-Math.floor(a/b)*b); 
}


// Modulo PI
function Mod2Pi(x)
{
	x = Mod(x, 2*pi);
	return(x);
}


function round100000(x) 
{ 
	return(Math.round(100000*x)/100000.); 
}

function round10000(x) 
{ 
	return(Math.round(10000*x)/10000.); 
}

function round1000(x) 
{ 
	return(Math.round(1000*x)/1000.); 
}

function round100(x) 
{ 
	return(Math.round(100*x)/100.); 
}

function round10(x) 
{ 
	return(Math.round(10*x)/10.); 
}

var empty = "--";

function HHMM(hh) 
{
	if (hh==0) 
		return(empty);
  
	m = frac(hh)*60;
	var h = Int(hh);
	if (m>=59.5) 
	{ 
		h++; 
		m -=60; 
	}
	m = Math.round(m);
	if (h<10) 
		h = "0"+h;
	h = h+":";
	if (m<10) 
		h = h+"0";
	h = h+m + " Uhr";
	//return(h+" = "+round1000(hh));
	return(h);
}

function HHMMSS(hh) 
{
	if (hh==0) 
		return(empty);
  
	m = frac(hh)*60;
	var h = Int(hh);
	s = frac(m)*60;
	m = Int(m);
	if (s>=59.5) 
	{ 
		m++; 
		s -=60; 
	}
	if (m>=60)   
	{ 
		h++; 
		m -=60; 
	}
	s = Math.round(s);
	if (h<10) 
		h = "0"+h;
	h = h+":";
	if (m<10) 
		h = h+"0";
	h = h+m+":";
	if (s<10) 
		h = h+"0";
	h = h+s;
	return(h+" = "+round10000(hh));
}

function Sign(lon)
{ 
	var signs= new Array("Widder", "Stier", "Zwillinge", "Krebs", "Löwe", "Jungfrau", 
		"Waage", "Skorpion", "Schütze", "Steinbock", "Wassermann", "Fische");
	return( signs[Math.floor(lon*RAD/30)] );
}

// Calculate Julian date: valid only from 1.3.1901 to 28.2.2100
function CalcJD(day,month,year)
{
	jd = 2415020.5-64; // 1.1.1900 - correction of algorithm
	if (month<=2) 
	{ 
		year--; 
		month += 12; 
	}
	jd += Int( (year-1900)*365.25 );
	jd += Int( 30.6001*(1+month) );
	return(jd + day);
}

// Julian Date to Greenwich Mean Sidereal Time
function GMST(JD)
{
	var UT = frac(JD-0.5)*24; // UT in hours
	JD = Math.floor(JD-0.5)+0.5;   // JD at 0 hours UT
	var T = (JD-2451545.0)/36525.0;
	T0 = 6.697374558 + T*(2400.051336 + T*0.000025862);
	return(Mod(T0+UT*1.002737909, 24));
}

// Convert Greenweek mean sidereal time to UT
function GMST2UT(JD, gmst)
{
	JD = Math.floor(JD-0.5)+0.5;   // JD at 0 hours UT
	var T = (JD-2451545.0)/36525.0;
	var T0 = Mod(6.697374558 + T*(2400.051336 + T*0.000025862), 24);
	//var UT = 0.9972695663*Mod((gmst-T0), 24.);
	var UT = 0.9972695663*((gmst-T0));
	return(UT);
}

// Local Mean Sidereal Time, geographical longitude in radians, East is positive
function GMST2LMST(gmst, lon)
{
	var lmst = Mod(gmst+RAD*lon/15, 24.);
	return( lmst );
}

// Transform ecliptical coordinates (lon/lat) to equatorial coordinates (RA/dec)
function Ecl2Equ(coor, TDT)
{
	T = (TDT-2451545.0)/36525.; // Epoch 2000 January 1.5
	eps = (23.+(26+21.45/60)/60 + T*(-46.815 +T*(-0.0006 + T*0.00181) )/3600 )*DEG;
	var coseps = Math.cos(eps);
	var sineps = Math.sin(eps);
  
	var sinlon = Math.sin(coor.lon);
	coor.ra  = Mod2Pi( Math.atan2( (sinlon*coseps-Math.tan(coor.lat)*sineps), Math.cos(coor.lon) ) );
	coor.dec = Math.asin( Math.sin(coor.lat)*coseps + Math.cos(coor.lat)*sineps*sinlon );
  
	return coor;
}

// Calculate coordinates for Sun
// Coordinates are accurate to about 10s (right ascension) 
// and a few minutes of arc (declination)
function SunPosition(TDT)
{
	var D = TDT-2447891.5;
  
	var eg = 279.403303*DEG;
	var wg = 282.768422*DEG;
	var e  = 0.016713;
	var a  = 149598500; // km
	var diameter0 = 0.533128*DEG; // angular diameter of Moon at a distance
  
	var MSun = 360*DEG/365.242191*D+eg-wg;
	var nu = MSun + 360.*DEG/pi*e*Math.sin(MSun);
  
	var sunCoor = new Object();
	sunCoor.lon =  Mod2Pi(nu+wg);
	sunCoor.lat = 0;
	sunCoor.anomalyMean = MSun;
  
	sunCoor.distance = (1-sqr(e))/(1+e*Math.cos(nu));   // distance in astronomical units
	sunCoor.diameter = diameter0/sunCoor.distance; // angular diameter in radians
	sunCoor.distance *= a;                                                  // distance in km
	sunCoor.parallax = 6378.137/sunCoor.distance;  // horizonal parallax

	sunCoor = Ecl2Equ(sunCoor, TDT);
  
	sunCoor.sign = Sign(sunCoor.lon);
	return sunCoor;
}

// Calculate data and coordinates for the Moon
// Coordinates are accurate to about 1/5 degree (in ecliptic coordinates)
function MoonPosition(sunCoor, TDT)
{
	var D = TDT-2447891.5;
  
	// Mean Moon orbit elements as of 1990.0
	var l0 = 318.351648*DEG;
	var P0 =  36.340410*DEG;
	var N0 = 318.510107*DEG;
	var i  = 5.145396*DEG;
	var e  = 0.054900;
	var a  = 384401; // km
	var diameter0 = 0.5181*DEG; // angular diameter of Moon at a distance
	var parallax0 = 0.9507*DEG; // parallax at distance a
  
	var l = 13.1763966*DEG*D+l0;
	var MMoon = l-0.1114041*DEG*D-P0; // Moon's mean anomaly M
	var N = N0-0.0529539*DEG*D;       // Moon's mean ascending node longitude
	var C = l-sunCoor.lon;
	var Ev = 1.2739*DEG*Math.sin(2*C-MMoon);
	var Ae = 0.1858*DEG*Math.sin(sunCoor.anomalyMean);
	var A3 = 0.37*DEG*Math.sin(sunCoor.anomalyMean);
	var MMoon2 = MMoon+Ev-Ae-A3;  // corrected Moon anomaly
	var Ec = 6.2886*DEG*Math.sin(MMoon2);  // equation of centre
	var A4 = 0.214*DEG*Math.sin(2*MMoon2);
	var l2 = l+Ev+Ec-Ae+A4; // corrected Moon's longitude
	var V = 0.6583*DEG*Math.sin(2*(l2-sunCoor.lon));
	var l3 = l2+V; // true orbital longitude;

	var N2 = N-0.16*DEG*Math.sin(sunCoor.anomalyMean);
  
	var moonCoor = new Object();  
	moonCoor.lon = Mod2Pi( N2 + Math.atan2( Math.sin(l3-N2)*Math.cos(i), Math.cos(l3-N2) ) );
	moonCoor.lat = Math.asin( Math.sin(l3-N2)*Math.sin(i) );
	moonCoor.orbitLon = l3;
  
	moonCoor = Ecl2Equ(moonCoor, TDT);
	// relative distance to semi mayor axis of lunar oribt
	moonCoor.distance = (1-sqr(e)) / (1+e*Math.cos(MMoon2+Ec) );
	moonCoor.diameter = diameter0/moonCoor.distance; // angular diameter in radians
	moonCoor.parallax = parallax0/moonCoor.distance; // horizontal parallax in radians
	moonCoor.distance *= a;       // distance in km

	// Age of Moon in radians since New Moon (0) - Full Moon (pi)
	moonCoor.moonAge = Mod2Pi(l3-sunCoor.lon);   
	moonCoor.phase   = 0.5*(1-Math.cos(moonCoor.moonAge)); // Moon phase, 0-1
  
	var phases = new Array("Neumond", "Zunehmende Sichel", "Erstes Viertel", "Zunnehmender Mond", 
		"Vollmond", "Abnehmender Mond", "Letztes Viertel", "Abnehmende Sichel", "Neumond");
	var bilder = new Array("/dpl-full-theme/images/dpl/layout/neumond.gif", "/dpl-full-theme/images/dpl/layout/mond_erstesviertel.gif", "/dpl-full-theme/images/dpl/layout/mond_erstesviertel.gif", "/dpl-full-theme/images/dpl/layout/mond_erstesviertel.gif", 
		"/dpl-full-theme/images/dpl/layout/vollmond.gif", "/dpl-full-theme/images/dpl/layout/mond_letztesviertel.gif", "/dpl-full-theme/images/dpl/layout/mond_letztesviertel.gif", "/dpl-full-theme/images/dpl/layout/mond_letztesviertel.gif", "/dpl-full-theme/images/dpl/layout/neumond.gif");
	var mainPhase = 1./29.53*360*DEG; // show 'Newmoon, 'Quarter' for +/-1 day arond the actual event
	var p = Mod(moonCoor.moonAge, 90.*DEG);
	if (p < mainPhase || p > 90*DEG-mainPhase) 
		p = 2*Math.round(moonCoor.moonAge / (90.*DEG));
	else 
		p = 2*Math.floor(moonCoor.moonAge / (90.*DEG))+1;
	moonCoor.moonPhase = phases[p];
	moonCoor.moonPic = bilder[p];
  
	moonCoor.sign = Sign(moonCoor.lon);

	return(moonCoor);
}

// returns Greenwich sidereal time (hours) of time of rise 
// and set of object with coordinates coor.ra/coor.dec
// at geographic position lon/lat (all values in radians)
// for mathematical horizon, without refraction and body radius correction
function GMSTRiseSet(coor, lon, lat)
{
	var riseset = new Object();	
	var tagbogen = Math.acos(-Math.tan(lat)*Math.tan(coor.dec));
 
	riseset.transit =     RAD/15*(         +coor.ra-lon);
	riseset.rise    = 24.+RAD/15*(-tagbogen+coor.ra-lon); // calculate GMST of rise of object
	riseset.set     =     RAD/15*(+tagbogen+coor.ra-lon); // calculate GMST of set of object

	riseset.transit = Mod(riseset.transit, 24);
	riseset.rise    = Mod(riseset.rise, 24);
	riseset.set     = Mod(riseset.set, 24);
 
	return(riseset);
}

// Find GMST of rise/set of object from the two calculates 
// (start)points (day 1 and 2) and at midnight UT(0)
function InterpolateGMST(gmst0, gmst1, gmst2, timefactor)
{
 	return( (timefactor*24.07*gmst1- gmst0*(gmst2-gmst1)) / (timefactor*24.07+gmst1-gmst2) );
}

// JD is the Julian Date of 0h UTC time (midnight)
function RiseSet(jd0UT, coor1, coor2, lon, lat, timeinterval)
{
	var rise1 = GMSTRiseSet(coor1, lon, lat);
	var rise2 = GMSTRiseSet(coor2, lon, lat);
  
	var rise = new Object();
  
	//alert( rise1.set  +"  "+ rise2.set );
	// unwrap GMST in case we move across 24h -> 0h
	if (rise1.transit > rise2.transit && Math.abs(rise1.transit-rise2.transit)>18) 
		rise2.transit += 24;
	if (rise1.rise    > rise2.rise    && Math.abs(rise1.rise   -rise2.rise)>18)    
		rise2.rise += 24;
	if (rise1.set     > rise2.set     && Math.abs(rise1.set    -rise2.set)>18)     
		rise2.set  += 24;
	var T0 = GMST(jd0UT);
	//  var T02 = T0-zone*1.002738; // Greenwich sidereal time at 0h time zone (zone: hours)

	// Greenwich sidereal time for 0h at selected longitude
	var T02 = T0-lon*RAD/15*1.002738; if (T02 < 0) T02 += 24; 

	if (rise1.transit < T02) 
	{ 
		rise1.transit += 24; 
		rise2.transit += 24; 
	}
	if (rise1.rise    < T02) 
	{ 
		rise1.rise    += 24; 
		rise2.rise    += 24; 
	}
	if (rise1.set     < T02) 
	{ 
		rise1.set     += 24; 
		rise2.set     += 24; 
	}
  
	//alert("after="+ rise1.set  +"  "+ rise2.set+ " T0="+ T0 );
 
	// Refraction and Parallax correction
  
	var decMean = 0.5*(coor1.dec+coor2.dec);
	var psi = Math.acos(Math.sin(lat)/Math.cos(decMean));
	// altitude of sun center: semi-diameter, horizontal parallax and (standard) refraction of 34'
	var alt = 0.5*coor1.diameter-coor1.parallax+34./60*DEG;
	var y = Math.asin(Math.sin(alt)/Math.sin(psi));
	var dt = 240*RAD*y/Math.cos(decMean)/3600; // time correction due to refraction, parallax

	rise.transit = GMST2UT( jd0UT, InterpolateGMST( T0, rise1.transit, rise2.transit, timeinterval) );
	rise.rise    = GMST2UT( jd0UT, InterpolateGMST( T0, rise1.rise,    rise2.rise,    timeinterval) -dt );
	rise.set     = GMST2UT( jd0UT, InterpolateGMST( T0, rise1.set,     rise2.set,     timeinterval) +dt );
  
	//rise.transit = Mod(rise.transit, 24.);
	//rise.rise    = Mod(rise.rise, 24.);
	//rise.set     = Mod(rise.set,  24.);
 
	return(rise);  
}

// Find (local) time of sunrise and sunset
// JD is the Julian Date of 0h local time (midnight)
// Accurate to about 1-2 minutes
// recursive: 1 - calculate rise/set in UTC
// recursive: 0 - find rise/set on the current local day (set could also be first)
function SunRise(JD, deltaT, lon, lat, zone, recursive)
{
	var jd0UT = Math.floor(JD-0.5)+0.5;   // JD at 0 hours UT
	//  alert("jd0UT="+jd0UT+"  JD="+JD);
	var coor1 = SunPosition(jd0UT+0*deltaT/24/3600);
	var coor2 = SunPosition(jd0UT+1+0*deltaT/24/3600); // calculations for next day's UTC midnight
  
	var risetemp = new Object();
	var rise = new Object();
	rise = RiseSet(jd0UT, coor1, coor2, lon, lat, 1 ); // rise/set time in UTC
	if (!recursive) 
	{ // check and adjust to have rise/set time on local calendar day
		if (zone>0) 
		{
			// rise time was yesterday local time -> calculate rise time for next UTC day
			if (rise.rise>=24-zone || rise.transit>=24-zone || rise.set>=24-zone) 
			{
				risetemp = SunRise(JD+1, deltaT, lon, lat, zone, 1);
				if (rise.rise>=24-zone) 
					rise.rise = risetemp.rise;
				if (rise.transit >=24-zone) 
					rise.transit = risetemp.transit;
				if (rise.set >=24-zone) 
					rise.set  = risetemp.set;
			}
		}
		else if (zone<0) 
		{
			// rise time was yesterday local time -> calculate rise time for next UTC day
			if (rise.rise<-zone || rise.transit<-zone || rise.set<-zone) 
			{
				risetemp = SunRise(JD-1, deltaT, lon, lat, zone, 1);
				if (rise.rise<-zone) 
					rise.rise = risetemp.rise;
				if (rise.transit<-zone) 
					rise.transit = risetemp.transit;
				if (rise.set <-zone) 
					rise.set  = risetemp.set;
			}
		}
		rise.transit = Mod(rise.transit+zone, 24);
		rise.rise    = Mod(rise.rise   +zone, 24);
		rise.set     = Mod(rise.set    +zone, 24);
	}
	return( rise );  
}

// Find local time of moonrise and moonset
// JD is the Julian Date of 0h local time (midnight)
// Accurate to about 5 minutes or better
// recursive: 1 - calculate rise/set in UTC
// recursive: 0 - find rise/set on the current local day (set could also be first)
// returns '' for moonrise/set does not occur on selected day
function MoonRise(JD, deltaT, lon, lat, zone, recursive)
{
	var timeinterval = 0.5;
  
	var jd0UT = Math.floor(JD-0.5)+0.5;   // JD at 0 hours UT
	var suncoor1 = SunPosition(jd0UT+0*deltaT/24/3600);
	var coor1 = MoonPosition(suncoor1, jd0UT+0*deltaT/24/3600);

	var suncoor2 = SunPosition(jd0UT +timeinterval +0*deltaT/24/3600); // calculations for noon
	// calculations for next day's midnight
	var coor2 = MoonPosition(suncoor2, jd0UT +timeinterval +0*deltaT/24/3600); 
  
	var risetemp = new Object();
	var rise = new Object();
  
	// rise/set time in UTC, time zone corrected later
	rise = RiseSet(jd0UT, coor1, coor2, lon, lat, timeinterval); 
  
	if (!recursive) 
	{ // check and adjust to have rise/set time on local calendar day
		if (zone>0) 
		{
			// recursive call to MoonRise returns events in UTC
			riseprev = MoonRise(JD-1, deltaT, lon, lat, zone, 1); 
			// recursive call to MoonRise returns events in UTC
			//risenext = MoonRise(JD+1, deltaT, lon, lat, zone, 1);
			//alert("yesterday="+riseprev.transit+"  today="+rise.transit+" tomorrow="+risenext.transit);
			//alert("yesterday="+riseprev.rise+"  today="+rise.rise+" tomorrow="+risenext.rise);
			//alert("yesterday="+riseprev.set+"  today="+rise.set+" tomorrow="+risenext.set);

			if (rise.transit >= 24-zone || rise.transit < -zone) 
			{ // transit time is tomorrow local time
				if (riseprev.transit < 24-zone) 
					rise.transit = ''; // there is no moontransit today
				else 
					rise.transit  = riseprev.transit;
			}

			if (rise.rise >= 24-zone || rise.rise < -zone) 
			{ // transit time is tomorrow local time
				if (riseprev.rise < 24-zone) 
					rise.rise = ''; // there is no moontransit today
				else 
					rise.rise  = riseprev.rise;
			}

			if (rise.set >= 24-zone || rise.set < -zone) 
			{ // transit time is tomorrow local time
				if (riseprev.set < 24-zone) 
					rise.set = ''; // there is no moontransit today
				else 
					rise.set  = riseprev.set;
			}
		}
		else if (zone<0) 
		{
			// rise/set time was tomorrow local time -> calculate rise time for former UTC day
			if (rise.rise<-zone || rise.set<-zone || rise.transit<-zone) 
			{ 
				risetemp = MoonRise(JD+1, deltaT, lon, lat, zone, 1);
                
				if (rise.rise < -zone) 
				{
					if (risetemp.rise > -zone) 
						rise.rise = ''; // there is no moonrise today
					else 
						rise.rise = risetemp.rise;
				}
                
				if (rise.transit < -zone)
				{
					if (risetemp.transit > -zone)  
						rise.transit = ''; // there is no moonset today
					else 
						rise.transit  = risetemp.transit;
				}
                
				if (rise.set < -zone)
				{
					if (risetemp.set > -zone)  
						rise.set = ''; // there is no moonset today
					else 
						rise.set  = risetemp.set;
				}
                
			}
		}
        
		if (rise.rise)    
			rise.rise = Mod(rise.rise+zone, 24);    // correct for time zone, if time is valid
		if (rise.transit) 
			rise.transit  = Mod(rise.transit +zone, 24); // correct for time zone, if time is valid
		if (rise.set)     
			rise.set  = Mod(rise.set +zone, 24);    // correct for time zone, if time is valid
	}
	return( rise );  
}

function Compute(form)
{

	if (eval(form.Year.value)<=1900 || eval(form.Year.value)>=2100 ) 
	{
		alert("Dies Script erlaubt nur Berechnungen"+
			"in der Zeitperiode 1901-2099. Angezeigte Resultat sind ungültig.");
		return;
	}
	JD0 = CalcJD( eval(form.Day.value), eval(form.Month.value), eval(form.Year.value) );
	JD  = JD0
		+( eval(form.Hour.value) -eval(form.Zone.value.replace(/,/,'.')) +eval(form.Minute.value)/60. 
		+ eval(form.Second.value.replace(/,/,'.'))/3600) /24.;
	TDT = JD+eval(form.DeltaT.value.replace(/,/,'.'))/24./3600.;

	lat      = eval(form.Lat.value.replace(/,/,'.'))*DEG;
	lon      = eval(form.Lon.value.replace(/,/,'.'))*DEG;

	sunCoor  = SunPosition(TDT);                    // Calculate data for the Sun at given time
	moonCoor = MoonPosition(sunCoor, TDT);  // Calculate data for the Moon at given time
        
	form.JD.value = round100000(JD);
	var gmst = GMST(JD);
	form.GMST.value = HHMMSS(gmst);
	form.LMST.value = HHMMSS(GMST2LMST(gmst, lon));
        
	if (eval(form.Minute.value)<10) 
		form.Minute.value = "0"+eval(form.Minute.value);
	if (eval(form.Month.value)<10) 
		form.Month.value = "0"+eval(form.Month.value);
        
	form.SunLon.value  = round1000(sunCoor.lon*RAD);
	form.SunRA.value   = HHMM(sunCoor.ra*RAD/15);
	form.SunDec.value  = round1000(sunCoor.dec*RAD);

	form.SunDistance.value = round10(sunCoor.distance);
	form.SunDiameter.value = round100(sunCoor.diameter*RAD*60); // angular diameter in arc seconds

	form.SunSign.value = sunCoor.sign;
        
	// JD0: JD of 0h UTC time
	sunRise = SunRise(JD0, eval(form.DeltaT.value.replace(/,/,'.')), lon, lat, eval(form.Zone.value.replace(/,/,'.')), 0);

	form.SunTransit.value = HHMM(sunRise.transit);
	form.SunRise.value    = HHMM(sunRise.rise);
	form.SunSet.value     = HHMM(sunRise.set);
        
	form.MoonLon.value = round1000(moonCoor.lon*RAD);
	form.MoonLat.value = round1000(moonCoor.lat*RAD);
	form.MoonRA.value  = HHMM(moonCoor.ra*RAD/15);
	form.MoonDec.value = round1000(moonCoor.dec*RAD);
	form.MoonAge.value = round1000(moonCoor.moonAge*RAD);
	form.MoonPhaseNumber.value = round1000(moonCoor.phase);
	form.MoonPhase.value = moonCoor.moonPhase;
        
	form.MoonDistance.value = round10(moonCoor.distance);
	form.MoonDiameter.value = round100(moonCoor.diameter*RAD*60); // angular diameter in arc seconds
        
	form.MoonSign.value = moonCoor.sign;
        
	moonRise = MoonRise(JD0, eval(form.DeltaT.value.replace(/,/,'.')), lon, lat, eval(form.Zone.value.replace(/,/,'.')), 0);

	form.MoonTransit.value = HHMM(moonRise.transit);
	form.MoonRise.value    = HHMM(moonRise.rise);
	form.MoonSet.value     = HHMM(moonRise.set);
}

var Lon_value = "8.366003";
var Lat_value = "51.654228";
var DeltaT_value = "65";
var Hour_value;
var Minute_value;
var Second_value;
var Day_value;
var Month_value;
var Year_value;
var Zone_value;

function InitPatriot(lon, lat)
{       
	var now=new Date();
	Hour_value    = now.getHours();
	Minute_value  = now.getMinutes();
	if (Minute_value<10) 
		Minute_value = "0"+Minute_value;
	Second_value  = now.getSeconds();
	if (Second_value<10) 
		Second_value = "0"+Second_value;
	Day_value     = now.getDate();
	Month_value   = now.getMonth()+1;
	if (Month_value<10) 
		Month_value = "0"+Month_value;
	if (now.getYear()<1900) 
		Year_value = now.getYear()+1900; // MSIE returns 2004, but NS years since 1900
	else 
		Year_value    = now.getYear();
	Zone_value    = -now.getTimezoneOffset()/60;

	DeltaT_value  = "65"; // deltaT - difference among 'earth center' versus 'observered' time (TDT-UT), in seconds
	Lon_value     = lon;  
	Lat_value     = lat;
	document.getElementById("SunRise").innerHTML = "";
	document.getElementById("SunSet").innerHTML = "";
	document.getElementById("MoonRise").innerHTML = "";
	document.getElementById("MoonSet").innerHTML = "";
	document.getElementById("MoonPhase").innerHTML = "";
}

function ComputePatriot()
{
	JD0 = CalcJD(eval(Day_value), eval(Month_value), eval(Year_value));
	JD  = JD0 + (eval(Hour_value) - eval(Zone_value) + eval(Minute_value)/60 + eval(Second_value)/3600) / 24.;

	TDT = JD + eval(DeltaT_value)/24./ 3600.;

	lat      = Lat_value*DEG;
	lon      = Lon_value*DEG;

	sunCoor  = SunPosition(TDT);                    // Calculate data for the Sun at given time
	moonCoor = MoonPosition(sunCoor, TDT);  // Calculate data for the Moon at given time
        
	//form.JD.value = round100000(JD);
	var gmst = GMST(JD);
	//form.GMST.value = HHMMSS(gmst);
	//form.LMST.value = HHMMSS(GMST2LMST(gmst, lon));
        
	//if (eval(form.Minute.value)<10) 
	//	form.Minute.value = "0"+eval(form.Minute.value);
	//if (eval(form.Month.value)<10) 
	//	form.Month.value = "0"+eval(form.Month.value);
        
	//form.SunLon.value  = round1000(sunCoor.lon*RAD);
	//form.SunRA.value   = HHMM(sunCoor.ra*RAD/15);
	//form.SunDec.value  = round1000(sunCoor.dec*RAD);

	//form.SunDistance.value = round10(sunCoor.distance);
	//form.SunDiameter.value = round100(sunCoor.diameter*RAD*60); // angular diameter in arc seconds

	//form.SunSign.value = sunCoor.sign;
        
	// JD0: JD of 0h UTC time
	sunRise = SunRise(JD0, DeltaT_value, lon, lat, Zone_value, 0);

	//form.SunTransit.value = HHMM(sunRise.transit);
	document.getElementById("SunRise").innerHTML = HHMM(sunRise.rise);
	document.getElementById("SunSet").innerHTML = HHMM(sunRise.set);
        
	//form.MoonLon.value = round1000(moonCoor.lon*RAD);
	//form.MoonLat.value = round1000(moonCoor.lat*RAD);
	//form.MoonRA.value  = HHMM(moonCoor.ra*RAD/15);
	//form.MoonDec.value = round1000(moonCoor.dec*RAD);
	//form.MoonAge.value = round1000(moonCoor.moonAge*RAD);
	//form.MoonPhaseNumber.value = round1000(moonCoor.phase);
	document.getElementById("MoonPhase").innerHTML = moonCoor.moonPhase;
	document.getElementById("MoonPhase_b_bild").src = moonCoor.moonPic;
	//form.MoonDistance.value = round10(moonCoor.distance);
	//form.MoonDiameter.value = round100(moonCoor.diameter*RAD*60); // angular diameter in arc seconds
        
	//form.MoonSign.value = moonCoor.sign;
        
	moonRise = MoonRise(JD0, DeltaT_value, lon, lat, Zone_value, 0);

	//form.MoonTransit.value = HHMM(moonRise.transit);
	 document.getElementById("MoonRise").innerHTML = HHMM(moonRise.rise);
	 document.getElementById("MoonSet").innerHTML = HHMM(moonRise.set);
}