2 calls to
CALC_P_RHO
SOLVE_EM
SOLVE_EM