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)
tb <- gsw_thermobaric(SA, CT, p)
stopifnot(all.equal(tb*1e11,
c(0.141342632944971, 0.142352284525832, 0.163216280125501,
0.226030772122855, 0.246185239871747, 0.261474794884197)))
Run the code above in your browser using DataLab