#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)