#variables = {'strongestwind':{}} ; variables['strongestwind'] = {'lv':2000, 'itime':itime, 'ftime':ftime, 'range':(thresh,27.5,1.25)}
#variables['PMSL'] = {'lv':2000,'pt':times,'plottype':'contour','smooth':5}
#variables['wind10'] = {'lv':2000,'pt':times}
#variables['buoyancy'] = {'lv':2000,'pt':times}
#variables['shear'] = {'top':3,'bottom':0,'pt':shear_times}
# variables = {'shear':{'pt':shear_times, 'top':3, 'bottom':0}}
# variables = {'thetae':{'lv':2000,'pt':times}, 'CAPE':{'pt':times}}
# 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:
#variables = {'strongestwind':{}} ; variables['strongestwind'] = {'lv':2000, 'itime':itime, 'ftime':ftime, 'range':(thresh,27.5,1.25)}
#variables['PMSL'] = {'lv':2000,'pt':times,'plottype':'contour','smooth':5}
#variables['wind10'] = {'lv':2000,'pt':times}
#variables['buoyancy'] = {'lv':2000,'pt':times}
#variables['shear'] = {'top':3,'bottom':0,'pt':shear_times}
# variables = {'shear':{'pt':shear_times, 'top':3, 'bottom':0}}
# variables = {'thetae':{'lv':2000,'pt':times}, 'CAPE':{'pt':times}}
# 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)
