function ROT = Rotationnel(II)

% ROT = Rotationnel(II)
%
%Calcul du rotationnel d'un champ de vecteurs. 

[N1,N2,N3] = size(II);
T1  = [2:N1 1]; T2 = [2:N2 1]; 
ROT = (II(:,T2,2)-II(:,:,2))- (II(T1,:,1)-II(:,:,1)); 
end
