function [hs1,hs2,hs3]=panel_day(timestamp,Vh,theta,beta,heights)

timestartdisp=timestamp(1);
start_vec=datevec(timestartdisp);
start_str=datestr(start_vec,'mmddyyyy');
daynum=datenum(start_vec(1), start_vec(2),start_vec(3),0,0,0);
date_time=datestr(timestamp(1),'yyyy-mm-dd');

hs1=figure('Position',[50 50 1000 300]);
p2=pcolor(timestamp-daynum,heights,Vh');
set(p2,'edgecolor','none')
colormap(gca,jet(128))
title(['Horizontal Wind Speed (m/s) for ',datestr(timestamp(1),'mm/dd/yy')],'FontSize',14,'FontWeight','bold')
caxis([0 25])
datetick('x',15)
ylabel('Heights AGL (m)','FontSize',14)
xlabel('Time of Day (UTC)','FontSize',14)
cbt=colorbar('YTick',[0,5,10,15,20,25]);

ylabel('Height AGL (m)','FontSize',14)
xlabel('Time (UTC) ','FontSize',14)
saveas(hs1,['Vh-',start_str,'.png'])
saveas(hs1,['Vh-',start_str,'.fig'])

%ht1=title(['Wind Profile for  ',date_time],'FontSize',16)

hs2=figure('Position',[50 50 1000 300]);
p2=pcolor(timestamp-daynum,heights,theta');
set(p2,'edgecolor','none')
title(['Wind Direction for ',datestr(timestamp(1),'mm/dd/yy')],'FontSize',14,'FontWeight','bold')
caxis([0 360])
datetick('x',15)
colormap(hsv(50))
ylabel('Heights AGL (m)','FontSize',14)
xlabel('Time of Day (UTC)','FontSize',14)
cbt=colorbar('YTick',[0,90,180,270,360]);
set(cbt,'YTickLabel',{'N', 'E', 'S', 'W', 'N'})
ylabel('Height AGL (m)','FontSize',14)
xlabel('Time (UTC) ','FontSize',14)


hs3=figure('Position',[50 50 1000 300]);
p3=pcolor(timestamp-daynum,heights,10*log10(beta'/1e-3));
colormap(gca,jet(256))
set(p3,'edgecolor','none')
title(['Backscatter Coefficient (dBm) for ',datestr(timestamp(1),'mm/dd/yy')],'FontSize',14,'FontWeight','bold')
caxis([-30 -5])
datetick('x',15)
colormap(jet(128))
colorbar
ylabel('Height AGL (m)','FontSize',14)
xlabel('Time (UTC) ','FontSize',14)

end