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