Added function to calculate magnetosonic speed
Break MHD speeds calculation into two : magnetosonic speed (=sqrt(cs2 + v_Alfven2)) and fast magnetosonic wave phase velocity (same, augmented with a projection factor). First is used in Dedner, second is used in v_sig