diff --git a/Project.toml b/Project.toml
index aa99e07..d36bbe6 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "UltraDark"
 uuid = "1c8d022d-dfc0-4b41-80ab-3fc7e88cdfea"
 authors = ["Nathan Musoke <n.musoke@auckland.ac.nz>"]
-version = "0.2.2"
+version = "0.3.0"
 
 [deps]
 AbstractFFTs = "621f4979-c628-5d54-868e-fcf4e3e8185c"
diff --git a/src/UltraDark.jl b/src/UltraDark.jl
index 28ac5c2..4c1d36c 100644
--- a/src/UltraDark.jl
+++ b/src/UltraDark.jl
@@ -25,6 +25,14 @@ import .Output: output_summary_row, output_summary_header, output_grids, output_
 import .Output: SummaryStatistics, SummaryStatisticsMeanMaxRms
 import .Config: SimulationConfig, constant_scale_factor, TimeStepOptions
 
+"""
+    outer_step!(Δt, grids, constants; a=1.0)
+
+Perform the "outer" time step in the symmetrized split-step Fourier method.
+
+This step only updates the phase of ψ applying accelerations due to gravity,
+the amplitude is not changed.
+"""
 function outer_step!(Δt, grids, constants; a=1.0)
     psi_whole_step!(Δt, grids, constants)
 end
@@ -53,6 +61,13 @@ function psi_whole_step!(Δt, grids, constants)
     end
 end
 
+"""
+    inner_step!(Δt, grids, constants; a=1.0)
+
+Perform the "inner" time step in the symmetrized split-step Fourier method.
+
+This step applies the diffusion terms and updates the gravitational potential.
+"""
 function inner_step!(Δt, grids, constants; a=1.0)
     phi_whole_step!(Δt, grids, constants; a=1.0)
 end
@@ -145,8 +160,6 @@ function take_steps!(grids, t_start, Δt, n, output_config, a, constants)
     outer_step!(Δt/2, grids, constants)
     t += Δt / 2
 
-    output_summary_row(grids, output_config, t, a(t), Δt)
-
     t
 end