2 calls to
THETA_AND_THETAM_LBC_ONLY
SOLVE_EM
SOLVE_EM