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_rho_second_derivatives(SA, CT, p)
stopifnot(all.equal(r$rho_SA_SA/1e-3, c(0.207364734477357, 0.207415414547223,
0.192903197286004, 0.135809142211237, 0.122627562106076,
0.114042431905783)))
stopifnot(all.equal(r$rho_SA_CT, c(-0.001832856561477, -0.001837354806146,
-0.001988065808078, -0.002560181494807, -0.002708939446458,
-0.002798484050141)))
stopifnot(all.equal(r$rho_CT_CT, c(-0.007241243828334, -0.007267807914635,
-0.007964270843331, -0.010008164822017, -0.010572200761984,
-0.010939294762200)))
all.equal(r$rho_SA_p, 1e-9*c(-0.617330965378778, -0.618403843947729,
-0.655302447133274, -0.764800777480716, -0.792168044875350,
-0.810125648949170))
all.equal(r$rho_CT_p, 1e-8*c(-0.116597992537549, -0.117744271236102,
-0.141712549466964, -0.214414626736539, -0.237704139801551,
-0.255296606034074))
Run the code above in your browser using DataLab