import sys
import datetime
sys.path.append('/home/jrlawson/gitprojects/')
from WEM.postWRF.postWRF.obs import StormReports
fpath = '/home/jrlawson/storm_data_search_results.csv'
itime = datetime.datetime(2013, 8, 15, 12, 0, 0)
ftime = datetime.datetime(2013, 8, 16, 12, 0, 0)
lims = {'Nlim': 42.0, 'Elim': -94.0, 'Slim': 31.0, 'Wlim': -103.4}
fname = 'KSOK13_wind'
outdir = '/home/jrlawson/public_html/bowecho'
SR = StormReports(fpath)
SR.plot('Wind', itime, ftime, fname, outdir, **lims)
import sys
import datetime
sys.path.append('/home/jrlawson/gitprojects/')
from WEM.postWRF.postWRF.obs import StormReports
fpath = '/home/jrlawson/storm_data_search_results.csv'
itime = datetime.datetime(2013,8,15,12,0,0)
ftime = datetime.datetime(2013,8,16,12,0,0)
lims = {'Nlim':42.0,'Elim':-94.0,'Slim':31.0,'Wlim':-103.4}
fname = 'KSOK13_wind'
outdir = '/home/jrlawson/public_html/bowecho'
SR = StormReports(fpath)
SR.plot('Wind',itime,ftime,fname,outdir,**lims)
# radartimes = [datetime.datetime(2011,4,19,22,0,0),
# datetime.datetime(2011,4,20,2,30,0),
# datetime.datetime(2011,4,20,5,0,0)]
radarthresh = 25
# radarthresh = 35
radartimes = datetime.datetime(2011,4,20,3,0,0)
lims = {'Nlim':45.0,'Elim':-80.0,'Slim':32.0,'Wlim':-95.0}
fname = 'ILIN11_reports_radar.png'
outdir = '/home/jrlawson/public_html/bowecho'
# FIGURE SETUP
fig, ax = plt.subplots(1,figsize=(5,6))
p = WRFEnviron()
# Composite radar
p.plot_radar(radartimes,datadir,outdir,compthresh=radarthresh,
composite=False,fig=fig,ax=ax,**lims)
# Storm reports
SR = StormReports(fpath)
SR.plot('Wind',itime,ftime,fname,outdir,ax=ax,fig=fig,ss=8,color='magenta',**lims)
# Annotate four bow echoes
# Save fig
fig.tight_layout()
fig.savefig(os.path.join(outdir,fname))
