I found a formula in your jupyter-notebook that make me puzzled.

[JLU-ICCMS-MaYuan]

$$
f_{1} f_{2}=\sqrt[3]{1+\left[\frac{\lambda}{2.46\left(1+3.8 \mu^{}\right)}\right]^{3 / 2}} \cdot\left[1-\frac{\lambda^{2}\left(1-\omega_{2} / \omega_{\log }\right)}{\lambda^{2}+3.312\left(1+6.3 \mu^{}\right)^{2}}\right] \nonumber
$$

I suspect that the denominator of f2 in the above formula is missing an item: (w2/wlog)^2.

The same problem occurs in the following code, which position is tc.py
f2 = 1 - (np.power(_lambda[pos], 2) * (1 - w2[pos] / wlog[pos])) / (
np.power(_lambda[pos], 2) + 3.312 * np.power(1 + 6.3 * mu, 2))