function ws = MultiSURE2(wc,L)
%
%%%%%%It's My Own Function!!!
%
% MultiSURE2 -- Apply Shrinkage to 2-d Wavelet Coefficients
%  Usage 
%    ws = MultiSURE2(wc,L)
%  Inputs 
%    wc    Wavelet Transform of noisy sequence with N(0,1) noise
%    L     low-frequency cutoff for Wavelet Transform
%  Outputs 
%    ws    result of applying SUREThresh to each dyadic block
%
	[n,J] = dyadlength(wc);
    ws = wc;
    for j=(J-1):-1:L
        [t1,t2]=dyad2HH(j);
        [ws(t1,t2),thresh]=SUREThresh2(wc(t1,t2));
        [t1,t2]=dyad2HL(j);
        [ws(t1,t2),thresh]=SUREThresh2(wc(t1,t2));
        [t1,t2]=dyad2LH(j);
        [ws(t1,t2),thresh]=SUREThresh2(wc(t1,t2));
    end