AddPlasma: Avoid a calculation resulting in NaN when num_ppc is 0 (#1787)

* avoid a NaN when num_ppc is 0

* fix style

1 files changed, 6 insertions, 3 deletions

diff --git a/Source/Particles/PhysicalParticleContainer.cpp b/Source/Particles/PhysicalParticleContainer.cpp
index f0320d00f..412ed53d4 100644
--- a/Source/Particles/PhysicalParticleContainer.cpp
+++ b/Source/Particles/PhysicalParticleContainer.cpp
@@ -541,12 +541,15 @@ PhysicalParticleContainer::AddPlasma (int lev, RealBox part_realbox)
     const auto dx = geom.CellSizeArray();
     const auto problo = geom.ProbLoArray();
 
-    Real scale_fac;
+    Real scale_fac = 0.0_rt;
+
+    if(num_ppc != 0){
 #if AMREX_SPACEDIM==3
-    scale_fac = dx[0]*dx[1]*dx[2]/num_ppc;
+        scale_fac = dx[0]*dx[1]*dx[2]/num_ppc;
 #elif AMREX_SPACEDIM==2
-    scale_fac = dx[0]*dx[1]/num_ppc;
+        scale_fac = dx[0]*dx[1]/num_ppc;
 #endif
+    }
 
     defineAllParticleTiles();