old_contents: import psycopg2
  import datetime
  
  import network
  nt = network.Table(('IA_ASOS','MO_ASOS','IL_ASOS', 'ND_ASOS', 'AWOS',
            'WI_ASOS','MN_ASOS', 'SD_ASOS', 'NE_ASOS', 'KS_ASOS',
            'IN_ASOS','KY_ASOS','OH_ASOS','MI_ASOS', 'WFO'))
  AFOS = psycopg2.connect(database='afos', host='iemdb', user='nobody')
  acursor = AFOS.cursor()
  acursor2 = AFOS.cursor()
  
  
  
  acursor.execute("""SELECT distinct pil from products where
   substr(pil,1,3) = 'CLI' and entered > 'YESTERDAY' """)
  
  lats = []
  lons = []
  vals = []
  
  for row in acursor:
      stid = row[0][-3:]
      if stid in ['RAP', 'DVN', 'FGF', 'OAX', 'MPX']:
          continue
      if not nt.sts.has_key(row[0][-3:]):
          continue
      if nt.sts[stid]['state'] in ['WY', 'CO']:
          continue
      acursor2.execute("""
      SELECT data from products WHERE pil = %s 
      and entered > 'YESTERDAY' ORDER by entered DESC
      LIMIT 1
      """, (row[0],))
      row2 = acursor2.fetchone()
      hasSnowfall = False
      for line in row2[0].split("\n"):
          if line.find("SNOWFALL (IN)") > -1:
              hasSnowfall = True
          if line.strip() == "":
              hasSnowfall = False
          if hasSnowfall and line.find("SINCE JUL 1") > 0:
              tokens = line.split()
              if tokens[3] not in ['MM','T'] and float(tokens[3]) >= 0:
                  #print row[0], tokens
                  lats.append( nt.sts[row[0][-3:]]['lat'])
                  lons.append( nt.sts[row[0][-3:]]['lon'])
                  vals.append( float(tokens[3]) )
              break
  
  from pyiem.plot import MapPlot
  
  yesterday = datetime.datetime.today() - datetime.timedelta(days=1)
  
  m = MapPlot(sector='midwest',
         title='NWS Total Snowfall (inches) thru %s' % (yesterday.strftime("%-d %B %Y"),),
         subtitle= '1 July 2013 - %s' % (datetime.datetime.today().strftime("%-d %B %Y"),))
  m.plot_values(lons, lats, vals, fmt='%.1f')
  pqstr = "data ac %s0000 summary/mw_season_snowfall.png mw_season_snowfall.png png" % (datetime.datetime.today().strftime("%Y%m%d"),)
  m.postprocess(pqstr=pqstr)
  
new_contents: import psycopg2
  import datetime
  
  import network
  nt = network.Table(('IA_ASOS','MO_ASOS','IL_ASOS', 'ND_ASOS', 'AWOS',
            'WI_ASOS','MN_ASOS', 'SD_ASOS', 'NE_ASOS', 'KS_ASOS',
            'IN_ASOS','KY_ASOS','OH_ASOS','MI_ASOS', 'WFO'))
  AFOS = psycopg2.connect(database='afos', host='iemdb', user='nobody')
  acursor = AFOS.cursor()
  acursor2 = AFOS.cursor()
  
  
  
  acursor.execute("""SELECT distinct pil from products where
   substr(pil,1,3) = 'CLI' and entered > 'YESTERDAY' """)
  
  lats = []
  lons = []
  vals = []
  
  for row in acursor:
      stid = row[0][-3:]
      if stid in ['RAP', 'DVN', 'FGF', 'OAX', 'MPX']:
          continue
      if not nt.sts.has_key(row[0][-3:]):
          continue
      if nt.sts[stid]['state'] in ['WY', 'CO']:
          continue
      acursor2.execute("""
      SELECT data from products WHERE pil = %s 
      and entered > 'YESTERDAY' ORDER by entered DESC
      LIMIT 1
      """, (row[0],))
      row2 = acursor2.fetchone()
      hasSnowfall = False
      for line in row2[0].split("\n"):
          if line.find("SNOWFALL (IN)") > -1:
              hasSnowfall = True
          if line.strip() == "":
              hasSnowfall = False
          if hasSnowfall and line.find("SINCE JUL 1") > 0:
              tokens = line.split()
              if tokens[3] not in ['MM','T'] and float(tokens[3]) >= 0:
                  #print row[0], tokens
                  lats.append( nt.sts[row[0][-3:]]['lat'])
                  lons.append( nt.sts[row[0][-3:]]['lon'])
                  vals.append( float(tokens[3]) )
              break
  
  from pyiem.plot import MapPlot
  
  yesterday = datetime.datetime.today() - datetime.timedelta(days=1)
  
  m = MapPlot(sector='midwest', axisbg='white',
         title='NWS Total Snowfall (inches) thru %s' % (yesterday.strftime("%-d %B %Y"),),
         subtitle= '1 July 2013 - %s' % (datetime.datetime.today().strftime("%-d %B %Y"),))
  m.plot_values(lons, lats, vals, fmt='%.1f')
  pqstr = "data ac %s0000 summary/mw_season_snowfall.png mw_season_snowfall.png png" % (datetime.datetime.today().strftime("%Y%m%d"),)
  m.postprocess(pqstr=pqstr)
  
subject: make axis background white