figure('color','white','units','normalized','outerposition',[0 0 1 1],'name','1');
hh = tight_subplot(2,2);

axes(hh(1))
quikpcolor(Umit');
%c = colorbar('southoutside');
caxis([-.2 .5])
title('U Ice Fixed Velocity 1 day avg')
%xlabel(c,'Velocity (m/s)')
set(get(gca,'title'),'fontsize',16);
%a = get(c,'xlabel');
%set(a,'fontsize',16);
set(gca,'xtick',[])
set(gca,'ytick',[])
%axis([0 8640 1500 6500])
%daspect([1 1 1])

axes(hh(2))
quikpcolor(Vmit');
%c = colorbar('southoutside');
caxis([-.2 .5])
title('V Ice Fixed Velocity 1 day avg')
%xlabel(c,'Velocity (m/s)')
set(get(gca,'title'),'fontsize',16);
%a = get(c,'xlabel');
%set(a,'fontsize',16);
set(gca,'xtick',[])
set(gca,'ytick',[])
%axis([0 8640 1500 6500])
%daspect([1 1 1])

load iceGeo;

axes(hh(3))
quikpcolor(uIce');
c = colorbar('southoutside');
caxis([-.2 .5])
title('U No Ice Fix Velocity 1 day avg')
xlabel(c,'Velocity (m/s)')
set(get(gca,'title'),'fontsize',16);
a = get(c,'xlabel');
set(a,'fontsize',16);
set(gca,'xtick',[])
set(gca,'ytick',[])
%axis([0 8640 1500 6500])
%daspect([1 1 1])

axes(hh(4))
quikpcolor(vIce');
c = colorbar('southoutside');
caxis([-.2 .5])
title('V No Ice Fix Velocity 1 day avg')
xlabel(c,'Velocity (m/s)')
set(get(gca,'title'),'fontsize',16);
a = get(c,'xlabel');
set(a,'fontsize',16);
set(gca,'xtick',[])
set(gca,'ytick',[])
%axis([0 8640 1500 6500])
%daspect([1 1 1])

figure('color','white','units','normalized','outerposition',[0 0 1 1],'name','1');
quikpcolor((sqrt(vIce.^2 + uIce.^2) - sqrt(Vmit.^2 + Umit.^2))');
c = colorbar('southoutside');
caxis([0 .5])
title('Diff Mag (IceNotFiced-IceFixed) Velocity 1 day avg')
xlabel(c,'Velocity (m/s)')
set(get(gca,'title'),'fontsize',16);
a = get(c,'xlabel');
set(a,'fontsize',16);
set(gca,'xtick',[])
set(gca,'ytick',[])