- TODO Fix  Ndens  argument treatment:  

Right  now  code  is oblivius  to  --Ndens
parameter, since eta ~1/Ndens sits  in global parameter section and already
defined. So it is not recalculated after the parameter assignment.

Ugly hack  would be to  make variable  eta_times_Ndens instead of  eta, and
recalculate  eta=eta_times_Ndens/Ndens  in the  ODE  part  around where  it
needed. It will have performance penalty but probably not too large.


- TODO Fix  Temperature  argument treatment:  

Same problem as above