for ibl=1:3

  xy(ibl).psi=zeros(mxbl(ibl)+1,mybl(ibl)+1);
end


xy(1).psi(2:mxbl(1),2:mybl(1)+1)=reshape(psi(1:(mxbl(1)-1)*mybl(1)), ...
                                         mxbl(1)-1,mybl(1));

b1e=(mxbl(1)-1)*mybl(1);					 
xy(2).psi(2:mxbl(2),2:mybl(2))=reshape(psi(b1e+1:b1e+(mybl(2)-1)*(mxbl(2)-1)), ...
                                         mxbl(2)-1,mybl(2)-1);
xy(2).psi(ixstart:ixstart+mxbl(1),1)=xy(1).psi(:,end);

b3b=b1e+(mybl(2)-1)*(mxbl(2)-1)+1;
xy(3).psi(2:mxbl(3),1:mybl(3))=reshape(psi(b3b:end),mxbl(3)-1,mybl(3));

xy(2).psi(ixstart:ixstart+mxbl(3),mybl(2)+1)=xy(3).psi(:,1);