[7]  1  %'griddata_uvmat': function griddata_uvmat(vec2_X,vec2_Y,vec2_U,vec_X,vec_Y,'linear')


 2  %adapt the input of the matlab function griddata to the appropriate version of Matlab


[399]  3  function ZI = griddata_uvmat(X,Y,Z,XI,YI)


 4  % if ~exist('rho','var') isequal(rho,0)


[427]  5  ZI=griddata(X,Y,Z,XI,YI,'linear');


 6  % txt=ver('MATLAB');


 7  % Release=txt.Release;


 8  % relnumb=str2num(Release(3:4));


 9  % if relnumb >= 20


 10  % ZI=griddata(double(X),double(Y),double(Z),double(XI),double(YI),'linear',{'QJ'});


 11  % elseif relnumb >=14


 12  % ZI=griddata(X,Y,Z,XI,YI,'linear',{'QJ'});


 13  % else


 14  % ZI=griddata(X,Y,Z,XI,YI,'linear');


 15  % end


[399]  16  % else %smooth with thin plate spline


 17  % [ZI,Z_diff]=patch_uvmat(X,Y,Z,XI,YI,rho);


 18  % diff_norm=mean(Z_diff.*Z_diff)


 19  % ind_good=find(abs(Z_diff)<5*diff_norm);


 20  % nb_remove=numel(Z_diff)numel(ind_good)


 21  % if nb_remove>0


 22  % X=X(ind_good);


 23  % Y=Y(ind_good);


 24  % Z=Z(ind_good);


 25  % [ZI,Z_diff]=patch_uvmat(X,Y,Z,XI,YI,rho);


 26  % diff_norm_new=mean(Z_diff.*Z_diff)


 27  % end


 28  % end 
