elim RHSfunc, replace with j_solver_alpha, p and g now fields