function [bb_diff_sct, drysmc_diff_sct, f11_diff_sct, maxsmc_diff_sct,...
    refsmc_diff_sct, satpsi_diff_sct, satdk_diff_sct, satdw_diff_sct, ...
    wltsmc_diff_sct, qtz_diff_sct]...
    = func_sct_only(bb_diff, drysmc_diff, f11_diff, maxsmc_diff,...
    refsmc_diff, satpsi_diff, satdk_diff, satdw_diff, wltsmc_diff, qtz_diff,sct_mat, sct)

    bb_diff_sct = bb_diff;
    bb_diff_sct(sct_mat~=sct) = nan;

    drysmc_diff_sct = drysmc_diff;
    drysmc_diff_sct(sct_mat~=sct) = nan;
    
    f11_diff_sct = f11_diff;
    f11_diff_sct(sct_mat~=sct) = nan;
    
    maxsmc_diff_sct = maxsmc_diff;
    maxsmc_diff_sct(sct_mat~=sct) = nan;
    
    refsmc_diff_sct = refsmc_diff;
    refsmc_diff_sct(sct_mat~=sct) = nan;
    
    satpsi_diff_sct = satpsi_diff;
    satpsi_diff_sct(sct_mat~=sct) = nan;
    
    satdk_diff_sct = satdk_diff;
    satdk_diff_sct(sct_mat~=sct) = nan;
    
    satdw_diff_sct = satdw_diff;
    satdw_diff_sct(sct_mat~=sct) = nan;
    
    wltsmc_diff_sct = wltsmc_diff;
    wltsmc_diff_sct(sct_mat~=sct) = nan;
    
    qtz_diff_sct = qtz_diff;
    qtz_diff_sct(sct_mat~=sct) = nan;

end

