
%yrbcolormap
colormap(1-.3*(0:.1:1)' * [0 1 1])
if Manifold==0
  % plot in computational plane
  daspect([1 1 1])
  axis([-3 1 -1 1])
  % If you want to see the full domain with both copies of sphere:
  % axis([-3 1 -3 1])
  view(2)
  return
  end

caxis([0 2])
daspect([1 1 1])
grid off
axis off
title(' ')
text(.8,0,2.2,['t = ' num2str(t)],'FontSize',15)

view([100 38])
return

% plot some grid lines:
hold on
[xeq,yeq,zeq] = mapc2m(zeros(1,100),linspace(-1,1,100));
plot3(xeq,yeq,zeq,'k')
[xeq,yeq,zeq] = mapc2m(-ones(1,100),linspace(-1,1,100));
plot3(xeq,yeq,zeq,'k')
[xeq,yeq,zeq] = mapc2m(ones(1,100),linspace(-1,1,100));
plot3(xeq,yeq,zeq,'k')
[xeq,yeq,zeq] = mapc2m(-2*ones(1,100),linspace(-1,1,100));
plot3(xeq,yeq,zeq,'k')
[xeq,yeq,zeq] = mapc2m(.5*ones(1,100),linspace(-1,1,100));
plot3(xeq,yeq,zeq,'k')
[xeq,yeq,zeq] = mapc2m(-.5*ones(1,100),linspace(-1,1,100));
plot3(xeq,yeq,zeq,'k')
[xeq,yeq,zeq] = mapc2m(-1.5*ones(1,100),linspace(-1,1,100));
plot3(xeq,yeq,zeq,'k')
[xeq,yeq,zeq] = mapc2m(-2.5*ones(1,100),linspace(-1,1,100));
plot3(xeq,yeq,zeq,'k')
[xeq,yeq,zeq] = mapc2m(-1.5*ones(1,100),linspace(-1,1,100));
plot3(xeq,yeq,zeq,'k')

[xeq,yeq,zeq] = mapc2m(linspace(-1,1,100),.5*ones(1,100));
plot3(xeq,yeq,zeq,'k')
[xeq,yeq,zeq] = mapc2m(linspace(-1,1,100),zeros(1,100));
plot3(xeq,yeq,zeq,'k')
[xeq,yeq,zeq] = mapc2m(linspace(-1,1,100),-.5*ones(1,100));
plot3(xeq,yeq,zeq,'k')
[xeq,yeq,zeq] = mapc2m(linspace(-1,1,100)-2,.5*ones(1,100));
plot3(xeq,yeq,zeq,'k')
[xeq,yeq,zeq] = mapc2m(linspace(-1,1,100)-2,-.5*ones(1,100));
plot3(xeq,yeq,zeq,'k')
[xeq,yeq,zeq] = mapc2m(linspace(-1,1,100)-2,zeros(1,100));
plot3(xeq,yeq,zeq,'k')

% plot the axis of rotation:
%plot3(xax,yax,zax,'b','LineWidth',2)

%showgridlines
hold off

