ylimits = [0, 2e8]
outdir = os.path.join(outroot, case)
p.all_error_growth(outdir, EXS, f_prefix='DTE', ylim=ylimits)
if delta_plot:
for t in times:
# nc1
IC = 'GEFSR2'
en = 'p09'
enstype = 'STCH'
ex = 's01'
MP = 'ICBC'
outdir, ncdir1 = get_folders(en, ex)
# nc2
IC = 'GEFSR2'
en = 'p09'
ex = 'ICBC'
enstype = 'ICBC'
MP = 'ICBC'
outdir, ncdir2 = get_folders(en, ex)
clvs = N.arange(-9, 10, 1)
p.plot_delta('wind10',
t,
ncdir1=ncdir1,
ncdir2=ncdir2,
outdir=outdir,
cb=True,
clvs=clvs)
if delta_plot:
for t in times:
# nc1
IC = 'GEFSR2'
en = 'p09'
enstype = 'STCH'
ex = 's01'
MP = 'ICBC'
outdir, ncdir1 = get_folders(en,ex)
# nc2
IC = 'GEFSR2'
en = 'p09'
ex = 'ICBC'
enstype = 'ICBC'
MP = 'ICBC'
outdir, ncdir2 = get_folders(en,ex)
clvs = N.arange(-9,10,1)
p.plot_delta('wind10',t,ncdir1=ncdir1,ncdir2=ncdir2,outdir=outdir,cb=True,clvs=clvs,)
if probability:
ensemble = create_ensemble()
outdir2, ncdir = get_folders(ensnames[0],experiments[0])
outdir = os.path.dirname(outdir2)
# outdir = '/home/jrlawson/public_html/bowecho/{0}/{1}/'.format(case,IC)
smooth = 'maxfilter'
thresh = 20
p.probability_threshold(ensemble,'wind10','over',thresh,iwind,fwind,smooth=smooth,outdir=outdir,level=False)
