diff --git a/nufi/parameters.h b/nufi/parameters.h
index 65fb282..e071530 100644
--- a/nufi/parameters.h
+++ b/nufi/parameters.h
@@ -34,7 +34,7 @@ constexpr unsigned int CONVERGENCE_ITERATIONS = 5000;
 constexpr double CONVERGENCE_LIMIT = 1e-8;
 
 // Adaptive refinement options
-constexpr unsigned int REFINE_FREQUENCY = 30;
+constexpr unsigned int REFINE_FREQUENCY = 10;
 constexpr double REFINEMENT_TOP_FRACTION = 0.8;
 constexpr double REFINEMENT_BOTTOM_FRACTION = 0.1;
 
@@ -51,7 +51,7 @@ constexpr double DT = 1. / 10.;
 constexpr unsigned int TMAX = 100;
 
 // Plotting options
-constexpr int PLOT_FREQUENCY = 10;
+constexpr int PLOT_FREQUENCY = 5;
 constexpr size_t PLOT_NX = 512;
 constexpr double PLOT_DX = LX / PLOT_NX;
 const std::string PLOT_DIR = "results/";
diff --git a/src/nufi_solver.cc b/src/nufi_solver.cc
index 8d0511d..55731d9 100644
--- a/src/nufi_solver.cc
+++ b/src/nufi_solver.cc
@@ -152,7 +152,7 @@ NuFISolver::eval_rho(unsigned int n, std::vector<double> &X,
 
   const double dv =
       (Parameters::V_DOMAIN_RIGHT - Parameters::V_DOMAIN_LEFT) / Nv;
-  const double v_min = Parameters::V_DOMAIN_LEFT;
+  const double v_min = Parameters::V_DOMAIN_LEFT + .5 * dv;
 
 #pragma omp parallel for
   for (unsigned int i = 0; i < Nv; ++i) {