Switch slice_thickness unit from um to mm

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]

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