# NOT RUN {
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_shallow <- c(10,      50,     125,     250,     600,    1000)
p_deep <- c(  110,     150,     225,     350,     700,    1100)
ed <- gsw_enthalpy_diff(SA, CT, p_shallow, p_deep)
expect_equal(ed/1e2, c(9.784180644568052, 9.780195056105020, 9.759587700515114,
                       9.727552719534447, 9.708223170174454, 9.687871289079633))
# }
Run the code above in your browser using DataLab