url='http://www.smast.umassd.edu:8080/thredds/dodsC/fvcom/hindcasts/30yr_gom3';
ax=[ -71.2591  -68.3224   41.3899   43.9527];
movie_name='april87_fix';
cax=[30 33]
tic
nc=ncgeodataset(url);
%%
lon=nc{'lon'}(:);
lat=nc{'lat'}(:);
nv=nc{'nv'}(:).';
jd=nj_time(nc,'salinity');
ii=date_index(jd,[1987 4 1 0 0 0],[1987 5 10 1 0 0]);
jj=0;
isub=3
vidobj=VideoWriter(movie_name); %creates AVI file, test.avi 
open(vidobj);
fig=figure(1);
for i=1:isub:length(ii);
  jj=jj+1;
  s=nc{'salinity'}(ii(i),1,:);
  trisurf(nv,lon,lat,zeros(size(s)),s);view(2);shading interp;colorbar;dasp(40);caxis(cax)
  title(['FVCOM GOM3 surface salinity: ' datestr(jd(ii(i)))])
  set(gca,'tickdir','out');
  axis(ax);
  F=getframe(fig);
  writeVideo(vidobj,F);
end
close(vidobj);
toc