cd ~dmenemen/llc_4320/regions/latlon
U=readbin('U/U_17280x8640x1.20111129T200000',[17280 8640]);
V=readbin('V/V_17280x8640x1.20111129T200000',[17280 8640]);
U=(U(1:end-1,1:end-1)+U(2:end,1:end-1))/2;
V=(V(1:end-1,1:end-1)+V(1:end-1,2:end))/2;
S=sqrt(U.^2+V.^2);
s=0*S;
s(S>0)=1;
s(S>.15)=2;
s(S>.5)=3;
orient landscape
clf
cm=cmap;
cm1=[cm(1,:); cm(51:124,:); cm(201:375,:); cm(463:512,:)];
colormap(cm1)
subplot(211)
mypcolor(S');
set(gca,'xtick',[],'ytick',[])
caxis([0 .6])
h=thincolorbar('horizontal');
set(h,'xtick',[0 .15 .5])
print -dtiff -r2000 ~/speed
