gsw_Turner_Rsubrho {gsw} | R Documentation |
This uses the 75-term density equation. The values of Turner Angle
Tu and density ratio Rrho are calculated at mid-point pressures, p_mid
.
gsw_Turner_Rsubrho(SA, CT, p)
SA |
Absolute Salinity [ g/kg ] |
CT |
Conservative Temperature [ degC ] |
p |
sea pressure [dbar], i.e. absolute pressure [dbar] minus 10.1325 dbar |
The present R function works with a wrapper to a C function contained within the GSW-C system (Version 3.05-4 dated 2017-08-07, available at https://github.com/TEOS-10/GSW-C, as git commit '5b4d959e54031f9e972f3e863f63e67fa4f5bfec'), which stems from the GSW-Fortran system (https://github.com/TEOS-10/GSW-Fortran) which in turn stems from the GSW-Matlab system (https://github.com/TEOS-10/GSW-Matlab). Consult http://www.teos-10.org to learn more about these software systems, their authorships, and the science behind it all.
List containing Tu
[ degrees ], Rsubrho
[ unitless ], and p_mid
[ dbar ]
http://www.teos-10.org/pubs/gsw/html/gsw_Turner_Rsubrho.html
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) r <- gsw_Turner_Rsubrho(SA, CT, p) expect_equal(r$Tu, c(-2.063858905281147, 41.758435216784427, 47.606966981687535, 53.710351151706369, 45.527063858211527)) expect_equal(r$Rsubrho, 100*c(-0.009304335069039, -0.176564834348709, 0.219627771740757, 0.065271424662002, 1.087044054679743)) expect_equal(r$p_mid, 100*c(0.300, 0.875, 1.875, 4.250, 8.000))