Exemplo n.º 1
0
# 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')
    
Exemplo n.º 2
0
# -*- 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)