file1='http://coast-enviro.er.usgs.gov/thredds/dodsC/wake/ucla/run002';
par1=load('ucla_10.mat');

file2='http://coast-enviro.er.usgs.gov/thredds/dodsC/wake/rutgers/run007';
par2=load('rutgers_10.mat');

N=10;
coef=1e4;
jdmat1=cf_time(file1,'zeta');
jdmat2=cf_time(file2,'zeta');
nt=length(jdmat1);
set(gcf,'color','white');
days1=jdmat-jdmat(1);
days2=jdmat-jdmat(1);
set(gcf,'pos', [30         50         986         900]);
h1=axes('pos',[.05 .05 .9 .45]);
h2=axes('pos',[.05 .525 .9 .45]);
nc1=netcdf(file1);
nc2=netcdf(file2);
for i=1:2
%for i=1:nt
  axes(h1);
  u=nc1{'u'}(i,N,:,:);
  v=nc1{'v'}(i,N,:,:);
  xi=coef.*par1.mask.*psi2rho(vorticity(u,v,par1.pm,par1.pn));
%pslice(lon,lat,xi.*mask,[-1 1]);dasp(mean(lat(:)));  % plot in lon/lat
 pslice2(par1.x,par1.y,xi,[-.25 .25]);dasp;
 title(sprintf('Vorticity UCLA: Day= %5.2f',days(i)));
 xlabel('km')
 grid2
 set(gca,'tickdir','out');

 axes(h2);
  u=nc2{'u'}(i,N,:,:);
  v=nc2{'v'}(i,N,:,:);
  xi=coef.*par2.mask.*psi2rho(vorticity(u,v,par2.pm,par2.pn));
%pslice(lon,lat,xi.*mask,[-1 1]);dasp(mean(lat(:)));  % plot in lon/lat
 pslice2(par2.x,par2.y,xi,[-.25 .25]);dasp;
 title(sprintf('Vorticity Rutgers: Day= %5.2f',days(i)));
 xlabel('km')
 grid2
 set(gca,'tickdir','out');
 anim_frame('wake_both',i);
end
close(nc1);
close(nc2);
anim_make