function zopsout=zops(dM,n,mc)
%        zopsout=zops(dM,n,mc);
%  maps  R^{n x n} --->  R^{n-1}xR^{mc}
%           matrix space n by n to vectors
global indsEc x Z
Zs=(Z*dM+dM'*Z)/2;    % adjoint multiplication by Z
%% adjoint of partial diagonal operator
zopsout=diag(Zs);
zopsout=zopsout(2:n)-zopsout(1);   % vector length n-1
%%
zopsout=[zopsout
         usvecEc(Zs,n,mc) ];  % final vector length n-1 + mc
