Enable `divE` field diagnostic with Ohm's law solver (#4178)

* use same finite difference algo as for Yee with hybrid-PIC to calculate divE

* add `divE` and `divB` to field diagnostic in picmi.py

1 files changed, 4 insertions, 2 deletions

diff --git a/Source/FieldSolver/FiniteDifferenceSolver/ComputeDivE.cpp b/Source/FieldSolver/FiniteDifferenceSolver/ComputeDivE.cpp
index 7f4016450..163792a27 100644
--- a/Source/FieldSolver/FiniteDifferenceSolver/ComputeDivE.cpp
+++ b/Source/FieldSolver/FiniteDifferenceSolver/ComputeDivE.cpp
@@ -46,7 +46,8 @@ void FiniteDifferenceSolver::ComputeDivE (
    // Select algorithm (The choice of algorithm is a runtime option,
    // but we compile code for each algorithm, using templates)
 #ifdef WARPX_DIM_RZ
-    if (m_fdtd_algo == ElectromagneticSolverAlgo::Yee){
+    if (m_fdtd_algo == ElectromagneticSolverAlgo::Yee ||
+        m_fdtd_algo == ElectromagneticSolverAlgo::HybridPIC){
 
         ComputeDivECylindrical <CylindricalYeeAlgorithm> ( Efield, divEfield );
 
@@ -55,7 +56,8 @@ void FiniteDifferenceSolver::ComputeDivE (
 
         ComputeDivECartesian <CartesianNodalAlgorithm> ( Efield, divEfield );
 
-    } else if (m_fdtd_algo == ElectromagneticSolverAlgo::Yee) {
+    } else if (m_fdtd_algo == ElectromagneticSolverAlgo::Yee ||
+               m_fdtd_algo == ElectromagneticSolverAlgo::HybridPIC) {
 
         ComputeDivECartesian <CartesianYeeAlgorithm> ( Efield, divEfield );