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;
clf
set(gcf,'color','white');
set(gcf,'pos', [30         50         986         900]);
h1=axes('pos',[.07 .05 .85 .45]);
h2=axes('pos',[.07 .525 .85 .45]);
nc1=netcdf(file1);
days1=nc1{'ocean_time'}(:)/3600/24;
nc2=netcdf(file2);
days2=nc2{'ocean_time'}(:)/3600/24;
nt=length(days1);
%for i=1:2
for i=1:nt
  axes(h1);
  u1=nc1{'u'}(i,N,:,:);
  v1=nc1{'v'}(i,N,:,:);
  xi1=coef.*par1.mask.*psi2rho(vorticity(u1,v1,par1.pm,par1.pn));
%pslice(lon,lat,xi.*mask,[-1 1]);dasp(mean(lat(:)));  % plot in lon/lat
 pslice2(par1.x,par1.y,xi1,[-.25 .25]);dasp;
% t1=text(-30,30,sprintf('Vorticity UCLA: Day= %5.2f',days1(i)));
 title(sprintf('Vorticity (1e4) UCLA: Day= %5.2f',days1(i)));
 xlabel('km')
 ylabel('km')
 grid2
 set(gca,'tickdir','out');

 axes(h2);
  u2=nc2{'u'}(i,N,:,:);
  v2=nc2{'v'}(i,N,:,:);
  xi2=coef.*par2.mask.*psi2rho(vorticity(u2,v2,par2.pm,par2.pn));
%pslice(lon,lat,xi.*mask,[-1 1]);dasp(mean(lat(:)));  % plot in lon/lat
 pslice2(par2.x,par2.y,xi2,[-.25 .25]);dasp;
% t2=text(-30,30,sprintf('Vorticity Rutgers: Day= %5.2f',days2(i)));
 t2=title(sprintf('Vorticity (1e4) Rutgers: Day= %5.2f',days2(i)));
 ylabel('km')
 grid2
 set(gca,'tickdir','out');
 pclegend([-.25 .25].',[.96 .2 .03 .6],gcf);
 anim_frame('wake_both',i);
end
close(nc1);
close(nc2);
anim_make