From cf8a956293a23708d00e4589aafde9a55dabf844 Mon Sep 17 00:00:00 2001 From: John Omotani Date: Mon, 25 Sep 2023 12:51:31 +0100 Subject: [PATCH] Restrict check for 1V in Krook collisions to moment-kinetic cases 2V Krook collisions are now implemented for full-f simulations. --- src/krook_collisions.jl | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/krook_collisions.jl b/src/krook_collisions.jl index 751d6caa1..1671326c5 100644 --- a/src/krook_collisions.jl +++ b/src/krook_collisions.jl @@ -42,8 +42,8 @@ Currently Krook collisions function krook_collisions!(pdf_out, fvec_in, moments, composition, collisions, vperp, vpa, dt) begin_s_r_z_region() - if vperp.n > 1 - error("Krook collisions not implemented for 2V case yet") + if vperp.n > 1 && (moments.evolve_density || moments.evolve_upar || moments.evolve_ppar) + error("Krook collisions not implemented for 2V moment-kinetic cases yet") end # Note: do not need 1/sqrt(pi) for the 'Maxwellian' term because the pdf is already