Moon Phase Script?

Since the Fantastic Weather plugins are no more I lost my moon phase device state info.
Before using the plugin I used an Applescript to determine the moon phase,
which of course is also defunct.

Hoping there might be a Python version of this script to return a numeric description of the current moon phase?

Code: Select all

tell application "IndigoServer"
   set theYear to year of (current date)
   set theMonth to month of (current date) as integer
   set theDay to day of (current date)
   
   if (theMonth = 1) then
      set theDay to theDay - 1
   else if (theMonth = 2) then
      set theDay to theDay + 30
   else
      set theDay to theDay + 28 + (theMonth - 2) * 3059 / 100
   end if
   set g to (theYear - 1900) mod 19 + 1
   set e to (11 * g + 18) mod 30
   if (((e = 25) and (g > 11)) or (e = 24)) then
      set e to e + 1
   end if
   set fullness to ((((e + theDay) * 6 + 11) mod 177) / 22) as integer
   if (fullness = 8) then set fullness to 0
   set value of variable "MoonPhase" to fullness
end tell


Thanks,

Carl

https://gist.github.com/miklb/ed145757971096565723

Looks great but how do I apply that to populating a variable?

Thanks,

Carl

Instead of the print statement, save the value to an Indigo variable. How to do that is in the Wiki.

Am I close? Cant get it to compile.

Thanks, Carl

Code: Select all

#!/usr/bin/env python

"""

moonphase.py - Calculate Lunar Phase

Author: Sean B. Palmer, inamidst.com

Cf. http://en.wikipedia.org/wiki/Lunar_phase#Lunar_phase_calculation

"""

import math, decimal, datetime

dec = decimal.Decimal


def position(now=None):

   if now is None:

      now = datetime.datetime.now()


   diff = now - datetime.datetime(2001, 1, 1)

   days = dec(diff.days) + (dec(diff.seconds) / dec(86400))

   lunations = dec("0.20439731") + (days * dec("0.03386319269"))


   return lunations % dec(1)


def phase(pos):

   index = (pos * dec(8)) + dec("0.5")

   index = math.floor(index)

   return {

      0: "New Moon",

      1: "Waxing Crescent",

      2: "First Quarter",

      3: "Waxing Gibbous",

      4: "Full Moon",

      5: "Waning Gibbous",

      6: "Last Quarter",

      7: "Waning Crescent"

   }[int(index) & 7]


def main():

   pos = position()

   phasename = phase(pos)


   roundedpos = round(float(pos), 3)

   indigo.variable.updateValue(indigo.variables[103086076], value= roundedpos)

   ##print "%s (%s)" % (phasename, roundedpos)


if __name__=="__main__":

   main()


What's with all the blank lines?

This works for me. You have to cast the value to a string before you store it in the Indigo variable. I also removed the "main" function because that doesn't work in an embedded script.

Code: Select all

import math, decimal, datetime

dec = decimal.Decimal

def position(now = None):
   if now is None:
      now = datetime.datetime.now()

   diff = now - datetime.datetime(2001, 1, 1)
   days = dec(diff.days) + (dec(diff.seconds) / dec(86400))
   lunations = dec("0.20439731") + (days * dec("0.03386319269"))
   return lunations % dec(1)

def phase(pos):
   index = (pos * dec(8)) + dec("0.5")
   index = math.floor(index)
   return {
      0: "New Moon",
      1: "Waxing Crescent",
      2: "First Quarter",
      3: "Waxing Gibbous",
      4: "Full Moon",
      5: "Waning Gibbous",
      6: "Last Quarter",
      7: "Waning Crescent"
   }[int(index) & 7]

pos = position()
phasename = phase(pos)
roundedpos = round(float(pos), 3)

indigo.variable.updateValue(indigo.variables[504685077], value = str(roundedpos))
#print(roundedpos)


Working great. Trying modify it to return a whole number without a decimal.

Got it working.

Thanks,

Carl

Code: Select all

roundedpos = int(round(float(pos), 3))