Astronomical clock

[admin]

Only analog (numeric) and binary values are supported by BACnet swrver.

[Erwin van der Zwart]

Hi! Great script, thanks!

I'm having an issue, trying to export the 250 byte strings for sunrise and sundown and calculated sunrise and sundown, but it won't let me export those objects.
Also I did some adjustment to the script and byte the same times as 3 byte date / time but I still can't export to BACnet. Any reason it won't let me?

Do you need this values in a SE (EBO) system? If yes, i have ported this function to Plain English so you have all values as BACnet datetime objects there..

[asgeirsb]

Hi! Great script, thanks!

I'm having an issue, trying to export the 250 byte strings for sunrise and sundown and calculated sunrise and sundown, but it won't let me export those objects.
Also I did some adjustment to the script and byte the same times as 3 byte date / time but I still can't export to BACnet. Any reason it won't let me?

Do you need this values in a SE (EBO) system? If yes, i have ported this function to Plain English so you have all values as BACnet datetime objects there..

Actually, yes! Would love that

[Erwin van der Zwart]

Here you go (:

Created for EBO 2.0 and higher 

  Sun Tracking & Astronomic Clock Calculation EZT.zip (Size: 17.86 KB / Downloads: 28)

   

[asgeirsb]

Here you go (:

Created for EBO 2.0 and higher 

Thank you very much!

[asgeirsb]

Hello again.

Is it possible to transfer time and/or dates through BACnet?

I need to send sunset and sunrise times to an atvise server through bacnet but I am struggeling

[admin]

Date/time BACnet type is not supported. A possible workaround is to split hours and minutes into separate numeric objects.

[Domoticatorino]

Hi Josep,


Make sure this function is inside common functions:

Code:

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465

-- sunrise / sunset calculation function rscalc(latitude, longitude, when)  local pi = math.pi  local doublepi = pi * 2  local rads = pi / 180.0  local TZ = function(when)    local ts = os.time(when)    local utcdate, localdate = os.date('!*t', ts), os.date('*t', ts)    localdate.isdst = false    local diff = os.time(localdate) - os.time(utcdate)    return math.floor(diff / 3600)  end  local range = function(x)    local a = x / doublepi    local b = doublepi * (a - math.floor(a))    return b < 0 and (doublepi + b) or b  end  when = when or os.date('*t')  local y2k = { year = 2000, month = 1, day = 1 }  local y2kdays = os.time(when) - os.time(y2k)  y2kdays = math.ceil(y2kdays / 86400)  local meanlongitude = range(280.461 * rads + 0.9856474 * rads * y2kdays)  local meananomaly = range(357.528 * rads + 0.9856003 * rads * y2kdays)  local lambda = range(meanlongitude + 1.915 * rads * math.sin(meananomaly) + rads / 50 * math.sin(2 * meananomaly))  local obliq = 23.439 * rads - y2kdays * rads / 2500000  local alpha = math.atan2(math.cos(obliq) * math.sin(lambda), math.cos(lambda))  local declination = math.asin(math.sin(obliq) * math.sin(lambda))  local LL = meanlongitude - alpha  if meanlongitude < pi then    LL = LL + doublepi  end  local dfo = pi / 216.45  if latitude < 0 then    dfo = -dfo  end  local fo = math.min(math.tan(declination + dfo) * math.tan(latitude * rads), 1)  local ha = 12 * math.asin(fo) / pi + 6  local timezone = TZ(when)  local equation = 12 + timezone + 24 * (1 - LL / doublepi) - longitude / 15  local sunrise, sunset = equation - ha, equation + ha  if sunrise > 24 then    sunrise = sunrise - 24  end  if sunset > 24 then    sunset = sunset - 24  end  return math.floor(sunrise * 60), math.ceil(sunset * 60) end

Then create a resident script with time 60 seconds and add this code:

Code:

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677787980818283

-- ***** Astronomic clock module Version 1.0 Created by Erwin van der Zwart ***** -- -- ********* Use this script as resident script and run every 60 seconds ******** -- -- ****** This script makes use of rscalc function inside common functions ****** -- -- ******************************** SET PARAMETERS ****************************** -- -- Set output address adressoutput = '1/1/1' -- 1 bit output for control object -- Set time shifting addresses addressshifttimesunrice = '1/1/2' -- 1 or 2 byte signed integer to be able to adjust/offset calculated sunrise time (set min/max on object) addressshifttimesundown = '1/1/3'-- 1 or 2 byte signed integer to be able to adjust/offset calculated down time (set min/max on object) -- Set feedback adresses addresssunrise = '1/1/4' -- 255 byte object for user feedback, shows calculated sunrise time addresssundown = '1/1/5' -- 255 byte object for user feedback, shows calculated sundown time addresssunriseshifted = '1/1/6' -- 255 byte object for user feedback, shows with offset calculated sunrise time addressundownshifted = '1/1/7' -- 255 byte object for user feedback, shows with offset calculatedsundown time -- Set location coordinates latitude = 52.52198 --Zwolle Netherlands, Look at 'http://www.mapcoordinates.net/en' for other locations longitude = 6.08213 --Zwolle Netherlands, Look at 'http://www.mapcoordinates.net/en' for other locations -- ******************************** END PARAMETERS ****************************** -- -- ********************* DON'T CHANGE ANYTHING UNDER THIS LINE ****************** -- -- Function to calculate switching times based and output value function astronomic(latitude,longitude,shiftedsunrise,shiftedsundown)  local now = os.date('*t')  local sunrise, sunset = rscalc(latitude, longitude, now) --Zwolle Netherlands, See http://www.mapcoordinates.net/en for other locations  sunriseorg = string.format('%d:%.2d', math.floor(sunrise / 60), sunrise % 60)  sunsetorg = string.format('%d:%.2d', math.floor(sunset / 60), sunset % 60)  sunrisenew = sunrise + shiftedsunrise  sunsetnew = sunset + shiftedsundown  sunrisetime = string.format('%d:%.2d', math.floor(sunrisenew / 60), sunrisenew % 60)  sunsettime = string.format('%d:%.2d', math.floor(sunsetnew / 60), sunsetnew % 60)  local minutes = now.hour * 60 + now.min  nowdark = minutes < sunrisenew or sunsetnew < minutes  local t = {sunricenormal = sunriseorg, sundownnormal = sunsetorg, sunriceshifted = sunrisetime, sundownshifted = sunsettime, output = nowdark}    return (t) end -- Get values for timeshifting shiftedsunrise = grp.getvalue(addressshifttimesunrice) or 0 -- Make minutes negative to switch earlyer than normal calculated sunrice and positive to switch later than normal calculated sunrice shiftedsundown = grp.getvalue(addressshifttimesundown) or 0 -- Make minutes negative to switch earlyer than normal calculated sundown and positive to switch later than normal calculated sundown -- Call function to calculate astroinfo astroinfo = astronomic(latitude,longitude,shiftedsunrise,shiftedsundown) -- Set new sunrice time (unshifted) currentvalue = grp.getvalue(addresssunrise) if astroinfo.sunricenormal ~= currentvalue then    grp.update(addresssunrise, astroinfo.sunricenormal) end -- Set new sundown time (unshifted) currentvalue = grp.getvalue(addresssundown) if astroinfo.sundownnormal ~= currentvalue then    grp.update(addresssundown, astroinfo.sundownnormal) end -- Set new sunrice time (shifted) currentvalue = grp.getvalue(addresssunriseshifted) if astroinfo.sunriceshifted ~= currentvalue then    grp.update(addresssunriseshifted, astroinfo.sunriceshifted) end -- Set new sundown time (shifted) currentvalue = grp.getvalue(addressundownshifted) if astroinfo.sundownshifted ~= currentvalue then  grp.update(addressundownshifted, astroinfo.sundownshifted) end -- Set output value to output object currentvalue = grp.getvalue(adressoutput) if astroinfo.output == true then  if currentvalue == false then    grp.write(adressoutput, true)  end else  if currentvalue == true then    grp.write(adressoutput, false)  end end

If you want to combine it with a scheduler then make a 'and' logic, put the output of this script to input 1, output of the schedular to input 2 and the control object to the output of your 'and' logic.

BR,

Erwin

Hi Erwin,
very easy question. 

adressoutput = '7/0/4' -- 1 bit output for control object

It gives day or night?

0 for day and 1 for night?

Thanks.

BR

[Erwin van der Zwart]

I have no idea anymore, this code is so old and not used anymore from the day that schedulers has the option to set on sunrise/sunset.

Given the variable is called “nowdark” i assume it’s true when dark but i strongly advise to use the build-in method on the schedulers..

BR,

Erwin