import sys
import numpy as N
sys.path.append('/home/jrlawson/gitprojects/WEM')
from WEM.postWRF.postWRF import WRFEnviron
# Initialise the environment
p = WRFEnviron()
# output directory
outdir = '/home/jrlawson/public_html/bowecho'
# netcdf directory.
ncdir = '/chinook2/jrlawson/bowecho/20130815/GEFSR2/c00/ICBC'
# initial and final times.
itime = (2013, 8, 15, 3, 0, 0)
ftime = (2013, 8, 15, 3, 0, 0)
hourly = 6
times = p.generate_times(itime, ftime, hourly * 60 * 60)
# Loop over all times.
# clvs lets you change the contour values
# The second argument is accumulation time, in hours
for t in times:
p.plot_accum_rain(t, hourly, ncdir, outdir, clvs=N.arange(5, 85, 1))
p = WRFEnviron()
outdir = ('/home/jrlawson/public_html/bowecho/'
'20110419/GEFSR2/p04/WDM6_Hail/test')
ncdir = ('/chinook2/jrlawson/bowecho/'
'20110419/GEFSR2/p04/WDM6_Hail')
plot2D = 1
streamlines = 0
# Create list of times
itime = (2013,8,16,1,0,0)
ftime = (2013,8,16,5,0,0)
hourly = 1
times = p.generate_times(itime,ftime,hourly*60*60)
levels = (2000,)
vrbls = ('cref',)
# import pdb; pdb.set_trace()
# Plot standard 2D plots
if plot2D:
for t in times:
for lv in levels:
for v in vrbls:
p.plot2D(v,t,lv,ncdir,outdir)
# Streamline 2D plots
if streamlines:
for t in times:
for lv in levels:
