cf0, cf1 = p.cold_pool_strength(t,wrf_sd=wrf_sd,out_sd=out_sd,
swath_width=130,fig=fig,axes=(ax0,ax1),dz=1)
plt.close(fig)
if spaghetti:
wrf_sds = []
for en in ensnames:
for ex in experiments:
out_sd, wrf_sd = get_folders(en,ex)
wrf_sds.append(wrf_sd)
lv = 2000
# Save to higher directory
out_d = os.path.dirname(out_sd)
for t in times:
p.spaghetti(t,lv,'cref',40,wrf_sds[:4],out_d)
if std:
wrf_sds = []
for en in ensnames:
for ex in experiments:
out_sd, wrf_sd = get_folders(en,ex)
wrf_sds.append(wrf_sd)
lv = 2000
# Save to higher directory
out_d = os.path.dirname(out_sd)
if enstype == 'ICBC':
out_d = os.path.dirname(out_d)
for t in times:
p.std(t,lv,'RH',wrf_sds,out_d,clvs=N.arange(0,26,1))
axes=(ax0, ax1),
dz=1)
plt.close(fig)
if spaghetti:
wrf_sds = []
for en in ensnames:
for ex in experiments:
out_sd, wrf_sd = get_folders(en, ex)
wrf_sds.append(wrf_sd)
lv = 2000
# Save to higher directory
out_d = os.path.dirname(out_sd)
for t in times:
p.spaghetti(t, lv, 'cref', 40, wrf_sds[:4], out_d)
if std:
wrf_sds = []
for en in ensnames:
for ex in experiments:
out_sd, wrf_sd = get_folders(en, ex)
wrf_sds.append(wrf_sd)
lv = 2000
# Save to higher directory
out_d = os.path.dirname(out_sd)
if enstype == 'ICBC':
out_d = os.path.dirname(out_d)
for t in times:
p.std(t, lv, 'RH', wrf_sds, out_d, clvs=N.arange(0, 26, 1))
