Daniel@0: function bigT = divide_by_table(bigT, bigdom, bigsz, smallT, smalldom, smallsz)
Daniel@0: % DIVIDE_BY_TABLE 
Daniel@0: % bigT = divide_by_table(bigT, bigdom, bigsz, smallT, smalldom, smallsz)
Daniel@0: %
Daniel@0: 
Daniel@0: 
Daniel@0: Ts = extend_domain_table(smallT, smalldom, smallsz, bigdom, bigsz);
Daniel@0: % Replace 0s by 1s before dividing. This is valid, Ts(i)=0 iff Tbig(i)=0.
Daniel@0: Ts = Ts + (Ts==0);
Daniel@0: %Tbig.T(:) = Tbig.T(:) ./ Ts(:);
Daniel@0: bigT(:) = bigT(:) ./ Ts(:);
Daniel@0: