PlotEspic2dgriddata.m
a2847b9464fageom_testpardiso

PlotEspic2dgriddata.m
View Options

			function PlotEspic2dgriddata(fig,M,fieldstep)
			%PlotEspic2dgriddata Plot the 2d data of espic2d at time step fieldstep
			sgtitle(fig,sprintf('step=%d t=%0.5e s',(fieldstep-1)*M.it1,M.t2d(fieldstep)))


			ax1=subplot(2,2,1,'Parent',fig);
			pcolor(ax1,M.zgrid,M.rgrid,M.N(:,:,fieldstep));
			xlim(ax1,[M.zgrid(1) M.zgrid(end)])
			ylim(ax1,[M.rgrid(1) M.rgrid(end)])
			xlabel(ax1,'Z [m]')
			ylabel(ax1,'R [m]')
			title(ax1,'Position')
			c = colorbar(ax1);
			c.Label.String= 'N';
			view(ax1,2)
			%set(ax1,'colorscale','log')


			ax2=subplot(2,2,2,'Parent',fig);
			contourf(ax2,M.zgrid,M.rgrid,M.pot(:,:,fieldstep)');
			%plot(ax2,M.zgrid,data.pot(:,5))
			xlabel(ax2,'Z [m]')
			ylabel(ax2,'R [m]')
			colormap(ax2,'jet')
			c = colorbar(ax2);
			c.Label.String= '\Phi [V]';
			title(ax2,'Es potential')
			grid(ax2, 'on')


			ax3=subplot(2,2,3,'Parent',fig);
			contourf(ax3,M.zgrid,M.rgrid,M.Er(:,:,fieldstep)')
			xlabel(ax3,'Z [m]')
			ylabel(ax3,'R [m]')
			c = colorbar(ax3);
			c.Label.String= 'E_r [V/m]';
			title(ax3,'Radial Electric field')
			grid(ax3, 'on')

			ax4=subplot(2,2,4,'Parent',fig);
			contourf(ax4,M.zgrid,M.rgrid,M.Ez(:,:,fieldstep)')
			xlabel(ax4,'Z [m]')
			ylabel(ax4,'R [m]')
			c = colorbar(ax4);
			c.Label.String= 'E_z [V/m]';
			title(ax4,'Axial Electric field')
			grid(ax4, 'on')
			end