locs = {'KTOP':(39.073,-95.626),'KOAX':(41.320,-96.366),
'KOUN':(35.244,-97.471)}
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:
out_sd, wrf_sd = get_folders(en,ex)
for time in times:
p.frontogenesis(time,850,wrf_sd=wrf_sd,out_sd=out_sd,
clvs=N.arange(-5.0,5.25,0.25)*10**-11
# clvs = N.arange(0,1.3,0.01)*10**-3
)
if upperlevel:
for en in ensnames:
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:
out_sd, wrf_sd = get_folders(en, ex)
for time in times:
p.frontogenesis(
time,
925,
nc_sd=wrf_sd,
nc_init=inittime,
