diff --git a/Source/Evolve.H b/Source/Evolve.H
index eca88b07..395634e7 100644
--- a/Source/Evolve.H
+++ b/Source/Evolve.H
@@ -28,8 +28,6 @@ namespace GIdx
 
 amrex::Real compute_dt(const amrex::Geometry& geom, const MultiFab& state, const FlavoredNeutrinoContainer& neutrinos, const TestParams* parms);
 
-void init_background_to_mesh(FlavoredNeutrinoContainer& neutrinos_rhs, const MultiFab& state, const Geometry& geom, const TestParams* parms);
-
 void deposit_to_mesh(const FlavoredNeutrinoContainer& neutrinos, amrex::MultiFab& state, const amrex::Geometry& geom);
 
 void interpolate_rhs_from_mesh(FlavoredNeutrinoContainer& neutrinos_rhs, const amrex::MultiFab& state, const amrex::Geometry& geom, const TestParams* parms);
diff --git a/Source/Evolve.cpp b/Source/Evolve.cpp
index 16cc23a5..2b6dcb73 100644
--- a/Source/Evolve.cpp
+++ b/Source/Evolve.cpp
@@ -103,31 +103,6 @@ Real compute_dt(const Geometry& geom, const MultiFab& state, const FlavoredNeutr
     return dt;
 }
 
-void init_background_to_mesh(FlavoredNeutrinoContainer& neutrinos_rhs, const MultiFab& state, const Geometry& geom, const TestParams* parms)
-{
-    const auto plo = geom.ProbLoArray();
-    const auto dxi = geom.InvCellSizeArray();
-
-    amrex::MeshToParticle(neutrinos_rhs, state, 0,
-    [=] AMREX_GPU_DEVICE (FlavoredNeutrinoContainer::ParticleType& p,
-                          amrex::Array4<const amrex::Real> const& sarr)
-    { 
-        double x,y,z;
-
-        for (int k = 0; k < parms->ncell[2]; ++k) {
-            for (int j = 0; j < parms->ncell[1]; ++j) {
-                for (int i = 0; i < parms->ncell[0]; ++i) {
-
-                    x =  plo[0] + dxi[0] * i + dxi[0] / 2.0;        
-                    y =  plo[1] + dxi[1] * j + dxi[1] / 2.0;         
-                    z =  plo[2] + dxi[2] * k + dxi[2] / 2.0;         
-
-                }
-            }
-        }
-    });
-}
-
 void deposit_to_mesh(const FlavoredNeutrinoContainer& neutrinos, MultiFab& state, const Geometry& geom)
 {
     const auto plo = geom.ProbLoArray();
diff --git a/Source/main.cpp b/Source/main.cpp
index c49639f9..d541c970 100644
--- a/Source/main.cpp
+++ b/Source/main.cpp
@@ -117,9 +117,7 @@ void evolve_flavor(const TestParams* parms)
 
     // Deposit particles to grid
     deposit_to_mesh(neutrinos_old, state, geom);
-    
-    init_background_to_mesh(neutrinos_old, state, geom, parms);
-    
+        
     // Write plotfile after initialization
     DataReducer rd;
     if (not parms->do_restart) {