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']])