Added Brunt-Vaisala frequency squared profile to MESA + example