sphdem runs now, but calculation of omega is not good. need to change back to paper version, since rho is not consistant