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