function []=check_HxHy(H2,Hx2,Hy2)
% Check that Hx2 and Hy2 don't violate basic rules
% This should *not* ever happen
%return; % spead up process
[Hxcc,Hycc]=gen_hxhy(H2);
ic=find(Hx2<Hxcc-0.001);
if prod(size(ic)) ~= 0
%[Hx2(ic)/1e3 Hxcc(ic)/1e3 Hx2(ic)-Hxcc(ic)]
 error('The X cross-sectional area turned out inconsistently')
end
jc=find(Hy2<Hycc-0.001);
if prod(size(jc)) ~= 0
%[Hy2(jc)/1e3 Hycc(jc)/1e3 Hy2(jc)-Hycc(jc)]
 error('The Y cross-sectional area turned out inconsistently')
end