def Plot_Station_Mean():

    lons = load('%s/lons_good_stations' % workspace)
    lats = load('%s/lats_good_stations' % workspace)

    avg = load('%s/6var_avg_good_station' % (workspace))
    vars = ['Pan', 'Tair', 'Rn', 'Wind', 'Humidity', 'VPD']
    units = ['(mm)', r'($^oC$)', '(MJ/m2/d)', '(m/s)', '(hpa)', '(hpa)']

    for i in xrange(0, 6):
        Plotting.Mapshow(avg[:, i], lons, lats, 30, 0.7, None, None, 'jet',
                         vars[i], None, units[i], None, None)
        savefig('%s/map_%s_avg_good_station.tif' % (figdir, vars[i]), dpi=200)
        plt.clf()
    return
def Plot_Selected_Station_Locations(j):
	lons = load('%s/lons%s' %(workspace, station_list_names[j]))
	lats = load('%s/lats%s' %(workspace, station_list_names[j]))
	Plotting.Mapshow(cols[j], lons, lats, 45, 0.5, None, None, None, 'Qualified pan evaporation met stations (%s)' %len(lons), 'Met Stations', '', figdir, '%s_stations_scattermap.png' %fignames[j])

	return