#!/mesonet/python/bin/python

import iemAccess, pg
mydb = pg.connect('mesosite')

dict = {}
rs = mydb.query("SELECT * from stations WHERE network = 'IA_COOP'").dictresult()
for i in range(len(rs)):
  dict[rs[i]["id"]] = rs[i]

rs = iemAccess.iemdb.query("SELECT station, snow, snowd from summary \
  WHERE network = 'IA_COOP' and day = '2005-01-05' and snow > 0").dictresult()

db = {}
for i in range(len(rs)):
  station = rs[i]["station"]
  snow = float( rs[i]["snow"] )
  snowd = float( rs[i]["snowd"] )
  db[station] = {'snow': snow, 'snowd': snowd }


rs = iemAccess.iemdb.query("SELECT station, snow, snowd from summary \
  WHERE network = 'IA_COOP' and day = '2005-01-06' and snow > 0").dictresult()

for i in range(len(rs)):
  station = rs[i]["station"]
  snow = float( rs[i]["snow"] )
  snowd = float( rs[i]["snowd"] )
  if (not db.has_key(station)):
    db[station] = {'snow': snow, 'snowd': snowd }
    continue
  db[station]['snow'] = db[station]['snow'] + snow
  db[station]['snowd'] = db[station]['snowd'] 

  id = station
  s = db[id]['snow']
  if (s < db[id]['snowd']): s = db[id]['snowd']
  print "%s,%s (%s),%s,%s,%s,%s" % ('COOP', dict[id]["name"], id, dict[id]["county"], s, dict[id]["latitude"], dict[id]["longitude"])