diff --git a/src/initial_conditions.jl b/src/initial_conditions.jl index 14a7758c6..1634afbb8 100644 --- a/src/initial_conditions.jl +++ b/src/initial_conditions.jl @@ -608,7 +608,7 @@ function init_charged_pdf_over_density!(pdf, spec, composition, vpa, vperp, z, # # Implemented by multiplying by a smooth 'notch' function # notch(v,u0,width) = 1 - exp(-(v-u0)^2/width) - width = 0.1 + width = 0.1 * vth[iz] inverse_width = 1.0 / width @. pdf[:,ivperp,iz] *= 1.0 - exp(-vpa.grid^2*inverse_width)