ebase_rho function

Seawater density calculation