SA <- c(34.7118, 34.8915, 35.0256, 34.8472, 34.7366, 34.7324)
CT <- c(28.7856, 28.4329, 22.8103, 10.2600, 6.8863, 4.4036)
p <- c( 10, 50, 125, 250, 600, 1000)
r <- gsw_specvol_second_derivatives(SA, CT, p)
stopifnot(all.equal(r$specvol_SA_SA/1e-8, c(0.080906777599140,
0.080915086639384, 0.084568844270812, 0.096725108896007,
0.099111765836648, 0.100302277946072)))
stopifnot(all.equal(r$specvol_SA_CT/1e-8, c(0.129965332117084,
0.130523053162130, 0.149555815430615, 0.217023290441810,
0.233892039070486, 0.243659989480325)))
stopifnot(all.equal(r$specvol_CT_CT/1e-7, c(0.071409582006642,
0.071582962051991, 0.077436153664104, 0.095329736274850,
0.100105336953738, 0.103044572835472)))
stopifnot(all.equal(r$specvol_SA_p/1e-14, c(0.116889015000936,
0.116897424150385, 0.121500614193893, 0.136008673596132,
0.139023051292893, 0.140581903529772)))
stopifnot(all.equal(r$specvol_CT_p/1e-14, c(0.085542828707964,
0.086723632576213, 0.112156562396990, 0.188269893599500,
0.211615556759369, 0.228609575049911)))
Run the code above in your browser using DataLab