fpe with WIND star accretion type
In that mode, the radial speed is initialized with:
vradial(i,h,j) = -(2*torqueCoeff/(1
+exp(-steepness*(1
-dead(i,h,j)
/userActiveDensity)))
/radii(i));});
But dead(i,h,j)
can have significantly big values (>150) with -steepness*(1-dead(i,h,j)/userActiveDensity)
reaching 839.119345 which is too big for an exponential.