IMPORTANT: Corrected bug in interp_boundary.m

[arango]

Corrected a bug in Matlab script interp_boundary.m. In line 201 we need to have:

        if (~isempty(ind)),
          if (RemoveNaN),
            B.(field)(ind) = FN(I.XR.(edge)(ind), I.YR.(edge)(ind));  % <===
            Rmin = min(Rmin, min(B.(field)(ind)));
            Rmax = max(Rmax, max(B.(field)(ind)));

            ind = find(isnan(B.(field)));
            if (~isempty(ind)),
              Ncount = length(ind);
            end       
          else
            Ncount = length(ind);
          end
        end

instead of

            B.(field)(ind) = FN(I.XR.(edge)(ind), I.YR.(field)(ind));

We need to have edge instead of field in I.YR.(?) structure. This bug is only relevant when removing interpolated variable NaNs values with a nearest neighbor interpolant. This is the second interpolation to remove outliers.

This is an important correction and users need to update this file if using this function to prepare boundary conditions.