vis3d/interp_12.m

   1 function r=interp_12(a,i1,i2)
   2 % horizontal interpolation: in the first two coordinates, to given i1 i2
   3     if isempty(a),
   4         r=[];
   5     else
   6         nv=size(a,3); % number of vertical layers
   7         r=zeros(nv,1);
   8         for i=1:nv
   9             b=squeeze(double(a(:,:,i)));
  10             r(i)=interp2(b,i1,i2);
  11         end
  12     end
  13 end
  14