Why doesn't Johnson-SU distribution give positive skewness in scipy.stats?

Question

The code below maps the statistical moments (mean, variance, skewness, excess kurtosis) generated by corresponding parameters of the Johnson-SU distribution (johnsonsu