PlotEspic2dpartdata.m
b18e08af36dcgeom_testmumps

PlotEspic2dpartdata.m
View Options

			function PlotEspic2dpartdata(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(3,2,1,'Parent',fig);
			plot(ax1,M.Z(:,fieldstep),M.R(:,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')
			grid(ax1, 'on')


			ax2=subplot(3,2,2,'Parent',fig);
			plot(ax2,M.VZ(:,fieldstep),M.VR(:,fieldstep),'.');
			xlabel(ax2,'VZ [m/s]')
			ylabel(ax2,'VR [m/s]')
			axis(ax2, 'equal')
			grid(ax2, 'on')


			ax3=subplot(3,2,3,'Parent',fig);
			plot(ax3,M.Z(:,fieldstep),M.VZ(:,fieldstep),'.');
			xlim(ax3,[M.zgrid(1) M.zgrid(end)])
			xlabel(ax3,'Z [m]')
			ylabel(ax3,'VZ [m/s]')
			grid(ax3, 'on')

			ax4=subplot(3,2,4,'Parent',fig);
			plot(ax4,M.VZ(:,fieldstep),M.R(:,fieldstep),'.')
			ylabel(ax4,'R [m]')
			xlabel(ax4,'VZ [m/s]')
			ylim(ax4,[M.rgrid(1) M.rgrid(end)])
			grid(ax4, 'on')

			ax5=subplot(3,2,5,'Parent',fig);
			plot(ax5,M.Z(:,fieldstep),M.VR(:,fieldstep),'.');
			xlim(ax5,[M.zgrid(1) M.zgrid(end)])
			xlabel(ax5,'Z [m]')
			ylabel(ax5,'VR [m/s]')
			grid(ax5, 'on')

			ax6=subplot(3,2,6,'Parent',fig);
			plot(ax6,M.VR(:,fieldstep),M.R(:,fieldstep),'.');
			ylim(ax6,[M.rgrid(1) M.rgrid(end)])
			ylabel(ax6,'R [m]')
			xlabel(ax6,'VR [m/s]')
			grid(ax6, 'on')
			end