Ejemplo n.º 1
0
import calendar
import time
import matplotlib.pyplot as plt

sys.path.append('/uufs/chpc.utah.edu/common/home/u0737349/gitprojects/')

from DKE_settings import Settings
from WEM.postWRF import WRFEnviron
import WEM.utils as utils

# Time script
scriptstart = time.time()

# Initialise settings and environment
config = Settings()
p = WRFEnviron(config)

# User settings
init_time = utils.string_from_time('dir', (2009, 9, 10, 0, 0, 0),
                                   strlen='hour')
rootdir = '/uufs/chpc.utah.edu/common/home/horel-group2/lawson2/'
outdir = '/uufs/chpc.utah.edu/common/home/u0737349/public_html/paper2/'

# Create lists of all fig, ax, cb objects
figs = {}
axes = {}
cbs = {}

ensnames = ['c00'] + ['p{0:02d}'.format(n) for n in range(1, 11)]

#for rundate in ('25','27','29'):
Ejemplo n.º 2
0
import time
import matplotlib.pyplot as plt

sys.path.append('/uufs/chpc.utah.edu/common/home/u0737349/gitprojects/')

from DKE_settings import Settings
from WEM.postWRF import WRFEnviron
import WEM.utils as utils


# Time script
scriptstart = time.time()

# Initialise settings and environment
config = Settings()
p = WRFEnviron(config)

# User settings
init_time = utils.string_from_time('dir',(2009,9,10,0,0,0),strlen='hour')
rootdir = '/uufs/chpc.utah.edu/common/home/horel-group2/lawson2/'
outdir = '/uufs/chpc.utah.edu/common/home/u0737349/public_html/paper2/'

# Create lists of all fig, ax, cb objects
figs = {}
axes = {}
cbs = {}

ensnames = ['c00'] + ['p{0:02d}'.format(n) for n in range(1,11)]

#for rundate in ('25','27','29'):
plot_all = 1
Ejemplo n.º 3
0
import pdb
import sys
import matplotlib as M
M.use('gtkagg')
import matplotlib.pyplot as plt
import numpy as N

sys.path.append('/home/jrlawson/gitprojects/')

from WEM.postWRF import WRFEnviron
from settings import Settings
import WEM.utils as utils
#from WEM.postWRF.postWRF.rucplot import RUCPlot

config = Settings()
p = WRFEnviron(config)

skewT = 0
plot2D = 0
streamlines = 0
rucplot = 0
coldpoolstrength = 0
spaghetti = 0
std = 0
profiles = 0
frontogenesis = 1
upperlevel = 0

# enstype = 'STCH'
# enstype = 'ICBC'
enstype = 'MXMP'
Ejemplo n.º 4
0
import os
import pdb
import sys
import matplotlib as M
M.use('gtkagg')
import matplotlib.pyplot as plt

sys.path.append('/home/jrlawson/gitprojects/')

from WEM.postWRF import WRFEnviron
from barents_settings import Settings
import WEM.utils as utils
#from WEM.postWRF.postWRF.rucplot import RUCPlot

config = Settings()
p = WRFEnviron(config)

skewT = 0
plot2D = 0
streamlines = 1
rucplot = 0
coldpoolstrength = 0

itime = (2007,7,20,0,0,0)
ftime = (2007,7,22,0,0,0)
# times = [(2013,8,15,18,0,0),]
levels = 2000

times = utils.generate_times(itime,ftime,2*60*60)

#shear_times = utils.generate_times(itime,ftime,3*60*60)
Ejemplo n.º 5
0
import pdb
import sys
import matplotlib as M
M.use('gtkagg')
import matplotlib.pyplot as plt
import numpy as N

sys.path.append('/home/jrlawson/gitprojects/')

from WEM.postWRF import WRFEnviron
from settings import Settings
import WEM.utils as utils
#from WEM.postWRF.postWRF.rucplot import RUCPlot

config = Settings()
p = WRFEnviron(config)

skewT = 0
plot2D = 0
streamlines = 0
rucplot = 0
coldpoolstrength = 0
spaghetti = 0
std = 0
profiles = 0
frontogenesis = 1
upperlevel = 0

# enstype = 'STCH'
# enstype = 'ICBC'
enstype = 'MXMP'
Ejemplo n.º 6
0
import sys
import matplotlib as M

M.use("gtkagg")
import matplotlib.pyplot as plt

sys.path.append("/home/jrlawson/gitprojects/")

from WEM.postWRF import WRFEnviron
from barents_settings import Settings
import WEM.utils as utils

# from WEM.postWRF.postWRF.rucplot import RUCPlot

config = Settings()
p = WRFEnviron(config)

skewT = 0
plot2D = 0
streamlines = 1
rucplot = 0
coldpoolstrength = 0

itime = (2007, 7, 20, 0, 0, 0)
ftime = (2007, 7, 22, 0, 0, 0)
# times = [(2013,8,15,18,0,0),]
levels = 2000

times = utils.generate_times(itime, ftime, 2 * 60 * 60)

# shear_times = utils.generate_times(itime,ftime,3*60*60)
Ejemplo n.º 7
0
Archivo: DKE.py Proyecto: qingu/WEM
IC = 'NAM'
experiment = 'MXMP'
ens = 'anl'
MP = 'ICBC'
# Time script
scriptstart = time.time()
stoch_names = ["s{0:02d}".format(n) for n in range(1,11)]
ens_names = ['c00',] + ["p{0:02d}".format(n) for n in range(1,11)]

MPs = ['ICBC','WSM6_Grau','WSM6_Hail','Kessler','Ferrier','WSM5',
		'WDM5','Lin','WDM6_Grau','WDM6_Hail',
		'Morrison_Grau','Morrison_Hail']

# Initialise settings and environment
config = Settings()
p = WRFEnviron(config)

if case[:4] == '2006':
    itime = (2006,5,26,0,0,0)
    ftime = (2006,5,27,12,0,0)
elif case[:4] == '2009':
    itime = (2009,9,10,0,0,0)
    ftime = (2009,9,11,18,0,0)
elif case[:4] == '2011':
    itime = (2011,4,19,0,0,0)
    ftime = (2011,4,20,15,0,0)
else:
    itime = (2013,8,15,0,0,0)
    ftime = (2013,8,16,15,0,0)

times = utils.generate_times(itime,ftime,2*3600)
Ejemplo n.º 8
0
Archivo: DKE.py Proyecto: wf835334/WEM
MP = 'ICBC'
# Time script
scriptstart = time.time()
stoch_names = ["s{0:02d}".format(n) for n in range(1, 11)]
ens_names = [
    'c00',
] + ["p{0:02d}".format(n) for n in range(1, 11)]

MPs = [
    'ICBC', 'WSM6_Grau', 'WSM6_Hail', 'Kessler', 'Ferrier', 'WSM5', 'WDM5',
    'Lin', 'WDM6_Grau', 'WDM6_Hail', 'Morrison_Grau', 'Morrison_Hail'
]

# Initialise settings and environment
config = Settings()
p = WRFEnviron(config)

if case[:4] == '2006':
    itime = (2006, 5, 26, 0, 0, 0)
    ftime = (2006, 5, 27, 12, 0, 0)
elif case[:4] == '2009':
    itime = (2009, 9, 10, 0, 0, 0)
    ftime = (2009, 9, 11, 18, 0, 0)
elif case[:4] == '2011':
    itime = (2011, 4, 19, 0, 0, 0)
    ftime = (2011, 4, 20, 15, 0, 0)
else:
    itime = (2013, 8, 15, 0, 0, 0)
    ftime = (2013, 8, 16, 15, 0, 0)

times = utils.generate_times(itime, ftime, 2 * 3600)
Ejemplo n.º 9
0
import sys
import numpy as N
import pdb
import matplotlib as M
import os
M.use('gtkagg')
import matplotlib.pyplot as plt
from mpl_toolkits.axes_grid1 import make_axes_locatable

sys.path.append('/home/jrlawson/gitprojects/')

from WEM.postWRF import WRFEnviron

p = WRFEnviron()

outdir = ('/home/jrlawson/public_html/')
ncdir = ('/chinook2/jrlawson/')
nct = (2014, 7, 7, 12, 0, 0)

plot2D = 1
streamlines = 0

# Create list of times
times = [
    (2014, 7, 7, 18, 0, 0),
]
levels = (2000, )
vrbls = ('RAINNC', )

# Plot standard 2D plots
if plot2D:
Ejemplo n.º 10
0
import os
import pdb
import sys
import matplotlib as M
M.use('gtkagg')
import matplotlib.pyplot as plt
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:
Ejemplo n.º 11
0
import sys
import numpy as N
import pdb
import matplotlib as M
import os

M.use("gtkagg")
import matplotlib.pyplot as plt
from mpl_toolkits.axes_grid1 import make_axes_locatable

sys.path.append("/home/jrlawson/gitprojects/")

from WEM.postWRF import WRFEnviron

p = WRFEnviron()

outdir = "/home/jrlawson/public_html/"
ncdir = "/chinook2/jrlawson/"
nct = (2014, 7, 7, 12, 0, 0)

plot2D = 1
streamlines = 0

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

# Plot standard 2D plots
if plot2D:
    for t in times:
Ejemplo n.º 12
0
import pdb
import sys
import matplotlib as M
M.use('gtkagg')
import matplotlib.pyplot as plt
import numpy as N

sys.path.append('/home/jrlawson/gitprojects/')

from WEM.postWRF import WRFEnviron
from settings import Settings
import WEM.utils as utils
#from WEM.postWRF.postWRF.rucplot import RUCPlot

config = Settings()
p = WRFEnviron(config)

skewT = 0
plot2D = 0
streamlines = 0
rucplot = 0
coldpoolstrength = 0
spaghetti = 0
std = 0
profiles = 0
frontogenesis = 1
upperlevel = 0

# enstype = 'STCH'
# enstype = 'ICBC'
enstype = 'MXMP'
Ejemplo n.º 13
0
import pdb
import sys
import matplotlib as M
M.use('gtkagg')
import matplotlib.pyplot as plt
import numpy as N

sys.path.append('/home/jrlawson/gitprojects/')

from WEM.postWRF import WRFEnviron
from settings import Settings
import WEM.utils as utils
#from WEM.postWRF.postWRF.rucplot import RUCPlot

config = Settings()
p = WRFEnviron(config)

skewT = 0
plot2D = 0
streamlines = 0
rucplot = 0
coldpoolstrength = 0
spaghetti = 0
std = 0
profiles = 0
frontogenesis = 1
upperlevel = 0

# enstype = 'STCH'
# enstype = 'ICBC'
enstype = 'MXMP'