コード例 #1
0
                        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')

    if upperlevel:
        for en in ensnames:
            for ex in experiments:
                out_sd, wrf_sd = get_folders(en, ex)
                for time in times:
                    p.upperlevel_W(time,
                                   850,
                                   wrf_sd=wrf_sd,
                                   out_sd=out_sd,
                                   clvs=N.arange(0, 1.0, 0.01))

    # if strongestwind:
    # for en in ensnames:
    # for ex in experiments:
    # outdir, ncdir = get_folders(en,ex)
    # p.plot_strongest_wind(iwind,fwind,2000,ncdir,outdir,clvs=windlvs)

    if accum_rain:
        for en in ensnames:
            for ex in experiments:
                for t in times:
                    outdir, ncdir = get_folders(en, ex)
                    p.plot_accum_rain(t,
コード例 #2
0
ファイル: bowecho_plot.py プロジェクト: NicWayand/WEM
            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'
                                )

if upperlevel:
    for en in ensnames:
        for ex in experiments:
            out_sd, wrf_sd = get_folders(en,ex)
            for time in times: 
                p.upperlevel_W(time,850,wrf_sd=wrf_sd,out_sd=out_sd,
                                clvs = N.arange(0,1.0,0.01)
                                )

windlvs = N.arange(10,31,1)
if strongestwind:
    for en in ensnames:
        for ex in experiments:
            outdir, ncdir = get_folders(en,ex)
            # p.plot_strongest_wind(iwind,fwind,2000,ncdir,outdir,clvs=windlvs)
            p.plot_strongest_wind(iwind,fwind,2000,ncdir=ncdir,nct=nct,outdir=outdir,clvs=windlvs,dom=1,
                    cmap='jet',cb=True)

if accum_rain:
    for en in ensnames:
        for ex in experiments:
            for t in times: