2 calls to
CALC_P_RHO_PHI
SOLVE_EM
SOLVE_EM