# variables = {'cref':{'lv':2000,'pt':times}, 'shear':{'pt':shear_times, 'top':3, 'bottom':0}}
#variables = {'PMSL':{'lv':2000,'pt':times,'plottype':'contour','smooth':5}}
#shear06 = {'shear':{'top':6,'bottom':0,'pt':shear_times}}
if skewT:
for en in ensnames:
for ex in experiments:
# Reload settings
p.C = Settings()
# Change paths to new location
p.C.output_root = os.path.join(config.output_root,case,IC,en,MP,ex)
p.C.wrfout_root = os.path.join(config.wrfout_root,case,IC,en,MP,ex)
p.C.pickledir = os.path.join(config.wrfout_root,case,IC,en,MP,ex)
p.plot_skewT(skewT_time,skewT_latlon,save_output=1)
#p.plot_skewT(skewT_time,skewT_latlon,composite=1)
if plot2D:
# p.plot_strongest_wind(itime,ftime,2000,wrf_sd=wrf_sd,out_sd=out_sd)
p.plot2D('olr',times,levels)
if streamlines:
for t in times:
p.plot_streamlines(2000,t)
if rucplot:
# RUC file is one-per-time so .nc file is specified beforehand
en = ensnames[0]
RC = Settings()
skewT_time = (2006, 5, 27, 0, 0, 0)
skewT_latlon = (36.73, -102.51) # Boise City, OK
if skewT:
for en in ensnames:
for ex in experiments:
# Reload settings
p.C = Settings()
# Change paths to new location
out_sd, wrf_sd = get_folders(en, ex)
p.plot_skewT(skewT_time,
skewT_latlon,
out_sd=out_sd,
wrf_sd=wrf_sd,
save_output=0)
#p.plot_skewT(skewT_time,skewT_latlon,composite=1)
if plot2D:
for en in ensnames:
for ex in experiments:
for t in times:
out_sd, wrf_sd = get_folders(en, ex)
# p.plot_strongest_wind(itime,ftime,2000,wrf_sd=wrf_sd,out_sd=out_sd)
# p.plot2D('Z',t,500,wrf_sd=wrf_sd,out_sd=out_sd,plottype='contour',smooth=10)
p.plot2D('cref', t, 2000, wrf_sd=wrf_sd, out_sd=out_sd)
if streamlines:
# variables = {'cref':{'lv':2000,'pt':times}, 'shear':{'pt':shear_times, 'top':3, 'bottom':0}}
#variables = {'PMSL':{'lv':2000,'pt':times,'plottype':'contour','smooth':5}}
#shear06 = {'shear':{'top':6,'bottom':0,'pt':shear_times}}
if skewT:
for en in ensnames:
for ex in experiments:
# Reload settings
p.C = Settings()
# Change paths to new location
p.C.output_root = os.path.join(config.output_root,case,IC,en,MP,ex)
p.C.wrfout_root = os.path.join(config.wrfout_root,case,IC,en,MP,ex)
p.C.pickledir = os.path.join(config.wrfout_root,case,IC,en,MP,ex)
p.plot_skewT(skewT_time,skewT_latlon,save_output=1)
#p.plot_skewT(skewT_time,skewT_latlon,composite=1)
if plot2D:
for en in ensnames:
for ex in experiments:
out_sd, wrf_sd = get_folders(en,ex)
# p.plot_strongest_wind(itime,ftime,2000,wrf_sd=wrf_sd,out_sd=out_sd)
p.plot2D('cref',times,levels,wrf_sd=wrf_sd,out_sd=out_sd)
if streamlines:
for en in ensnames:
# Reload settings
p.C = Settings()
#shear06 = {'shear':{'top':6,'bottom':0,'pt':shear_times}}
skewT_time = (2006,5,27,0,0,0)
skewT_latlon = (36.73,-102.51) # Boise City, OK
if skewT:
for en in ensnames:
for ex in experiments:
# Reload settings
p.C = Settings()
# Change paths to new location
out_sd, wrf_sd = get_folders(en,ex)
p.plot_skewT(skewT_time,skewT_latlon,out_sd=out_sd,wrf_sd=wrf_sd,save_output=0)
#p.plot_skewT(skewT_time,skewT_latlon,composite=1)
if plot2D:
for en in ensnames:
for ex in experiments:
for t in times:
out_sd, wrf_sd = get_folders(en,ex)
# p.plot_strongest_wind(itime,ftime,2000,wrf_sd=wrf_sd,out_sd=out_sd)
# p.plot2D('Z',t,500,wrf_sd=wrf_sd,out_sd=out_sd,plottype='contour',smooth=10)
p.plot2D('cref',t,2000,wrf_sd=wrf_sd,out_sd=out_sd)
if streamlines:
for en in ensnames:
for ex in experiments:
