def pickle_to_cappi(gp_pickle, ber_pickle, **kwargs):
path=kwargs.get('path', '/bm/gkeep/scollis/deal_ber/')
debug=kwargs.get('debug', False)
if debug: print "Loading Pickles"
gp=pickle_zip.load(path+gp_pickle)
ber=pickle_zip.load(path+ber_pickle)
ber_loc=[-12.457, 130.925]
gp_loc=[-12.2492, 131.0444]
displace=mathematics.corner_to_point(gp_loc, ber_loc)
ldict={'lat_0':gp_loc[0], 'lon_0':gp_loc[1],'llcrnrlat':-13.0, 'llcrnrlon':130.2, 'urcrnrlat':-12.0 , 'urcrnrlon':131.2, 'lat_ts':gp_loc[0]}
levs=linspace(500,10000, 30)
xar=linspace(-50.,50., 100)*1000.0
yar=linspace(-50.,50., 100)*1000.0
gp_cube=radar_to_cart.make_cube(gp, xar, yar, levs)
ber_cube=radar_to_cart.make_cube(ber, xar-displace[0], yar-displace[1], levs)
print gp_cube['CZ'].shape
netcdf_utis.save_data_cube(ber_cube, gp_cube, '/bm/gdata/scollis/cube_data/'+std_datestr(gp[0]['date'])+'_deal.nc', gp_loc)
def dump_pickle(pickle_name):
locs = {
"Berrimah": [-12.457, 130.925],
"C-POL": [-12.2492, 131.0444],
"Berrimah_deal": [-12.457, 130.925],
"C-POL_deal": [-12.2492, 131.0444],
}
basedir = "/bm/gkeep/scollis/deal_ber/"
radar = pickle_zip.load(basedir + pickle_name)
pres.dump_radar(radar, locs[radar[0]["radar_name"]])
def dump_pickle(pickle_name):
locs={'Berrimah':[-12.457, 130.925], 'C-POL':[-12.2492, 131.0444], 'Berrimah_deal':[-12.457, 130.925], 'C-POL_deal':[-12.2492, 131.0444]}
basedir='/bm/gdata/scollis/radar_shelves/'
radar=pickle_zip.load(basedir+pickle_name)
pres.dump_radar(radar, locs[radar[0]['radar_name']])
