Calc_phi_s() corrected