function [] = plot_property_vs_winltimes(axn, midwinpoints, time, lEnergy)
global fntsize res xsize ysize lmarg bmarg widthv heightv xstr titlestr fstr pstr
axes(axn);
for i=1:length(midwinpoints);
a(i,1:3) = quantile(lEnergy{i},[0.25, 0.5, 0.75]);
end
h1 = plot(time(midwinpoints),a);
set(h1,'Marker','o', 'MarkerSize', 5, 'LineStyle','-');
datetick(axn,'x','hh:00');
xlabel('Universal Time (hh:mm)', 'fontsize', fntsize, 'fontname', 'helvetica'); % label the axes
ylabel('Lightning Energy (J)','fontsize', fntsize, 'fontname', 'helvetica');
title('Quantiles of Matched Lightning Energies over Time','fontsize', fntsize, 'fontname', 'helvetica');
end