year = 2015
doy = 76
instr_nums = [5] # In N-to-S order


fig = figure(figsize=(4,2))

# Setup basic stuff based on a working instr
instr_ref = 'minime05'
tZ, _, _, _, _, _, _, _, _ = get_data_level1(instr_ref, year, doy, 'Zenith')
tref = datetime(tZ[-1].year, tZ[-1].month, tZ[-1].day) # midnight
t0 = tref - timedelta(hours=1)
t1 = tref + timedelta(hours=12)

successes = []
for i in range(len(instr_nums)):
    instr_name = 'minime%02i' % instr_nums[i]

    tN, wN, weN, IN, qN, _, _, _, _  = get_data_level1(instr_name, year, doy, 'North')
    tS, wS, weS, IS, qE, _, _, _, _  = get_data_level1(instr_name, year, doy, 'South')
    tE, wE, weE, IE, qS, _, _, _, _  = get_data_level1(instr_name, year, doy, 'East')
    tW, wW, weW, IW, qW, _, _, _, _  = get_data_level1(instr_name, year, doy, 'West')
    tZ, wZ, weZ, IZ, qZ, _, _, _, _  = get_data_level1(instr_name, year, doy, 'Zenith')
    t_ave, ulo, umid, uhi, vlo,vmid, vhi = average_winds(instr_name, year, doy)

    ############## MERID + VERTICAL WINDS ##############
    plot(tN,wN,'g.-',label='Merid North-look')
    plot(tS,wS,'r.-',label='Merid South-look')  
    #plot(tZ,wZ,'k.-',label='Vertical')  
    fill_between(t_ave, vlo, vhi, facecolor='gray', edgecolor='none')
    plot([], [], color='gray', linewidth=6, label='Merid Variability') #to show on legend    
    legend(loc='best',ncol=3, prop={'size':5})
    plt.gca().xaxis.set_major_formatter(matplotlib.dates.DateFormatter('%H'))
    grid()
    xlim((t0,t1))
    ylim((-250,250)) 
    ylabel('Meridional Wind [m/s]')
    xlabel('UT [hr]')
fig.savefig('/home/bhardin2/ipython_figures/agu_2015/meridional_%s_%i_%03i.pdf' % (instr_name, year, doy))

years =      [2015,  2015,  2015,      2015,      2015,      2015,    2013,    2014, ]
doys  =      [76,    191,   227,       250,       257,       263,     190,     143, ]
instr_nums = [[5,6], [5,6], [8,5,9,6], [8,5,9,6], [8,5,9,6], [8,5,6], [8,5,6], [5,9,6]]

delta = False

for j in range(len(years)):
    fig = hodochrone(instr_nums[j], years[j], doys[j], delta=delta)
    

Difference from quiet behavior

years =      [2015,  2015,  2015,      2015,      2015,      2015,    2013,    2014, ]
doys  =      [76,    191,   227,       250,       257,       263,     190,     143, ]
instr_nums = [[5,6], [5,6], [8,5,9,6], [8,5,9,6], [8,5,9,6], [8,5,6], [8,5,6], [5,9,6]]

delta = True

for j in range(len(years)):
    fig = hodochrone(instr_nums[j], years[j], doys[j], delta=delta)
    

years =      [2015,  2015,  2015,      2015,      2015,      2015,    2013,    2014, ]
doys  =      [76,    191,   227,       250,       257,       263,     190,     143, ]
instr_nums = [[5,6], [5,6], [8,5,9,6], [8,5,9,6], [8,5,9,6], [8,5,6], [8,5,6], [5,9,6]]

for j in range(len(years)):
    fig = hodochrone_brightness(instr_nums[j], years[j], doys[j])
    

import subprocess
import shutil
nb_name = 'AGU_2015_Presentation_1.ipynb' # There's probably a way to get this programmatically
cmd = ['jupyter', 'nbconvert', '--to=slides', nb_name, \
       '--reveal-prefix', '"http://cdn.bootcss.com/reveal.js/3.1.0"',\
      '--template','output_toggle']
# The --template command points to a custom template which opens and closes the input code cells with a click
subprocess.call(cmd)
html_name = '%s.slides.html' % (nb_name.split('.')[0])
shutil.copy(html_name, '/home/bhardin2/public_html/slideshows/')