diff --git a/src/PointwiseFunctions/AnalyticSolutions/RelativisticEuler/FishboneMoncriefDisk.cpp b/src/PointwiseFunctions/AnalyticSolutions/RelativisticEuler/FishboneMoncriefDisk.cpp
index 41361ed7b53e7..4f997412f5f33 100644
--- a/src/PointwiseFunctions/AnalyticSolutions/RelativisticEuler/FishboneMoncriefDisk.cpp
+++ b/src/PointwiseFunctions/AnalyticSolutions/RelativisticEuler/FishboneMoncriefDisk.cpp
@@ -154,13 +154,20 @@ FishboneMoncriefDisk::variables(
   const auto specific_enthalpy = get<hydro::Tags::SpecificEnthalpy<DataType>>(
       variables(x, tmpl::list<hydro::Tags::SpecificEnthalpy<DataType>>{}, vars,
                 index));
-  auto rest_mass_density = make_with_value<Scalar<DataType>>(x, 0.0);
+  auto rest_mass_density = make_with_value<Scalar<DataType>>(x, 1e-15);
   variables_impl(vars, [&rest_mass_density, &specific_enthalpy, this](
                            const size_t s, const double /*potential_at_s*/) {
     get_element(get(rest_mass_density), s) =
         get(equation_of_state_.rest_mass_density_from_enthalpy(
             Scalar<double>{get_element(get(specific_enthalpy), s)}));
   });
+  //  if constexpr (std::is_same_v<DataType, DataVector>){
+  //  get(rest_mass_density) = max(1.0e-15, get(rest_mass_density));
+  //} else {
+  //  get(rest_mass_density) =
+  //  max(DataVector{get(rest_mass_density).size(), 1.0e-15},
+  //  get(rest_mass_density));
+  //}
   return {std::move(rest_mass_density)};
 }