- changed milestone to 8.0.1
Incompatible options in collisions
Issue #42
resolved
If ei_coll_only
is true then this effectively ignores the const_v=.true.
case.
I think the solution is to replace
if (ei_coll_only) then
vnew(ik,ie,is) = spec(is)%vnewk/energy(ie)**1.5 &
*zeff*0.5*tunits(ik)
with
if (ei_coll_only) then
if (const_v) then
vnew(ik,ie,is) = spec(is)%vnewk &
*zeff*0.5*tunits(ik)
else
vnew(ik,ie,is) = spec(is)%vnewk/energy(ie)**1.5 &
*zeff*0.5*tunits(ik)
endif
I have a small patch that implements this that I plan to push in the near future.
Does this seem sensible @colinmroach ?
Comments (2)
-
reporter -
- changed status to resolved
Closed by PR #55
- Log in to comment