Пример #1
0
            if radarplot:
                verifdir = '/chinook2/jrlawson/bowecho/{0}/VERIF'.format(case)
                p.plot_radar(t,
                             verifdir,
                             outdir=outdir,
                             ncdir=ncroot,
                             nct=nct,
                             dom=dom)

    if axesofdilatation:
        outdir = os.path.join(outroot, 'hires', 'd0{0}'.format(dom))
        for t in times:
            p.plot_axes_of_dilatation(utc=t,
                                      ncdir=ncroot,
                                      nct=nct,
                                      outdir=outdir,
                                      cb=True,
                                      dom=dom)

    if radarcomp:
        en = ensnames[0]
        ex = experiments[0]
        out_sd, wrf_sd = get_folders(en, ex)
        outdir, datadir = get_verif_dirs()
        # p.plot_radar(compt,datadir,outdir,ncdir=wrf_sd,composite=True)
        p.plot_radar(
            compt,
            datadir,
            outdir,
            composite=True,
            # Nlim=40.1,Elim=-94.9,Slim=34.3,Wlim=-100.8)
Пример #2
0
                # p.plot2D('lyapunov',t,700,ncdir=ncdir,nct=nct,outdir=outdir,cb=True,dom=dom,clvs=N.arange(-7.5,8.0,0.5)*10**-3,cmap='bwr')
                # print(ncdir)
                # p.plot2D('WSPD10MAX',t,ncdir=ncdir,nct=nct,outdir=outdir,cb=True,dom=dom,clvs=N.arange(10,31,1))
                p.plot2D('cref',t,ncdir=ncdir,nct=nct,outdir=outdir,cb=True,dom=dom,**lims)
                # p.plot2D('REFL_comp',t,ncdir=ncdir,nct=nct,outdir=outdir,cb=True,dom=dom)
                # p.plot2D('shear',t,ncdir=ncdir,nct=nct,outdir=outdir,cb=True,dom=dom)
                # p.plot2D('wind10',t,ncdir=ncdir,outdir=outdir,locations=locs,cb=True,clvs=N.arange(5,32,2))

            if radarplot:
                verifdir = '/chinook2/jrlawson/bowecho/{0}/VERIF'.format(case)
                p.plot_radar(t,verifdir,outdir=outdir,ncdir=ncroot,nct=nct,dom=dom)

    if axesofdilatation:
        outdir = os.path.join(outroot,'hires','d0{0}'.format(dom))
        for t in times:
            p.plot_axes_of_dilatation(utc=t,ncdir=ncroot,nct=nct,outdir=outdir,cb=True,dom=dom)

    if radarcomp:
        en = ensnames[0]
        ex = experiments[0]
        out_sd, wrf_sd = get_folders(en,ex)
        outdir, datadir = get_verif_dirs()
        # p.plot_radar(compt,datadir,outdir,ncdir=wrf_sd,composite=True)
        p.plot_radar(compt,datadir,outdir,composite=True,
                        # Nlim=40.1,Elim=-94.9,Slim=34.3,Wlim=-100.8)
                        Nlim=42.7,Elim=-94.9,Slim=37.0,Wlim=-101.8)

    if streamlines:
        if ens:
            outdir = os.path.join(outroot,'d0{0}'.format(dom),ens)
            ncdir = os.path.join(ncroot,ens)