Remove a normalisation by vth in the krook operator that shouldn't be…

… there (grid is already normalised when evolve_ppar = true
mabarnes · May 22, 2024 · e2a6c33 · e2a6c33
1 parent 7650b7d
commit e2a6c33
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/moment_kinetics/src/krook_collisions.jl b/moment_kinetics/src/krook_collisions.jl
@@ -115,8 +115,8 @@ function krook_collisions!(pdf_out, fvec_in, moments, composition, collisions, v
             @loop_vperp_vpa ivperp ivpa begin
                 pdf_out[ivpa,ivperp,iz,ir,is] -= dt * nu_ii *
                     (fvec_in.pdf[ivpa,ivperp,iz,ir,is]
-                     - exp(-((vpa.grid[ivpa] - fvec_in.upar[iz,ir,is])/vth)^2
-                           - (vperp.grid[ivperp]/vth)^2))
+                     - exp(-((vpa.grid[ivpa] - fvec_in.upar[iz,ir,is]))^2
+                           - (vperp.grid[ivperp])^2))
             end
         end
     elseif moments.evolve_upar