In [1]:

    

from matplotlib import pyplot as plt
import numpy as np
import pyproj
from datetime import datetime as dt
from matplotlib.ticker import MultipleLocator

from samuraiAnalysis import samImport,samPlt,getFLpathData

getFLpathData = getFLpathData.getFLpathData
samImport = samImport.samImport

%matplotlib inline


    

    




/Users/danstechman/anaconda3/envs/py35/lib/python3.5/site-packages/cmocean/tools.py:76: MatplotlibDeprecationWarning: The is_string_like function was deprecated in version 2.1.
  if not mpl.cbook.is_string_like(rgbin[0]):

In [2]:

    

flFile = '/Users/danstechman/GoogleDrive/PECAN-Data/FlightLevelData/Processed/20150620_FltLvl_Processed.nc'


    

In [3]:

    

startT2 = '20150620-050350'
endT2 = '20150620-051303'

lat2_0 = 44.530
lon2_0 = -102.3462

# xsAngleD2 = -12.5653
xsAngleD2 = 0

xsAngle2 = np.deg2rad(xsAngleD2)

savePath2 = '/Users/danstechman/GoogleDrive/PECAN-Data/samurai/20150620/s2_udx/figs/xs_fltTrx_ref/'

samFileMaster2 = '/Users/danstechman/GoogleDrive/PECAN-Data/samurai/20150620/s2_udx/output/samurai_XYZ_analysis.nc'
samDataMaster2 = samImport(samFileMaster2)

lat1d2 = samDataMaster2['lat']
lon1d2 = samDataMaster2['lon']
x1d2 = samDataMaster2['x']
y1d2 = samDataMaster2['y']
alt1d2 = samDataMaster2['alt']


flS2 = dt.strptime(startT2,'%Y%m%d-%H%M%S')
flE2 = dt.strptime(endT2,'%Y%m%d-%H%M%S')

flData2 = getFLpathData(flFile,pathStrt=flS2,pathEnd=flE2)

flLat2 = flData2['lat']
flLon2 = flData2['lon']
flAlt2_msl = flData2['alt']

flAlt2_elevAdj = 790 # Elevation at spiral location
flAlt2 = flAlt2_msl - flAlt2_elevAdj

p2 = pyproj.Proj(proj='laea',ellps='WGS84',lat_0=lat2_0,lon_0=lon2_0)

fl_x2,fl_y2 = p2(flLon2,flLat2)

flX_rot2 = ((fl_x2 * np.cos(xsAngle2)) + (fl_y2 * np.sin(xsAngle2)))/1000
flY_rot2 = ((fl_y2 * np.cos(xsAngle2)) - (fl_x2 * np.sin(xsAngle2)))/1000

fig = plt.figure(figsize=(9.2,6.4))
ax = fig.add_subplot(111)
plt.plot(flX_rot2,flAlt2/1000,color='k',linewidth=3)
plt.plot(flX_rot2,flAlt2/1000,color='w',linewidth=1)
ax.set_xlim((-50,50))
ax.set_ylim((0,16))
ax.xaxis.set_major_locator(MultipleLocator(5))
ax.xaxis.set_minor_locator(MultipleLocator(1))
ax.tick_params(labelsize=12)
# saveStr2 = '{}/s2_diagXS_fltTrk.pdf'.format(savePath2)
saveStr2 = '{}/s2_weXS_fltTrk.pdf'.format(savePath2)
fig.savefig(saveStr2, transparent=True,bbox_inches='tight')


    

In [4]:

    

startT5 = '20150620-060437'
endT5 = '20150620-062248'

lat5_0 = 44.2180
lon5_0 = -100.9570

# xsAngleD5 = -44.3883
xsAngleD5 = 0

xsAngle5 = np.deg2rad(xsAngleD5)

savePath5 = '/Users/danstechman/GoogleDrive/PECAN-Data/samurai/20150620/s5_udxabr/figs/xs_fltTrx_ref/'

samFileMaster5 = '/Users/danstechman/GoogleDrive/PECAN-Data/samurai/20150620/s5_udxabr/output/samurai_XYZ_analysis.nc'
samDataMaster5 = samImport(samFileMaster5)

lat1d5 = samDataMaster5['lat']
lon1d5 = samDataMaster5['lon']
x1d5 = samDataMaster5['x']
y1d5 = samDataMaster5['y']
alt1d5 = samDataMaster5['alt']


flS5 = dt.strptime(startT5,'%Y%m%d-%H%M%S')
flE5 = dt.strptime(endT5,'%Y%m%d-%H%M%S')

flData5 = getFLpathData(flFile,pathStrt=flS5,pathEnd=flE5)

flLat5 = flData5['lat']
flLon5 = flData5['lon']
flAlt5_msl = flData5['alt']

flAlt5_elevAdj = 630 # Elevation at spiral location
flAlt5 = flAlt5_msl - flAlt5_elevAdj


p5 = pyproj.Proj(proj='laea',ellps='WGS84',lat_0=lat5_0,lon_0=lon5_0)

fl_x5,fl_y5 = p5(flLon5,flLat5)

flX_rot5 = ((fl_x5 * np.cos(xsAngle5)) + (fl_y5 * np.sin(xsAngle5)))/1000
flY_rot5 = ((fl_y5 * np.cos(xsAngle5)) - (fl_x5 * np.sin(xsAngle5)))/1000

fig = plt.figure(figsize=(9.2,6.4))
ax = fig.add_subplot(111)
plt.plot(flX_rot5,flAlt5/1000,color='k',linewidth=3)
plt.plot(flX_rot5,flAlt5/1000,color='w',linewidth=1)
ax.set_xlim((-50,50))
ax.set_ylim((0,16))
ax.xaxis.set_major_locator(MultipleLocator(5))
ax.xaxis.set_minor_locator(MultipleLocator(1))
ax.tick_params(labelsize=12)
# saveStr5 = '{}/s5_diagXS_fltTrk.pdf'.format(savePath5)
saveStr5 = '{}/s5_weXS_fltTrk.pdf'.format(savePath5)
fig.savefig(saveStr5, transparent=True,bbox_inches='tight')


    

In [5]:

    

startT7 = '20150620-074547'
endT7 = '20150620-075822'

lat7_0 = 43.8385
lon7_0 = -98.2041

# xsAngleD7 = -28.5112
xsAngleD7 = 0

xsAngle7 = np.deg2rad(xsAngleD7)

savePath7 = '/Users/danstechman/GoogleDrive/PECAN-Data/samurai/20150620/s7_fsdabr/figs/xs_fltTrx_ref/'

samFileMaster7 = '/Users/danstechman/GoogleDrive/PECAN-Data/samurai/20150620/s7_fsdabr/output/samurai_XYZ_analysis.nc'
samDataMaster7 = samImport(samFileMaster7)

lat1d7 = samDataMaster7['lat']
lon1d7 = samDataMaster7['lon']
x1d7 = samDataMaster7['x']
y1d7 = samDataMaster7['y']
alt1d7 = samDataMaster7['alt']


flS7 = dt.strptime(startT7,'%Y%m%d-%H%M%S')
flE7 = dt.strptime(endT7,'%Y%m%d-%H%M%S')

flData7 = getFLpathData(flFile,pathStrt=flS7,pathEnd=flE7)

flLat7 = flData7['lat']
flLon7 = flData7['lon']
flAlt7 = flData7['alt']
flAlt7_msl = flData7['alt']

flAlt7_elevAdj = 400 # Elevation at spiral location
flAlt7 = flAlt7_msl - flAlt7_elevAdj


p7 = pyproj.Proj(proj='laea',ellps='WGS84',lat_0=lat7_0,lon_0=lon7_0)

fl_x7,fl_y7 = p7(flLon7,flLat7)

flX_rot7 = ((fl_x7 * np.cos(xsAngle7)) + (fl_y7 * np.sin(xsAngle7)))/1000
flY_rot7 = ((fl_y7 * np.cos(xsAngle7)) - (fl_x7 * np.sin(xsAngle7)))/1000

fig = plt.figure(figsize=(9.2,6.4))
ax = fig.add_subplot(111)
plt.plot(flX_rot7,flAlt7/1000,color='k',linewidth=3)
plt.plot(flX_rot7,flAlt7/1000,color='w',linewidth=1)
ax.set_xlim((-50,50))
ax.set_ylim((0,16))
ax.xaxis.set_major_locator(MultipleLocator(5))
ax.xaxis.set_minor_locator(MultipleLocator(1))
ax.tick_params(labelsize=12)
# saveStr7 = '{}/s7_diagXS_fltTrk.pdf'.format(savePath7)
saveStr7 = '{}/s7_weXS_fltTrk.pdf'.format(savePath7)
fig.savefig(saveStr7, transparent=True,bbox_inches='tight')