예제 #1
0
double CSunRiseSet::calcSunsetGMT(int iJulDay, double dLatitude, double dLongitude) 
{ 
		// First calculates sunrise and approx length of day 
 
		double dGamma = CalcGamma(iJulDay + 1); 
		double eqTime = CalcEqofTime(dGamma); 
		double solarDec = CalcSolarDec(dGamma); 
		double hourAngle = CalcHourAngle(dLatitude, solarDec, 0); 
		double delta = dLongitude - dRadToDeg(hourAngle); 
		double timeDiff = 4 * delta; 
		double setTimeGMT = 720 + timeDiff - eqTime; 
 
		// first pass used to include fractional day in gamma calc 
 
		double gamma_sunset = CalcGamma2(iJulDay, setTimeGMT/60); 
		eqTime = CalcEqofTime(gamma_sunset); 
 
 
		solarDec = CalcSolarDec(gamma_sunset); 
		 
		hourAngle = CalcHourAngle(dLatitude, solarDec, false); 
		delta = dLongitude - dRadToDeg(hourAngle); 
		timeDiff = 4 * delta; 
		setTimeGMT = 720 + timeDiff - eqTime; // in minutes 
 
		return setTimeGMT; 
	} 
예제 #2
0
파일: positioner.c 프로젝트: nvertigo/vdr
void cPositioner::StartMovementTimer(int Longitude)
{
  if (Setup.PositionerSpeed <= 0)
     return;
  cMutexLock MutexLock(&mutex);
  lastLongitude = CurrentLongitude(); // in case the dish was already in motion
  targetLongitude = Longitude;
  lastHourAngle = CalcHourAngle(lastLongitude);
  targetHourAngle = CalcHourAngle(targetLongitude);
  swingTime = abs(targetHourAngle - lastHourAngle) * 1000 / Setup.PositionerSpeed; // time (ms) it takes to move the dish from lastHourAngle to targetHourAngle
  movementStart.Set();
  Setup.PositionerLastLon = targetLongitude;
}
예제 #3
0
파일: positioner.c 프로젝트: nvertigo/vdr
cPositioner::cPositioner(void)
{
  capabilities = pcCanNothing;
  frontend = -1;
  targetLongitude = lastLongitude = Setup.PositionerLastLon;
  targetHourAngle = lastHourAngle = CalcHourAngle(lastLongitude);
  swingTime = 0;
  delete positioner;
  positioner = this;
}
예제 #4
0
double CSunRiseSet::calcSunriseGMT(int iJulDay, double dLatitude, double dLongitude) 
{ 
	// *** First pass to approximate sunrise 
 
	double gamma = CalcGamma(iJulDay); 
	double eqTime = CalcEqofTime(gamma); 
	double solarDec = CalcSolarDec(gamma); 
	double hourAngle = CalcHourAngle(dLatitude, solarDec, true); 
	double delta = dLongitude - dRadToDeg(hourAngle); 
	double timeDiff = 4 * delta; 
	double timeGMT = 720 + timeDiff - eqTime; 
 
	// *** Second pass includes fractional jday in gamma calc 
 
	double gamma_sunrise = CalcGamma2(iJulDay, timeGMT/60); 
	eqTime = CalcEqofTime(gamma_sunrise); 
	solarDec = CalcSolarDec(gamma_sunrise); 
	hourAngle = CalcHourAngle(dLatitude, solarDec, 1); 
	delta = dLongitude - dRadToDeg(hourAngle); 
	timeDiff = 4 * delta; 
	timeGMT = 720 + timeDiff - eqTime; // in minutes 
 
	return timeGMT; 
}