In [1]:
import os
import pyart
import shutil
import warnings
from datetime import datetime as dt
import matplotlib as mpl
from matplotlib import pyplot as plt
import numpy as np
import numpy.ma as ma
import yaml
from glob import glob
import argparse
from samuraiAnalysis import samImport_masked,samImport,samPlt
samImport_masked = samImport_masked.samImport_masked
plotContour = samPlt.plotContour
plotVec = samPlt.plotVec
plotXS = samPlt.plotXS
warnings.filterwarnings('ignore', 'invalid value encountered in less')
warnings.filterwarnings('ignore', 'tight_layout : falling back to Agg renderer')
%matplotlib inline
In [2]:
m1 = '/Users/danstechman/GoogleDrive/PECAN-Data/samurai/20150620/s7_fsdabr/output/samurai_XYZ_analysis.nc'
s1 = '/Users/danstechman/GoogleDrive/PECAN-Data/samurai/20150620/s7_only/output/samurai_XYZ_analysis.nc'
s2 = '/Users/danstechman/GoogleDrive/PECAN-Data/samurai/20150620/s7sub_abr_only/output/samurai_XYZ_analysis.nc'
s3 = '/Users/danstechman/GoogleDrive/PECAN-Data/samurai/20150620/s7sub_fsd_only/output/samurai_XYZ_analysis.nc'
samDict = samImport_masked(m1,s1,s2,s3)
In [3]:
x1d = samDict['x']
y1d = samDict['y']
alt1d = samDict['alt']
dbz3d = samDict['dbz']
u3d = samDict['u']
v3d = samDict['v']
w3d = samDict['w']
In [4]:
xsStrt = (-75,10)
xsEnd = (99.9,10)
x1,y1 = xsStrt
x2,y2 = xsEnd
x_XY = np.array([x1,x2])
y_XY = np.array([y1,y2])
In [5]:
xStrd = 7
yStrd = 7
plotContour(dbz3d,4,x1d,y1d,'xy','dbz')
plt.plot(x_XY,y_XY)
plt.barbs(x1d[1::xStrd],y1d[1::yStrd],u3d[4,1::xStrd,1::yStrd],v3d[4,1::xStrd,1::yStrd],length=6.5)
Out[5]:
In [7]:
plotXS(dbz3d,x1d,y1d,alt1d,u3d,v3d,w3d,xsStrt,xsEnd,'dbz',figsize=(20,5),barbPlt=True)
Out[7]: