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