mirror of
https://github.com/FlightControl-Master/MOOSE.git
synced 2025-10-29 16:58:06 +00:00
Astro Update
**COORDINATE** * Added functions to get sunrise and sunset times. **UTILS** * Added functions to calculate sunset and sunrise * Added function to replace illegal characters. **AIRBOSS v1.1.5** * Fixed wrong function name when recovery is resumed. * Replace illegcal characters in player names when writing trapsheet.
This commit is contained in:
parent
21652de804
commit
61bb59d8b3
@ -1969,6 +1969,240 @@ do -- COORDINATE
|
|||||||
return InSphere
|
return InSphere
|
||||||
end
|
end
|
||||||
|
|
||||||
|
--- Get sun rise time for a specific date at the coordinate.
|
||||||
|
-- @param #COORDINATE self
|
||||||
|
-- @param #number Day The day.
|
||||||
|
-- @param #number Month The month.
|
||||||
|
-- @param #number Year The year.
|
||||||
|
-- @param #boolean InSeconds If true, return the sun rise time in seconds.
|
||||||
|
-- @return #string Sunrise time, e.g. "05:41".
|
||||||
|
function COORDINATE:GetSunriseAtDate(Day, Month, Year, InSeconds)
|
||||||
|
|
||||||
|
-- Day of the year.
|
||||||
|
local DayOfYear=UTILS.GetDayOfYear(Year, Month, Day)
|
||||||
|
|
||||||
|
local Latitude, Longitude=self:GetLLDDM()
|
||||||
|
|
||||||
|
local Tdiff=UTILS.GMTToLocalTimeDifference()
|
||||||
|
|
||||||
|
local sunrise=UTILS.GetSunRiseAndSet(DayOfYear, Latitude, Longitude, true, Tdiff)
|
||||||
|
|
||||||
|
if InSeconds then
|
||||||
|
return sunrise
|
||||||
|
else
|
||||||
|
return UTILS.SecondsToClock(sunrise, true)
|
||||||
|
end
|
||||||
|
|
||||||
|
end
|
||||||
|
|
||||||
|
--- Get todays sun rise time.
|
||||||
|
-- @param #COORDINATE self
|
||||||
|
-- @param #boolean InSeconds If true, return the sun rise time in seconds.
|
||||||
|
-- @return #string Sunrise time, e.g. "05:41".
|
||||||
|
function COORDINATE:GetSunrise(InSeconds)
|
||||||
|
|
||||||
|
-- Get current day of the year.
|
||||||
|
local DayOfYear=UTILS.GetMissionDayOfYear()
|
||||||
|
|
||||||
|
-- Lat and long at this point.
|
||||||
|
local Latitude, Longitude=self:GetLLDDM()
|
||||||
|
|
||||||
|
-- GMT time diff.
|
||||||
|
local Tdiff=UTILS.GMTToLocalTimeDifference()
|
||||||
|
|
||||||
|
-- Sunrise in seconds of the day.
|
||||||
|
local sunrise=UTILS.GetSunRiseAndSet(DayOfYear, Latitude, Longitude, true, Tdiff)
|
||||||
|
|
||||||
|
-- Debug output.
|
||||||
|
--self:I(string.format("Sun rise at lat=%.3f long=%.3f on %s + %d days (DayOfYear=%d): %s (GMT %d)", Latitude, Longitude, date, x, DayOfYear, UTILS.SecondsToClock(sunrise), Tdiff))
|
||||||
|
|
||||||
|
if InSeconds then
|
||||||
|
return sunrise
|
||||||
|
else
|
||||||
|
return UTILS.SecondsToClock(sunrise, true)
|
||||||
|
end
|
||||||
|
|
||||||
|
end
|
||||||
|
|
||||||
|
--- Get minutes until the next sun rise at this coordinate.
|
||||||
|
-- @param #COORDINATE self
|
||||||
|
-- @param OnlyToday If true, only calculate the sun rise of today. If sun has already risen, the time in negative minutes since sunrise is reported.
|
||||||
|
-- @return #number Minutes to the next sunrise.
|
||||||
|
function COORDINATE:GetMinutesToSunrise(OnlyToday)
|
||||||
|
|
||||||
|
-- Seconds of today
|
||||||
|
local time=UTILS.SecondsOfToday()
|
||||||
|
|
||||||
|
-- Next Sunrise in seconds.
|
||||||
|
local sunrise=nil
|
||||||
|
|
||||||
|
-- Time to sunrise.
|
||||||
|
local delta=nil
|
||||||
|
|
||||||
|
if OnlyToday then
|
||||||
|
|
||||||
|
---
|
||||||
|
-- Sunrise of today
|
||||||
|
---
|
||||||
|
|
||||||
|
sunrise=self:GetSunrise(true)
|
||||||
|
|
||||||
|
delta=sunrise-time
|
||||||
|
|
||||||
|
else
|
||||||
|
|
||||||
|
---
|
||||||
|
-- Sunrise of tomorrow
|
||||||
|
---
|
||||||
|
|
||||||
|
-- Tomorrows day of the year.
|
||||||
|
local DayOfYear=UTILS.GetMissionDayOfYear()+1
|
||||||
|
|
||||||
|
local Latitude, Longitude=self:GetLLDDM()
|
||||||
|
|
||||||
|
local Tdiff=UTILS.GMTToLocalTimeDifference()
|
||||||
|
|
||||||
|
sunrise=UTILS.GetSunRiseAndSet(DayOfYear, Latitude, Longitude, true, Tdiff)
|
||||||
|
|
||||||
|
delta=sunrise+UTILS.SecondsToMidnight()
|
||||||
|
|
||||||
|
end
|
||||||
|
|
||||||
|
return delta/60
|
||||||
|
end
|
||||||
|
|
||||||
|
--- Check if it is day, i.e. if the sun has risen about the horizon at this coordinate.
|
||||||
|
-- @param #COORDINATE self
|
||||||
|
-- @return #boolean If true, it is day. If false, it is night time.
|
||||||
|
function COORDINATE:IsDay()
|
||||||
|
|
||||||
|
-- Todays sun rise in sec.
|
||||||
|
local sunrise=self:GetSunrise(true)
|
||||||
|
|
||||||
|
-- Todays sun set in sec.
|
||||||
|
local sunset=self:GetSunset(true)
|
||||||
|
|
||||||
|
-- Seconds passed since midnight.
|
||||||
|
local time=UTILS.SecondsOfToday()
|
||||||
|
|
||||||
|
-- Check if time is between sunrise and sunset.
|
||||||
|
if time>sunrise and time<=sunset then
|
||||||
|
return true
|
||||||
|
else
|
||||||
|
return false
|
||||||
|
end
|
||||||
|
|
||||||
|
end
|
||||||
|
|
||||||
|
--- Check if it is night, i.e. if the sun has set below the horizon at this coordinate.
|
||||||
|
-- @param #COORDINATE self
|
||||||
|
-- @return #boolean If true, it is night. If false, it is day time.
|
||||||
|
function COORDINATE:IsNight()
|
||||||
|
return not self:IsDay()
|
||||||
|
end
|
||||||
|
|
||||||
|
--- Get sun set time for a specific date at the coordinate.
|
||||||
|
-- @param #COORDINATE self
|
||||||
|
-- @param #number Day The day.
|
||||||
|
-- @param #number Month The month.
|
||||||
|
-- @param #number Year The year.
|
||||||
|
-- @param #boolean InSeconds If true, return the sun rise time in seconds.
|
||||||
|
-- @return #string Sunset time, e.g. "20:41".
|
||||||
|
function COORDINATE:GetSunsetAtDate(Day, Month, Year, InSeconds)
|
||||||
|
|
||||||
|
-- Day of the year.
|
||||||
|
local DayOfYear=UTILS.GetDayOfYear(Year, Month, Day)
|
||||||
|
|
||||||
|
local Latitude, Longitude=self:GetLLDDM()
|
||||||
|
|
||||||
|
local Tdiff=UTILS.GMTToLocalTimeDifference()
|
||||||
|
|
||||||
|
local sunset=UTILS.GetSunRiseAndSet(DayOfYear, Latitude, Longitude, false, Tdiff)
|
||||||
|
|
||||||
|
if InSeconds then
|
||||||
|
return sunset
|
||||||
|
else
|
||||||
|
return UTILS.SecondsToClock(sunset, true)
|
||||||
|
end
|
||||||
|
|
||||||
|
end
|
||||||
|
|
||||||
|
--- Get todays sun set time.
|
||||||
|
-- @param #COORDINATE self
|
||||||
|
-- @param #boolean InSeconds If true, return the sun set time in seconds.
|
||||||
|
-- @return #string Sunrise time, e.g. "20:41".
|
||||||
|
function COORDINATE:GetSunset(InSeconds)
|
||||||
|
|
||||||
|
-- Get current day of the year.
|
||||||
|
local DayOfYear=UTILS.GetMissionDayOfYear()
|
||||||
|
|
||||||
|
-- Lat and long at this point.
|
||||||
|
local Latitude, Longitude=self:GetLLDDM()
|
||||||
|
|
||||||
|
-- GMT time diff.
|
||||||
|
local Tdiff=UTILS.GMTToLocalTimeDifference()
|
||||||
|
|
||||||
|
-- Sunrise in seconds of the day.
|
||||||
|
local sunrise=UTILS.GetSunRiseAndSet(DayOfYear, Latitude, Longitude, false, Tdiff)
|
||||||
|
|
||||||
|
-- Debug output.
|
||||||
|
--self:I(string.format("Sun rise at lat=%.3f long=%.3f on %s + %d days (DayOfYear=%d): %s (GMT %d)", Latitude, Longitude, date, x, DayOfYear, UTILS.SecondsToClock(sunrise), Tdiff))
|
||||||
|
|
||||||
|
if InSeconds then
|
||||||
|
return sunrise
|
||||||
|
else
|
||||||
|
return UTILS.SecondsToClock(sunrise, true)
|
||||||
|
end
|
||||||
|
|
||||||
|
end
|
||||||
|
|
||||||
|
--- Get minutes until the next sun set at this coordinate.
|
||||||
|
-- @param #COORDINATE self
|
||||||
|
-- @param OnlyToday If true, only calculate the sun set of today. If sun has already set, the time in negative minutes since sunset is reported.
|
||||||
|
-- @return #number Minutes to the next sunrise.
|
||||||
|
function COORDINATE:GetMinutesToSunset(OnlyToday)
|
||||||
|
|
||||||
|
-- Seconds of today
|
||||||
|
local time=UTILS.SecondsOfToday()
|
||||||
|
|
||||||
|
-- Next Sunset in seconds.
|
||||||
|
local sunset=nil
|
||||||
|
|
||||||
|
-- Time to sunrise.
|
||||||
|
local delta=nil
|
||||||
|
|
||||||
|
if OnlyToday then
|
||||||
|
|
||||||
|
---
|
||||||
|
-- Sunset of today
|
||||||
|
---
|
||||||
|
|
||||||
|
sunset=self:GetSunset(true)
|
||||||
|
|
||||||
|
delta=sunset-time
|
||||||
|
|
||||||
|
else
|
||||||
|
|
||||||
|
---
|
||||||
|
-- Sunset of tomorrow
|
||||||
|
---
|
||||||
|
|
||||||
|
-- Tomorrows day of the year.
|
||||||
|
local DayOfYear=UTILS.GetMissionDayOfYear()+1
|
||||||
|
|
||||||
|
local Latitude, Longitude=self:GetLLDDM()
|
||||||
|
|
||||||
|
local Tdiff=UTILS.GMTToLocalTimeDifference()
|
||||||
|
|
||||||
|
sunset=UTILS.GetSunRiseAndSet(DayOfYear, Latitude, Longitude, false, Tdiff)
|
||||||
|
|
||||||
|
delta=sunset+UTILS.SecondsToMidnight()
|
||||||
|
|
||||||
|
end
|
||||||
|
|
||||||
|
return delta/60
|
||||||
|
end
|
||||||
|
|
||||||
|
|
||||||
--- Return a BR string from a COORDINATE to the COORDINATE.
|
--- Return a BR string from a COORDINATE to the COORDINATE.
|
||||||
-- @param #COORDINATE self
|
-- @param #COORDINATE self
|
||||||
@ -2036,6 +2270,14 @@ do -- COORDINATE
|
|||||||
return ""
|
return ""
|
||||||
end
|
end
|
||||||
|
|
||||||
|
--- Get Latitude and Longitude in Degrees Decimal Minutes (DDM).
|
||||||
|
-- @param #COORDINATE self
|
||||||
|
-- @return #number Latitude in DDM.
|
||||||
|
-- @return #number Lontitude in DDM.
|
||||||
|
function COORDINATE:GetLLDDM()
|
||||||
|
return coord.LOtoLL( self:GetVec3() )
|
||||||
|
end
|
||||||
|
|
||||||
--- Provides a Lat Lon string in Degree Minute Second format.
|
--- Provides a Lat Lon string in Degree Minute Second format.
|
||||||
-- @param #COORDINATE self
|
-- @param #COORDINATE self
|
||||||
-- @param Core.Settings#SETTINGS Settings (optional) The settings. Can be nil, and in this case the default settings are used. If you want to specify your own settings, use the _SETTINGS object.
|
-- @param Core.Settings#SETTINGS Settings (optional) The settings. Can be nil, and in this case the default settings are used. If you want to specify your own settings, use the _SETTINGS object.
|
||||||
|
|||||||
@ -1697,7 +1697,7 @@ AIRBOSS.MenuF10Root=nil
|
|||||||
|
|
||||||
--- Airboss class version.
|
--- Airboss class version.
|
||||||
-- @field #string version
|
-- @field #string version
|
||||||
AIRBOSS.version="1.1.4"
|
AIRBOSS.version="1.1.5"
|
||||||
|
|
||||||
-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
|
-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
|
||||||
-- TODO list
|
-- TODO list
|
||||||
@ -4094,7 +4094,7 @@ function AIRBOSS:onafterRecoveryUnpause(From, Event, To)
|
|||||||
self:T(self.lid..string.format("Unpausing aircraft recovery."))
|
self:T(self.lid..string.format("Unpausing aircraft recovery."))
|
||||||
|
|
||||||
-- Resume recovery.
|
-- Resume recovery.
|
||||||
self:_MarshalCallRecoveryResume()
|
self:_MarshalCallResumeRecovery()
|
||||||
|
|
||||||
end
|
end
|
||||||
|
|
||||||
@ -17641,7 +17641,8 @@ function AIRBOSS:_SaveTrapSheet(playerData, grade)
|
|||||||
if self.trapprefix then
|
if self.trapprefix then
|
||||||
filename=string.format("%s_%s-%04d.csv", self.trapprefix, playerData.actype, i)
|
filename=string.format("%s_%s-%04d.csv", self.trapprefix, playerData.actype, i)
|
||||||
else
|
else
|
||||||
filename=string.format("AIRBOSS-%s_Trapsheet-%s_%s-%04d.csv", self.alias, playerData.name, playerData.actype, i)
|
local name=UTILS.ReplaceIllegalCharacters(playerData.name, "_")
|
||||||
|
filename=string.format("AIRBOSS-%s_Trapsheet-%s_%s-%04d.csv", self.alias, name, playerData.actype, i)
|
||||||
end
|
end
|
||||||
|
|
||||||
-- Set path.
|
-- Set path.
|
||||||
|
|||||||
@ -718,7 +718,8 @@ function UTILS.SecondsToClock(seconds, short)
|
|||||||
local clock=hours..":"..mins..":"..secs.."+"..days
|
local clock=hours..":"..mins..":"..secs.."+"..days
|
||||||
if short then
|
if short then
|
||||||
if hours=="00" then
|
if hours=="00" then
|
||||||
clock=mins..":"..secs
|
--clock=mins..":"..secs
|
||||||
|
clock=hours..":"..mins..":"..secs
|
||||||
else
|
else
|
||||||
clock=hours..":"..mins..":"..secs
|
clock=hours..":"..mins..":"..secs
|
||||||
end
|
end
|
||||||
@ -727,6 +728,26 @@ function UTILS.SecondsToClock(seconds, short)
|
|||||||
end
|
end
|
||||||
end
|
end
|
||||||
|
|
||||||
|
--- Seconds of today.
|
||||||
|
-- @return #number Seconds passed since last midnight.
|
||||||
|
function UTILS.SecondsOfToday()
|
||||||
|
|
||||||
|
-- Time in seconds.
|
||||||
|
local time=timer.getAbsTime()
|
||||||
|
|
||||||
|
-- Short format without days since mission start.
|
||||||
|
local clock=UTILS.SecondsToClock(time, true)
|
||||||
|
|
||||||
|
-- Time is now the seconds passed since last midnight.
|
||||||
|
return UTILS.ClockToSeconds(clock)
|
||||||
|
end
|
||||||
|
|
||||||
|
--- Cound seconds until next midnight.
|
||||||
|
-- @return #number Seconds to midnight.
|
||||||
|
function UTILS.SecondsToMidnight()
|
||||||
|
return 24*60*60-UTILS.SecondsOfToday()
|
||||||
|
end
|
||||||
|
|
||||||
--- Convert clock time from hours, minutes and seconds to seconds.
|
--- Convert clock time from hours, minutes and seconds to seconds.
|
||||||
-- @param #string clock String of clock time. E.g., "06:12:35" or "5:1:30+1". Format is (H)H:(M)M:((S)S)(+D) H=Hours, M=Minutes, S=Seconds, D=Days.
|
-- @param #string clock String of clock time. E.g., "06:12:35" or "5:1:30+1". Format is (H)H:(M)M:((S)S)(+D) H=Hours, M=Minutes, S=Seconds, D=Days.
|
||||||
-- @return #number Seconds. Corresponds to what you cet from timer.getAbsTime() function.
|
-- @return #number Seconds. Corresponds to what you cet from timer.getAbsTime() function.
|
||||||
@ -783,6 +804,15 @@ function UTILS.DisplayMissionTime(duration)
|
|||||||
MESSAGE:New(text, duration):ToAll()
|
MESSAGE:New(text, duration):ToAll()
|
||||||
end
|
end
|
||||||
|
|
||||||
|
--- Replace illegal characters [<>|/?*:\\] in a string.
|
||||||
|
-- @param #string Text Input text.
|
||||||
|
-- @param #string ReplaceBy Replace illegal characters by this character or string. Default underscore "_".
|
||||||
|
-- @return #string The input text with illegal chars replaced.
|
||||||
|
function UTILS.ReplaceIllegalCharacters(Text, ReplaceBy)
|
||||||
|
ReplaceBy=ReplaceBy or "_"
|
||||||
|
local text=Text:gsub("[<>|/?*:\\]", ReplaceBy)
|
||||||
|
return text
|
||||||
|
end
|
||||||
|
|
||||||
--- Generate a Gaussian pseudo-random number.
|
--- Generate a Gaussian pseudo-random number.
|
||||||
-- @param #number x0 Expectation value of distribution.
|
-- @param #number x0 Expectation value of distribution.
|
||||||
@ -1002,15 +1032,62 @@ function UTILS.GetDCSMap()
|
|||||||
return env.mission.theatre
|
return env.mission.theatre
|
||||||
end
|
end
|
||||||
|
|
||||||
--- Returns the mission date. This is the date the mission started.
|
--- Returns the mission date. This is the date the mission **started**.
|
||||||
-- @return #string Mission date in yyyy/mm/dd format.
|
-- @return #string Mission date in yyyy/mm/dd format.
|
||||||
|
-- @return #number The year anno domini.
|
||||||
|
-- @return #number The month.
|
||||||
|
-- @return #number The day.
|
||||||
function UTILS.GetDCSMissionDate()
|
function UTILS.GetDCSMissionDate()
|
||||||
local year=tostring(env.mission.date.Year)
|
local year=tostring(env.mission.date.Year)
|
||||||
local month=tostring(env.mission.date.Month)
|
local month=tostring(env.mission.date.Month)
|
||||||
local day=tostring(env.mission.date.Day)
|
local day=tostring(env.mission.date.Day)
|
||||||
return string.format("%s/%s/%s", year, month, day)
|
return string.format("%s/%s/%s", year, month, day), tonumber(year), tonumber(month), tonumber(day)
|
||||||
end
|
end
|
||||||
|
|
||||||
|
--- Returns the day of the mission.
|
||||||
|
-- @return #number Day of the mission. Mission starts on day 0.
|
||||||
|
function UTILS.GetMissionDay()
|
||||||
|
|
||||||
|
local time=timer.getAbsTime()
|
||||||
|
|
||||||
|
local clock=UTILS.SecondsToClock(time, false)
|
||||||
|
|
||||||
|
local x=tonumber(UTILS.Split(clock, "+")[2])
|
||||||
|
|
||||||
|
return x
|
||||||
|
end
|
||||||
|
|
||||||
|
--- Returns the current day of the year of the mission.
|
||||||
|
-- @return #number Current day of year of the mission. For example, January 1st returns 0, January 2nd returns 1 etc.
|
||||||
|
function UTILS.GetMissionDayOfYear()
|
||||||
|
|
||||||
|
local Date, Year, Month, Day=UTILS.GetDCSMissionDate()
|
||||||
|
|
||||||
|
local d=UTILS.GetMissionDay()
|
||||||
|
|
||||||
|
return UTILS.GetDayOfYear(Year, Month, Day)+d
|
||||||
|
|
||||||
|
end
|
||||||
|
|
||||||
|
--- Returns the current date.
|
||||||
|
-- @return #string Mission date in yyyy/mm/dd format.
|
||||||
|
-- @return #number The year anno domini.
|
||||||
|
-- @return #number The month.
|
||||||
|
-- @return #number The day.
|
||||||
|
function UTILS.GetDate()
|
||||||
|
|
||||||
|
-- Mission start date
|
||||||
|
local date, year, month, day=UTILS.GetDCSMissionDate()
|
||||||
|
|
||||||
|
local time=timer.getAbsTime()
|
||||||
|
|
||||||
|
local clock=UTILS.SecondsToClock(time, false)
|
||||||
|
|
||||||
|
local x=tonumber(UTILS.Split(clock, "+")[2])
|
||||||
|
|
||||||
|
local day=day+x
|
||||||
|
|
||||||
|
end
|
||||||
|
|
||||||
--- Returns the magnetic declination of the map.
|
--- Returns the magnetic declination of the map.
|
||||||
-- Returned values for the current maps are:
|
-- Returned values for the current maps are:
|
||||||
@ -1144,3 +1221,175 @@ function UTILS.GetCallsignName(Callsign)
|
|||||||
|
|
||||||
return "Ghostrider"
|
return "Ghostrider"
|
||||||
end
|
end
|
||||||
|
|
||||||
|
--- Get the time difference between GMT and local time.
|
||||||
|
-- @return #number Local time difference in hours compared to GMT. E.g. Dubai is GMT+4 ==> +4 is returned.
|
||||||
|
function UTILS.GMTToLocalTimeDifference()
|
||||||
|
|
||||||
|
local theatre=UTILS.GetDCSMap()
|
||||||
|
|
||||||
|
if theatre==DCSMAP.Caucasus then
|
||||||
|
return 4 -- Caucasus UTC+4 hours
|
||||||
|
elseif theatre==DCSMAP.PersianGulf then
|
||||||
|
return 4 -- Abu Dhabi UTC+4 hours
|
||||||
|
elseif theatre==DCSMAP.NTTR then
|
||||||
|
return -7 -- Las Vegas UTC-7 hours
|
||||||
|
elseif theatre==DCSMAP.Normandy then
|
||||||
|
return 1 -- Calais UTC+1 hour
|
||||||
|
else
|
||||||
|
BASE:E(string.format("ERROR: Unknown Map %s in UTILS.GMTToLocal function. Returning 0", tostring(theatre)))
|
||||||
|
return 0
|
||||||
|
end
|
||||||
|
|
||||||
|
end
|
||||||
|
|
||||||
|
|
||||||
|
--- Get the day of the year. Counting starts on 1st of January.
|
||||||
|
-- @param #number Year The year.
|
||||||
|
-- @param #number Month The month.
|
||||||
|
-- @param #number Day The day.
|
||||||
|
-- @return #number The day of the year.
|
||||||
|
function UTILS.GetDayOfYear(Year, Month, Day)
|
||||||
|
|
||||||
|
local floor = math.floor
|
||||||
|
|
||||||
|
local n1 = floor(275 * Month / 9)
|
||||||
|
local n2 = floor((Month + 9) / 12)
|
||||||
|
local n3 = (1 + floor((Year - 4 * floor(Year / 4) + 2) / 3))
|
||||||
|
|
||||||
|
return n1 - (n2 * n3) + Day - 30
|
||||||
|
end
|
||||||
|
|
||||||
|
--- Get sunrise or sun set of a specific day of the year at a specific location.
|
||||||
|
-- @param #number DayOfYear The day of the year.
|
||||||
|
-- @param #number Latitude Latitude.
|
||||||
|
-- @param #number Longitude Longitude.
|
||||||
|
-- @param #boolean Rising If true, calc sun rise, or sun set otherwise.
|
||||||
|
-- @param #number Tlocal Local time offset in hours. E.g. +4 for a location which has GMT+4.
|
||||||
|
-- @return #number Sun rise/set in seconds of the day.
|
||||||
|
function UTILS.GetSunRiseAndSet(DayOfYear, Latitude, Longitude, Rising, Tlocal)
|
||||||
|
|
||||||
|
-- Defaults
|
||||||
|
local zenith=90.83
|
||||||
|
local latitude=Latitude
|
||||||
|
local longitude=Longitude
|
||||||
|
local rising=Rising
|
||||||
|
local n=DayOfYear
|
||||||
|
Tlocal=Tlocal or 0
|
||||||
|
|
||||||
|
|
||||||
|
-- Short cuts.
|
||||||
|
local rad = math.rad
|
||||||
|
local deg = math.deg
|
||||||
|
local floor = math.floor
|
||||||
|
local frac = function(n) return n - floor(n) end
|
||||||
|
local cos = function(d) return math.cos(rad(d)) end
|
||||||
|
local acos = function(d) return deg(math.acos(d)) end
|
||||||
|
local sin = function(d) return math.sin(rad(d)) end
|
||||||
|
local asin = function(d) return deg(math.asin(d)) end
|
||||||
|
local tan = function(d) return math.tan(rad(d)) end
|
||||||
|
local atan = function(d) return deg(math.atan(d)) end
|
||||||
|
|
||||||
|
local function fit_into_range(val, min, max)
|
||||||
|
local range = max - min
|
||||||
|
local count
|
||||||
|
if val < min then
|
||||||
|
count = floor((min - val) / range) + 1
|
||||||
|
return val + count * range
|
||||||
|
elseif val >= max then
|
||||||
|
count = floor((val - max) / range) + 1
|
||||||
|
return val - count * range
|
||||||
|
else
|
||||||
|
return val
|
||||||
|
end
|
||||||
|
end
|
||||||
|
|
||||||
|
-- Convert the longitude to hour value and calculate an approximate time
|
||||||
|
local lng_hour = longitude / 15
|
||||||
|
|
||||||
|
local t
|
||||||
|
if rising then -- Rising time is desired
|
||||||
|
t = n + ((6 - lng_hour) / 24)
|
||||||
|
else -- Setting time is desired
|
||||||
|
t = n + ((18 - lng_hour) / 24)
|
||||||
|
end
|
||||||
|
|
||||||
|
-- Calculate the Sun's mean anomaly
|
||||||
|
local M = (0.9856 * t) - 3.289
|
||||||
|
|
||||||
|
-- Calculate the Sun's true longitude
|
||||||
|
local L = fit_into_range(M + (1.916 * sin(M)) + (0.020 * sin(2 * M)) + 282.634, 0, 360)
|
||||||
|
|
||||||
|
-- Calculate the Sun's right ascension
|
||||||
|
local RA = fit_into_range(atan(0.91764 * tan(L)), 0, 360)
|
||||||
|
|
||||||
|
-- Right ascension value needs to be in the same quadrant as L
|
||||||
|
local Lquadrant = floor(L / 90) * 90
|
||||||
|
local RAquadrant = floor(RA / 90) * 90
|
||||||
|
RA = RA + Lquadrant - RAquadrant
|
||||||
|
|
||||||
|
-- Right ascension value needs to be converted into hours
|
||||||
|
RA = RA / 15
|
||||||
|
|
||||||
|
-- Calculate the Sun's declination
|
||||||
|
local sinDec = 0.39782 * sin(L)
|
||||||
|
local cosDec = cos(asin(sinDec))
|
||||||
|
|
||||||
|
-- Calculate the Sun's local hour angle
|
||||||
|
local cosH = (cos(zenith) - (sinDec * sin(latitude))) / (cosDec * cos(latitude))
|
||||||
|
|
||||||
|
if rising and cosH > 1 then
|
||||||
|
return "N/R" -- The sun never rises on this location on the specified date
|
||||||
|
elseif cosH < -1 then
|
||||||
|
return "N/S" -- The sun never sets on this location on the specified date
|
||||||
|
end
|
||||||
|
|
||||||
|
-- Finish calculating H and convert into hours
|
||||||
|
local H
|
||||||
|
if rising then
|
||||||
|
H = 360 - acos(cosH)
|
||||||
|
else
|
||||||
|
H = acos(cosH)
|
||||||
|
end
|
||||||
|
H = H / 15
|
||||||
|
|
||||||
|
-- Calculate local mean time of rising/setting
|
||||||
|
local T = H + RA - (0.06571 * t) - 6.622
|
||||||
|
|
||||||
|
-- Adjust back to UTC
|
||||||
|
local UT = fit_into_range(T - lng_hour, 0, 24)
|
||||||
|
|
||||||
|
return floor(UT)*60*60+frac(UT)*60*60+Tlocal*60*60
|
||||||
|
end
|
||||||
|
|
||||||
|
--- Get sun rise of a specific day of the year at a specific location.
|
||||||
|
-- @param #number Day Day of the year.
|
||||||
|
-- @param #number Month Month of the year.
|
||||||
|
-- @param #number Year Year.
|
||||||
|
-- @param #number Latitude Latitude.
|
||||||
|
-- @param #number Longitude Longitude.
|
||||||
|
-- @param #boolean Rising If true, calc sun rise, or sun set otherwise.
|
||||||
|
-- @param #number Tlocal Local time offset in hours. E.g. +4 for a location which has GMT+4. Default 0.
|
||||||
|
-- @return #number Sun rise in seconds of the day.
|
||||||
|
function UTILS.GetSunrise(Day, Month, Year, Latitude, Longitude, Tlocal)
|
||||||
|
|
||||||
|
local DayOfYear=UTILS.GetDayOfYear(Year, Month, Day)
|
||||||
|
|
||||||
|
return UTILS.GetSunRiseAndSet(DayOfYear, Latitude, Longitude, true, Tlocal)
|
||||||
|
end
|
||||||
|
|
||||||
|
--- Get sun set of a specific day of the year at a specific location.
|
||||||
|
-- @param #number Day Day of the year.
|
||||||
|
-- @param #number Month Month of the year.
|
||||||
|
-- @param #number Year Year.
|
||||||
|
-- @param #number Latitude Latitude.
|
||||||
|
-- @param #number Longitude Longitude.
|
||||||
|
-- @param #boolean Rising If true, calc sun rise, or sun set otherwise.
|
||||||
|
-- @param #number Tlocal Local time offset in hours. E.g. +4 for a location which has GMT+4. Default 0.
|
||||||
|
-- @return #number Sun rise in seconds of the day.
|
||||||
|
function UTILS.GetSunset(Day, Month, Year, Latitude, Longitude, Tlocal)
|
||||||
|
|
||||||
|
local DayOfYear=UTILS.GetDayOfYear(Year, Month, Day)
|
||||||
|
|
||||||
|
return UTILS.GetSunRiseAndSet(DayOfYear, Latitude, Longitude, false, Tlocal)
|
||||||
|
end
|
||||||
Loading…
x
Reference in New Issue
Block a user