diff --git a/src/standalone/model/host_memory.cc b/src/standalone/model/host_memory.cc
index 4697690..94dbb5a 100644
--- a/src/standalone/model/host_memory.cc
+++ b/src/standalone/model/host_memory.cc
@@ -71,7 +71,7 @@ acmesh_create(const AcMeshInfo& mesh_info)
     return mesh;
 }
 
-static void
+void
 vertex_buffer_set(const VertexBufferHandle& key, const AcReal& val, AcMesh* mesh)
 {
     const int n = acVertexBufferSize(mesh->info);
diff --git a/src/standalone/model/host_memory.h b/src/standalone/model/host_memory.h
index e2958cd..c724e5e 100644
--- a/src/standalone/model/host_memory.h
+++ b/src/standalone/model/host_memory.h
@@ -48,6 +48,8 @@ extern const char* init_type_names[]; // Defined in host_memory.cc
 
 AcMesh* acmesh_create(const AcMeshInfo& mesh_info);
 
+void vertex_buffer_set(const VertexBufferHandle& key, const AcReal& val, AcMesh* mesh);
+
 void acmesh_clear(AcMesh* mesh);
 
 void acmesh_init_to(const InitType& type, AcMesh* mesh);
diff --git a/src/standalone/simulation.cc b/src/standalone/simulation.cc
index f4f7b36..8159084 100644
--- a/src/standalone/simulation.cc
+++ b/src/standalone/simulation.cc
@@ -96,7 +96,6 @@ write_mesh_info(const AcMeshInfo* config)
     fprintf(infotxt, "real AC_sink_pos_z %e \n", (double)config->real_params[AC_sink_pos_z]);
     fprintf(infotxt, "real AC_M_sink     %e \n", (double)config->real_params[AC_M_sink]);
     fprintf(infotxt, "real AC_soft     %e \n", (double)config->real_params[AC_soft]);
- 
     fprintf(infotxt, "real AC_G_const     %e \n", (double)config->real_params[AC_G_const]);
     
     fclose(infotxt);
@@ -196,6 +195,8 @@ run_simulation(void)
     // TODO: This need to be possible to define in astaroth.conf
     acmesh_init_to(INIT_TYPE_GAUSSIAN_RADIAL_EXPL, mesh);
 
+    vertex_buffer_set(VTXBUF_ACCRETION, 0.0, mesh);
+
     acInit(mesh_info);
     acLoad(*mesh);