diff --git a/src/overlapping/trapezoid_gradients.jl b/src/overlapping/trapezoid_gradients.jl
index 04b3d04b22f8aadf9ac34e2abc09db80eb02f6bd..39fb32eaa8a3351cc26c7b19da391e8dec4d9269 100644
--- a/src/overlapping/trapezoid_gradients.jl
+++ b/src/overlapping/trapezoid_gradients.jl
@@ -37,7 +37,7 @@ If not set, they will be determined during the sequence optimisation.
 - [`qval`](@ref)/[`qvec`](@ref): Spatial scale on which spins will be dephased due to this pulsed gradient in rad/um (given by `δ` * `gradient_strength`).
 ### Pulse variables
 Any variables defined for the specific pulse added. Also:
-- [`slice_thickness`](@ref): vector with the slice thickness
+- [`slice_thickness`](@ref): vector with the slice thickness in mm.
 
 The [`bvalue`](@ref) can be constrained for multiple gradient pulses.
 """
@@ -130,7 +130,7 @@ slew_rate(g::TrapezoidGradient) = g.slew_rate
 δ(g::TrapezoidGradient) = rise_time(g) + flat_time(g)
 duration(g::TrapezoidGradient) = 2 * rise_time(g) + flat_time(g)
 qvec(g::TrapezoidGradient, ::Nothing, ::Nothing) = δ(g) .* gradient_strength(g) .* 2π
-inverse_slice_thickness(g::TrapezoidGradient) = isnothing(g.pulse) ? nothing : inverse_bandwidth(g) .* gradient_strength(g)
+inverse_slice_thickness(g::TrapezoidGradient) = isnothing(g.pulse) ? nothing : inverse_bandwidth(g) .* gradient_strength(g) .* 1000
 
 function variables(tg::TrapezoidGradient) 
     list = [qvec, δ, gradient_strength, duration, rise_time, flat_time]
diff --git a/src/variables.jl b/src/variables.jl
index 92025d45c9ee3c323dc954fca9293693fd5badd2..7539c5b02f67d996f6d62aaa1d3e1351c8be5fcb 100644
--- a/src/variables.jl
+++ b/src/variables.jl
@@ -19,8 +19,8 @@ all_variables_symbols = [
         :inverse_bandwidth => "Inverse of the [`bandwidth`](@ref) of the RF pulse in ms",
         :N_left => "The number of zero crossings of the RF pulse before the main peak",
         :N_right => "The number of zero crossings of the RF pulse after the main peak",
-        :slice_thickness => "Slice thickness of an RF pulse that is active during a gradient.",
-        :inverse_slice_thickness => "Inverse of the [`slice_thickness`](@ref) in 1/um.",
+        :slice_thickness => "Slice thickness of an RF pulse that is active during a gradient in mm.",
+        :inverse_slice_thickness => "Inverse of the [`slice_thickness`](@ref) in 1/mm.",
     ],
 
     # gradients