public PositionVector GetApparentPosition(double tjd)
{
CatEntry star = new CatEntry();
PositionVector positionVector = new PositionVector();
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[] pos = 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[] pos2_5 = new double[4];
double[] vel = new double[4];
double[] pos2_6 = new double[4];
if (!(this.m_rav & this.m_decv))
{
throw new ASCOM.Astrometry.Exceptions.ValueNotSetException("Star.GetApparentPosition RA or DEC 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_6[0] = 0.0;
pos2_6[1] = 0.0;
pos2_6[2] = 0.0;
throw new NOVASFunctionException("Star.GetApparentPosition", "get_earth", this.hr);
}
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 pos, ref vel);
NOVAS2.ProperMotion(2451545.0, pos, vel, tdb, ref pos2_1);
double lighttime = 0;
NOVAS2.BaryToGeo(pos2_1, bary_earthp, ref pos2_2, ref lighttime);
int num1 = (int)NOVAS2.SunField(pos2_2, helio_earthp, ref pos2_3);
int num2 = (int)NOVAS2.Aberration(pos2_3, bary_earthv, lighttime, ref pos2_4);
NOVAS2.Precession(2451545.0, pos2_4, tdb, ref pos2_5);
int num3 = (int)NOVAS2.Nutate(tdb, NutationDirection.MeanToTrue, pos2_5, ref pos2_6);
positionVector.x = pos2_6[0];
positionVector.y = pos2_6[1];
positionVector.z = pos2_6[2];
return(positionVector);
}
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 GetLocalPosition(double tjd, Site site)
{
CatEntry star = new CatEntry();
PositionVector positionVector = new PositionVector();
SiteInfo locale = new SiteInfo();
double[] pos2_1 = new double[4];
double[] pos2_2 = new double[4];
double[] earthvector1 = new double[4];
double[] vel1 = new double[4];
double[] earthvector2 = new double[4];
double[] numArray = new double[4];
double[] pos = new double[4];
double[] vel2 = 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[] 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_7 = new double[4];
if (!(this.m_rav & this.m_decv))
{
throw new ASCOM.Astrometry.Exceptions.ValueNotSetException("Star.GetLocalPosition RA or DEC not available");
}
double jd_high = !this.m_bDTValid ? tjd - DeltatCode.DeltaTCalc(tjd) / 86400.0 : tjd - this.m_deltat;
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;
try
{
locale.Latitude = site.Latitude;
}
catch (Exception ex)
{
//ProjectData.SetProjectError(ex);
throw new ValueNotAvailableException("Star:GetLocalPosition Site.Latitude is not available");
}
try
{
locale.Longitude = site.Longitude;
}
catch (Exception ex)
{
//ProjectData.SetProjectError(ex);
throw new ValueNotAvailableException("Star:GetLocalPosition Site.Longitude is not available");
}
try
{
locale.Height = site.Height;
}
catch (Exception ex)
{
//ProjectData.SetProjectError(ex);
throw new ValueNotAvailableException("Star:GetLocalPosition 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_7[0] = 0.0;
pos2_7[1] = 0.0;
pos2_7[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 locale, gst, ref pos, ref vel2);
int num1 = (int)NOVAS2.Nutate(tdb, NutationDirection.TrueToMean, pos, ref pos2_3);
NOVAS2.Precession(tdb, pos2_3, 2451545.0, ref pos2_1);
int num2 = (int)NOVAS2.Nutate(tdb, NutationDirection.TrueToMean, vel2, ref pos2_4);
NOVAS2.Precession(tdb, pos2_4, 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];
numArray[index] = helio_earthv[index] + pos2_2[index];
checked { ++index; }
}while (index <= 2);
NOVAS2.StarVectors(star, ref pos, ref vel2);
NOVAS2.ProperMotion(2451545.0, pos, vel2, tdb, ref pos2_3);
double lighttime = 0;
NOVAS2.BaryToGeo(pos2_3, earthvector1, ref pos2_5, ref lighttime);
int num3 = (int)NOVAS2.SunField(pos2_5, earthvector2, ref pos2_6);
int num4 = (int)NOVAS2.Aberration(pos2_6, vel1, lighttime, ref pos2_7);
positionVector.x = pos2_7[0];
positionVector.y = pos2_7[1];
positionVector.z = pos2_7[2];
return(positionVector);
}
public short GetEarth(double tjd, ref BodyDescription earth, ref double tdb, ref double[] bary_earthp, ref double[] bary_earthv, ref double[] helio_earthp, ref double[] helio_earthv)
{
return(NOVAS2.GetEarth(tjd, ref earth, ref tdb, ref bary_earthp, ref bary_earthv, ref helio_earthp, ref helio_earthv));
}