#ftime = (2013,8,16,11,30,0) times = utils.generate_times(itime,ftime,60*60) shear_times = utils.generate_times(itime,ftime,3*60*60) thresh = 10 variables = {'cref':{}, 'wind10':{}, 'buoyancy':{},'shear':{}, 'strongestwind':{}} variables['cref'] = {'lv':2000,'pt':times} variables['wind10'] = {'lv':2000,'pt':times} variables['buoyancy'] = {'lv':2000,'pt':times} variables['shear'] = {'top':3,'bottom':0,'pt':shear_times} variables['strongestwind'] = {'lv':2000, 'itime':itime, 'ftime':ftime, 'range':(thresh,27.5,1.25)} # 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 = {'wind10':{'lv':2000,'pt':times}} #shear06 = {'shear':{'top':6,'bottom':0,'pt':shear_times}} for en in ensnames: # pdb.set_trace() # Reload settings p.C = Settings() # Change paths to new location p.C.output_root = os.path.join(config.output_root,case,IC,en,experiment.keys()[0]) p.C.wrfout_root = os.path.join(config.wrfout_root,case,IC,en,experiment.keys()[0]) p.plot_2D(variables)
# ... or set just one time skewT_time = (2013, 8, 16, 3, 0, 0) # NUMBER 4: SET THE LAT/LON skewT_latlon = (35.2435, -97.4708) # NUMBER 5: # This is a loop over different folders # For one plot, just use the (last line) plot_skewT command # For Skew Ts 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) # save_output saves pickle files for later use p.plot_skewT(skewT_time, skewT_latlon, save_output=1) # For composite profiles path_to_wrfouts = [] for ens in ICens: for ex in experiments: