%
% output_test.m
% requires having matlab function grdread.mexmaci loaded in the toolbox
%
clear

[lon,lat,topo]=grdread('topo_masked.grd');
x=grdread('x.grd');
y=grdread('y.grd');
z=grdread('z.grd');


figure(1),clf
subplot(2,2,1),imagesc(lon,lat,topo,[-200,200]),axis xy,colorbar,title('input topo')
subplot(2,2,2),imagesc(lon,lat,x),axis xy,colorbar,title('output x')
subplot(2,2,3),imagesc(lon,lat,y),axis xy,colorbar,title('output y')
subplot(2,2,4),imagesc(lon,lat,z),axis xy,colorbar,title('output z')
  axis([-117,-115,32,34])
