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)
latitude <- c(4, 4, 4, 4, 4, 4)
longitude <- c(188, 188, 188, 188, 188, 188)
O2sol <- gsw_O2sol(SA,CT,p,longitude,latitude)
stopifnot(all.equal(O2sol/100, c(1.949651126384804, 1.958728907684003,
2.148922307892045, 2.738656506758550, 2.955109771828408,
3.133584919106894)))
Run the code above in your browser using DataLab