function sh = smooth(h)
%smooth applies a 1-2-1 smoother to the input field

sh=(circshift(h,1)+circshift(h,-1)+2*h)/4;
end

