if PlotType == 1
  % Assume presssure is being plotted.
  rybcolormap;
  caxis([0.2 0.6]);
  daspect([1 1 1]);
  setsurfalpha(0.5);
  if (Frame == 0)
    camlight;
  end;
end

if PlotType==4
   axis([0 3.5 0 1]);
   % scatter plot.  Also plot rad1d solution for comparison
   [rad1d,tref] = readamrdata(1,Frame,'./rad1d/');
   if (abs(tref  - t) > 1e-8)
     error('Times for reference solution and 3d solution not compatible');
   end
   hold on;
   [qref,xref,p] = plotframe1ez(rad1d,mq,'r-',UserVariableFile);
   set(p,'LineWidth',2.5);
   h = getlegendinfo;
   legend([h,p],'3d solution','1d solution');

   hold off

 end;

 clear afterframe;
