eta = (A0(flow) / tau0) * mtotal**(-5./3)