vrbl = 'RH'
xlim = [0, 110, 10]
# vrbl = 'wind'; xlim=[0,50,5]
# Save to higher directory
ml = -2
out_d = os.path.dirname(out_sd)
if enstype == 'ICBC':
out_d = os.path.dirname(out_d)
ml = -3
for t in times:
for ln, ll in locs.items():
p.twopanel_profile(vrbl,
t,
wrf_sds,
out_d,
two_panel=1,
xlim=xlim,
ylim=[500, 1000, 50],
latlon=ll,
locname=ln,
ml=ml)
if frontogenesis:
for en in ensnames:
for ex in experiments:
outdir, ncdir = get_folders(en, ex)
for t in times:
lv = 2000
p.frontogenesis(
t,
lv,
ncdir=ncdir,
# locs = {'KTOP':(39.073,-95.626),'KOAX':(41.320,-96.366),'KOUN':(35.244,-97.471)}
locs = {'KAMA':(35.2202,-101.7173)}
lv = 2000
vrbl = 'RH'; xlim=[0,110,10]
# vrbl = 'wind'; xlim=[0,50,5]
# Save to higher directory
ml = -2
out_d = os.path.dirname(out_sd)
if enstype == 'ICBC':
out_d = os.path.dirname(out_d)
ml = -3
for t in times:
for ln,ll in locs.iteritems():
p.twopanel_profile(vrbl,t,wrf_sds,out_d,two_panel=1,
xlim=xlim,ylim=[500,1000,50],
latlon=ll,locname=ln,ml=ml)
if frontogenesis:
for en in ensnames:
for ex in experiments:
outdir, ncdir = get_folders(en,ex)
for t in times:
lv = 2000
p.frontogenesis(t,lv,ncdir=ncdir,outdir=outdir,
clvs=N.arange(-2.0,2.125,0.125)*10**-7,
# clvs = N.arange(-500,510,10)
smooth=3, cmap='bwr',cb='only'
)
