Add code to solve for q^(1) and check J_R conservation.