function h = fill_plot(time,data_mat,color)
% function h = fill_plot(time, data_mat, color)
%
%

x = [time fliplr(time)];

sem = nanstd(data_mat)./sqrt(length(data_mat(:,1)));

y = [(nanmean(data_mat)-sem) fliplr(nanmean(data_mat)+sem)];
h = fill(x,y,color,'EdgeColor','none');


