#!/mesonet/python/bin/python

import pg
mydb = pg.connect('postgis')

lines = open('snowfall.csv', 'r').readlines()

print lines[0]
for line in lines[1:]:
  t = line.split(",")
  city = t[1]
  lon = t[5]
  lat = t[4]
  if (len(t[4]) == 0):
    rs = mydb.query("SELECT x(the_geom) as lon, y(the_geom) as lat from \
     cities WHERE name = '%s' " % (city,)).dictresult()
    if len(rs) > 0:
      lon = rs[0]["lon"]
      lat = rs[0]["lat"]

  print "%s,%s,%s,%s,%s,%s" % (t[0], city, t[2], t[3], lat, lon)