locs = {'Norman': (35.2, -97.4)}
if plot2D or radarplot or strongestwind:
if ens:
outdir = os.path.join(outroot, 'd0{0}'.format(dom), ens)
ncdir = os.path.join(ncroot, ens)
else:
ncdir = ncroot
outdir = os.path.join(outroot, 'd0{0}'.format(dom))
if strongestwind:
p.plot_strongest_wind(iwind,
fwind,
2000,
ncdir=ncdir,
nct=nct,
outdir=outdir,
clvs=windlvs,
dom=dom,
cmap='jet',
cb=True)
for t in times:
if plot2D:
# p.plot2D('Z',t,500,wrf_sd=wrf_sd,out_sd=out_sd,plottype='contour',smooth=10)
# p.plot2D('Td2',t,ncdir=ncdir,outdir=outdir,nct=t,match_nc=matchnc,clvs=N.arange(260,291,1))
# p.plot2D('Q2',t,ncdir=ncdir,outdir=outdir,nct=t,match_nc=matchnc,clvs=N.arange(1,20.5,0.5)*10**-3)
# p.plot2D('RAINNC',t,ncdir=wrf_sd,outdir=out_sd,locations=locs,clvs=N.arange(1,100,2))
# p.plot2D('fluidtrapping',t,ncdir=ncroot,nct=nct,outdir=outdir,cb=True,dom=dom,clvs=N.arange(-1,1,0.1)*10**-6)
# p.plot2D('lyapunov',t,700,ncdir=ncdir,nct=nct,outdir=outdir,cb=True,dom=dom,clvs=N.arange(-7.5,8.0,0.5)*10**-3,cmap='bwr')
# print(ncdir)
# p.plot2D('WSPD10MAX',t,ncdir=ncdir,nct=nct,outdir=outdir,cb=True,dom=dom,clvs=N.arange(10,31,1))
if upperlevel:
for en in ensnames:
for ex in experiments:
out_sd, wrf_sd = get_folders(en,ex)
for time in times:
p.upperlevel_W(time,850,wrf_sd=wrf_sd,out_sd=out_sd,
clvs = N.arange(0,1.0,0.01)
)
windlvs = N.arange(10,31,1)
if strongestwind:
for en in ensnames:
for ex in experiments:
outdir, ncdir = get_folders(en,ex)
# p.plot_strongest_wind(iwind,fwind,2000,ncdir,outdir,clvs=windlvs)
p.plot_strongest_wind(iwind,fwind,2000,ncdir=ncdir,nct=nct,outdir=outdir,clvs=windlvs,dom=1,
cmap='jet',cb=True)
if accum_rain:
for en in ensnames:
for ex in experiments:
for t in times:
outdir, ncdir = get_folders(en,ex)
p.plot_accum_rain(t,6,ncdir,outdir,
Nlim=42.7,Elim=-94.9,Slim=37.0,Wlim=-101.8
)
if compute_dte or plot_3D_dte or plot_1D_dte or powerspectrum:
pfname = 'DTE_' + enstype
ofname = enstype
pickledir,outdir = get_pickle_dirs(ensnames[0])
path_to_wrfouts = []
