Update SO.py

pho critical = 2.7846602822750852e-08  #Unit [10^10 M_sun/h] / [(kpc/h)^3] at redshift 0