public void Equ2Hor(double tjd, double deltat, double x, double y, ref SiteInfo location, double ra, double dec, RefractionOption ref_option, ref double zd, ref double az, ref double rar, ref double decr)
{
NOVAS2.Equ2Hor(tjd, deltat, x, y, ref location, ra, dec, ref_option, ref zd, ref az, ref rar, ref decr);
}
public PositionVector GetTopocentricPosition(double tjd, Site site, bool Refract)
{
CatEntry star = new CatEntry();
SiteInfo siteInfo = new SiteInfo();
double[] earthvector1 = new double[4];
double[] pos2_1 = new double[4];
double[] vel1 = new double[4];
double[] pos2_2 = new double[4];
double[] earthvector2 = new double[4];
double[] pos1 = new double[4];
double[] pos2_3 = new double[4];
double[] pos2_4 = new double[4];
double[] pos2_5 = new double[4];
double[] pos2_6 = new double[4];
double[] pos2_7 = new double[4];
double[] vel2 = new double[4];
double[] pos2_8 = new double[4];
double[] bary_earthp = new double[4];
double[] bary_earthv = new double[4];
double[] helio_earthp = new double[4];
double[] helio_earthv = new double[4];
double[] pos2 = new double[4];
if (!(this.m_rav & this.m_decv))
{
throw new ASCOM.Astrometry.Exceptions.ValueNotSetException("Star.GetTopocentricPosition RA or DEC not available");
}
double jd_high = !this.m_bDTValid ? tjd - DeltatCode.DeltaTCalc(tjd) / 86400.0 : tjd - this.m_deltat;
try
{
siteInfo.Latitude = site.Latitude;
}
catch (Exception ex)
{
//ProjectData.SetProjectError(ex);
throw new ValueNotAvailableException("Star:GetTopocentricPosition Site.Latitude is not available");
}
try
{
siteInfo.Longitude = site.Longitude;
}
catch (Exception ex)
{
//ProjectData.SetProjectError(ex);
throw new ValueNotAvailableException("Star:GetTopocentricPosition Site.Longitude is not available");
}
try
{
siteInfo.Height = site.Height;
}
catch (Exception ex)
{
//ProjectData.SetProjectError(ex);
throw new ValueNotAvailableException("Star:GetTopocentricPosition Site.Height is not available");
}
double tdb = 0;
this.hr = NOVAS2.GetEarth(tjd, ref this.m_earth, ref tdb, ref bary_earthp, ref bary_earthv, ref helio_earthp, ref helio_earthv);
if ((int)this.hr > 0)
{
pos2[0] = 0.0;
pos2[1] = 0.0;
pos2[2] = 0.0;
throw new NOVASFunctionException("Star.GetApparentPosition", "get_earth", this.hr);
}
double mobl = 0;
double tobl = 0;
double eq = 0;
double dpsi = 0;
double deps = 0;
NOVAS2.EarthTilt(tdb, ref mobl, ref tobl, ref eq, ref dpsi, ref deps);
double gst = 0;
NOVAS2.SiderealTime(jd_high, 0.0, eq, ref gst);
NOVAS2.Terra(ref siteInfo, gst, ref pos1, ref vel2);
int num1 = (int)NOVAS2.Nutate(tdb, NutationDirection.TrueToMean, pos1, ref pos2_3);
NOVAS2.Precession(tdb, pos2_3, 2451545.0, ref pos2_1);
int num2 = (int)NOVAS2.Nutate(tdb, NutationDirection.TrueToMean, vel2, ref pos2_8);
NOVAS2.Precession(tdb, pos2_8, 2451545.0, ref pos2_2);
int index = 0;
do
{
earthvector1[index] = bary_earthp[index] + pos2_1[index];
vel1[index] = bary_earthv[index] + pos2_2[index];
earthvector2[index] = helio_earthp[index] + pos2_1[index];
checked { ++index; }
}while (index <= 2);
star.RA = this.m_ra;
star.Dec = this.m_dec;
star.ProMoRA = this.m_pmra;
star.ProMoDec = this.m_pmdec;
star.Parallax = this.m_plx;
star.RadialVelocity = this.m_rv;
NOVAS2.StarVectors(star, ref pos1, ref vel2);
NOVAS2.ProperMotion(2451545.0, pos1, vel2, tdb, ref pos2_3);
double lighttime = 0;
NOVAS2.BaryToGeo(pos2_3, earthvector1, ref pos2_4, ref lighttime);
int num3 = (int)NOVAS2.SunField(pos2_4, earthvector2, ref pos2_5);
int num4 = (int)NOVAS2.Aberration(pos2_5, vel1, lighttime, ref pos2_6);
NOVAS2.Precession(2451545.0, pos2_6, tdb, ref pos2_7);
int num5 = (int)NOVAS2.Nutate(tdb, NutationDirection.MeanToTrue, pos2_7, ref pos2);
double ra = 0;
double dec = 0;
int num6 = (int)NOVAS2.Vector2RADec(pos2, ref ra, ref dec);
double num7 = Math.Sqrt(Math.Pow(pos2[0], 2.0) + Math.Pow(pos2[1], 2.0) + Math.Pow(pos2[2], 2.0));
RefractionOption ref_option = RefractionOption.NoRefraction;
if (Refract)
{
bool flag = true;
try
{
siteInfo.Temperature = site.Temperature;
}
catch (Exception ex)
{
//ProjectData.SetProjectError(ex);
flag = false;
//ProjectData.ClearProjectError();
}
try
{
siteInfo.Pressure = site.Pressure;
}
catch (Exception ex)
{
//ProjectData.SetProjectError(ex);
flag = false;
//ProjectData.ClearProjectError();
}
ref_option = !flag ? RefractionOption.StandardRefraction : RefractionOption.LocationRefraction;
}
double zd = 0;
double az = 0;
double rar = 0;
double decr = 0;
if (this.m_bDTValid)
{
NOVAS2.Equ2Hor(tjd, this.m_deltat, 0.0, 0.0, ref siteInfo, ra, dec, ref_option, ref zd, ref az, ref rar, ref decr);
}
else
{
NOVAS2.Equ2Hor(tjd, DeltatCode.DeltaTCalc(tjd), 0.0, 0.0, ref siteInfo, ra, dec, ref_option, ref zd, ref az, ref rar, ref decr);
}
if (ref_option > RefractionOption.NoRefraction)
{
NOVAS2.RADec2Vector(rar, decr, num7, ref pos2);
}
return(new PositionVector(pos2[0], pos2[1], pos2[2], rar, decr, num7, num7 / 173.14463348, az, 90.0 - zd));
}
public PositionVector GetTopocentricPosition(double tjd, Site site, bool Refract)
{
double[] pos1 = new double[4];
double[] pos2_1 = new double[4];
double[] pos2_2 = new double[4];
double[] pos2_3 = new double[4];
double[] pos2_4 = new double[4];
double[] vel1 = new double[4];
double[] pos2_5 = new double[4];
double[] pos2_6 = new double[4];
double[] pos2_7 = new double[4];
double[] earthvector1 = new double[4];
double[] pos2 = new double[4];
double[] vel2 = new double[4];
double[] earthvector2 = new double[4];
double[] peb = new double[4];
double[] veb = new double[4];
double[] pes = new double[4];
double[] ves = new double[4];
Object3 SsBody = new Object3();
OnSurface onSurface = new OnSurface();
if (!this.m_bDTValid)
{
this.m_deltat = DeltatCode.DeltaTCalc(tjd);
}
double num1 = tjd - this.m_deltat / 86400.0;
SiteInfo siteInfo = new SiteInfo();
try
{
siteInfo.Latitude = site.Latitude;
}
catch (Exception ex)
{
//ProjectData.SetProjectError(ex);
throw new ValueNotAvailableException("Star:GetTopocentricPosition Site.Latitude is not available");
}
try
{
siteInfo.Longitude = site.Longitude;
}
catch (Exception ex)
{
//ProjectData.SetProjectError(ex);
throw new ValueNotAvailableException("Star:GetTopocentricPosition Site.Longitude is not available");
}
try
{
siteInfo.Height = site.Height;
}
catch (Exception ex)
{
//ProjectData.SetProjectError(ex);
throw new ValueNotAvailableException("Star:GetTopocentricPosition Site.Height is not available");
}
PositionVector positionVector;
if (this.m_type == BodyType.Moon & (this.m_number == 10 | this.m_number == 11))
{
onSurface.Height = site.Height;
onSurface.Latitude = site.Latitude;
onSurface.Longitude = site.Longitude;
onSurface.Pressure = site.Pressure;
onSurface.Temperature = site.Temperature;
SsBody.Number = CommonCode.NumberToBody(this.m_number);
SsBody.Type = ASCOM.Astrometry.ObjectType.MajorPlanetSunOrMoon;
RefractionOption RefOption = !Refract ? RefractionOption.NoRefraction : RefractionOption.LocationRefraction;
double Ra = 0;
double Dec = 0;
double Dis = 0;
this.Nov31.TopoPlanet(tjd, SsBody, this.m_deltat, onSurface, Accuracy.Full, ref Ra, ref Dec, ref Dis);
double Zd = 0;
double Az = 0;
double RaR = 0;
double DecR = 0;
this.Nov31.Equ2Hor(num1, this.m_deltat, Accuracy.Full, 0.0, 0.0, onSurface, Ra, Dec, RefOption, ref Zd, ref Az, ref RaR, ref DecR);
this.Nov31.RaDec2Vector(RaR, DecR, Dis, ref pos2);
positionVector = new PositionVector(pos2[0], pos2[1], pos2[2], RaR, DecR, Dis, Dis / 173.14463348, Az, 90.0 - Zd);
}
else
{
double tdb = 0;
EphemerisCode.get_earth_nov(ref this.m_earthephobj, tjd, ref tdb, ref peb, ref veb, ref pes, ref ves);
double mobl = 0;
double tobl = 0;
double eq = 0;
double dpsi = 0;
double deps = 0;
NOVAS2.EarthTilt(tdb, ref mobl, ref tobl, ref eq, ref dpsi, ref deps);
double gst = 0;
NOVAS2.SiderealTime(num1, 0.0, eq, ref gst);
NOVAS2.Terra(ref siteInfo, gst, ref pos1, ref vel1);
int num2 = (int)NOVAS2.Nutate(tdb, NutationDirection.TrueToMean, pos1, ref pos2_1);
NOVAS2.Precession(tdb, pos2_1, 2451545.0, ref pos2_6);
int num3 = (int)NOVAS2.Nutate(tdb, NutationDirection.TrueToMean, vel1, ref pos2_5);
NOVAS2.Precession(tdb, pos2_5, 2451545.0, ref pos2_7);
short num4 = 0;
do
{
earthvector1[(int)num4] = peb[(int)num4] + pos2_6[(int)num4];
vel2[(int)num4] = veb[(int)num4] + pos2_7[(int)num4];
earthvector2[(int)num4] = pes[(int)num4] + pos2_6[(int)num4];
++num4;
}while ((int)num4 <= 2);
EphemerisCode.ephemeris_nov(ref this.m_ephobj, tdb, this.m_type, this.m_number, this.m_name, Origin.Barycentric, ref pos1, ref vel1);
double lighttime = 0;
NOVAS2.BaryToGeo(pos1, earthvector1, ref pos2_1, ref lighttime);
double num5 = tdb - lighttime;
int num6 = 0;
double tjd1;
do
{
tjd1 = num5;
EphemerisCode.ephemeris_nov(ref this.m_ephobj, tjd1, this.m_type, this.m_number, this.m_name, Origin.Barycentric, ref pos1, ref vel1);
NOVAS2.BaryToGeo(pos1, earthvector1, ref pos2_1, ref lighttime);
num5 = tdb - lighttime;
checked { ++num6; }
}while (Math.Abs(num5 - tjd1) > 1E-06 & num6 < 100);
if (num6 >= 100)
{
throw new HelperException("Planet:GetTopocentricPoition ephemeris_nov did not converge in 100 iterations");
}
int num7 = (int)NOVAS2.SunField(pos2_1, earthvector2, ref pos2_2);
int num8 = (int)NOVAS2.Aberration(pos2_2, vel2, lighttime, ref pos2_3);
NOVAS2.Precession(2451545.0, pos2_3, tdb, ref pos2_4);
int num9 = (int)NOVAS2.Nutate(tdb, NutationDirection.MeanToTrue, pos2_4, ref pos2);
double ra = 0;
double dec = 0;
int num10 = (int)NOVAS2.Vector2RADec(pos2, ref ra, ref dec);
double num11 = Math.Sqrt(Math.Pow(pos2[0], 2.0) + Math.Pow(pos2[1], 2.0) + Math.Pow(pos2[2], 2.0));
RefractionOption ref_option = RefractionOption.NoRefraction;
if (Refract)
{
bool flag = true;
try
{
siteInfo.Temperature = site.Temperature;
}
catch (Exception ex)
{
//ProjectData.SetProjectError(ex);
flag = false;
//ProjectData.ClearProjectError();
}
try
{
siteInfo.Pressure = site.Pressure;
}
catch (Exception ex)
{
//ProjectData.SetProjectError(ex);
flag = false;
//ProjectData.ClearProjectError();
}
ref_option = !flag ? RefractionOption.StandardRefraction : RefractionOption.LocationRefraction;
}
double zd = 0;
double az = 0;
double rar = 0;
double decr = 0;
if (this.m_bDTValid)
{
NOVAS2.Equ2Hor(tjd, this.m_deltat, 0.0, 0.0, ref siteInfo, ra, dec, ref_option, ref zd, ref az, ref rar, ref decr);
}
else
{
NOVAS2.Equ2Hor(tjd, DeltatCode.DeltaTCalc(tjd), 0.0, 0.0, ref siteInfo, ra, dec, ref_option, ref zd, ref az, ref rar, ref decr);
}
if (ref_option != RefractionOption.NoRefraction)
{
NOVAS2.RADec2Vector(rar, decr, num11, ref pos2);
}
positionVector = new PositionVector(pos2[0], pos2[1], pos2[2], rar, decr, num11, num11 / 173.14463348, az, 90.0 - zd);
}
return(positionVector);
}