p.C.output_root = os.path.join(config.output_root, case, IC, ICens[0], MP)
for ens in ICens:
for ex in experiments:
fpath = os.path.join(config.wrfout_root, case, IC, ens, MP, ex)
path_to_wrfouts.append(p.wrfout_files_in(fpath, dom=1)[0])
elif ensembletype == 'MXMP':
experiments = ('ICBC', 'WSM6_Grau', 'WSM6_Hail', 'Kessler', 'Ferrier',
'WSM5', 'WDM5', 'Lin', 'WDM6_Grau', 'WDM6_Hail',
'Morrison_Grau', 'Morrison_Hail')
p.C.output_root = os.path.join(config.output_root, case, IC, ICens[0])
for ens in ICens:
for ex in experiments:
fpath = os.path.join(config.wrfout_root, case, IC, ens, ex)
path_to_wrfouts.append(p.wrfout_files_in(fpath, dom=1)[0])
ylim = 0
va = 'theta'
xlim = (298, 320, 4)
ylim = (1000, 500, 50)
#va = 'RH' ; xlim = (0,100,5)
#va = 'wind' ; xlim = (0,35,5)
p.composite_profile(va,
skewT_time,
skewT_latlon,
path_to_wrfouts,
mean=1,
std=1,
xlim=xlim,
ylim=ylim)
# save_output saves pickle files for later use
p.plot_skewT(skewT_time,skewT_latlon,save_output=1)
"""
path_to_wrfouts = []
if ensembletype == 'STCH':
experiments = ['s'+"%02d" %n for n in range(1,11)]
p.C.output_root = os.path.join(config.output_root,case,IC,ICens[0],MP)
for ens in ICens:
for ex in experiments:
fpath = os.path.join(config.wrfout_root,case,IC,ens,MP,ex)
path_to_wrfouts.append(p.wrfout_files_in(fpath,dom=1)[0])
elif ensembletype == 'MXMP':
experiments = ('ICBC','WSM6_Grau','WSM6_Hail','Kessler','Ferrier','WSM5',
'WDM5','Lin','WDM6_Grau','WDM6_Hail',
'Morrison_Grau','Morrison_Hail')
p.C.output_root = os.path.join(config.output_root,case,IC,ICens[0])
for ens in ICens:
for ex in experiments:
fpath = os.path.join(config.wrfout_root,case,IC,ens,ex)
path_to_wrfouts.append(p.wrfout_files_in(fpath,dom=1)[0])
ylim = 0
va = 'theta' ; xlim = (298,320,4) ; ylim=(1000,500,50)
#va = 'RH' ; xlim = (0,100,5)
#va = 'wind' ; xlim = (0,35,5)
p.composite_profile(va,skewT_time,skewT_latlon,path_to_wrfouts,mean=1,std=1,xlim=xlim,ylim=ylim)
