% lookkk at   ice differences
close all
ts=140160
nx =  2160;


dyy=ts2dte(ts,45,2011,1,17,30);
landV = quikreadpcolor_llc('~dmenemen/llc_2160/grid/hFacS.data',nx);


fnam = ['/nobackupp8/awinetee/test/Eta_seaIceRemoved_30hLPF_hamming_daily/' 'Eta_seaIceRemoved_' int2str((nx)) 'x' int2str((nx*13)) '.' dyy];
etaIce = quikreadpcolor_llc(fnam,nx);

fnam = ['/nobackupp8/awinetee/test/Eta_seaIceRemoved_30hLPF_hamming_daily/' 'Eta_seaIceACTUALLYRemoved_' int2str((nx)) 'x' int2str((nx*13)) '.' dyy];
etaNoIce = quikreadpcolor_llc(fnam,nx);


etaIce(landV==0) = nan;
etaNoIce(landV==0) = nan;
close all;

figure('color','white');
quikpcolor(etaIce');
set(gca,'xticklabel','')
set(gca,'yticklabel','')
grid on
c = colorbar
ylabel(c,'Sea Surface Height (m)');
colormap([160/255 160/255 160/255 ; brewermap(256,'RdYlBu')])
caxis([-3.5 2])

print('-dpng','-r400','~awinetee/Desktop/thesis/seaSurfaceHeight_withIce_gray.png');
eta(landV==0) = nan;
h = figure('color','white');

quikpcolor(eta');
set(gca,'xticklabel','')
set(gca,'yticklabel','')
c = colorbar
ylabel(c,'Sea Surface Height (m)');
grid on
colormap([160/255 160/255 160/255 ; brewermap(256,'RdYlBu')])
caxis([-3.5 2])
print('-dpng','-r400','~awinetee/Desktop/thesis/seaSurfaceHeight_NoIce_filtered_gray.png');
