% Compare extract fields to my (alex's) code

efU = readbin('~dmenemen/llc_2160/regions/latlon/V/V_8640x4320x1.20120123T070000',[8640,4320,1]);

efL = readbin('~dmenemen/llc_2160/regions/latlon/grid/YC_8640x4320',[8640,4320,1]);
nx = 2160;
fnamaU = '~dmenemen/llc_2160/MITgcm/run/U.00092160';
fnamaV = '~dmenemen/llc_2160/MITgcm/run/V.00092160';
[aU, aV] = quikreadpcolor_uv_llc(finu,finv,nx);
lats = quikreadpcolor_llc('~dmenemen/llc_2160/grid/YC.data',nx);

aaU = aV;
aaU(lats>57) = nan;
aaU(lats<-70) = nan;
U = aaU(:,1441:5760);


figure
quikpcolor(aV');
colorbar
caxis([-.3,.3]);
%axis([5000,5500,2000,2500])
print('-dpng','-r600','alexV.png');
figure
quikpcolor(foutv');
colorbar
caxis([-.3,.3]);
%axis([5000,5500,2000,2500])
print('-dpng','-r600','dimitrisV.png');