22 double x =
MIN2( 0.01, 1. - r_star/rad );
23 double v_r = v_star + (v_inf - v_0) * sqrt( Beta1 * x + (1.-Beta1) *
pow(x, Beta2) );
27 double density = Mdot / ( PI4 * ATOMIC_MASS_UNIT * mu *
pow2(rad) * v_r);
double density(const genericState &gs)
double pow(double x, int i)
double dense_parametric_wind(double rad)