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)