% draw W map
% Create W from C for +/- r vs z
% Plot in Figure 2

figure(2);clf
set(figure(2), 'position', [10 100 1000 700])

maxC             = max(max(C));
temp             = C(2:Nz-1, 2:Nr-1)*exp(-mua*c*time)/Co;
Nw = 2*(Nr-2);
W(:,Nw/2+1:Nw)   = temp;
W(:,Nw/2:-1:1)   = temp;


imagesc(rw,zw,log10(W),[-12 0])
%imagesc(rw,zw,W,[0.999 1.001])
%imagesc(rw,zw,W, [0 1])


colorbar
colormap(c2)
axis('square')
set(gca,'fontsize',24)
set(colorbar,'fontsize',24)
text(rmax,-zmax*0.05,'log_1_0(C/U_0)','fontsize',24)
xlabel('r [cm]')
ylabel('z [cm]')
axis([-rmax rmax 0 zmax])

if time<1000e-12
	title(sprintf('time = %5.0f [ps]',time*1e12),'fontsize',36)
else
	title(sprintf('time = %5.3f [ns]',time*1e9),'fontsize',36)
end

hold on
