vis3d/meshmovie.m

   1 function meshmovie(q,name)
   2 % input
   3 % q     structure with fields from wrfout
   4 % name  name of the field to show
   5
   6 v = q.(name);
   7 minv = min(v(:));
   8 maxv = max(v(:));
   9 for i=1:size(q.times,2)
  10     mesh(q.xlat(:,:,i),q.xlong(:,:,i),v(:,:,i))
  11     t = char(q.times(:,i)')
  12     zlim([minv,maxv])
  13     title(t,interpreter='none')
  14     drawnow
  15 end