flx_fill_in.m

来自「Petroleum Reservoir Simulation」· M 代码 · 共 31 行

function flx=flx_fill_in(flxy1,flxy2,flxy3,nr)
for ii=1:nr
   if ii==1
      flx(ii,1)=flxy2(ii);
      flx(ii,2)=flxy3(ii);
   end
   
   if ii~=1&ii~=nr
      for jj=(ii-1):(ii+1)
         if ii>jj
            flx(ii,jj)=flxy1(ii-1);
         end
         if ii==jj
            flx(ii,jj)=flxy2(ii);
         end
         if ii<jj
            flx(ii,jj)=flxy3(ii-1);
         end
      end
   end
   
   if ii==nr
      flx(ii,nr-1)=flxy1(nr-1);
      flx(ii,nr)=flxy2(nr);
   end
   
end

%disp('next step is fl_fill_in:fl=fl_fill_in(fwp,fws,nr)');