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