using PhysicalOceanography
isapprox(f(45),1.03125867182e-4)