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