function Zopout=zop(dd,dx,n,mc)
%        Zopout=zop(dd,dx,n,mc)
%  maps (dd,dx) in R^{n-1}xR^{mc} ---> R^{n,n}
%         vectors   to  matrix space n by n
global Z 
Zopout=usMatEc(dx,n,mc);
Zopout(1:n+1:end)=[-sum(dd);dd];
Zopout=Z*Zopout;
%Zopout=Z*(  spdiags( [-sum(dd);dd],0,n,n)+usMatEc(dx,n,mc)  ); % old way
