data = readdump_all('dump.lammpstrj');
%dump 1 all custom 10 dump.lammpstrj id type x y vx vy
t = data.timestep;
N = data.Natoms;
nt = length(t);

vx=squeeze(data.atom_data(:,5,:));
vy=squeeze(data.atom_data(:,6,:));
vxsq=mean(vx.^2);
vysq=mean(vy.^2);
T=(vxsq+vysq)/2;
figure(1), plot(t,vxsq,'-r',t,vysq,'-b',t,T,'-k')
xlabel('t'), ylabel('T')
%%
out=load('log.strip');
t=out(:,1);
Ep=out(:,2);
Ek=out(:,3);
T=out(:,4);
P=out(:,5);
figure(2), plot(t,T)
xlabel('t'), ylabel('T')
figure(3), plot(t,Ep+Ek)
xlabel('t'), ylabel('E_{tot}')
