clear all
clc
omega='0.13';
level=1;
filename=strcat(omega,'rads\Vel_data_level',num2str(level),'.nc');
[x, info] = readNetCDFVariable(filename,'x');
[y, info] = readNetCDFVariable(filename,'y');
[t, info] = readNetCDFVariable(filename,'t');
[u, info] = readNetCDFVariable(filename,'u');
[v, info] = readNetCDFVariable(filename,'v');

for snapshot=1:1:100
[x2d,y2d]=meshgrid(x,y);
u_snapshot=u(:,:,snapshot);
v_snapshot=v(:,:,snapshot);
[vort,cav]=curl(x2d,y2d,u_snapshot,v_snapshot);
figure(1)
pcolor(x2d,y2d,vort)
shading interp
colormap(jet)
caxis([-0.6 0.6])
pause(0.05)
end