Plot diagrams indicating when each volume of ground-based radars occured relative to the periods of TDR data collected during PDDs and spirals during PECAN (20150706).
In [ ]:
from datetime import datetime as dt
import matplotlib.pyplot as plt
import numpy as np
from matplotlib.pyplot import cm
import matplotlib.dates as mdates
import os
%matplotlib inline
dfmt = mdates.DateFormatter('%H:%M') # Date/time format for plot tick labels
In [ ]:
def plotGantt(strtTs,endTs,c,lbl,labs,tickloc,col,radNum):
for i in range(len(strtTs)):
plt.hlines(radNum,strtTs[i],endTs[i],color=c,linewidth=40)
labs.append(lbl)
tickloc.append(radNum)
col.append(c)
radNum -= 1
return labs,tickloc,col,radNum
In [ ]:
flight = '20150706'
savePath = '/Users/danstechman/GoogleDrive/PECAN-Data/radar-data/' + flight + '/RadarVolumeTimePlots/'
fType = 'png'
# Make the output figure directory if it doesn't exist
if not os.path.exists(savePath):
os.makedirs(savePath)
In [ ]:
TDR_pdds_s = ['20150706-015628','20150706-025109','20150706-031216','20150706-034610','20150706-040115',
'20150706-041740','20150706-045600','20150706-051818','20150706-061517','20150706-070004']
TDR_pdds_e = ['20150706-021900','20150706-031035','20150706-032103','20150706-035612','20150706-041620',
'20150706-042340','20150706-051457','20150706-054001','20150706-062230','20150706-073103']
TDR_pdds_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in TDR_pdds_s])
TDR_pdds_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in TDR_pdds_e])
TDR_sprls_s = ['20150706-031947','20150706-033053','20150706-042344','20150706-044005','20150706-054019',
'20150706-055624','20150706-062230','20150706-063759']
TDR_sprls_e = ['20150706-033052','20150706-034534','20150706-044004','20150706-045335','20150706-055549',
'20150706-061313','20150706-063758','20150706-065435']
TDR_sprls_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in TDR_sprls_s])
TDR_sprls_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in TDR_sprls_e])
In [ ]:
#
tNum = 1
totalRad = 8
radNum = totalRad
############# Start/End Times for Radar Volumes #############
D6_p01_s = ['20150706-015707','20150706-020119','20150706-020132','20150706-020239','20150706-020459',
'20150706-020959','20150706-021459']
D6_p01_e = ['20150706-020109','20150706-020126','20150706-020222','20150706-020426','20150706-020919',
'20150706-021420','20150706-021920']
D6_p01_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_p01_s])
D6_p01_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_p01_e])
D7_p01_s = ['20150706-015500','20150706-020000','20150706-020500','20150706-021000','20150706-021502']
D7_p01_e = ['20150706-015921','20150706-020420','20150706-020921','20150706-021421','20150706-021921']
D7_p01_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_p01_s])
D7_p01_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_p01_e])
D8_p01_s = ['20150706-015501','20150706-015959','20150706-020500','20150706-021001','20150706-021500']
D8_p01_e = ['20150706-015946','20150706-020453','20150706-020922','20150706-021454','20150706-021954']
D8_p01_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_p01_s])
D8_p01_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_p01_e])
SR1_p01_s = ['20150706-015508','20150706-020008','20150706-020508','20150706-021007','20150706-021508']
SR1_p01_e = ['20150706-015914','20150706-020413','20150706-020916','20150706-021411','20150706-021912']
SR1_p01_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_p01_s])
SR1_p01_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_p01_e])
SR2_p01_s = ['20150706-015504','20150706-020004','20150706-020504','20150706-021003','20150706-021503']
SR2_p01_e = ['20150706-015913','20150706-020412','20150706-020912','20150706-021411','20150706-021912']
SR2_p01_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_p01_s])
SR2_p01_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_p01_e])
# RX_p01_s = ['20150706-']
# RX_p01_e = ['20150706-']
# RX_p01_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_p01_s])
# RX_p01_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_p01_e])
# XP_p01_s = ['20150706-']
# XP_p01_e = ['20150706-']
# XP_p01_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_p01_s])
# XP_p01_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_p01_e])
FSD_p01_s = ['20150706-015506','20150706-015915','20150706-020339','20150706-020749','20150706-021159',
'20150706-021610']
FSD_p01_e = ['20150706-015806','20150706-020228','20150706-020638','20150706-021048','20150706-021459',
'20150706-021923']
FSD_p01_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_p01_s])
FSD_p01_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_p01_e])
ABR_p01_s = ['20150706-015733','20150706-020200']
ABR_p01_e = ['20150706-020045','20150706-020512']
ABR_p01_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_p01_s])
ABR_p01_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_p01_e])
############# Plotting #############
labs, tickloc, col = [], [], []
# create color iterator for multi-color lines in gantt chart
color=iter(cm.Dark2(np.linspace(0,1,radNum)))
fig, ax = plt.subplots(1,figsize=(15,10))
plt.title('6 July 2015 - PDD ' + str(tNum) + ' - Radar Volume Overlaps',fontsize=28)
c = next(color)
plt.hlines(radNum,TDR_pdds_s_DT[tNum-1],TDR_pdds_e_DT[tNum-1],color=c,linewidth=40)
labs.append('TDR')
tickloc.append(radNum)
col.append(c)
radNum -= 1
c = next(color)
labs,tickloc,col,radNum = plotGantt(D6_p01_s_DT,D6_p01_e_DT,c,'DOW6',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D7_p01_s_DT,D7_p01_e_DT,c,'DOW7',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D8_p01_s_DT,D8_p01_e_DT,c,'DOW8',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR1_p01_s_DT,SR1_p01_e_DT,c,'SR1',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR2_p01_s_DT,SR2_p01_e_DT,c,'SR2',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(FSD_p01_s_DT,FSD_p01_e_DT,c,'KFSD',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(ABR_p01_s_DT,ABR_p01_e_DT,c,'KABR',labs,tickloc,col,radNum)
plt.ylim(-1,totalRad+2)
plt.yticks(tickloc, labs)
plt.grid(which='both',axis='x')
plt.xlabel('Time (UTC)',fontsize=20)
plt.xlim((ax.get_xlim()[0]-0.0007),(ax.get_xlim()[1]+0.0007))
ytl = plt.gca().get_yticklabels()
for i in range(len(ytl)):
ytl[i].set_color(col[i])
plt.tick_params(axis='both', which='major', labelsize=20)
fig.autofmt_xdate()
ax.xaxis.set_major_locator(mdates.MinuteLocator(interval=2))
ax.xaxis.set_minor_locator(mdates.MinuteLocator(interval=1))
ax.xaxis.set_major_formatter(dfmt)
saveStr = '{}{}_PDD{:02d}_RadarVolTimes.{}'.format(savePath,flight,tNum,fType)
fig.savefig(saveStr,bbox_inches='tight')
In [ ]:
#
tNum = 2
totalRad = 10
radNum = totalRad
############# Start/End Times for Radar Volumes #############
D6_p02_s = ['20150706-025001','20150706-025501','20150706-030000','20150706-030500','20150706-031001']
D6_p02_e = ['20150706-025421','20150706-025921','20150706-030420','20150706-030921','20150706-031047']
D6_p02_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_p02_s])
D6_p02_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_p02_e])
D7_p02_s = ['20150706-025001','20150706-025501','20150706-030000','20150706-030501','20150706-031001']
D7_p02_e = ['20150706-025421','20150706-025921','20150706-030420','20150706-030921','20150706-031002']
D7_p02_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_p02_s])
D7_p02_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_p02_e])
D8_p02_s = ['20150706-025001','20150706-025501','20150706-030000','20150706-030500','20150706-031000']
D8_p02_e = ['20150706-025422','20150706-025955','20150706-030448','20150706-030955','20150706-031107']
D8_p02_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_p02_s])
D8_p02_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_p02_e])
SR1_p02_s = ['20150706-025013','20150706-025507','20150706-030007','20150706-030508','20150706-031008']
SR1_p02_e = ['20150706-025418','20150706-025911','20150706-030411','20150706-030912','20150706-031412']
SR1_p02_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_p02_s])
SR1_p02_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_p02_e])
SR2_p02_s = ['20150706-025020','20150706-025504','20150706-030004','20150706-030504','20150706-031003']
SR2_p02_e = ['20150706-025428','20150706-025913','20150706-030412','20150706-030912','20150706-031412']
SR2_p02_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_p02_s])
SR2_p02_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_p02_e])
RX_p02_s = ['20150706-024950','20150706-025452','20150706-030011','20150706-030448','20150706-031010']
RX_p02_e = ['20150706-025425','20150706-025859','20150706-030423','20150706-030902','20150706-031416']
RX_p02_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_p02_s])
RX_p02_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_p02_e])
XP_p02_s = ['20150706-025006','20150706-025506','20150706-030006','20150706-030506','20150706-031006']
XP_p02_e = ['20150706-025402','20150706-025902','20150706-030403','20150706-030902','20150706-031403']
XP_p02_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_p02_s])
XP_p02_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_p02_e])
FSD_p02_s = ['20150706-024844','20150706-025335','20150706-025826','20150706-030318','20150706-030809']
FSD_p02_e = ['20150706-025224','20150706-025715','20150706-030207','20150706-030658','20150706-031149']
FSD_p02_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_p02_s])
FSD_p02_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_p02_e])
ABR_p02_s = ['20150706-024859','20150706-025351','20150706-025846','20150706-030340','20150706-030834']
ABR_p02_e = ['20150706-025238','20150706-025731','20150706-030225','20150706-030719','20150706-031214']
ABR_p02_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_p02_s])
ABR_p02_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_p02_e])
############# Plotting #############
labs, tickloc, col = [], [], []
color=iter(cm.Dark2(np.linspace(0,1,radNum)))
fig, ax = plt.subplots(1,figsize=(15,10))
plt.title('6 July 2015 - PDD ' + str(tNum) + ' - Radar Volume Overlaps',fontsize=28)
c = next(color)
plt.hlines(radNum,TDR_pdds_s_DT[tNum-1],TDR_pdds_e_DT[tNum-1],color=c,linewidth=40)
labs.append('TDR')
tickloc.append(radNum)
col.append(c)
radNum -= 1
c = next(color)
labs,tickloc,col,radNum = plotGantt(D6_p02_s_DT,D6_p02_e_DT,c,'DOW6',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D7_p02_s_DT,D7_p02_e_DT,c,'DOW7',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D8_p02_s_DT,D8_p02_e_DT,c,'DOW8',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR1_p02_s_DT,SR1_p02_e_DT,c,'SR1',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR2_p02_s_DT,SR2_p02_e_DT,c,'SR2',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(RX_p02_s_DT,RX_p02_e_DT,c,'RAXPOL',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(XP_p02_s_DT,XP_p02_e_DT,c,'NOXP',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(FSD_p02_s_DT,FSD_p02_e_DT,c,'KFSD',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(ABR_p02_s_DT,ABR_p02_e_DT,c,'KABR',labs,tickloc,col,radNum)
plt.ylim(-0.25,totalRad+1.25)
plt.yticks(tickloc, labs)
plt.grid(which='both',axis='x')
plt.xlabel('Time (UTC)',fontsize=20)
plt.xlim((ax.get_xlim()[0]-0.0007),(ax.get_xlim()[1]+0.0007))
ytl = plt.gca().get_yticklabels()
for i in range(len(ytl)):
ytl[i].set_color(col[i])
plt.tick_params(axis='both', which='major', labelsize=20)
fig.autofmt_xdate()
ax.xaxis.set_major_locator(mdates.MinuteLocator(interval=2))
ax.xaxis.set_minor_locator(mdates.MinuteLocator(interval=1))
ax.xaxis.set_major_formatter(dfmt)
saveStr = '{}{}_PDD{:02d}_RadarVolTimes.{}'.format(savePath,flight,tNum,fType)
fig.savefig(saveStr,bbox_inches='tight')
In [ ]:
#
tNum = 3
totalRad = 10
radNum = totalRad
############# Start/End Times for Radar Volumes #############
D6_p03_s = ['20150706-031125','20150706-031529','20150706-031625','20150706-032034']
D6_p03_e = ['20150706-031520','20150706-031616','20150706-032024','20150706-032049']
D6_p03_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_p03_s])
D6_p03_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_p03_e])
D7_p03_s = ['20150706-031117','20150706-031523','20150706-031619','20150706-032032']
D7_p03_e = ['20150706-031513','20150706-031610','20150706-032018','20150706-032043']
D7_p03_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_p03_s])
D7_p03_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_p03_e])
D8_p03_s = ['20150706-031151','20150706-031601','20150706-031657','20150706-031959']
D8_p03_e = ['20150706-031552','20150706-031649','20150706-031946','20150706-032406']
D8_p03_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_p03_s])
D8_p03_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_p03_e])
SR1_p03_s = ['20150706-031008','20150706-031507','20150706-032007']
SR1_p03_e = ['20150706-031412','20150706-031912','20150706-032411']
SR1_p03_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_p03_s])
SR1_p03_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_p03_e])
SR2_p03_s = ['20150706-031003','20150706-031504','20150706-032004']
SR2_p03_e = ['20150706-031412','20150706-031912','20150706-032412']
SR2_p03_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_p03_s])
SR2_p03_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_p03_e])
RX_p03_s = ['20150706-031010','20150706-031442','20150706-031940']
RX_p03_e = ['20150706-031416','20150706-031915','20150706-032413']
RX_p03_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_p03_s])
RX_p03_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_p03_e])
XP_p03_s = ['20150706-031006','20150706-031506','20150706-032006']
XP_p03_e = ['20150706-031403','20150706-031902','20150706-032403']
XP_p03_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_p03_s])
XP_p03_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_p03_e])
FSD_p03_s = ['20150706-031300','20150706-031751']
FSD_p03_e = ['20150706-031640','20150706-032132']
FSD_p03_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_p03_s])
FSD_p03_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_p03_e])
ABR_p03_s = ['20150706-031328','20150706-031835']
ABR_p03_e = ['20150706-031720','20150706-032214']
ABR_p03_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_p03_s])
ABR_p03_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_p03_e])
############# Plotting #############
labs, tickloc, col = [], [], []
color=iter(cm.Dark2(np.linspace(0,1,radNum)))
fig, ax = plt.subplots(1,figsize=(15,10))
plt.title('6 July 2015 - PDD ' + str(tNum) + ' - Radar Volume Overlaps',fontsize=28)
c = next(color)
plt.hlines(radNum,TDR_pdds_s_DT[tNum-1],TDR_pdds_e_DT[tNum-1],color=c,linewidth=40)
labs.append('TDR')
tickloc.append(radNum)
col.append(c)
radNum -= 1
c = next(color)
labs,tickloc,col,radNum = plotGantt(D6_p03_s_DT,D6_p03_e_DT,c,'DOW6',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D7_p03_s_DT,D7_p03_e_DT,c,'DOW7',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D8_p03_s_DT,D8_p03_e_DT,c,'DOW8',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR1_p03_s_DT,SR1_p03_e_DT,c,'SR1',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR2_p03_s_DT,SR2_p03_e_DT,c,'SR2',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(RX_p03_s_DT,RX_p03_e_DT,c,'RAXPOL',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(XP_p03_s_DT,XP_p03_e_DT,c,'NOXP',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(FSD_p03_s_DT,FSD_p03_e_DT,c,'KFSD',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(ABR_p03_s_DT,ABR_p03_e_DT,c,'KABR',labs,tickloc,col,radNum)
plt.ylim(-0.25,totalRad+1.25)
plt.yticks(tickloc, labs)
plt.grid(which='both',axis='x')
plt.xlabel('Time (UTC)',fontsize=20)
plt.xlim((ax.get_xlim()[0]-0.0007),(ax.get_xlim()[1]+0.0007))
ytl = plt.gca().get_yticklabels()
for i in range(len(ytl)):
ytl[i].set_color(col[i])
plt.tick_params(axis='both', which='major', labelsize=20)
fig.autofmt_xdate()
ax.xaxis.set_major_locator(mdates.MinuteLocator(interval=2))
ax.xaxis.set_minor_locator(mdates.MinuteLocator(interval=1))
ax.xaxis.set_major_formatter(dfmt)
saveStr = '{}{}_PDD{:02d}_RadarVolTimes.{}'.format(savePath,flight,tNum,fType)
fig.savefig(saveStr,bbox_inches='tight')
In [ ]:
#
tNum = 4
totalRad = 10
radNum = totalRad
############# Start/End Times for Radar Volumes #############
D6_p04_s = ['20150706-034507','20150706-034916','20150706-035000','20150706-035412','20150706-035508']
D6_p04_e = ['20150706-034906','20150706-034931','20150706-035402','20150706-035459','20150706-035907']
D6_p04_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_p04_s])
D6_p04_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_p04_e])
D7_p04_s = ['20150706-034507','20150706-034921','20150706-034952','20150706-035001','20150706-035412',
'20150706-035508']
D7_p04_e = ['20150706-034907','20150706-034931','20150706-034953','20150706-035402','20150706-035459',
'20150706-035907']
D7_p04_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_p04_s])
D7_p04_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_p04_e])
D8_p04_s = ['20150706-034504','20150706-034914','20150706-035000','20150706-035409','20150706-035505']
D8_p04_e = ['20150706-034905','20150706-034949','20150706-035400','20150706-035457','20150706-035905']
D8_p04_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_p04_s])
D8_p04_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_p04_e])
SR1_p04_s = ['20150706-034507','20150706-035009','20150706-035509']
SR1_p04_e = ['20150706-034911','20150706-035415','20150706-035927']
SR1_p04_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_p04_s])
SR1_p04_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_p04_e])
SR2_p04_s = ['20150706-034505','20150706-035005','20150706-035522']
SR2_p04_e = ['20150706-034924','20150706-035423','20150706-035946']
SR2_p04_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_p04_s])
SR2_p04_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_p04_e])
RX_p04_s = ['20150706-034432','20150706-034930','20150706-035428']
RX_p04_e = ['20150706-034905','20150706-035403','20150706-035901']
RX_p04_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_p04_s])
RX_p04_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_p04_e])
XP_p04_s = ['20150706-034506','20150706-035021','20150706-035506']
XP_p04_e = ['20150706-034902','20150706-035436','20150706-035921']
XP_p04_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_p04_s])
XP_p04_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_p04_e])
FSD_p04_s = ['20150706-034657','20150706-035148']
FSD_p04_e = ['20150706-035038','20150706-035529']
FSD_p04_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_p04_s])
FSD_p04_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_p04_e])
ABR_p04_s = ['20150706-034346','20150706-034853','20150706-035400']
ABR_p04_e = ['20150706-034739','20150706-035246','20150706-035753']
ABR_p04_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_p04_s])
ABR_p04_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_p04_e])
############# Plotting #############
labs, tickloc, col = [], [], []
color=iter(cm.Dark2(np.linspace(0,1,radNum)))
fig, ax = plt.subplots(1,figsize=(15,10))
plt.title('6 July 2015 - PDD ' + str(tNum) + ' - Radar Volume Overlaps',fontsize=28)
c = next(color)
plt.hlines(radNum,TDR_pdds_s_DT[tNum-1],TDR_pdds_e_DT[tNum-1],color=c,linewidth=40)
labs.append('TDR')
tickloc.append(radNum)
col.append(c)
radNum -= 1
c = next(color)
labs,tickloc,col,radNum = plotGantt(D6_p04_s_DT,D6_p04_e_DT,c,'DOW6',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D7_p04_s_DT,D7_p04_e_DT,c,'DOW7',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D8_p04_s_DT,D8_p04_e_DT,c,'DOW8',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR1_p04_s_DT,SR1_p04_e_DT,c,'SR1',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR2_p04_s_DT,SR2_p04_e_DT,c,'SR2',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(RX_p04_s_DT,RX_p04_e_DT,c,'RAXPOL',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(XP_p04_s_DT,XP_p04_e_DT,c,'NOXP',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(FSD_p04_s_DT,FSD_p04_e_DT,c,'KFSD',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(ABR_p04_s_DT,ABR_p04_e_DT,c,'KABR',labs,tickloc,col,radNum)
plt.ylim(-0.25,totalRad+1.25)
plt.yticks(tickloc, labs)
plt.grid(which='both',axis='x')
plt.xlabel('Time (UTC)',fontsize=20)
plt.xlim((ax.get_xlim()[0]-0.0007),(ax.get_xlim()[1]+0.0007))
ytl = plt.gca().get_yticklabels()
for i in range(len(ytl)):
ytl[i].set_color(col[i])
plt.tick_params(axis='both', which='major', labelsize=20)
fig.autofmt_xdate()
ax.xaxis.set_major_locator(mdates.MinuteLocator(interval=2))
ax.xaxis.set_minor_locator(mdates.MinuteLocator(interval=1))
ax.xaxis.set_major_formatter(dfmt)
saveStr = '{}{}_PDD{:02d}_RadarVolTimes.{}'.format(savePath,flight,tNum,fType)
fig.savefig(saveStr,bbox_inches='tight')
In [ ]:
#
tNum = 5
totalRad = 10
radNum = totalRad
############# Start/End Times for Radar Volumes #############
D6_p05_s = ['20150706-040000','20150706-040411','20150706-040507','20150706-040916','20150706-041000',
'20150706-041412','20150706-041508']
D6_p05_e = ['20150706-040401','20150706-040458','20150706-040906','20150706-040931','20150706-041402',
'20150706-041459','20150706-041907']
D6_p05_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_p05_s])
D6_p05_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_p05_e])
D7_p05_s = ['20150706-040000','20150706-040411','20150706-040507','20150706-040916','20150706-041001',
'20150706-041412','20150706-041508']
D7_p05_e = ['20150706-040402','20150706-040458','20150706-040907','20150706-040932','20150706-041402',
'20150706-041459','20150706-041907']
D7_p05_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_p05_s])
D7_p05_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_p05_e])
D8_p05_s = ['20150706-035959','20150706-040409','20150706-040504','20150706-040914','20150706-041000',
'20150706-041409','20150706-041505']
D8_p05_e = ['20150706-040400','20150706-040456','20150706-040905','20150706-040949','20150706-041400',
'20150706-041457','20150706-041905']
D8_p05_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_p05_s])
D8_p05_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_p05_e])
SR1_p05_s = ['20150706-040016','20150706-040520','20150706-041016','20150706-041522']
SR1_p05_e = ['20150706-040434','20150706-040930','20150706-041430','20150706-041931']
SR1_p05_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_p05_s])
SR1_p05_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_p05_e])
SR2_p05_s = ['20150706-040011','20150706-040511','20150706-041011','20150706-041511']
SR2_p05_e = ['20150706-040435','20150706-040935','20150706-041435','20150706-041934']
SR2_p05_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_p05_s])
SR2_p05_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_p05_e])
RX_p05_s = ['20150706-035927','20150706-040425','20150706-040925','20150706-041422']
RX_p05_e = ['20150706-040400','20150706-040858','20150706-041358','20150706-041855']
RX_p05_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_p05_s])
RX_p05_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_p05_e])
XP_p05_s = ['20150706-040006','20150706-040511','20150706-041010','20150706-041511']
XP_p05_e = ['20150706-040437','20150706-040942','20150706-041440','20150706-041800']
XP_p05_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_p05_s])
XP_p05_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_p05_e])
FSD_p05_s = ['20150706-040130','20150706-040621','20150706-041112']
FSD_p05_e = ['20150706-040510','20150706-041001','20150706-041452']
FSD_p05_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_p05_s])
FSD_p05_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_p05_e])
ABR_p05_s = ['20150706-035908','20150706-040416','20150706-040925','20150706-041434']
ABR_p05_e = ['20150706-040301','20150706-040811','20150706-041320','20150706-041828']
ABR_p05_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_p05_s])
ABR_p05_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_p05_e])
############# Plotting #############
labs, tickloc, col = [], [], []
color=iter(cm.Dark2(np.linspace(0,1,radNum)))
fig, ax = plt.subplots(1,figsize=(15,10))
plt.title('6 July 2015 - PDD ' + str(tNum) + ' - Radar Volume Overlaps',fontsize=28)
c = next(color)
plt.hlines(radNum,TDR_pdds_s_DT[tNum-1],TDR_pdds_e_DT[tNum-1],color=c,linewidth=40)
labs.append('TDR')
tickloc.append(radNum)
col.append(c)
radNum -= 1
c = next(color)
labs,tickloc,col,radNum = plotGantt(D6_p05_s_DT,D6_p05_e_DT,c,'DOW6',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D7_p05_s_DT,D7_p05_e_DT,c,'DOW7',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D8_p05_s_DT,D8_p05_e_DT,c,'DOW8',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR1_p05_s_DT,SR1_p05_e_DT,c,'SR1',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR2_p05_s_DT,SR2_p05_e_DT,c,'SR2',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(RX_p05_s_DT,RX_p05_e_DT,c,'RAXPOL',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(XP_p05_s_DT,XP_p05_e_DT,c,'NOXP',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(FSD_p05_s_DT,FSD_p05_e_DT,c,'KFSD',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(ABR_p05_s_DT,ABR_p05_e_DT,c,'KABR',labs,tickloc,col,radNum)
plt.ylim(-0.25,totalRad+1.25)
plt.yticks(tickloc, labs)
plt.grid(which='both',axis='x')
plt.xlabel('Time (UTC)',fontsize=20)
plt.xlim((ax.get_xlim()[0]-0.0007),(ax.get_xlim()[1]+0.0007))
ytl = plt.gca().get_yticklabels()
for i in range(len(ytl)):
ytl[i].set_color(col[i])
plt.tick_params(axis='both', which='major', labelsize=20)
fig.autofmt_xdate()
ax.xaxis.set_major_locator(mdates.MinuteLocator(interval=2))
ax.xaxis.set_minor_locator(mdates.MinuteLocator(interval=1))
ax.xaxis.set_major_formatter(dfmt)
saveStr = '{}{}_PDD{:02d}_RadarVolTimes.{}'.format(savePath,flight,tNum,fType)
fig.savefig(saveStr,bbox_inches='tight')
In [ ]:
#
tNum = 6
totalRad = 10
radNum = totalRad
############# Start/End Times for Radar Volumes #############
D6_p06_s = ['20150706-041508','20150706-041917','20150706-042000']
D6_p06_e = ['20150706-041907','20150706-041932','20150706-042402']
D6_p06_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_p06_s])
D6_p06_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_p06_e])
D7_p06_s = ['20150706-041508','20150706-041922','20150706-042000']
D7_p06_e = ['20150706-041907','20150706-041932','20150706-042401']
D7_p06_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_p06_s])
D7_p06_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_p06_e])
D8_p06_s = ['20150706-041505','20150706-041914','20150706-042000']
D8_p06_e = ['20150706-041905','20150706-041938','20150706-042400']
D8_p06_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_p06_s])
D8_p06_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_p06_e])
SR1_p06_s = ['20150706-041522','20150706-042016']
SR1_p06_e = ['20150706-041931','20150706-042428']
SR1_p06_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_p06_s])
SR1_p06_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_p06_e])
SR2_p06_s = ['20150706-041511','20150706-042011']
SR2_p06_e = ['20150706-041934','20150706-042435']
SR2_p06_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_p06_s])
SR2_p06_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_p06_e])
RX_p06_s = ['20150706-041422','20150706-041921']
RX_p06_e = ['20150706-041855','20150706-042441']
RX_p06_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_p06_s])
RX_p06_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_p06_e])
XP_p06_s = ['20150706-041511','20150706-042011']
XP_p06_e = ['20150706-041942','20150706-042440']
XP_p06_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_p06_s])
XP_p06_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_p06_e])
FSD_p06_s = ['20150706-041603','20150706-042053']
FSD_p06_e = ['20150706-041943','20150706-042434']
FSD_p06_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_p06_s])
FSD_p06_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_p06_e])
ABR_p06_s = ['20150706-041434','20150706-041942']
ABR_p06_e = ['20150706-041828','20150706-042335']
ABR_p06_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_p06_s])
ABR_p06_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_p06_e])
############# Plotting #############
labs, tickloc, col = [], [], []
color=iter(cm.Dark2(np.linspace(0,1,radNum)))
fig, ax = plt.subplots(1,figsize=(15,10))
plt.title('6 July 2015 - PDD ' + str(tNum) + ' - Radar Volume Overlaps',fontsize=28)
c = next(color)
plt.hlines(radNum,TDR_pdds_s_DT[tNum-1],TDR_pdds_e_DT[tNum-1],color=c,linewidth=40)
labs.append('TDR')
tickloc.append(radNum)
col.append(c)
radNum -= 1
c = next(color)
labs,tickloc,col,radNum = plotGantt(D6_p06_s_DT,D6_p06_e_DT,c,'DOW6',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D7_p06_s_DT,D7_p06_e_DT,c,'DOW7',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D8_p06_s_DT,D8_p06_e_DT,c,'DOW8',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR1_p06_s_DT,SR1_p06_e_DT,c,'SR1',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR2_p06_s_DT,SR2_p06_e_DT,c,'SR2',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(RX_p06_s_DT,RX_p06_e_DT,c,'RAXPOL',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(XP_p06_s_DT,XP_p06_e_DT,c,'NOXP',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(FSD_p06_s_DT,FSD_p06_e_DT,c,'KFSD',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(ABR_p06_s_DT,ABR_p06_e_DT,c,'KABR',labs,tickloc,col,radNum)
plt.ylim(-0.25,totalRad+1.25)
plt.yticks(tickloc, labs)
plt.grid(which='both',axis='x')
plt.xlabel('Time (UTC)',fontsize=20)
plt.xlim((ax.get_xlim()[0]-0.0007),(ax.get_xlim()[1]+0.0007))
ytl = plt.gca().get_yticklabels()
for i in range(len(ytl)):
ytl[i].set_color(col[i])
plt.tick_params(axis='both', which='major', labelsize=20)
fig.autofmt_xdate()
ax.xaxis.set_major_locator(mdates.MinuteLocator(interval=2))
ax.xaxis.set_minor_locator(mdates.MinuteLocator(interval=1))
ax.xaxis.set_major_formatter(dfmt)
saveStr = '{}{}_PDD{:02d}_RadarVolTimes.{}'.format(savePath,flight,tNum,fType)
fig.savefig(saveStr,bbox_inches='tight')
In [ ]:
#
tNum = 7
totalRad = 10
radNum = totalRad
############# Start/End Times for Radar Volumes #############
D6_p07_s = ['20150706-045508','20150706-045917','20150706-050000','20150706-050411','20150706-050508',
'20150706-050916','20150706-051001','20150706-051412']
D6_p07_e = ['20150706-045907','20150706-045932','20150706-050401','20150706-050458','20150706-050906',
'20150706-050931','20150706-051402','20150706-051459']
D6_p07_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_p07_s])
D6_p07_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_p07_e])
D7_p07_s = ['20150706-051004','20150706-051414']
D7_p07_e = ['20150706-051404','20150706-051501']
D7_p07_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_p07_s])
D7_p07_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_p07_e])
D8_p07_s = ['20150706-045506','20150706-045914','20150706-050001','20150706-050409','20150706-050505',
'20150706-050914','20150706-051000','20150706-051409']
D8_p07_e = ['20150706-045905','20150706-045938','20150706-050400','20150706-050457','20150706-050905',
'20150706-050957','20150706-051400','20150706-051457']
D8_p07_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_p07_s])
D8_p07_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_p07_e])
SR1_p07_s = ['20150706-045516','20150706-050010','20150706-050515','20150706-051015']
SR1_p07_e = ['20150706-045940','20150706-050439','20150706-050943','20150706-051443']
SR1_p07_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_p07_s])
SR1_p07_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_p07_e])
SR2_p07_s = ['20150706-045526','20150706-050018','20150706-050511','20150706-051011']
SR2_p07_e = ['20150706-045935','20150706-050442','20150706-050935','20150706-051435']
SR2_p07_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_p07_s])
SR2_p07_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_p07_e])
RX_p07_s = ['20150706-045516','20150706-050037','20150706-050531','20150706-051033']
RX_p07_e = ['20150706-050015','20150706-050522','20150706-051025','20150706-051415']
RX_p07_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_p07_s])
RX_p07_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_p07_e])
XP_p07_s = ['20150706-045512','20150706-050011','20150706-050511','20150706-051010']
XP_p07_e = ['20150706-045943','20150706-050441','20150706-050942','20150706-051440']
XP_p07_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_p07_s])
XP_p07_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_p07_e])
FSD_p07_s = ['20150706-045450','20150706-045940','20150706-050431','20150706-050923','20150706-051414']
FSD_p07_e = ['20150706-045830','20150706-050321','20150706-050812','20150706-051303','20150706-051755']
FSD_p07_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_p07_s])
FSD_p07_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_p07_e])
ABR_p07_s = ['20150706-045535','20150706-050042','20150706-050550','20150706-051058']
ABR_p07_e = ['20150706-045928','20150706-050435','20150706-050943','20150706-051451']
ABR_p07_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_p07_s])
ABR_p07_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_p07_e])
############# Plotting #############
labs, tickloc, col = [], [], []
color=iter(cm.Dark2(np.linspace(0,1,radNum)))
fig, ax = plt.subplots(1,figsize=(15,10))
plt.title('6 July 2015 - PDD ' + str(tNum) + ' - Radar Volume Overlaps',fontsize=28)
c = next(color)
plt.hlines(radNum,TDR_pdds_s_DT[tNum-1],TDR_pdds_e_DT[tNum-1],color=c,linewidth=40)
labs.append('TDR')
tickloc.append(radNum)
col.append(c)
radNum -= 1
c = next(color)
labs,tickloc,col,radNum = plotGantt(D6_p07_s_DT,D6_p07_e_DT,c,'DOW6',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D7_p07_s_DT,D7_p07_e_DT,c,'DOW7',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D8_p07_s_DT,D8_p07_e_DT,c,'DOW8',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR1_p07_s_DT,SR1_p07_e_DT,c,'SR1',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR2_p07_s_DT,SR2_p07_e_DT,c,'SR2',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(RX_p07_s_DT,RX_p07_e_DT,c,'RAXPOL',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(XP_p07_s_DT,XP_p07_e_DT,c,'NOXP',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(FSD_p07_s_DT,FSD_p07_e_DT,c,'KFSD',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(ABR_p07_s_DT,ABR_p07_e_DT,c,'KABR',labs,tickloc,col,radNum)
plt.ylim(-0.25,totalRad+1.25)
plt.yticks(tickloc, labs)
plt.grid(which='both',axis='x')
plt.xlabel('Time (UTC)',fontsize=20)
plt.xlim((ax.get_xlim()[0]-0.0007),(ax.get_xlim()[1]+0.0007))
ytl = plt.gca().get_yticklabels()
for i in range(len(ytl)):
ytl[i].set_color(col[i])
plt.tick_params(axis='both', which='major', labelsize=20)
fig.autofmt_xdate()
ax.xaxis.set_major_locator(mdates.MinuteLocator(interval=2))
ax.xaxis.set_minor_locator(mdates.MinuteLocator(interval=1))
ax.xaxis.set_major_formatter(dfmt)
saveStr = '{}{}_PDD{:02d}_RadarVolTimes.{}'.format(savePath,flight,tNum,fType)
fig.savefig(saveStr,bbox_inches='tight')
In [ ]:
#
tNum = 8
totalRad = 10
radNum = totalRad
############# Start/End Times for Radar Volumes #############
D6_p08_s = ['20150706-051508','20150706-051917','20150706-052000','20150706-052411','20150706-052507',
'20150706-052916','20150706-053000','20150706-053412','20150706-053508','20150706-053917']
D6_p08_e = ['20150706-051907','20150706-051932','20150706-052401','20150706-052458','20150706-052906',
'20150706-052931','20150706-053402','20150706-053459','20150706-053907','20150706-053932']
D6_p08_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_p08_s])
D6_p08_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_p08_e])
D7_p08_s = ['20150706-051510','20150706-052211','20150706-052618','20150706-052714','20150706-053000',
'20150706-053412','20150706-053508']
D7_p08_e = ['20150706-051909','20150706-052608','20150706-052705','20150706-052944','20150706-053402',
'20150706-053459','20150706-053907']
D7_p08_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_p08_s])
D7_p08_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_p08_e])
D8_p08_s = ['20150706-051505','20150706-051914','20150706-051959','20150706-052409','20150706-052504',
'20150706-052914','20150706-053000','20150706-053409','20150706-053505','20150706-053914']
D8_p08_e = ['20150706-051905','20150706-051950','20150706-052400','20150706-052456','20150706-052905',
'20150706-052957','20150706-053400','20150706-053457','20150706-053905','20150706-053938']
D8_p08_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_p08_s])
D8_p08_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_p08_e])
SR1_p08_s = ['20150706-051518','20150706-052013','20150706-052517','20150706-053022','20150706-053526']
SR1_p08_e = ['20150706-051942','20150706-052437','20150706-052940','20150706-053450','20150706-053955']
SR1_p08_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_p08_s])
SR1_p08_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_p08_e])
SR2_p08_s = ['20150706-051510','20150706-052011','20150706-052511','20150706-053011','20150706-053511']
SR2_p08_e = ['20150706-051934','20150706-052434','20150706-052934','20150706-053434','20150706-053934']
SR2_p08_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_p08_s])
SR2_p08_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_p08_e])
RX_p08_s = ['20150706-051508','20150706-052000','20150706-052459','20150706-052959','20150706-053458']
RX_p08_e = ['20150706-051951','20150706-052451','20150706-052950','20150706-053450','20150706-053949']
RX_p08_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_p08_s])
RX_p08_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_p08_e])
XP_p08_s = ['20150706-051511','20150706-052010','20150706-052510','20150706-053011','20150706-053510']
XP_p08_e = ['20150706-051940','20150706-052440','20150706-052941','20150706-053440','20150706-053940']
XP_p08_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_p08_s])
XP_p08_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_p08_e])
FSD_p08_s = ['20150706-051905','20150706-052357','20150706-052847','20150706-053338','20150706-053829']
FSD_p08_e = ['20150706-052246','20150706-052737','20150706-053228','20150706-053718','20150706-054210']
FSD_p08_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_p08_s])
FSD_p08_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_p08_e])
ABR_p08_s = ['20150706-051605','20150706-052113','20150706-052622','20150706-053115','20150706-053609']
ABR_p08_e = ['20150706-051959','20150706-052507','20150706-053001','20150706-053455','20150706-053948']
ABR_p08_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_p08_s])
ABR_p08_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_p08_e])
############# Plotting #############
labs, tickloc, col = [], [], []
color=iter(cm.Dark2(np.linspace(0,1,radNum)))
fig, ax = plt.subplots(1,figsize=(15,10))
plt.title('6 July 2015 - PDD ' + str(tNum) + ' - Radar Volume Overlaps',fontsize=28)
c = next(color)
plt.hlines(radNum,TDR_pdds_s_DT[tNum-1],TDR_pdds_e_DT[tNum-1],color=c,linewidth=40)
labs.append('TDR')
tickloc.append(radNum)
col.append(c)
radNum -= 1
c = next(color)
labs,tickloc,col,radNum = plotGantt(D6_p08_s_DT,D6_p08_e_DT,c,'DOW6',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D7_p08_s_DT,D7_p08_e_DT,c,'DOW7',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D8_p08_s_DT,D8_p08_e_DT,c,'DOW8',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR1_p08_s_DT,SR1_p08_e_DT,c,'SR1',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR2_p08_s_DT,SR2_p08_e_DT,c,'SR2',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(RX_p08_s_DT,RX_p08_e_DT,c,'RAXPOL',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(XP_p08_s_DT,XP_p08_e_DT,c,'NOXP',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(FSD_p08_s_DT,FSD_p08_e_DT,c,'KFSD',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(ABR_p08_s_DT,ABR_p08_e_DT,c,'KABR',labs,tickloc,col,radNum)
plt.ylim(-0.25,totalRad+1.25)
plt.yticks(tickloc, labs)
plt.grid(which='both',axis='x')
plt.xlabel('Time (UTC)',fontsize=20)
plt.xlim((ax.get_xlim()[0]-0.0007),(ax.get_xlim()[1]+0.0007))
ytl = plt.gca().get_yticklabels()
for i in range(len(ytl)):
ytl[i].set_color(col[i])
plt.tick_params(axis='both', which='major', labelsize=20)
fig.autofmt_xdate()
ax.xaxis.set_major_locator(mdates.MinuteLocator(interval=2))
ax.xaxis.set_minor_locator(mdates.MinuteLocator(interval=1))
ax.xaxis.set_major_formatter(dfmt)
saveStr = '{}{}_PDD{:02d}_RadarVolTimes.{}'.format(savePath,flight,tNum,fType)
fig.savefig(saveStr,bbox_inches='tight')
In [ ]:
#
tNum = 9
totalRad = 10
radNum = totalRad
############# Start/End Times for Radar Volumes #############
D6_p09_s = ['20150706-061508','20150706-061917','20150706-062000']
D6_p09_e = ['20150706-061907','20150706-061932','20150706-062402']
D6_p09_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_p09_s])
D6_p09_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_p09_e])
D7_p09_s = ['20150706-061413','20150706-061750','20150706-062003']
D7_p09_e = ['20150706-061737','20150706-061857','20150706-062355']
D7_p09_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_p09_s])
D7_p09_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_p09_e])
D8_p09_s = ['20150706-061505','20150706-061914','20150706-062000']
D8_p09_e = ['20150706-061905','20150706-061946','20150706-062400']
D8_p09_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_p09_s])
D8_p09_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_p09_e])
SR1_p09_s = ['20150706-061512','20150706-062013']
SR1_p09_e = ['20150706-061932','20150706-062440']
SR1_p09_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_p09_s])
SR1_p09_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_p09_e])
SR2_p09_s = ['20150706-061510','20150706-062010']
SR2_p09_e = ['20150706-061934','20150706-062434']
SR2_p09_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_p09_s])
SR2_p09_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_p09_e])
RX_p09_s = ['20150706-061452','20150706-061952']
RX_p09_e = ['20150706-061943','20150706-062443']
RX_p09_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_p09_s])
RX_p09_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_p09_e])
XP_p09_s = ['20150706-061510','20150706-062011']
XP_p09_e = ['20150706-061940','20150706-062442']
XP_p09_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_p09_s])
XP_p09_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_p09_e])
FSD_p09_s = ['20150706-061225','20150706-061716','20150706-062206']
FSD_p09_e = ['20150706-061606','20150706-062056','20150706-062546']
FSD_p09_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_p09_s])
FSD_p09_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_p09_e])
ABR_p09_s = ['20150706-061403','20150706-061830']
ABR_p09_e = ['20150706-061715','20150706-062142']
ABR_p09_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_p09_s])
ABR_p09_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_p09_e])
############# Plotting #############
labs, tickloc, col = [], [], []
color=iter(cm.Dark2(np.linspace(0,1,radNum)))
fig, ax = plt.subplots(1,figsize=(15,10))
plt.title('6 July 2015 - PDD ' + str(tNum) + ' - Radar Volume Overlaps',fontsize=28)
c = next(color)
plt.hlines(radNum,TDR_pdds_s_DT[tNum-1],TDR_pdds_e_DT[tNum-1],color=c,linewidth=40)
labs.append('TDR')
tickloc.append(radNum)
col.append(c)
radNum -= 1
c = next(color)
labs,tickloc,col,radNum = plotGantt(D6_p09_s_DT,D6_p09_e_DT,c,'DOW6',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D7_p09_s_DT,D7_p09_e_DT,c,'DOW7',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D8_p09_s_DT,D8_p09_e_DT,c,'DOW8',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR1_p09_s_DT,SR1_p09_e_DT,c,'SR1',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR2_p09_s_DT,SR2_p09_e_DT,c,'SR2',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(RX_p09_s_DT,RX_p09_e_DT,c,'RAXPOL',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(XP_p09_s_DT,XP_p09_e_DT,c,'NOXP',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(FSD_p09_s_DT,FSD_p09_e_DT,c,'KFSD',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(ABR_p09_s_DT,ABR_p09_e_DT,c,'KABR',labs,tickloc,col,radNum)
plt.ylim(-0.25,totalRad+1.25)
plt.yticks(tickloc, labs)
plt.grid(which='both',axis='x')
plt.xlabel('Time (UTC)',fontsize=20)
plt.xlim((ax.get_xlim()[0]-0.0007),(ax.get_xlim()[1]+0.0007))
ytl = plt.gca().get_yticklabels()
for i in range(len(ytl)):
ytl[i].set_color(col[i])
plt.tick_params(axis='both', which='major', labelsize=20)
fig.autofmt_xdate()
ax.xaxis.set_major_locator(mdates.MinuteLocator(interval=2))
ax.xaxis.set_minor_locator(mdates.MinuteLocator(interval=1))
ax.xaxis.set_major_formatter(dfmt)
saveStr = '{}{}_PDD{:02d}_RadarVolTimes.{}'.format(savePath,flight,tNum,fType)
fig.savefig(saveStr,bbox_inches='tight')
In [ ]:
#
tNum = 10
totalRad = 9
radNum = totalRad
############# Start/End Times for Radar Volumes #############
D6_p10_s = ['20150706-070000','20150706-070411','20150706-070507','20150706-070916','20150706-071001',
'20150706-071412','20150706-071508','20150706-071917','20150706-072000','20150706-072411',
'20150706-072507','20150706-072916','20150706-073001']
D6_p10_e = ['20150706-070401','20150706-070458','20150706-070906','20150706-070931','20150706-071402',
'20150706-071459','20150706-071907','20150706-071932','20150706-072401','20150706-072458',
'20150706-072906','20150706-072931','20150706-073402']
D6_p10_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_p10_s])
D6_p10_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_p10_e])
D7_p10_s = ['20150706-070000','20150706-070411','20150706-070508','20150706-070916','20150706-071000',
'20150706-071412','20150706-071508','20150706-071917','20150706-072000','20150706-072411',
'20150706-072507','20150706-072916','20150706-073001']
D7_p10_e = ['20150706-070402','20150706-070458','20150706-070907','20150706-070931','20150706-071402',
'20150706-071459','20150706-071907','20150706-071932','20150706-072402','20150706-072458',
'20150706-072907','20150706-072932','20150706-073402']
D7_p10_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_p10_s])
D7_p10_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_p10_e])
D8_p10_s = ['20150706-070000','20150706-070408','20150706-070504','20150706-070913','20150706-071000',
'20150706-071409','20150706-071505','20150706-071914','20150706-072000','20150706-072409',
'20150706-072504','20150706-072914','20150706-073000']
D8_p10_e = ['20150706-070400','20150706-070456','20150706-070905','20150706-071954','20150706-072400',
'20150706-071457','20150706-071905','20150706-071954','20150706-072400','20150706-072456',
'20150706-072905','20150706-072957','20150706-073401']
D8_p10_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_p10_s])
D8_p10_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_p10_e])
SR1_p10_s = ['20150706-070016','20150706-070507','20150706-071007','20150706-071508','20150706-072008',
'20150706-072507','20150706-073007']
SR1_p10_e = ['20150706-070438','20150706-070930','20150706-071431','20150706-071931','20150706-072430',
'20150706-072930','20150706-073430']
SR1_p10_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_p10_s])
SR1_p10_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_p10_e])
# SR2_p10_s = ['20150706-']
# SR2_p10_e = ['20150706-']
# SR2_p10_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_p10_s])
# SR2_p10_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_p10_e])
RX_p10_s = ['20150706-065947','20150706-070446','20150706-070945','20150706-071444','20150706-071943',
'20150706-072442','20150706-073000']
RX_p10_e = ['20150706-070438','20150706-070936','20150706-071435','20150706-071934','20150706-072433',
'20150706-072952','20150706-073450']
RX_p10_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_p10_s])
RX_p10_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_p10_e])
XP_p10_s = ['20150706-070010','20150706-070512','20150706-071010','20150706-071511','20150706-072010',
'20150706-072510','20150706-073011']
XP_p10_e = ['20150706-070441','20150706-070942','20150706-071442','20150706-071941','20150706-072440',
'20150706-072941','20150706-073442']
XP_p10_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_p10_s])
XP_p10_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_p10_e])
FSD_p10_s = ['20150706-070052','20150706-070545','20150706-071037','20150706-071528','20150706-072019',
'20150706-072511','20150706-073002']
FSD_p10_e = ['20150706-070434','20150706-070926','20150706-071417','20150706-071908','20150706-072400',
'20150706-072851','20150706-073342']
FSD_p10_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_p10_s])
FSD_p10_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_p10_e])
ABR_p10_s = ['20150706-070117','20150706-070557','20150706-071040','20150706-071534','20150706-072028',
'20150706-072746']
ABR_p10_e = ['20150706-070442','20150706-070923','20150706-071420','20150706-071914','20150706-072421',
'20150706-073139']
ABR_p10_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_p10_s])
ABR_p10_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_p10_e])
############# Plotting #############
labs, tickloc, col = [], [], []
color=iter(cm.Dark2(np.linspace(0,1,radNum)))
fig, ax = plt.subplots(1,figsize=(15,10))
plt.title('6 July 2015 - PDD ' + str(tNum) + ' - Radar Volume Overlaps',fontsize=28)
c = next(color)
plt.hlines(radNum,TDR_pdds_s_DT[tNum-1],TDR_pdds_e_DT[tNum-1],color=c,linewidth=40)
labs.append('TDR')
tickloc.append(radNum)
col.append(c)
radNum -= 1
c = next(color)
labs,tickloc,col,radNum = plotGantt(D6_p10_s_DT,D6_p10_e_DT,c,'DOW6',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D7_p10_s_DT,D7_p10_e_DT,c,'DOW7',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D8_p10_s_DT,D8_p10_e_DT,c,'DOW8',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR1_p10_s_DT,SR1_p10_e_DT,c,'SR1',labs,tickloc,col,radNum)
# c = next(color)
# labs,tickloc,col,radNum = plotGantt(SR2_p10_s_DT,SR2_p10_e_DT,c,'SR2',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(RX_p10_s_DT,RX_p10_e_DT,c,'RAXPOL',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(XP_p10_s_DT,XP_p10_e_DT,c,'NOXP',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(FSD_p10_s_DT,FSD_p10_e_DT,c,'KFSD',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(ABR_p10_s_DT,ABR_p10_e_DT,c,'KABR',labs,tickloc,col,radNum)
plt.ylim(-0.25,totalRad+1.25)
plt.yticks(tickloc, labs)
plt.grid(which='both',axis='x')
plt.xlabel('Time (UTC)',fontsize=20)
plt.xlim((ax.get_xlim()[0]-0.0007),(ax.get_xlim()[1]+0.0007))
ytl = plt.gca().get_yticklabels()
for i in range(len(ytl)):
ytl[i].set_color(col[i])
plt.tick_params(axis='both', which='major', labelsize=20)
fig.autofmt_xdate()
ax.xaxis.set_major_locator(mdates.MinuteLocator(interval=2))
ax.xaxis.set_minor_locator(mdates.MinuteLocator(interval=1))
ax.xaxis.set_major_formatter(dfmt)
saveStr = '{}{}_PDD{:02d}_RadarVolTimes.{}'.format(savePath,flight,tNum,fType)
fig.savefig(saveStr,bbox_inches='tight')
In [ ]:
#
tNum = 1
totalRad = 10
radNum = totalRad
############# Start/End Times for Radar Volumes #############
D6_s01_s = ['20150706-031625','20150706-032034','20150706-032303','20150706-032714','20150706-032810',
'20150706-032958']
D6_s01_e = ['20150706-032024','20150706-032049','20150706-032705','20150706-032801','20150706-032929',
'20150706-033400']
D6_s01_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_s01_s])
D6_s01_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_s01_e])
D7_s01_s = ['20150706-031619','20150706-032032','20150706-032107','20150706-032508','20150706-032604',
'20150706-032959']
D7_s01_e = ['20150706-032018','20150706-032043','20150706-032459','20150706-032555','20150706-032946',
'20150706-033357']
D7_s01_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_s01_s])
D7_s01_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_s01_e])
D8_s01_s = ['20150706-031959','20150706-032410','20150706-032506','20150706-032914','20150706-033000']
D8_s01_e = ['20150706-032406','20150706-032457','20150706-032906','20150706-032957','20150706-033401']
D8_s01_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_s01_s])
D8_s01_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_s01_e])
SR1_s01_s = ['20150706-032007','20150706-032508','20150706-033008']
SR1_s01_e = ['20150706-032411','20150706-032911','20150706-033412']
SR1_s01_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_s01_s])
SR1_s01_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_s01_e])
SR2_s01_s = ['20150706-032004','20150706-032504','20150706-033016']
SR2_s01_e = ['20150706-032412','20150706-032913','20150706-033434']
SR2_s01_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_s01_s])
SR2_s01_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_s01_e])
RX_s01_s = ['20150706-031940','20150706-032438','20150706-032936']
RX_s01_e = ['20150706-032413','20150706-032911','20150706-033409']
RX_s01_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_s01_s])
RX_s01_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_s01_e])
XP_s01_s = ['20150706-032006','20150706-032509','20150706-033007']
XP_s01_e = ['20150706-032403','20150706-032905','20150706-033404']
XP_s01_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_s01_s])
XP_s01_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_s01_e])
FSD_s01_s = ['20150706-031751','20150706-032242','20150706-032732']
FSD_s01_e = ['20150706-032132','20150706-032622','20150706-033113']
FSD_s01_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_s01_s])
FSD_s01_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_s01_e])
ABR_s01_s = ['20150706-031835','20150706-032327','20150706-032821']
ABR_s01_e = ['20150706-032214','20150706-032707','20150706-033214']
ABR_s01_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_s01_s])
ABR_s01_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_s01_e])
############# Plotting #############
labs, tickloc, col = [], [], []
color=iter(cm.Dark2(np.linspace(0,1,radNum)))
fig, ax = plt.subplots(1,figsize=(15,10))
plt.title('6 July 2015 - Spiral ' + str(tNum) + ' - Radar Volume Overlaps',fontsize=28)
c = next(color)
plt.hlines(radNum,TDR_sprls_s_DT[tNum-1],TDR_sprls_e_DT[tNum-1],color=c,linewidth=40)
labs.append('TDR')
tickloc.append(radNum)
col.append(c)
radNum -= 1
c = next(color)
labs,tickloc,col,radNum = plotGantt(D6_s01_s_DT,D6_s01_e_DT,c,'DOW6',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D7_s01_s_DT,D7_s01_e_DT,c,'DOW7',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D8_s01_s_DT,D8_s01_e_DT,c,'DOW8',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR1_s01_s_DT,SR1_s01_e_DT,c,'SR1',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR2_s01_s_DT,SR2_s01_e_DT,c,'SR2',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(RX_s01_s_DT,RX_s01_e_DT,c,'RAXPOL',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(XP_s01_s_DT,XP_s01_e_DT,c,'NOXP',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(FSD_s01_s_DT,FSD_s01_e_DT,c,'KFSD',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(ABR_s01_s_DT,ABR_s01_e_DT,c,'KABR',labs,tickloc,col,radNum)
plt.ylim(-0.25,totalRad+1.25)
plt.yticks(tickloc, labs)
plt.grid(which='both',axis='x')
plt.xlabel('Time (UTC)',fontsize=20)
plt.xlim((ax.get_xlim()[0]-0.0007),(ax.get_xlim()[1]+0.0007))
ytl = plt.gca().get_yticklabels()
for i in range(len(ytl)):
ytl[i].set_color(col[i])
plt.tick_params(axis='both', which='major', labelsize=20)
fig.autofmt_xdate()
ax.xaxis.set_major_locator(mdates.MinuteLocator(interval=2))
ax.xaxis.set_minor_locator(mdates.MinuteLocator(interval=1))
ax.xaxis.set_major_formatter(dfmt)
saveStr = '{}{}_Sprl{:02d}_RadarVolTimes.{}'.format(savePath,flight,tNum,fType)
fig.savefig(saveStr,bbox_inches='tight')
In [ ]:
#
tNum = 2
totalRad = 10
radNum = totalRad
############# Start/End Times for Radar Volumes #############
D6_s02_s = ['20150706-032958','20150706-033409','20150706-033505','20150706-033914','20150706-034000',
'20150706-034411','20150706-034507']
D6_s02_e = ['20150706-033400','20150706-033456','20150706-033930','20150706-033930','20150706-034401',
'20150706-034458','20150706-034906']
D6_s02_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_s02_s])
D6_s02_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_s02_e])
D7_s02_s = ['20150706-032959','20150706-033407','20150706-033503','20150706-033912','20150706-034000',
'20150706-034411','20150706-034507']
D7_s02_e = ['20150706-033357','20150706-033454','20150706-033902','20150706-033927','20150706-034402',
'20150706-034458','20150706-034907']
D7_s02_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_s02_s])
D7_s02_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_s02_e])
D8_s02_s = ['20150706-033000','20150706-033409','20150706-033505','20150706-033914','20150706-033959',
'20150706-034409','20150706-034504']
D8_s02_e = ['20150706-033401','20150706-033457','20150706-033905','20150706-033950','20150706-034400',
'20150706-034457','20150706-034905']
D8_s02_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_s02_s])
D8_s02_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_s02_e])
SR1_s02_s = ['20150706-033008','20150706-033508','20150706-034007','20150706-034507']
SR1_s02_e = ['20150706-033412','20150706-033912','20150706-033413','20150706-034911']
SR1_s02_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_s02_s])
SR1_s02_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_s02_e])
SR2_s02_s = ['20150706-033016','20150706-033505','20150706-034005','20150706-034505']
SR2_s02_e = ['20150706-033434','20150706-033923','20150706-034423','20150706-034924']
SR2_s02_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_s02_s])
SR2_s02_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_s02_e])
RX_s02_s = ['20150706-032936','20150706-033435','20150706-033933','20150706-034432']
RX_s02_e = ['20150706-033409','20150706-033908','20150706-034406','20150706-034905']
RX_s02_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_s02_s])
RX_s02_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_s02_e])
XP_s02_s = ['20150706-033007','20150706-033506','20150706-034006','20150706-034506']
XP_s02_e = ['20150706-033404','20150706-033902','20150706-034403','20150706-034902']
XP_s02_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_s02_s])
XP_s02_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_s02_e])
FSD_s02_s = ['20150706-032732','20150706-033224','20150706-033715','20150706-034205']
FSD_s02_e = ['20150706-033113','20150706-033605','20150706-034055','20150706-034546']
FSD_s02_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_s02_s])
FSD_s02_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_s02_e])
ABR_s02_s = ['20150706-032821','20150706-033329','20150706-033838','20150706-034346']
ABR_s02_e = ['20150706-033214','20150706-033723','20150706-034231','20150706-034739']
ABR_s02_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_s02_s])
ABR_s02_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_s02_e])
############# Plotting #############
labs, tickloc, col = [], [], []
color=iter(cm.Dark2(np.linspace(0,1,radNum)))
fig, ax = plt.subplots(1,figsize=(15,10))
plt.title('6 July 2015 - Spiral ' + str(tNum) + ' - Radar Volume Overlaps',fontsize=28)
c = next(color)
plt.hlines(radNum,TDR_sprls_s_DT[tNum-1],TDR_sprls_e_DT[tNum-1],color=c,linewidth=40)
labs.append('TDR')
tickloc.append(radNum)
col.append(c)
radNum -= 1
c = next(color)
labs,tickloc,col,radNum = plotGantt(D6_s02_s_DT,D6_s02_e_DT,c,'DOW6',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D7_s02_s_DT,D7_s02_e_DT,c,'DOW7',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D8_s02_s_DT,D8_s02_e_DT,c,'DOW8',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR1_s02_s_DT,SR1_s02_e_DT,c,'SR1',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR2_s02_s_DT,SR2_s02_e_DT,c,'SR2',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(RX_s02_s_DT,RX_s02_e_DT,c,'RAXPOL',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(XP_s02_s_DT,XP_s02_e_DT,c,'NOXP',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(FSD_s02_s_DT,FSD_s02_e_DT,c,'KFSD',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(ABR_s02_s_DT,ABR_s02_e_DT,c,'KABR',labs,tickloc,col,radNum)
plt.ylim(-0.25,totalRad+1.25)
plt.yticks(tickloc, labs)
plt.grid(which='both',axis='x')
plt.xlabel('Time (UTC)',fontsize=20)
plt.xlim((ax.get_xlim()[0]-0.0007),(ax.get_xlim()[1]+0.0007))
ytl = plt.gca().get_yticklabels()
for i in range(len(ytl)):
ytl[i].set_color(col[i])
plt.tick_params(axis='both', which='major', labelsize=20)
fig.autofmt_xdate()
ax.xaxis.set_major_locator(mdates.MinuteLocator(interval=2))
ax.xaxis.set_minor_locator(mdates.MinuteLocator(interval=1))
ax.xaxis.set_major_formatter(dfmt)
saveStr = '{}{}_Sprl{:02d}_RadarVolTimes.{}'.format(savePath,flight,tNum,fType)
fig.savefig(saveStr,bbox_inches='tight')
In [ ]:
#
tNum = 3
totalRad = 10
radNum = totalRad
############# Start/End Times for Radar Volumes #############
D6_s03_s = ['20150706-042000','20150706-042411','20150706-042507','20150706-042916','20150706-043001',
'20150706-043412','20150706-043508','20150706-043917']
D6_s03_e = ['20150706-042402','20150706-042458','20150706-042906','20150706-042931','20150706-043402',
'20150706-043459','20150706-043907','20150706-043932']
D6_s03_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_s03_s])
D6_s03_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_s03_e])
D7_s03_s = ['20150706-042000','20150706-042411','20150706-042507','20150706-042916','20150706-043001',
'20150706-043412','20150706-043508','20150706-043917']
D7_s03_e = ['20150706-042401','20150706-042458','20150706-042906','20150706-042931','20150706-043402',
'20150706-043459','20150706-043907','20150706-043932']
D7_s03_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_s03_s])
D7_s03_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_s03_e])
D8_s03_s = ['20150706-042000','20150706-042409','20150706-042504','20150706-042914','20150706-043000',
'20150706-043409','20150706-043505','20150706-043914']
D8_s03_e = ['20150706-042400','20150706-042456','20150706-042905','20150706-042949','20150706-043400',
'20150706-043457','20150706-043905','20150706-043950']
D8_s03_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_s03_s])
D8_s03_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_s03_e])
SR1_s03_s = ['20150706-042016','20150706-042511','20150706-043011','20150706-043510']
SR1_s03_e = ['20150706-042428','20150706-042928','20150706-043431','20150706-043930']
SR1_s03_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_s03_s])
SR1_s03_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_s03_e])
SR2_s03_s = ['20150706-042011','20150706-042511','20150706-043011','20150706-043511']
SR2_s03_e = ['20150706-042435','20150706-042935','20150706-043435','20150706-043935']
SR2_s03_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_s03_s])
SR2_s03_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_s03_e])
RX_s03_s = ['20150706-042017','20150706-042457','20150706-043000','20150706-043507']
RX_s03_e = ['20150706-042441','20150706-042944','20150706-043417','20150706-044011']
RX_s03_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_s03_s])
RX_s03_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_s03_e])
XP_s03_s = ['20150706-042011','20150706-042511','20150706-043010','20150706-043511']
XP_s03_e = ['20150706-042440','20150706-042943','20150706-043441','20150706-043941']
XP_s03_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_s03_s])
XP_s03_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_s03_e])
FSD_s03_s = ['20150706-042053','20150706-042545','20150706-043035','20150706-043526']
FSD_s03_e = ['20150706-042434','20150706-042925','20150706-043416','20150706-043907']
FSD_s03_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_s03_s])
FSD_s03_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_s03_e])
ABR_s03_s = ['20150706-042449','20150706-042957','20150706-043504']
ABR_s03_e = ['20150706-042842','20150706-043350','20150706-043857']
ABR_s03_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_s03_s])
ABR_s03_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_s03_e])
############# Plotting #############
labs, tickloc, col = [], [], []
color=iter(cm.Dark2(np.linspace(0,1,radNum)))
fig, ax = plt.subplots(1,figsize=(15,10))
plt.title('6 July 2015 - Spiral ' + str(tNum) + ' - Radar Volume Overlaps',fontsize=28)
c = next(color)
plt.hlines(radNum,TDR_sprls_s_DT[tNum-1],TDR_sprls_e_DT[tNum-1],color=c,linewidth=40)
labs.append('TDR')
tickloc.append(radNum)
col.append(c)
radNum -= 1
c = next(color)
labs,tickloc,col,radNum = plotGantt(D6_s03_s_DT,D6_s03_e_DT,c,'DOW6',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D7_s03_s_DT,D7_s03_e_DT,c,'DOW7',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D8_s03_s_DT,D8_s03_e_DT,c,'DOW8',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR1_s03_s_DT,SR1_s03_e_DT,c,'SR1',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR2_s03_s_DT,SR2_s03_e_DT,c,'SR2',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(RX_s03_s_DT,RX_s03_e_DT,c,'RAXPOL',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(XP_s03_s_DT,XP_s03_e_DT,c,'NOXP',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(FSD_s03_s_DT,FSD_s03_e_DT,c,'KFSD',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(ABR_s03_s_DT,ABR_s03_e_DT,c,'KABR',labs,tickloc,col,radNum)
plt.ylim(-0.25,totalRad+1.25)
plt.yticks(tickloc, labs)
plt.grid(which='both',axis='x')
plt.xlabel('Time (UTC)',fontsize=20)
plt.xlim((ax.get_xlim()[0]-0.0007),(ax.get_xlim()[1]+0.0007))
ytl = plt.gca().get_yticklabels()
for i in range(len(ytl)):
ytl[i].set_color(col[i])
plt.tick_params(axis='both', which='major', labelsize=20)
fig.autofmt_xdate()
ax.xaxis.set_major_locator(mdates.MinuteLocator(interval=2))
ax.xaxis.set_minor_locator(mdates.MinuteLocator(interval=1))
ax.xaxis.set_major_formatter(dfmt)
saveStr = '{}{}_Sprl{:02d}_RadarVolTimes.{}'.format(savePath,flight,tNum,fType)
fig.savefig(saveStr,bbox_inches='tight')
In [ ]:
#
tNum = 4
totalRad = 10
radNum = totalRad
############# Start/End Times for Radar Volumes #############
D6_s04_s = ['20150706-044000','20150706-044411','20150706-044507','20150706-044916','20150706-045001']
D6_s04_e = ['20150706-044402','20150706-044458','20150706-044907','20150706-044931','20150706-045402']
D6_s04_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_s04_s])
D6_s04_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_s04_e])
D7_s04_s = ['20150706-044000','20150706-044411','20150706-044508','20150706-044916','20150706-045000']
D7_s04_e = ['20150706-044401','20150706-044458','20150706-044906','20150706-044931','20150706-045402']
D7_s04_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_s04_s])
D7_s04_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_s04_e])
D8_s04_s = ['20150706-043959','20150706-044409','20150706-044504','20150706-044914','20150706-045001']
D8_s04_e = ['20150706-044400','20150706-044456','20150706-044905','20150706-044937','20150706-045400']
D8_s04_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_s04_s])
D8_s04_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_s04_e])
SR1_s04_s = ['20150706-044011','20150706-044511','20150706-045011']
SR1_s04_e = ['20150706-044435','20150706-044936','20150706-045437']
SR1_s04_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_s04_s])
SR1_s04_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_s04_e])
SR2_s04_s = ['20150706-044010','20150706-044510','20150706-045011']
SR2_s04_e = ['20150706-044434','20150706-044934','20150706-045435']
SR2_s04_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_s04_s])
SR2_s04_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_s04_e])
RX_s04_s = ['20150706-044020','20150706-044507','20150706-045008']
RX_s04_e = ['20150706-044439','20150706-044959','20150706-045508']
RX_s04_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_s04_s])
RX_s04_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_s04_e])
XP_s04_s = ['20150706-044010','20150706-044511','20150706-045011']
XP_s04_e = ['20150706-044440','20150706-044942','20150706-045442']
XP_s04_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_s04_s])
XP_s04_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_s04_e])
FSD_s04_s = ['20150706-044018','20150706-044509','20150706-044959']
FSD_s04_e = ['20150706-044358','20150706-044849','20150706-045339']
FSD_s04_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_s04_s])
FSD_s04_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_s04_e])
ABR_s04_s = ['20150706-044012','20150706-044520','20150706-045027']
ABR_s04_e = ['20150706-044405','20150706-044913','20150706-045420']
ABR_s04_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_s04_s])
ABR_s04_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_s04_e])
############# Plotting #############
labs, tickloc, col = [], [], []
color=iter(cm.Dark2(np.linspace(0,1,radNum)))
fig, ax = plt.subplots(1,figsize=(15,10))
plt.title('6 July 2015 - Spiral ' + str(tNum) + ' - Radar Volume Overlaps',fontsize=28)
c = next(color)
plt.hlines(radNum,TDR_sprls_s_DT[tNum-1],TDR_sprls_e_DT[tNum-1],color=c,linewidth=40)
labs.append('TDR')
tickloc.append(radNum)
col.append(c)
radNum -= 1
c = next(color)
labs,tickloc,col,radNum = plotGantt(D6_s04_s_DT,D6_s04_e_DT,c,'DOW6',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D7_s04_s_DT,D7_s04_e_DT,c,'DOW7',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D8_s04_s_DT,D8_s04_e_DT,c,'DOW8',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR1_s04_s_DT,SR1_s04_e_DT,c,'SR1',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR2_s04_s_DT,SR2_s04_e_DT,c,'SR2',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(RX_s04_s_DT,RX_s04_e_DT,c,'RAXPOL',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(XP_s04_s_DT,XP_s04_e_DT,c,'NOXP',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(FSD_s04_s_DT,FSD_s04_e_DT,c,'KFSD',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(ABR_s04_s_DT,ABR_s04_e_DT,c,'KABR',labs,tickloc,col,radNum)
plt.ylim(-0.25,totalRad+1.25)
plt.yticks(tickloc, labs)
plt.grid(which='both',axis='x')
plt.xlabel('Time (UTC)',fontsize=20)
plt.xlim((ax.get_xlim()[0]-0.0007),(ax.get_xlim()[1]+0.0007))
ytl = plt.gca().get_yticklabels()
for i in range(len(ytl)):
ytl[i].set_color(col[i])
plt.tick_params(axis='both', which='major', labelsize=20)
fig.autofmt_xdate()
ax.xaxis.set_major_locator(mdates.MinuteLocator(interval=2))
ax.xaxis.set_minor_locator(mdates.MinuteLocator(interval=1))
ax.xaxis.set_major_formatter(dfmt)
saveStr = '{}{}_Sprl{:02d}_RadarVolTimes.{}'.format(savePath,flight,tNum,fType)
fig.savefig(saveStr,bbox_inches='tight')
In [ ]:
#
tNum = 5
totalRad = 10
radNum = totalRad
############# Start/End Times for Radar Volumes #############
D6_s05_s = ['20150706-054000','20150706-054411','20150706-054507','20150706-054916','20150706-055000',
'20150706-055412','20150706-055508']
D6_s05_e = ['20150706-054401','20150706-054458','20150706-054906','20150706-054931','20150706-055402',
'20150706-055459','20150706-055907']
D6_s05_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_s05_s])
D6_s05_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_s05_e])
D7_s05_s = ['20150706-054015','20150706-054425','20150706-054521','20150706-055001','20150706-055201',
'20150706-055535']
D7_s05_e = ['20150706-054415','20150706-054512','20150706-054920','20150706-055128','20150706-055525',
'20150706-055557']
D7_s05_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_s05_s])
D7_s05_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_s05_e])
D8_s05_s = ['20150706-054000','20150706-054409','20150706-054504','20150706-054914','20150706-055001',
'20150706-055409','20150706-055506']
D8_s05_e = ['20150706-054400','20150706-054456','20150706-054905','20150706-054938','20150706-055401',
'20150706-055458','20150706-055905']
D8_s05_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_s05_s])
D8_s05_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_s05_e])
SR1_s05_s = ['20150706-054026','20150706-054511','20150706-055011','20150706-055516']
SR1_s05_e = ['20150706-054438','20150706-054936','20150706-055435','20150706-055939']
SR1_s05_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_s05_s])
SR1_s05_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_s05_e])
SR2_s05_s = ['20150706-054011','20150706-054511','20150706-055011','20150706-055511']
SR2_s05_e = ['20150706-054435','20150706-054935','20150706-055435','20150706-055935']
SR2_s05_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_s05_s])
SR2_s05_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_s05_e])
RX_s05_s = ['20150706-053958','20150706-054457','20150706-054956','20150706-055455']
RX_s05_e = ['20150706-054449','20150706-054948','20150706-055447','20150706-055946']
RX_s05_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_s05_s])
RX_s05_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_s05_e])
XP_s05_s = ['20150706-054010','20150706-054511','20150706-055010','20150706-055510']
XP_s05_e = ['20150706-054441','20150706-054942','20150706-055440','20150706-055942']
XP_s05_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_s05_s])
XP_s05_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_s05_e])
FSD_s05_s = ['20150706-053829','20150706-054321','20150706-054812','20150706-055302']
FSD_s05_e = ['20150706-054210','20150706-054701','20150706-055152','20150706-055642']
FSD_s05_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_s05_s])
FSD_s05_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_s05_e])
ABR_s05_s = ['20150706-054103','20150706-054557','20150706-055038','20150706-055519']
ABR_s05_e = ['20150706-054444','20150706-054923','20150706-055404','20150706-055844']
ABR_s05_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_s05_s])
ABR_s05_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_s05_e])
############# Plotting #############
labs, tickloc, col = [], [], []
color=iter(cm.Dark2(np.linspace(0,1,radNum)))
fig, ax = plt.subplots(1,figsize=(15,10))
plt.title('6 July 2015 - Spiral ' + str(tNum) + ' - Radar Volume Overlaps',fontsize=28)
c = next(color)
plt.hlines(radNum,TDR_sprls_s_DT[tNum-1],TDR_sprls_e_DT[tNum-1],color=c,linewidth=40)
labs.append('TDR')
tickloc.append(radNum)
col.append(c)
radNum -= 1
c = next(color)
labs,tickloc,col,radNum = plotGantt(D6_s05_s_DT,D6_s05_e_DT,c,'DOW6',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D7_s05_s_DT,D7_s05_e_DT,c,'DOW7',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D8_s05_s_DT,D8_s05_e_DT,c,'DOW8',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR1_s05_s_DT,SR1_s05_e_DT,c,'SR1',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR2_s05_s_DT,SR2_s05_e_DT,c,'SR2',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(RX_s05_s_DT,RX_s05_e_DT,c,'RAXPOL',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(XP_s05_s_DT,XP_s05_e_DT,c,'NOXP',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(FSD_s05_s_DT,FSD_s05_e_DT,c,'KFSD',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(ABR_s05_s_DT,ABR_s05_e_DT,c,'KABR',labs,tickloc,col,radNum)
plt.ylim(-0.25,totalRad+1.25)
plt.yticks(tickloc, labs)
plt.grid(which='both',axis='x')
plt.xlabel('Time (UTC)',fontsize=20)
plt.xlim((ax.get_xlim()[0]-0.0007),(ax.get_xlim()[1]+0.0007))
ytl = plt.gca().get_yticklabels()
for i in range(len(ytl)):
ytl[i].set_color(col[i])
plt.tick_params(axis='both', which='major', labelsize=20)
fig.autofmt_xdate()
ax.xaxis.set_major_locator(mdates.MinuteLocator(interval=2))
ax.xaxis.set_minor_locator(mdates.MinuteLocator(interval=1))
ax.xaxis.set_major_formatter(dfmt)
saveStr = '{}{}_Sprl{:02d}_RadarVolTimes.{}'.format(savePath,flight,tNum,fType)
fig.savefig(saveStr,bbox_inches='tight')
In [ ]:
#
tNum = 6
totalRad = 10
radNum = totalRad
############# Start/End Times for Radar Volumes #############
D6_s06_s = ['20150706-055508','20150706-055917','20150706-060000','20150706-060411','20150706-060508',
'20150706-060916','20150706-061000']
D6_s06_e = ['20150706-055907','20150706-055932','20150706-060401','20150706-060458','20150706-060906',
'20150706-060931','20150706-061402']
D6_s06_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_s06_s])
D6_s06_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_s06_e])
D7_s06_s = ['20150706-055616','20150706-060003','20150706-060332','20150706-060410','20150706-060746',
'20150706-060823','20150706-061047']
D7_s06_e = ['20150706-055938','20150706-060319','20150706-060352','20150706-060733','20150706-060805',
'20150706-061035','20150706-061403']
D7_s06_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_s06_s])
D7_s06_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_s06_e])
D8_s06_s = ['20150706-055506','20150706-055914','20150706-060000','20150706-060409','20150706-060504',
'20150706-060914','20150706-061000']
D8_s06_e = ['20150706-055905','20150706-055950','20150706-060400','20150706-060456','20150706-060905',
'20150706-060957','20150706-061401']
D8_s06_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_s06_s])
D8_s06_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_s06_e])
SR1_s06_s = ['20150706-055516','20150706-060012','20150706-060510','20150706-061011']
SR1_s06_e = ['20150706-055939','20150706-060436','20150706-060937','20150706-061437']
SR1_s06_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_s06_s])
SR1_s06_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_s06_e])
SR2_s06_s = ['20150706-055511','20150706-060014','20150706-060511','20150706-061010']
SR2_s06_e = ['20150706-055935','20150706-060438','20150706-060935','20150706-061434']
SR2_s06_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_s06_s])
SR2_s06_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_s06_e])
RX_s06_s = ['20150706-055455','20150706-055955','20150706-060454','20150706-060953']
RX_s06_e = ['20150706-055946','20150706-060446','20150706-060945','20150706-061444']
RX_s06_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_s06_s])
RX_s06_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_s06_e])
XP_s06_s = ['20150706-055510','20150706-060012','20150706-060510','20150706-061010']
XP_s06_e = ['20150706-055942','20150706-060442','20150706-060942','20150706-061441']
XP_s06_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_s06_s])
XP_s06_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_s06_e])
FSD_s06_s = ['20150706-055302','20150706-055752','20150706-060244','20150706-060735','20150706-061225']
FSD_s06_e = ['20150706-055642','20150706-060133','20150706-060624','20150706-061115','20150706-061606']
FSD_s06_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_s06_s])
FSD_s06_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_s06_e])
ABR_s06_s = ['20150706-055519','20150706-060000','20150706-060441','20150706-060921']
ABR_s06_e = ['20150706-055844','20150706-060326','20150706-060807','20150706-061247']
ABR_s06_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_s06_s])
ABR_s06_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_s06_e])
############# Plotting #############
labs, tickloc, col = [], [], []
color=iter(cm.Dark2(np.linspace(0,1,radNum)))
fig, ax = plt.subplots(1,figsize=(15,10))
plt.title('6 July 2015 - Spiral ' + str(tNum) + ' - Radar Volume Overlaps',fontsize=28)
c = next(color)
plt.hlines(radNum,TDR_sprls_s_DT[tNum-1],TDR_sprls_e_DT[tNum-1],color=c,linewidth=40)
labs.append('TDR')
tickloc.append(radNum)
col.append(c)
radNum -= 1
c = next(color)
labs,tickloc,col,radNum = plotGantt(D6_s06_s_DT,D6_s06_e_DT,c,'DOW6',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D7_s06_s_DT,D7_s06_e_DT,c,'DOW7',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D8_s06_s_DT,D8_s06_e_DT,c,'DOW8',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR1_s06_s_DT,SR1_s06_e_DT,c,'SR1',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR2_s06_s_DT,SR2_s06_e_DT,c,'SR2',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(RX_s06_s_DT,RX_s06_e_DT,c,'RAXPOL',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(XP_s06_s_DT,XP_s06_e_DT,c,'NOXP',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(FSD_s06_s_DT,FSD_s06_e_DT,c,'KFSD',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(ABR_s06_s_DT,ABR_s06_e_DT,c,'KABR',labs,tickloc,col,radNum)
plt.ylim(-0.25,totalRad+1.25)
plt.yticks(tickloc, labs)
plt.grid(which='both',axis='x')
plt.xlabel('Time (UTC)',fontsize=20)
plt.xlim((ax.get_xlim()[0]-0.0007),(ax.get_xlim()[1]+0.0007))
ytl = plt.gca().get_yticklabels()
for i in range(len(ytl)):
ytl[i].set_color(col[i])
plt.tick_params(axis='both', which='major', labelsize=20)
fig.autofmt_xdate()
ax.xaxis.set_major_locator(mdates.MinuteLocator(interval=2))
ax.xaxis.set_minor_locator(mdates.MinuteLocator(interval=1))
ax.xaxis.set_major_formatter(dfmt)
saveStr = '{}{}_Sprl{:02d}_RadarVolTimes.{}'.format(savePath,flight,tNum,fType)
fig.savefig(saveStr,bbox_inches='tight')
In [ ]:
#
tNum = 7
totalRad = 10
radNum = totalRad
############# Start/End Times for Radar Volumes #############
D6_s07_s = ['20150706-062000','20150706-062411','20150706-062507','20150706-062916','20150706-063000',
'20150706-063412','20150706-063508']
D6_s07_e = ['20150706-062402','20150706-062458','20150706-062906','20150706-062931','20150706-063402',
'20150706-063459','20150706-063907']
D6_s07_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_s07_s])
D6_s07_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_s07_e])
D7_s07_s = ['20150706-062003','20150706-062404','20150706-062502','20150706-062909','20150706-063000',
'20150706-063412','20150706-063508']
D7_s07_e = ['20150706-062355','20150706-062451','20150706-062900','20150706-062925','20150706-063402',
'20150706-063459','20150706-063907']
D7_s07_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_s07_s])
D7_s07_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_s07_e])
D8_s07_s = ['20150706-062000','20150706-062409','20150706-062504','20150706-062914','20150706-063001',
'20150706-063409','20150706-063505']
D8_s07_e = ['20150706-062400','20150706-062456','20150706-062905','20150706-062938','20150706-063401',
'20150706-063458','20150706-063905']
D8_s07_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_s07_s])
D8_s07_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_s07_e])
SR1_s07_s = ['20150706-062013','20150706-062516','20150706-063012','20150706-063517']
SR1_s07_e = ['20150706-062440','20150706-062944','20150706-063438','20150706-063945']
SR1_s07_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_s07_s])
SR1_s07_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_s07_e])
SR2_s07_s = ['20150706-062010','20150706-062510','20150706-063010','20150706-063511']
SR2_s07_e = ['20150706-062434','20150706-062934','20150706-063434','20150706-063935']
SR2_s07_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_s07_s])
SR2_s07_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_s07_e])
RX_s07_s = ['20150706-061952','20150706-062451','20150706-062951','20150706-063450']
RX_s07_e = ['20150706-062443','20150706-062942','20150706-063442','20150706-063941']
RX_s07_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_s07_s])
RX_s07_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_s07_e])
XP_s07_s = ['20150706-062011','20150706-062511','20150706-063010','20150706-063510']
XP_s07_e = ['20150706-062442','20150706-062941','20150706-063441','20150706-063941']
XP_s07_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_s07_s])
XP_s07_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_s07_e])
FSD_s07_s = ['20150706-062206','20150706-062656','20150706-063147','20150706-063637']
FSD_s07_e = ['20150706-062546','20150706-063036','20150706-063527','20150706-064018']
FSD_s07_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_s07_s])
FSD_s07_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_s07_e])
ABR_s07_s = ['20150706-062257','20150706-062725','20150706-063152','20150706-063659']
ABR_s07_e = ['20150706-062610','20150706-063037','20150706-063544','20150706-064025']
ABR_s07_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_s07_s])
ABR_s07_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_s07_e])
############# Plotting #############
labs, tickloc, col = [], [], []
color=iter(cm.Dark2(np.linspace(0,1,radNum)))
fig, ax = plt.subplots(1,figsize=(15,10))
plt.title('6 July 2015 - Spiral ' + str(tNum) + ' - Radar Volume Overlaps',fontsize=28)
c = next(color)
plt.hlines(radNum,TDR_sprls_s_DT[tNum-1],TDR_sprls_e_DT[tNum-1],color=c,linewidth=40)
labs.append('TDR')
tickloc.append(radNum)
col.append(c)
radNum -= 1
c = next(color)
labs,tickloc,col,radNum = plotGantt(D6_s07_s_DT,D6_s07_e_DT,c,'DOW6',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D7_s07_s_DT,D7_s07_e_DT,c,'DOW7',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D8_s07_s_DT,D8_s07_e_DT,c,'DOW8',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR1_s07_s_DT,SR1_s07_e_DT,c,'SR1',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR2_s07_s_DT,SR2_s07_e_DT,c,'SR2',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(RX_s07_s_DT,RX_s07_e_DT,c,'RAXPOL',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(XP_s07_s_DT,XP_s07_e_DT,c,'NOXP',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(FSD_s07_s_DT,FSD_s07_e_DT,c,'KFSD',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(ABR_s07_s_DT,ABR_s07_e_DT,c,'KABR',labs,tickloc,col,radNum)
plt.ylim(-0.25,totalRad+1.25)
plt.yticks(tickloc, labs)
plt.grid(which='both',axis='x')
plt.xlabel('Time (UTC)',fontsize=20)
plt.xlim((ax.get_xlim()[0]-0.0007),(ax.get_xlim()[1]+0.0007))
ytl = plt.gca().get_yticklabels()
for i in range(len(ytl)):
ytl[i].set_color(col[i])
plt.tick_params(axis='both', which='major', labelsize=20)
fig.autofmt_xdate()
ax.xaxis.set_major_locator(mdates.MinuteLocator(interval=2))
ax.xaxis.set_minor_locator(mdates.MinuteLocator(interval=1))
ax.xaxis.set_major_formatter(dfmt)
saveStr = '{}{}_Sprl{:02d}_RadarVolTimes.{}'.format(savePath,flight,tNum,fType)
fig.savefig(saveStr,bbox_inches='tight')
In [ ]:
#
tNum = 8
totalRad = 10
radNum = totalRad
############# Start/End Times for Radar Volumes #############
D6_s08_s = ['20150706-063508','20150706-063917','20150706-064000','20150706-064411','20150706-064507',
'20150706-064916','20150706-065001','20150706-065412']
D6_s08_e = ['20150706-063907','20150706-063932','20150706-064402','20150706-064458','20150706-064907',
'20150706-064931','20150706-065402','20150706-065459']
D6_s08_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_s08_s])
D6_s08_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D6_s08_e])
D7_s08_s = ['20150706-063508','20150706-063917','20150706-064000','20150706-064411','20150706-064507',
'20150706-064916','20150706-065000','20150706-065412']
D7_s08_e = ['20150706-063907','20150706-063932','20150706-064402','20150706-064458','20150706-064907',
'20150706-064931','20150706-065402','20150706-065459']
D7_s08_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_s08_s])
D7_s08_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D7_s08_e])
D8_s08_s = ['20150706-063505','20150706-063914','20150706-064000','20150706-064408','20150706-064504',
'20150706-064913','20150706-065000','20150706-065409']
D8_s08_e = ['20150706-063905','20150706-063938','20150706-064400','20150706-064456','20150706-064904',
'20150706-064945','20150706-065401','20150706-065457']
D8_s08_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_s08_s])
D8_s08_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in D8_s08_e])
SR1_s08_s = ['20150706-063517','20150706-064020','20150706-064513','20150706-065018']
SR1_s08_e = ['20150706-063945','20150706-064449','20150706-064947','20150706-065443']
SR1_s08_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_s08_s])
SR1_s08_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR1_s08_e])
SR2_s08_s = ['20150706-063511','20150706-064011','20150706-064511','20150706-065011']
SR2_s08_e = ['20150706-063935','20150706-064435','20150706-064934','20150706-065434']
SR2_s08_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_s08_s])
SR2_s08_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in SR2_s08_e])
RX_s08_s = ['20150706-063450','20150706-063950','20150706-064450','20150706-064949']
RX_s08_e = ['20150706-063941','20150706-064441','20150706-064949','20150706-065440']
RX_s08_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_s08_s])
RX_s08_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in RX_s08_e])
XP_s08_s = ['20150706-063510','20150706-064011','20150706-064511','20150706-065011']
XP_s08_e = ['20150706-063941','20150706-064441','20150706-064942','20150706-065442']
XP_s08_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_s08_s])
XP_s08_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in XP_s08_e])
FSD_s08_s = ['20150706-063637','20150706-064128','20150706-064619','20150706-065109']
FSD_s08_e = ['20150706-064018','20150706-064508','20150706-064959','20150706-065449']
FSD_s08_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_s08_s])
FSD_s08_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in FSD_s08_e])
ABR_s08_s = ['20150706-063659','20150706-064140','20150706-064701','20150706-065142']
ABR_s08_e = ['20150706-064025','20150706-064506','20150706-065027','20150706-065509']
ABR_s08_s_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_s08_s])
ABR_s08_e_DT = np.asarray([dt.strptime(fDate,'%Y%m%d-%H%M%S') for fDate in ABR_s08_e])
############# Plotting #############
labs, tickloc, col = [], [], []
color=iter(cm.Dark2(np.linspace(0,1,radNum)))
fig, ax = plt.subplots(1,figsize=(15,10))
plt.title('6 July 2015 - Spiral ' + str(tNum) + ' - Radar Volume Overlaps',fontsize=28)
c = next(color)
plt.hlines(radNum,TDR_sprls_s_DT[tNum-1],TDR_sprls_e_DT[tNum-1],color=c,linewidth=40)
labs.append('TDR')
tickloc.append(radNum)
col.append(c)
radNum -= 1
c = next(color)
labs,tickloc,col,radNum = plotGantt(D6_s08_s_DT,D6_s08_e_DT,c,'DOW6',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D7_s08_s_DT,D7_s08_e_DT,c,'DOW7',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(D8_s08_s_DT,D8_s08_e_DT,c,'DOW8',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR1_s08_s_DT,SR1_s08_e_DT,c,'SR1',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(SR2_s08_s_DT,SR2_s08_e_DT,c,'SR2',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(RX_s08_s_DT,RX_s08_e_DT,c,'RAXPOL',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(XP_s08_s_DT,XP_s08_e_DT,c,'NOXP',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(FSD_s08_s_DT,FSD_s08_e_DT,c,'KFSD',labs,tickloc,col,radNum)
c = next(color)
labs,tickloc,col,radNum = plotGantt(ABR_s08_s_DT,ABR_s08_e_DT,c,'KABR',labs,tickloc,col,radNum)
plt.ylim(-0.25,totalRad+1.25)
plt.yticks(tickloc, labs)
plt.grid(which='both',axis='x')
plt.xlabel('Time (UTC)',fontsize=20)
plt.xlim((ax.get_xlim()[0]-0.0007),(ax.get_xlim()[1]+0.0007))
ytl = plt.gca().get_yticklabels()
for i in range(len(ytl)):
ytl[i].set_color(col[i])
plt.tick_params(axis='both', which='major', labelsize=20)
fig.autofmt_xdate()
ax.xaxis.set_major_locator(mdates.MinuteLocator(interval=2))
ax.xaxis.set_minor_locator(mdates.MinuteLocator(interval=1))
ax.xaxis.set_major_formatter(dfmt)
saveStr = '{}{}_Sprl{:02d}_RadarVolTimes.{}'.format(savePath,flight,tNum,fType)
fig.savefig(saveStr,bbox_inches='tight')