diff --git a/Source/Microphysics/Kessler/ERF_Kessler.cpp b/Source/Microphysics/Kessler/ERF_Kessler.cpp
index 5268f369a..958eb6440 100644
--- a/Source/Microphysics/Kessler/ERF_Kessler.cpp
+++ b/Source/Microphysics/Kessler/ERF_Kessler.cpp
@@ -177,7 +177,7 @@ void Kessler::AdvanceKessler (const SolverChoice &solverChoice)
                                        {
                                            return { ma_fz_arr[box_no](i,j,k) };
                                        });
-        wt_max = get<0>(max) + 1.0e-16;
+        wt_max = get<0>(max) + std::numeric_limits<Real>::epsilon();
         n_substep = int( std::ceil(wt_max * coef / CFL_MAX) );
         AMREX_ALWAYS_ASSERT(n_substep >= 1);
         coef /= Real(n_substep);
diff --git a/Source/Microphysics/SAM/ERF_IceFall.cpp b/Source/Microphysics/SAM/ERF_IceFall.cpp
index 31f80a195..99bca7f52 100644
--- a/Source/Microphysics/SAM/ERF_IceFall.cpp
+++ b/Source/Microphysics/SAM/ERF_IceFall.cpp
@@ -80,7 +80,7 @@ void SAM::IceFall (const SolverChoice& sc) {
                          {
                              return { ma_fz_arr[box_no](i,j,k) };
                          });
-    wt_max = get<0>(max) + 1.0e-16;
+    wt_max = get<0>(max) + std::numeric_limits<Real>::epsilon();
     n_substep = int( std::ceil(wt_max * coef / CFL_MAX) );
     AMREX_ALWAYS_ASSERT(n_substep >= 1);
     coef /= Real(n_substep);
diff --git a/Source/Microphysics/SAM/ERF_PrecipFall.cpp b/Source/Microphysics/SAM/ERF_PrecipFall.cpp
index 0590d8990..07daef1d8 100644
--- a/Source/Microphysics/SAM/ERF_PrecipFall.cpp
+++ b/Source/Microphysics/SAM/ERF_PrecipFall.cpp
@@ -124,7 +124,7 @@ SAM::PrecipFall (const SolverChoice& sc)
                          {
                              return { ma_fz_arr[box_no](i,j,k) };
                          });
-    wt_max = get<0>(max) + 1.0e-16;
+    wt_max = get<0>(max) + std::numeric_limits<Real>::epsilon();
     n_substep = int( std::ceil(wt_max * coef / CFL_MAX) );
     AMREX_ALWAYS_ASSERT(n_substep >= 1);
     coef /= Real(n_substep);