|
|
|
|
|
using System;
|
|
|
|
|
|
|
|
|
|
|
|
namespace UnityEngine.AzureSky
|
|
|
|
|
|
{
|
|
|
|
|
|
[Serializable]
|
|
|
|
|
|
public class AzureSkyTimeOfDayComponent
|
|
|
|
|
|
{
|
|
|
|
|
|
//Used by Calendar.
|
|
|
|
|
|
public int selectableDayInt = 15;
|
|
|
|
|
|
public int dayOfYear = 1;
|
|
|
|
|
|
public string[] selectableDayList = new string[]
|
|
|
|
|
|
{
|
|
|
|
|
|
"0", "1", "2", "3", "4" , "5" , "6", "7" , "8", "9", "10",
|
|
|
|
|
|
"11", "12", "13", "14" , "15" , "16", "17" , "18", "19", "20",
|
|
|
|
|
|
"21", "22", "23", "24" , "25" , "26", "27" , "28", "29", "30",
|
|
|
|
|
|
"31", "32", "33", "34", "35", "36", "37", "38", "39", "40", "41"
|
|
|
|
|
|
};
|
|
|
|
|
|
|
|
|
|
|
|
//Used by Calendar and Time of Day tab.
|
|
|
|
|
|
public float hour = 6.0f;
|
|
|
|
|
|
public int day = 15;
|
|
|
|
|
|
public int month = 9;
|
|
|
|
|
|
public int year = 1903;
|
|
|
|
|
|
public int dayInMonth = 30;
|
|
|
|
|
|
public float utc = 0;
|
|
|
|
|
|
public float latitude = 0;
|
|
|
|
|
|
public float longitude = 0;
|
|
|
|
|
|
|
|
|
|
|
|
//Used by Time of Day tab.
|
|
|
|
|
|
public int timeMode = 0;
|
|
|
|
|
|
public float dayCycle = 10.0f;//In minutes
|
|
|
|
|
|
public bool setTimeByCurve = false;
|
|
|
|
|
|
public float hourByCurve = 6.0f;
|
|
|
|
|
|
public AnimationCurve dayCycleCurve = AnimationCurve.Linear(0, 0, 24, 24);
|
|
|
|
|
|
public float lst;
|
|
|
|
|
|
private float m_radians, m_radLatitude, m_sinLatitude, m_cosLatitude;
|
|
|
|
|
|
|
|
|
|
|
|
//Internal uses.
|
|
|
|
|
|
private Vector2 m_hourAndMinutes;
|
|
|
|
|
|
private Vector3 m_date;
|
|
|
|
|
|
private DateTime m_dateTime;
|
|
|
|
|
|
private int m_daysInMonth = 30;
|
|
|
|
|
|
private int m_dayOfWeek = 0;
|
|
|
|
|
|
private int m_fixLeapYear = 0;
|
|
|
|
|
|
|
|
|
|
|
|
/// <summary>
|
|
|
|
|
|
/// Used by AzureSkyController and AzureSkyControllerEditor scripts to correctly adjusts the day buttons and day profiles on the calendar.
|
|
|
|
|
|
/// </summary>
|
|
|
|
|
|
/// <param name="dayList">List containing 42 strings with the day numbers to rename the calendar buttons.</param>
|
|
|
|
|
|
public void UpdateCalendar (string[] dayList)
|
|
|
|
|
|
{
|
|
|
|
|
|
m_daysInMonth = DateTime.DaysInMonth(year, month);
|
|
|
|
|
|
if (day > m_daysInMonth) { day = m_daysInMonth; }
|
|
|
|
|
|
if (day < 1) { day = 1; }
|
|
|
|
|
|
|
|
|
|
|
|
m_dateTime = new DateTime(year, month, 1);
|
|
|
|
|
|
m_dayOfWeek = (int)m_dateTime.DayOfWeek;
|
|
|
|
|
|
selectableDayInt = day - 1 + m_dayOfWeek;
|
|
|
|
|
|
|
|
|
|
|
|
m_dateTime = new DateTime(year, month, day);
|
|
|
|
|
|
for (int i = 0; i < dayList.Length; i++)
|
|
|
|
|
|
{
|
|
|
|
|
|
if (i < m_dayOfWeek || i >= (m_dayOfWeek + m_daysInMonth))
|
|
|
|
|
|
{
|
|
|
|
|
|
dayList[i] = "";
|
|
|
|
|
|
continue;
|
|
|
|
|
|
}
|
|
|
|
|
|
m_dateTime = new DateTime(year, month, (i - m_dayOfWeek) + 1);
|
|
|
|
|
|
dayList[i] = m_dateTime.Day.ToString();
|
|
|
|
|
|
}
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
/// <summary>
|
|
|
|
|
|
/// Get the current day of the week and return an integer between 0 and 6.
|
|
|
|
|
|
/// </summary>
|
|
|
|
|
|
/// <returns></returns>
|
|
|
|
|
|
public int GetDayOfWeek ()
|
|
|
|
|
|
{
|
|
|
|
|
|
m_dateTime = new DateTime(year, month, day);
|
|
|
|
|
|
return (int)m_dateTime.DayOfWeek;
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
/// <summary>
|
|
|
|
|
|
/// Used by AzureSkyControllerEditor script to get the day of the week from a specific day of the month. Returns an integer between 0 and 6.
|
|
|
|
|
|
/// </summary>
|
|
|
|
|
|
/// <param name="dayNumber">Number of the day in the month.</param>
|
|
|
|
|
|
/// <returns></returns>
|
|
|
|
|
|
public int GetDayOfWeek (int dayNumber)
|
|
|
|
|
|
{
|
|
|
|
|
|
m_dateTime = new DateTime(year, month, dayNumber);
|
|
|
|
|
|
return (int)m_dateTime.DayOfWeek;
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
/// <summary>
|
|
|
|
|
|
/// Returns the current day number of a 366-day year.
|
|
|
|
|
|
/// The return value is fix when it is not leap year, this method is ideal for accessing, adding, or removing day profiles from the calendar.
|
|
|
|
|
|
/// If you want to get the correct number for the day in the year, then use System.DateTime.DayOfYear instead.
|
|
|
|
|
|
/// </summary>
|
|
|
|
|
|
/// <returns></returns>
|
|
|
|
|
|
public int GetDayOfYear ()
|
|
|
|
|
|
{
|
|
|
|
|
|
m_dateTime = new DateTime(year, month, day);
|
|
|
|
|
|
dayOfYear = m_dateTime.DayOfYear;
|
|
|
|
|
|
if ((dayOfYear > 59 && !IsLeapYear()))
|
|
|
|
|
|
{
|
|
|
|
|
|
m_fixLeapYear = 1;
|
|
|
|
|
|
}
|
|
|
|
|
|
else { m_fixLeapYear = 0; }
|
|
|
|
|
|
return m_dateTime.DayOfYear - 1 + m_fixLeapYear;
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
/// <summary>
|
|
|
|
|
|
/// Return true if the current year is a leap year.
|
|
|
|
|
|
/// </summary>
|
|
|
|
|
|
/// <returns></returns>
|
|
|
|
|
|
public bool IsLeapYear ()
|
|
|
|
|
|
{
|
|
|
|
|
|
return DateTime.IsLeapYear(year);
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
/// <summary>
|
|
|
|
|
|
/// Converts the timeline in hours and minutes and returns as a Vector2(hours, minutes).
|
|
|
|
|
|
/// </summary>
|
|
|
|
|
|
/// <returns></returns>
|
|
|
|
|
|
public Vector2 GetTime ()
|
|
|
|
|
|
{
|
|
|
|
|
|
if (setTimeByCurve)
|
|
|
|
|
|
{
|
|
|
|
|
|
m_hourAndMinutes.x = Mathf.Floor(hourByCurve);
|
|
|
|
|
|
m_hourAndMinutes.y = 60.0f * (hourByCurve - m_hourAndMinutes.x);
|
|
|
|
|
|
m_hourAndMinutes.y = Mathf.Floor(m_hourAndMinutes.y);
|
|
|
|
|
|
}
|
|
|
|
|
|
else
|
|
|
|
|
|
{
|
|
|
|
|
|
m_hourAndMinutes.x = Mathf.Floor(hour);
|
|
|
|
|
|
m_hourAndMinutes.y = 60.0f * (hour - m_hourAndMinutes.x);
|
|
|
|
|
|
m_hourAndMinutes.y = Mathf.Floor(m_hourAndMinutes.y);
|
|
|
|
|
|
}
|
|
|
|
|
|
return m_hourAndMinutes;
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
/// <summary>
|
|
|
|
|
|
/// Get the current date and returns as a Vector3(month, day, year)
|
|
|
|
|
|
/// </summary>
|
|
|
|
|
|
/// <returns></returns>
|
|
|
|
|
|
public Vector3 GetDate ()
|
|
|
|
|
|
{
|
|
|
|
|
|
m_date.x = month;
|
|
|
|
|
|
m_date.y = day;
|
|
|
|
|
|
m_date.z = year;
|
|
|
|
|
|
return m_date;
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
/// <summary>
|
|
|
|
|
|
/// Get the current system time and returns as a Float.
|
|
|
|
|
|
/// </summary>
|
|
|
|
|
|
/// <returns></returns>
|
|
|
|
|
|
public float GetSystemTime ()
|
|
|
|
|
|
{
|
|
|
|
|
|
return DateTime.Now.Hour + ((1.0f / 60.0f) * DateTime.Now.Minute);
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
/// <summary>
|
|
|
|
|
|
/// Used by AzureSkyController script to convert current system time into a float and apply to timeline.
|
|
|
|
|
|
/// </summary>
|
|
|
|
|
|
public void ApplySystemTime ()
|
|
|
|
|
|
{
|
|
|
|
|
|
hour = DateTime.Now.Hour + ((1.0f / 60.0f) * DateTime.Now.Minute);
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
/// <summary>
|
|
|
|
|
|
/// Used by AzureSkyController script to apply the current system date.
|
|
|
|
|
|
/// </summary>
|
|
|
|
|
|
public void ApplySystemDate ()
|
|
|
|
|
|
{
|
|
|
|
|
|
month = DateTime.Now.Month;
|
|
|
|
|
|
day = DateTime.Now.Day;
|
|
|
|
|
|
year = DateTime.Now.Year;
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
/// <summary>
|
|
|
|
|
|
/// Starts the next day and resets the time based on the Repeat Mode set in the Inspector. This method is commonly called at midnight by AzureSkyController script.
|
|
|
|
|
|
/// </summary>
|
|
|
|
|
|
/// <param name="repeatMode">Place here the Repeat Mode setted on the Options tab.</param>
|
|
|
|
|
|
public void StartNextDay (int repeatMode)
|
|
|
|
|
|
{
|
|
|
|
|
|
if(repeatMode != 1)day += 1;
|
|
|
|
|
|
hour = 0;
|
|
|
|
|
|
if(day > m_daysInMonth)
|
|
|
|
|
|
{
|
|
|
|
|
|
day = 1;
|
|
|
|
|
|
if (repeatMode != 2)
|
|
|
|
|
|
{
|
|
|
|
|
|
month += 1;
|
|
|
|
|
|
if (month > 12)
|
|
|
|
|
|
{
|
|
|
|
|
|
month = 1;
|
|
|
|
|
|
if (repeatMode != 3)
|
|
|
|
|
|
{
|
|
|
|
|
|
year += 1;
|
|
|
|
|
|
if (year > 9999)
|
|
|
|
|
|
{
|
|
|
|
|
|
year = 1;
|
|
|
|
|
|
}
|
|
|
|
|
|
}
|
|
|
|
|
|
}
|
|
|
|
|
|
}
|
|
|
|
|
|
}
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
/// <summary>
|
|
|
|
|
|
/// Skips to the next day and ignores the Repeat Mode set in the Inspector.
|
|
|
|
|
|
/// </summary>
|
|
|
|
|
|
/// <param name="keepTime">Keep the current time?</param>
|
|
|
|
|
|
public void JumpToNextDay ( bool keepTime = true)
|
|
|
|
|
|
{
|
|
|
|
|
|
day += 1;
|
|
|
|
|
|
if (!keepTime)
|
|
|
|
|
|
hour = 0;
|
|
|
|
|
|
if (day > m_daysInMonth)
|
|
|
|
|
|
{
|
|
|
|
|
|
day = 1;
|
|
|
|
|
|
month += 1;
|
|
|
|
|
|
if (month > 12)
|
|
|
|
|
|
{
|
|
|
|
|
|
month = 1;
|
|
|
|
|
|
year += 1;
|
|
|
|
|
|
if (year > 9999)
|
|
|
|
|
|
{
|
|
|
|
|
|
year = 1;
|
|
|
|
|
|
}
|
|
|
|
|
|
}
|
|
|
|
|
|
}
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
/// <summary>
|
|
|
|
|
|
/// Go to a custom date.
|
|
|
|
|
|
/// </summary>
|
|
|
|
|
|
/// <param name="month">The number of the month you want to go.</param>
|
|
|
|
|
|
/// <param name="day">The number of the day you want to go.</param>
|
|
|
|
|
|
/// <param name="year">The number of the year you want to go.</param>
|
|
|
|
|
|
public void GotoDate (int month, int day, int year)
|
|
|
|
|
|
{
|
|
|
|
|
|
this.month = month;
|
|
|
|
|
|
this.day = day;
|
|
|
|
|
|
this.year = year;
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
/// <summary>
|
|
|
|
|
|
/// Go to a custom time.
|
|
|
|
|
|
/// </summary>
|
|
|
|
|
|
/// <param name="hours">The number of the hour you want to go.</param>
|
|
|
|
|
|
/// <param name="minutes">The number of the minute you want to go.</param>
|
|
|
|
|
|
public void GotoTime (int hours, int minutes)
|
|
|
|
|
|
{
|
|
|
|
|
|
this.hour = hours + ((1.0f / 60.0f) * minutes);
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
/// <summary>
|
|
|
|
|
|
/// Used by AzureSkyController script to apply the time progression.
|
|
|
|
|
|
/// </summary>
|
|
|
|
|
|
/// <returns></returns>
|
|
|
|
|
|
public float GetDayLength ()
|
|
|
|
|
|
{
|
|
|
|
|
|
if (dayCycle > 0.0f)
|
|
|
|
|
|
return (24.0f / 60.0f) / dayCycle;
|
|
|
|
|
|
else
|
|
|
|
|
|
return 0.0f;
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
/// <summary>
|
|
|
|
|
|
/// Used by AzureSkyController script to calculate the time of day based on the curve time.
|
|
|
|
|
|
/// </summary>
|
|
|
|
|
|
public void CalculateTimeByCurve ()
|
|
|
|
|
|
{
|
|
|
|
|
|
hourByCurve = dayCycleCurve.Evaluate(hour);
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
/// <summary>
|
|
|
|
|
|
/// Used by AzureSkyController script to set the simple sun position.
|
|
|
|
|
|
/// </summary>
|
|
|
|
|
|
/// <returns></returns>
|
|
|
|
|
|
public float SetSimpleSunPosition ()
|
|
|
|
|
|
{
|
|
|
|
|
|
if (setTimeByCurve)
|
|
|
|
|
|
{
|
|
|
|
|
|
return ((hourByCurve + utc) * 360.0f / 24.0f) - 90.0f;
|
|
|
|
|
|
}
|
|
|
|
|
|
else
|
|
|
|
|
|
{
|
|
|
|
|
|
return ((hour + utc) * 360.0f / 24.0f) - 90.0f;
|
|
|
|
|
|
}
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
/// <summary>
|
|
|
|
|
|
/// Used by AzureSkyController script to set the simple moon position.
|
|
|
|
|
|
/// </summary>
|
|
|
|
|
|
/// <returns></returns>
|
|
|
|
|
|
public Quaternion SetSimpleMoonPosition ( Transform sun)
|
|
|
|
|
|
{
|
|
|
|
|
|
return sun.transform.rotation * Quaternion.Euler(0, -180, 0);
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
/// <summary>
|
|
|
|
|
|
/// Used by AzureSkyController script to set the realistic sun position based on Time, Date and Location.
|
|
|
|
|
|
/// </summary>
|
|
|
|
|
|
/// <returns></returns>
|
|
|
|
|
|
public Vector3 SetRealisticSunPosition ()
|
|
|
|
|
|
{
|
|
|
|
|
|
m_radians = (Mathf.PI * 2.0f) / 360.0f;//Used to convert degress to radians.
|
|
|
|
|
|
m_radLatitude = m_radians * latitude;
|
|
|
|
|
|
m_sinLatitude = Mathf.Sin(m_radLatitude);
|
|
|
|
|
|
m_cosLatitude = Mathf.Cos(m_radLatitude);
|
|
|
|
|
|
|
|
|
|
|
|
float hour = this.hour - utc;
|
|
|
|
|
|
//Time Scale.
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
//d = 367*y - 7 * ( y + (m+9)/12 ) / 4 + 275*m/9 + D - 730530
|
|
|
|
|
|
//d = d + UT/24.0
|
|
|
|
|
|
float d = 367 * year - 7 * (year + (month + 9) / 12) / 4 + 275 * month / 9 + day - 730530;
|
|
|
|
|
|
d = d + hour / 24.0f;
|
|
|
|
|
|
|
|
|
|
|
|
//Tilt of earth's axis.
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
//obliquity of the ecliptic.
|
|
|
|
|
|
float ecliptic = 23.4393f - 3.563E-7f * d;
|
|
|
|
|
|
//Need convert to radians before apply sine and cosine.
|
|
|
|
|
|
float radEcliptic = m_radians * ecliptic;
|
|
|
|
|
|
float sinEcliptic = Mathf.Sin(radEcliptic);
|
|
|
|
|
|
float cosEcliptic = Mathf.Cos(radEcliptic);
|
|
|
|
|
|
|
|
|
|
|
|
//Orbital elements of the Sun.
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
//float N = 0.0;
|
|
|
|
|
|
//float i = 0.0;
|
|
|
|
|
|
float w = 282.9404f + 4.70935E-5f * d;
|
|
|
|
|
|
//float a = 1.000000f;
|
|
|
|
|
|
float e = 0.016709f - 1.151E-9f * d;
|
|
|
|
|
|
float M = 356.0470f + 0.9856002585f * d;
|
|
|
|
|
|
|
|
|
|
|
|
//Eccentric anomaly.
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
//E = M + e*(180/pi) * sin(M) * ( 1.0 + e * cos(M) ) in degress.
|
|
|
|
|
|
//E = M + e * sin(M) * ( 1.0 + e * cos(M) ) in radians.
|
|
|
|
|
|
//Need convert to radians before apply sine and cosine.
|
|
|
|
|
|
float radM = m_radians * M;
|
|
|
|
|
|
float sinM = Mathf.Sin(radM);
|
|
|
|
|
|
float cosM = Mathf.Cos(radM);
|
|
|
|
|
|
|
|
|
|
|
|
//Need convert to radians before apply sine and cosine.
|
|
|
|
|
|
float radE = radM + e * sinM * (1.0f + e * cosM);
|
|
|
|
|
|
float sinE = Mathf.Sin(radE);
|
|
|
|
|
|
float cosE = Mathf.Cos(radE);
|
|
|
|
|
|
|
|
|
|
|
|
//Sun's distance (r) and its true anomaly (v).
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
//Xv = r * cos (v) = cos (E) - e
|
|
|
|
|
|
//Yv = r * sen (v) = sqrt (1,0 - e * e) * sen (E)
|
|
|
|
|
|
float xv = cosE - e;
|
|
|
|
|
|
float yv = Mathf.Sqrt(1.0f - e * e) * sinE;
|
|
|
|
|
|
|
|
|
|
|
|
//V = atan2 (yv, xv)
|
|
|
|
|
|
//R = sqrt (xv * xv + yv * yv)
|
|
|
|
|
|
float v = Mathf.Rad2Deg * Mathf.Atan2(yv, xv);
|
|
|
|
|
|
float r = Mathf.Sqrt(xv * xv + yv * yv);
|
|
|
|
|
|
|
|
|
|
|
|
//Sun's true longitude.
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
float radLongitude = m_radians * (v + w);
|
|
|
|
|
|
float sinLongitude = Mathf.Sin(radLongitude);
|
|
|
|
|
|
float cosLongitude = Mathf.Cos(radLongitude);
|
|
|
|
|
|
|
|
|
|
|
|
float xs = r * cosLongitude;
|
|
|
|
|
|
float ys = r * sinLongitude;
|
|
|
|
|
|
|
|
|
|
|
|
//Equatorial coordinates.
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
float xe = xs;
|
|
|
|
|
|
float ye = ys * cosEcliptic;
|
|
|
|
|
|
float ze = ys * sinEcliptic;
|
|
|
|
|
|
|
|
|
|
|
|
//Sun's Right Ascension(RA) and Declination(Dec).
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
float RA = Mathf.Atan2(ye, xe);
|
|
|
|
|
|
float Dec = Mathf.Atan2(ze, Mathf.Sqrt(xe * xe + ye * ye));
|
|
|
|
|
|
float sinDec = Mathf.Sin(Dec);
|
|
|
|
|
|
float cosDec = Mathf.Cos(Dec);
|
|
|
|
|
|
|
|
|
|
|
|
//The Sidereal Time.
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
float Ls = v + w;
|
|
|
|
|
|
|
|
|
|
|
|
float GMST0 = Ls + 180.0f;
|
|
|
|
|
|
float UT = 15.0f * hour;//Universal Time.
|
|
|
|
|
|
float GMST = GMST0 + UT;
|
|
|
|
|
|
float LST = m_radians * (GMST + longitude);
|
|
|
|
|
|
|
|
|
|
|
|
//Store local sideral time.
|
|
|
|
|
|
lst = LST;
|
|
|
|
|
|
|
|
|
|
|
|
//Azimuthal coordinates.
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
float HA = LST - RA;
|
|
|
|
|
|
float sinHA = Mathf.Sin(HA);
|
|
|
|
|
|
float cosHA = Mathf.Cos(HA);
|
|
|
|
|
|
|
|
|
|
|
|
float x = cosHA * cosDec;
|
|
|
|
|
|
float y = sinHA * cosDec;
|
|
|
|
|
|
float z = sinDec;
|
|
|
|
|
|
|
|
|
|
|
|
float xhor = x * m_sinLatitude - z * m_cosLatitude;
|
|
|
|
|
|
float yhor = y;
|
|
|
|
|
|
float zhor = x * m_cosLatitude + z * m_sinLatitude;
|
|
|
|
|
|
|
|
|
|
|
|
//az = atan2( yhor, xhor ) + 180_degrees
|
|
|
|
|
|
//alt = asin( zhor ) = atan2( zhor, sqrt(xhor*xhor+yhor*yhor) )
|
|
|
|
|
|
float azimuth = Mathf.Atan2(yhor, xhor) + m_radians * 180.0f;
|
|
|
|
|
|
float altitude = Mathf.Asin(zhor);
|
|
|
|
|
|
|
|
|
|
|
|
//Zenith angle.
|
|
|
|
|
|
//Zenith=90°−α Where α is the elevation angle.
|
|
|
|
|
|
float zenith = 90.0f * m_radians - altitude;
|
|
|
|
|
|
|
|
|
|
|
|
//Converts from Spherical(radius r, zenith-inclination θ, azimuth φ) to Cartesian(x,y,z) coordinates.
|
|
|
|
|
|
//https://en.wikipedia.org/wiki/Spherical_coordinate_system
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
//x = r sin(θ)cos(φ)
|
|
|
|
|
|
//y = r sin(θ)sin(φ)
|
|
|
|
|
|
//z = r cos(θ)
|
|
|
|
|
|
Vector3 ret;
|
|
|
|
|
|
|
|
|
|
|
|
//radius = 1
|
|
|
|
|
|
ret.z = Mathf.Sin(zenith) * Mathf.Cos(azimuth);
|
|
|
|
|
|
ret.x = Mathf.Sin(zenith) * Mathf.Sin(azimuth);
|
|
|
|
|
|
ret.y = Mathf.Cos(zenith);
|
|
|
|
|
|
|
|
|
|
|
|
return ret * -1.0f;
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
/// <summary>
|
|
|
|
|
|
/// Used by AzureSkyController script to set the realistic moon position based on Time, Date and Location.
|
|
|
|
|
|
/// </summary>
|
|
|
|
|
|
/// <returns></returns>
|
|
|
|
|
|
public Vector3 SetRealisticMoonPosition ()
|
|
|
|
|
|
{
|
|
|
|
|
|
float hour = this.hour - utc;
|
|
|
|
|
|
|
|
|
|
|
|
//Time Scale.
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
//d = 367*y - 7 * ( y + (m+9)/12 ) / 4 + 275*m/9 + D - 730530
|
|
|
|
|
|
//d = d + UT/24.0
|
|
|
|
|
|
float d = 367 * year - 7 * (year + (month + 9) / 12) / 4 + 275 * month / 9 + day - 730530;
|
|
|
|
|
|
d = d + hour / 24.0f;
|
|
|
|
|
|
|
|
|
|
|
|
//Tilt of earth's axis.
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
//obliquity of the ecliptic.
|
|
|
|
|
|
float ecliptic = 23.4393f - 3.563E-7f * d;
|
|
|
|
|
|
//Need convert to radians before apply sine and cosine.
|
|
|
|
|
|
float radEcliptic = m_radians * ecliptic;
|
|
|
|
|
|
float sinEcliptic = Mathf.Sin(radEcliptic);
|
|
|
|
|
|
float cosEcliptic = Mathf.Cos(radEcliptic);
|
|
|
|
|
|
|
|
|
|
|
|
//Orbital elements of the Moon.
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
float N = 125.1228f - 0.0529538083f * d;
|
|
|
|
|
|
float i = 5.1454f;
|
|
|
|
|
|
float w = 318.0634f + 0.1643573223f * d;
|
|
|
|
|
|
float a = 60.2666f;
|
|
|
|
|
|
float e = 0.054900f;
|
|
|
|
|
|
float M = 115.3654f + 13.0649929509f * d;
|
|
|
|
|
|
|
|
|
|
|
|
//Eccentric anomaly.
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
//E = M + e*(180/pi) * sin(M) * ( 1.0 + e * cos(M) )
|
|
|
|
|
|
float radM = m_radians * M;
|
|
|
|
|
|
float E = radM + e * Mathf.Sin(radM) * (1f + e * Mathf.Cos(radM));
|
|
|
|
|
|
|
|
|
|
|
|
//Planet's distance and true anomaly.
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
//xv = r * cos(v) = a * ( cos(E) - e )
|
|
|
|
|
|
//yv = r * sin(v) = a * ( sqrt(1.0 - e*e) * sin(E) )
|
|
|
|
|
|
float xv = a * (Mathf.Cos(E) - e);
|
|
|
|
|
|
float yv = a * (Mathf.Sqrt(1f - e * e) * Mathf.Sin(E));
|
|
|
|
|
|
//V = atan2 (yv, xv)
|
|
|
|
|
|
//R = sqrt (xv * xv + yv * yv)
|
|
|
|
|
|
float v = Mathf.Rad2Deg * Mathf.Atan2(yv, xv);
|
|
|
|
|
|
float r = Mathf.Sqrt(xv * xv + yv * yv);
|
|
|
|
|
|
|
|
|
|
|
|
//Moon position in 3D space.
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
float radLongitude = m_radians * (v + w);
|
|
|
|
|
|
float sinLongitude = Mathf.Sin(radLongitude);
|
|
|
|
|
|
float cosLongitude = Mathf.Cos(radLongitude);
|
|
|
|
|
|
|
|
|
|
|
|
//Geocentric (Earth-centered) coordinates.
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
//xh = r * ( cos(N) * cos(v+w) - sin(N) * sin(v+w) * cos(i) )
|
|
|
|
|
|
//yh = r * ( sin(N) * cos(v+w) + cos(N) * sin(v+w) * cos(i) )
|
|
|
|
|
|
//zh = r * ( sin(v+w) * sin(i) )
|
|
|
|
|
|
float radN = m_radians * N;
|
|
|
|
|
|
float radI = m_radians * i;
|
|
|
|
|
|
|
|
|
|
|
|
float xh = r * (Mathf.Cos(radN) * cosLongitude - Mathf.Sin(radN) * sinLongitude * Mathf.Cos(radI));
|
|
|
|
|
|
float yh = r * (Mathf.Sin(radN) * cosLongitude + Mathf.Cos(radN) * sinLongitude * Mathf.Cos(radI));
|
|
|
|
|
|
float zh = r * (sinLongitude * Mathf.Sin(radI));
|
|
|
|
|
|
|
|
|
|
|
|
//float xg = xh; //No needed to the moon.
|
|
|
|
|
|
//float yg = yh;
|
|
|
|
|
|
//float zg = zh;
|
|
|
|
|
|
|
|
|
|
|
|
//Equatorial coordinates.
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
float xe = xh;
|
|
|
|
|
|
float ye = yh * cosEcliptic - zh * sinEcliptic;
|
|
|
|
|
|
float ze = yh * sinEcliptic + zh * cosEcliptic;
|
|
|
|
|
|
|
|
|
|
|
|
//Planet's Right Ascension (RA) and Declination (Dec).
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
float RA = Mathf.Atan2(ye, xe);
|
|
|
|
|
|
float Dec = Mathf.Atan2(ze, Mathf.Sqrt(xe * xe + ye * ye));
|
|
|
|
|
|
|
|
|
|
|
|
//The Sidereal Time.
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
//It is already calculated for the sun and stored in the lst, it is not necessary to calculate again for the moon.
|
|
|
|
|
|
//float Ls = ls;
|
|
|
|
|
|
|
|
|
|
|
|
//float GMST0 = Ls + 180.0f;
|
|
|
|
|
|
//float UT = 15.0f * hour;
|
|
|
|
|
|
//float GMST = GMST0 + UT;
|
|
|
|
|
|
//float LST = radians * (GMST + Azure_Longitude);
|
|
|
|
|
|
|
|
|
|
|
|
//Azimuthal coordinates.
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
float HA = lst - RA;
|
|
|
|
|
|
|
|
|
|
|
|
float x = Mathf.Cos(HA) * Mathf.Cos(Dec);
|
|
|
|
|
|
float y = Mathf.Sin(HA) * Mathf.Cos(Dec);
|
|
|
|
|
|
float z = Mathf.Sin(Dec);
|
|
|
|
|
|
|
|
|
|
|
|
float xhor = x * m_sinLatitude - z * m_cosLatitude;
|
|
|
|
|
|
float yhor = y;
|
|
|
|
|
|
float zhor = x * m_cosLatitude + z * m_sinLatitude;
|
|
|
|
|
|
|
|
|
|
|
|
//az = atan2( yhor, xhor ) + 180_degrees
|
|
|
|
|
|
//alt = asin( zhor ) = atan2( zhor, sqrt(xhor*xhor+yhor*yhor) )
|
|
|
|
|
|
float azimuth = Mathf.Atan2(yhor, xhor) + m_radians * 180.0f;
|
|
|
|
|
|
float altitude = Mathf.Asin(zhor);
|
|
|
|
|
|
|
|
|
|
|
|
//Zenith angle.
|
|
|
|
|
|
//Zenith = 90°−α where α is the elevation angle.
|
|
|
|
|
|
float zenith = 90.0f * m_radians - altitude;
|
|
|
|
|
|
|
|
|
|
|
|
//Converts from Spherical(radius r, zenith-inclination θ, azimuth φ) to Cartesian(x,y,z) coordinates.
|
|
|
|
|
|
//https://en.wikipedia.org/wiki/Spherical_coordinate_system
|
|
|
|
|
|
//---------------------------------------------------------------------------------------------------
|
|
|
|
|
|
//x = r sin(θ)cos(φ)
|
|
|
|
|
|
//y = r sin(θ)sin(φ)
|
|
|
|
|
|
//z = r cos(θ)
|
|
|
|
|
|
Vector3 ret;
|
|
|
|
|
|
|
|
|
|
|
|
//radius = 1
|
|
|
|
|
|
ret.z = Mathf.Sin(zenith) * Mathf.Cos(azimuth);
|
|
|
|
|
|
ret.x = Mathf.Sin(zenith) * Mathf.Sin(azimuth);
|
|
|
|
|
|
ret.y = Mathf.Cos(zenith);
|
|
|
|
|
|
|
|
|
|
|
|
return ret * -1.0f;
|
|
|
|
|
|
}
|
|
|
|
|
|
}
|
|
|
|
|
|
}
|
