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)