Wavelet 2D decomposition - how can I input a total varitation constraint to be used in this algorithm?

I asked a very general question about the software last time and i have been using it successfully since thanks to your help :D Many thanks again! However, I have now hit a specific problem. I am trying to use the Wavelet 2D software with a total variation constraint for 3D data. I note that the documentation says

"The first consists of taking the wavelet expansion of the signal and keeping the largest absolute value coefficients. In this case, you can set a global threshold, a compression performance, or a relative square norm recovery performance. Thus, only a single parameter needs to be selected."

I'm assuming therefore that I need to set the relative square norm recovery to be whatever I want my total variation constraint to be. However, it is not clear how I input this constraint into matlab

[xd,cxd,lxd,perf0,perfl2] = wdencmp(opt,c,l,w,n,thr,sorh,keepapp); I can see that it is outputted by perf0 but can this be manipulated into a constraint to be used in the algorithm for compression?

Many thanks in advance guys :)