Пример #1
0
                                   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')

    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))
Пример #2
0
    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'
                                )

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: