# Select locations for plots
loclist = ['loc'+str(n) for n in range(1,10)]
# Grid labels for domains
gl = (r'12\,km',r'4\,km',r'1.3\,km')
gl3 = (r'3\,km')
# Dictionaries for cross-section:
XS = {}
XS['zonalwind'] = {'startll': (40.81, -96.7), 'endll': (40.81, -93.65),'vars':('uwind','uwind_sr','W','pertpres','qvapor','theta'),'storm_motion':(14,270)} # Lincoln to S of DMX
# ^^ Storm motion as if it's wind ^^
if domains_and_nests:
# Send cross-section data to plot on domains?
mapWRF.mapDomains(gl3,terrain=False,drawxs=XS)
#mapWRF.mapTerrain(includePMSL=True)
if xz_cross_sections:
# xzWRF needs time as a six-element tuple:
# (year,month,day,hour,min,sec)
#for t in xs_time: # wind has: normal, parallel, magnitude options.
xzWRF_MCS.xzPlot(nest=1, datetuple=xs_time, XS=XS)
if skew_T:
skewT.skewTPlot(nest=1, timetuple=skewT_time)
if tz_cross_sections:
tzWRF.tzCloudPlot(nest=3) # JRL: need to fix WRF time step issue
if time_series:
for l in loclist:
ts.timeseries(nest=3,loc=l)
if xz_animation:
print "Creating loop."
os.system('convert -delay 50 -loop 0 @'+outdir+'xz_timelist.txt '+outdir+naming+'loop.gif')
# -*- coding:utf-8 -*- """ @author Geir Arne Waagbø @see http://code.google.com/p/pywrfplot/ """ import mapWRF import xzWRF import tzWRF import skewT mapWRF.mapDomains(includePMSL=True) xzWRF.xzCloudPlot(nest=4, time=45, plotTemp=True, plotRH=False) skewT.skewTPlot(nest=4, time=45) tzWRF.tzCloudPlot(nest=4)