# NOT RUN {
SA <- c(34.7118, 34.8915, 35.0256, 34.8472, 34.7366, 34.7324)
CT <- c(28.8099, 28.4392, 22.7862, 10.2262, 6.8272, 4.3236)
p <- c( 10, 50, 125, 250, 600, 1000)
rho <- gsw_rho(SA,CT,p)
expect_equal(rho/1e3, c(1.021839935738108, 1.022262457966867, 1.024427195413316,
1.027790152759127, 1.029837779000189, 1.032002453224572))
# }
Run the code above in your browser using DataLab