Incompatible options in collisions

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)