add calcHsAndRho_s()