span class="keyword">for</span> ch1=1:Ch;
0036      J1=squeeze(J(:,:,ch1));
0037      C(1:nf,ch1,ch1)=1;
0038      phi(1:nf,ch1,ch1)=0;
0039 <span class="comment">%      if errtype==2;</span>
0040 <span class="comment">%           phierr(1:nf,ch1,ch1)=0;</span>
0041 <span class="comment">%           Cerr(1:2,1:nf,ch1,ch1)=0;</span>
0042 <span class="comment">%      elseif errtype==1</span>
0043 <span class="comment">%            phierr(1:2,1:nf,ch1,ch1)=0;</span>
0044 <span class="comment">%      e