Python WRFEnviron.plot_streamlines 예제들

예제 #1

파일: bowecho_plot.py 프로젝트: nishadhka/WEM

        for ex in experiments:
            out_sd, wrf_sd = get_folders(en,ex)
            
            # p.plot_strongest_wind(itime,ftime,2000,wrf_sd=wrf_sd,out_sd=out_sd)
            p.plot2D('cref',times,levels,wrf_sd=wrf_sd,out_sd=out_sd)

if streamlines:
    for en in ensnames:
        # Reload settings
        p.C = Settings()
        # Change paths to new location
        p.C.output_root = os.path.join(config.output_root,case,IC,en,experiment)
        p.C.wrfout_root = os.path.join(config.wrfout_root,case,IC,en,experiment)
        #p.plot_2D(variables)
        print p.C.wrfout_root
        p.plot_streamlines(2000,sl_times)

if rucplot:
    # RUC file is one-per-time so .nc file is specified beforehand
    en = ensnames[0]
    RC = Settings()
    RC.output_root = os.path.join(config.output_root,case,IC,en,experiment)
    RC.path_to_RUC = os.path.join(config.RUC_root,case,IC,en,experiment)
    WRF_dir = os.path.join(config.wrfout_root,case,'NAM',en,'ICBC')
    
    variables = ['streamlines',]
    level = 2000
    
    for t in sl_times:
        RUC = RUCPlot(RC,t,wrfdir=WRF_dir)
        #limits = RUC.colocate_WRF_map(WRF_dir)

예제 #2

            # Change paths to new location
            p.C.output_root = os.path.join(config.output_root,case,IC,en,MP,ex)
            p.C.wrfout_root = os.path.join(config.wrfout_root,case,IC,en,MP,ex)
            p.C.pickledir = os.path.join(config.wrfout_root,case,IC,en,MP,ex)
    
            p.plot_skewT(skewT_time,skewT_latlon,save_output=1)
    
    #p.plot_skewT(skewT_time,skewT_latlon,composite=1)

if plot2D:
    # p.plot_strongest_wind(itime,ftime,2000,wrf_sd=wrf_sd,out_sd=out_sd)
    p.plot2D('olr',times,levels)

if streamlines:
    for t in times:
        p.plot_streamlines(2000,t)

if rucplot:
    # RUC file is one-per-time so .nc file is specified beforehand
    en = ensnames[0]
    RC = Settings()
    RC.output_root = os.path.join(config.output_root,case,IC,en,experiment)
    RC.path_to_RUC = os.path.join(config.RUC_root,case,IC,en,experiment)
    WRF_dir = os.path.join(config.wrfout_root,case,'NAM',en,'ICBC')
    
    variables = ['streamlines',]
    level = 2000
    
    for t in sl_times:
        RUC = RUCPlot(RC,t,wrfdir=WRF_dir)
        #limits = RUC.colocate_WRF_map(WRF_dir)

예제 #3

파일: bowecho_plot.py 프로젝트: wf835334/WEM

if plot2D:
    for en in ensnames:
        for ex in experiments:
            for t in times:
                out_sd, wrf_sd = get_folders(en, ex)
                # p.plot_strongest_wind(itime,ftime,2000,wrf_sd=wrf_sd,out_sd=out_sd)
                # p.plot2D('Z',t,500,wrf_sd=wrf_sd,out_sd=out_sd,plottype='contour',smooth=10)
                p.plot2D('cref', t, 2000, wrf_sd=wrf_sd, out_sd=out_sd)

if streamlines:
    for en in ensnames:
        for ex in experiments:
            for t in times:
                out_sd, wrf_sd = get_folders(en, ex)
                p.plot_streamlines(t, 2000, out_sd=out_sd, wrf_sd=wrf_sd)

if rucplot:
    # RUC file is one-per-time so .nc file is specified beforehand
    en = ensnames[0]
    RC = Settings()
    RC.output_root = os.path.join(config.output_root, case, IC, en, experiment)
    RC.path_to_RUC = os.path.join(config.RUC_root, case, IC, en, experiment)
    WRF_dir = os.path.join(config.wrfout_root, case, 'NAM', en, 'ICBC')

    variables = [
        'streamlines',
    ]
    level = 2000

    for t in sl_times:

예제 #4

파일: barents_plot.py 프로젝트: NicWayand/WEM

            # Change paths to new location
            p.C.output_root = os.path.join(config.output_root, case, IC, en, MP, ex)
            p.C.wrfout_root = os.path.join(config.wrfout_root, case, IC, en, MP, ex)
            p.C.pickledir = os.path.join(config.wrfout_root, case, IC, en, MP, ex)

            p.plot_skewT(skewT_time, skewT_latlon, save_output=1)

    # p.plot_skewT(skewT_time,skewT_latlon,composite=1)

if plot2D:
    # p.plot_strongest_wind(itime,ftime,2000,wrf_sd=wrf_sd,out_sd=out_sd)
    p.plot2D("olr", times, levels)

if streamlines:
    for t in times:
        p.plot_streamlines(2000, t)

if rucplot:
    # RUC file is one-per-time so .nc file is specified beforehand
    en = ensnames[0]
    RC = Settings()
    RC.output_root = os.path.join(config.output_root, case, IC, en, experiment)
    RC.path_to_RUC = os.path.join(config.RUC_root, case, IC, en, experiment)
    WRF_dir = os.path.join(config.wrfout_root, case, "NAM", en, "ICBC")

    variables = ["streamlines"]
    level = 2000

    for t in sl_times:
        RUC = RUCPlot(RC, t, wrfdir=WRF_dir)
        # limits = RUC.colocate_WRF_map(WRF_dir)

예제 #5

파일: nonsquare.py 프로젝트: johnrobertlawson/evac

    for t in times:
        for lv in levels:
            for v in vrbls:
                p.plot2D(v,
                         t,
                         lv,
                         ncdir,
                         outdir,
                         nct=nct,
                         clvs=N.arange(5, 150, 5))

# Streamline 2D plots
if streamlines:
    for t in times:
        for lv in levels:
            p.plot_streamlines(t, lv, ncdir, outdir, smooth=5)

# Skew Ts
# skewT_time = (2013,8,16,3,0,0)
# skewT_latlon = (35.2435,-97.4708)
# p.plot_skewT(skewT_time,skewT_latlon)

# DKE
# p.compute_diff_energy('sum_z','total',path_to_wrfouts,times,
#   d_save=picklefolder, d_return=0,d_fname=pfname)

# latA, lonA, latB, lonB = (40.1,-89.0,38.1,-86.0)
# xstime = (2011,4,20,3,0,0)
# clvs = N.arange(-25,27.5,2.5)
#clvs = N.arange(0.0001,0.001,0.0001)
# p.plot_xs('parawind',xstime,latA=latA,lonA=lonA,latB=latB,lonB=lonB,

예제 #6

import numpy as N

sys.path.append('/path/to/WEM/')

from WEM.postWRF import WRFEnviron

p = WRFEnviron()

# Time of initialisation
nct = (2006,5,10,0,0,0)
# Or name of wrfout file
ncf = = 'wrfout_d01...'

# Directories for output and netCDF files
outdir = '/absolute/path/to/figures/'
ncdir = '/absolute/path/to/wrfoutdata/'

itime = (2006,5,10,12,0,0)
ftime = (2006,5,11,12,0,0)
hourly = 3
times = p.generate_times(itime,ftime,hourly*60*60)
level = 2000


for t in times:
    p.plot2D('cref',t,level,outdir=outdir,ncdir=ncdir,ncf=ncf,
                nct=nct,legend=True)
    p.plot2D('Z',t,(300,850),outdir=outdir,ncdir=ncdir,ncf=ncf,
                nct=nct,legend=True)
    p.plot_streamlines()

예제 #7

파일: nonsquare.py 프로젝트: NicWayand/WEM

times = [(2014, 7, 7, 18, 0, 0)]
levels = (2000,)
vrbls = ("RAINNC",)

# Plot standard 2D plots
if plot2D:
    for t in times:
        for lv in levels:
            for v in vrbls:
                p.plot2D(v, t, lv, ncdir, outdir, nct=nct, clvs=N.arange(5, 150, 5))

# Streamline 2D plots
if streamlines:
    for t in times:
        for lv in levels:
            p.plot_streamlines(t, lv, ncdir, outdir, smooth=5)

# Skew Ts
# skewT_time = (2013,8,16,3,0,0)
# skewT_latlon = (35.2435,-97.4708)
# p.plot_skewT(skewT_time,skewT_latlon)

# DKE
# p.compute_diff_energy('sum_z','total',path_to_wrfouts,times,
#   d_save=picklefolder, d_return=0,d_fname=pfname)

# latA, lonA, latB, lonB = (40.1,-89.0,38.1,-86.0)
# xstime = (2011,4,20,3,0,0)
# clvs = N.arange(-25,27.5,2.5)
# clvs = N.arange(0.0001,0.001,0.0001)
# p.plot_xs('parawind',xstime,latA=latA,lonA=lonA,latB=latB,lonB=lonB,

예제 #8

            # p.plot_strongest_wind(itime,ftime,2000,wrf_sd=wrf_sd,out_sd=out_sd)
            p.plot2D('cref', times, levels, wrf_sd=wrf_sd, out_sd=out_sd)

if streamlines:
    for en in ensnames:
        # Reload settings
        p.C = Settings()
        # Change paths to new location
        p.C.output_root = os.path.join(config.output_root, case, IC, en,
                                       experiment)
        p.C.wrfout_root = os.path.join(config.wrfout_root, case, IC, en,
                                       experiment)
        #p.plot_2D(variables)
        print p.C.wrfout_root
        p.plot_streamlines(2000, sl_times)

if rucplot:
    # RUC file is one-per-time so .nc file is specified beforehand
    en = ensnames[0]
    RC = Settings()
    RC.output_root = os.path.join(config.output_root, case, IC, en, experiment)
    RC.path_to_RUC = os.path.join(config.RUC_root, case, IC, en, experiment)
    WRF_dir = os.path.join(config.wrfout_root, case, 'NAM', en, 'ICBC')

    variables = [
        'streamlines',
    ]
    level = 2000

    for t in sl_times:

예제 #9

파일: bowecho_plot.py 프로젝트: qingu/WEM

if plot2D:
    for en in ensnames:
        for ex in experiments:
            for t in times:
                out_sd, wrf_sd = get_folders(en,ex)
                # p.plot_strongest_wind(itime,ftime,2000,wrf_sd=wrf_sd,out_sd=out_sd)
                # p.plot2D('Z',t,500,wrf_sd=wrf_sd,out_sd=out_sd,plottype='contour',smooth=10)
                p.plot2D('cref',t,2000,wrf_sd=wrf_sd,out_sd=out_sd)

if streamlines:
    for en in ensnames:
        for ex in experiments:
            for t in times:
                out_sd, wrf_sd = get_folders(en,ex)
                p.plot_streamlines(t,2000,out_sd=out_sd,wrf_sd=wrf_sd)

if rucplot:
    # RUC file is one-per-time so .nc file is specified beforehand
    en = ensnames[0]
    RC = Settings()
    RC.output_root = os.path.join(config.output_root,case,IC,en,experiment)
    RC.path_to_RUC = os.path.join(config.RUC_root,case,IC,en,experiment)
    WRF_dir = os.path.join(config.wrfout_root,case,'NAM',en,'ICBC')
    
    variables = ['streamlines',]
    level = 2000
    
    for t in sl_times:
        RUC = RUCPlot(RC,t,wrfdir=WRF_dir)
        #limits = RUC.colocate_WRF_map(WRF_dir)