case = '20130815'
#ensembletype = 'STCH'
ensembletype = 'MXMP'
IC = 'GEFSR2'
ICens = [
'p09',
]
#MP = 'Morrison_Hail'
MP = 'ICBC'
# NUMBER 3: SET THE TIME(s) YOU WANT PLOT
# Loop over times...
itime = (2013, 8, 15, 18, 0, 0)
ftime = (2013, 8, 16, 12, 0, 0)
times = utils.generate_times(itime, ftime, 60 * 60)
# ... 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 en in ensnames:
for ex in experiments:
# Reload settings
p.C = Settings()
#ensnames = ['c00'] + ['p'+"%02d" %n for n in range(1,11)]
ensnames = ['p'+"%02d" %n for n in range(8,11)]
#ensnames = ['p04']
#ensnames = ['anl']
experiment = {'ICBC':'CTRL'}
itime = (2006,5,26,22,0,0)
#itime = (2009,9,10,23,0,0)
#itime = (2011,4,19,18,0,0)
#itime = (2013,8,15,18,0,0)
ftime = (2006,5,27,11,0,0)
#ftime = (2009,9,11,14,0,0)
#ftime = (2011,4,20,10,30,0)
#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}}
if wy == 2006 or wy == 2013:
radar_view = "cent_plains"
elif wy == 2009:
radar_view = "north_plains"
elif wy == 2011:
radar_view = "upper_missvly"
elif plot == "shear03":
va_str = "shear"
lv_str = "0to3"
elif plot == "shear06":
va_str = "shear"
lv_str = "0to6"
times = utils.generate_times(itime, ftime, 60 * 60 * 1)
elif plot == "maxwind":
va_str = "strongestwind"
lv_str = "sfc"
elif plot == "streamlines":
va_str = "streamlines"
lv_str = "sfc"
if plot != "maxwind":
time_strs = get_times(itime, ftime, int_hr)
else:
time_strs = ["range"]
if plot == "streamlines":
p = WRFEnviron(config)
if case[:4] == '2006':
itime = (2006, 5, 26, 0, 0, 0)
ftime = (2006, 5, 27, 15, 0, 0)
elif case[:4] == '2009':
itime = (2009, 9, 10, 0, 0, 0)
ftime = (2009, 9, 11, 18, 0, 0)
elif case[:4] == '2011':
itime = (2011, 4, 19, 0, 0, 0)
ftime = (2011, 4, 20, 15, 0, 0)
else:
itime = (2013, 8, 15, 0, 0, 0)
ftime = (2013, 8, 16, 15, 0, 0)
times = utils.generate_times(itime, ftime, 3 * 3600)
if experiment == 'STCH':
picklefolder = os.path.join(config.wrfout_root, case, IC, ens, MP)
p.C.output_root = os.path.join(config.output_root, case, IC, ens, MP)
path_to_wrfouts = []
for s in stoch_names:
fpath = os.path.join(config.wrfout_root, case, IC, ens, MP, s)
path_to_wrfouts.append(p.wrfout_files_in(fpath, dom=1)[0])
sensitivity = 0
elif experiment == 'ICBC':
picklefolder = os.path.join(config.wrfout_root, case, IC)
p.C.output_root = os.path.join(config.output_root, case, IC)
sensitivity = ens_names
#ensnames = ['c00'] + ['p'+"%02d" %n for n in range(1,11)]
ensnames = ['p' + "%02d" % n for n in range(8, 11)]
#ensnames = ['p04']
#ensnames = ['anl']
experiment = {'ICBC': 'CTRL'}
itime = (2006, 5, 26, 22, 0, 0)
#itime = (2009,9,10,23,0,0)
#itime = (2011,4,19,18,0,0)
#itime = (2013,8,15,18,0,0)
ftime = (2006, 5, 27, 11, 0, 0)
#ftime = (2009,9,11,14,0,0)
#ftime = (2011,4,20,10,30,0)
#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'] = {
p = WRFEnviron(config)
if case[:4] == '2006':
itime = (2006,5,26,0,0,0)
ftime = (2006,5,27,15,0,0)
elif case[:4] == '2009':
itime = (2009,9,10,0,0,0)
ftime = (2009,9,11,18,0,0)
elif case[:4] == '2011':
itime = (2011,4,19,0,0,0)
ftime = (2011,4,20,15,0,0)
else:
itime = (2013,8,15,0,0,0)
ftime = (2013,8,16,15,0,0)
times = utils.generate_times(itime,ftime,3*3600)
if experiment=='STCH':
picklefolder = os.path.join(config.wrfout_root,case,IC,ens,MP)
p.C.output_root = os.path.join(config.output_root,case,IC,ens,MP)
path_to_wrfouts = []
for s in stoch_names:
fpath = os.path.join(config.wrfout_root,case,IC,ens,MP,s)
path_to_wrfouts.append(p.wrfout_files_in(fpath,dom=1)[0])
sensitivity=0
elif experiment=='ICBC':
picklefolder = os.path.join(config.wrfout_root,case,IC)
p.C.output_root = os.path.join(config.output_root,case,IC)
sensitivity=ens_names
