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