function [Xout] = scrop(X, mt, mb, nl, nr);
%%  remove first mt and last mb rows
%%  remove first nl and last nr columns 
    %%
[m, n, panels] =size(X);
Xout = X(mt+1:m-mb,nl+1:n-nr,:);



