sim_geometry.asv

function [torque, mass, j_rotor, r_phase] = sim_geometry(g, theta, id, iq, jd, jq)
    
    init_geometry_2(g, theta, id, iq, jd, jq);
    torque = calc_torque(g);
    i_norm = norm([id, iq]);
    j_norm = norm([jd, jq]);
    [mass, m_stator, m_rotor, j_rotor, r_phase, p_ir] = calc_phys_props(g, i_norm, j_norm);
    r_phase
    p_ir
end