condition for 0 uth

RevathiJambunathan · Mar 28, 2024 · 88a2f12 · 88a2f12
1 parent b5ad2bb
commit 88a2f12
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/Source/Particles/ParticleBoundaries_K.H b/Source/Particles/ParticleBoundaries_K.H
@@ -85,7 +85,7 @@ namespace ApplyParticleBoundaries {
     {
         u_tang1 = (uth > 0._rt) ? PhysConst::c * amrex::RandomNormal(0._rt, uth, engine) : 0._rt;
         u_tang2 = (uth > 0._rt) ? PhysConst::c * amrex::RandomNormal(0._rt, uth, engine) : 0._rt;
-        u_norm = std::copysign(1._prt, -u_norm) * PhysConst::c * generateGaussianFluxDist(0._rt, uth, engine);
+        u_norm = (uth > 0._rt) ? std::copysign(1._prt, -u_norm) * PhysConst::c * generateGaussianFluxDist(0._rt, uth, engine) : 0._rt;
     }