load si_sim.oum
theta=si_sim(:,1);	%crankangle
p    =si_sim(:,2);	%pressure (bar)
tu   =si_sim(:,3);	%unburned gas temp (k)
tb   =si_sim(:,4);	%bunred gas temp (k)
m_int=si_sim(:,5);	%cumulative mass inducted (g)
m_exh=si_sim(:,6);	%cumulative mass exhausted (g)
viv  =si_sim(:,7);	%velocity at intake valve (m/s)
vexh =si_sim(:,8);	%velocity at exhaust valve (m/s)
xb   =si_sim(:,9);	%mass fraction burned
x_net=si_sim(:,10);	%net ht rel / LHV*fuelmass
qdot_h=si_sim(:,11);	%heat transfer (J/deg)
w_pist=si_sim(:,12);	%work transfer to piston (Kj)
clear si_sim