Merge pull request #29 from burlen/sensei_missing_field

expose rho and F to sensei

Author: RemiLehe

Source/WarpXIO.cpp

@@ -474,6 +474,8 @@ WarpX::UpdateInSitu () const
         + static_cast<int>(plot_proc_number)
         + static_cast<int>(plot_divb)
         + static_cast<int>(plot_dive)
+        + static_cast<int>(plot_rho)
+        + static_cast<int>(plot_F)
         + static_cast<int>(plot_finepatch)*6
         + static_cast<int>(plot_crsepatch)*6
         + static_cast<int>(costs[0] != nullptr);
@@ -622,6 +624,26 @@ WarpX::UpdateInSitu () const
             dcomp += 1;
         }
 
+        if (plot_rho)
+        {
+            amrex::average_node_to_cellcenter(mf[lev], dcomp, *rho_fp[lev], 0, 1);
+            if (lev == 0)
+            {
+                varnames.push_back("rho");
+            }
+            dcomp += 1;
+        }
+
+        if (plot_F)
+        {
+            amrex::average_node_to_cellcenter(mf[lev], dcomp, *F_fp[lev], 0, 1);
+            if (lev == 0)
+            {
+                varnames.push_back("F");
+            }
+            dcomp += 1;
+        }
+
         if (plot_finepatch)
         {
             PackPlotDataPtrs(srcmf, Efield_fp[lev]);