issues/0039.R

library(gsw)
df <- structure(list(SA = c(36.855503939848, 36.4971650766637, 36.4977700026719,
36.2529246997542, 36.2534270655287, 36.002043110245, 36.0002357585083,
35.8278109198182, 35.8281142727647, 35.7015734367856, 35.7006709862616,
35.7848317547004, 35.7858367750849, 35.9496877557911, 35.9520999643311,
36.0623842502388, 36.0861726265337, 36.14217237462, 36.2275535176993,
36.2408756031207, 36.2187323903269, 36.1707742062984, 36.0918852617197,
35.9813688721453), CT = c(21.4947473754157, 16.6878140953016,
16.6748589933191, 14.5354547312802, 14.5417247209637, 13.1075384801307,
13.1059005757833, 11.9533561570296, 11.9365831694816, 10.834696340731,
10.8354205711472, 10.2515778760042, 10.2534330502052, 10.1594789517675,
10.1628546551997, 10.1319944946218, 10.000801742997, 9.97602409076545,
10.0160820710144, 9.89429410566537, 9.64473876660855, 9.23710853886353,
8.83733409331888, 8.23858450492611), p = c(11, 102.1, 103.1,
202.4, 202.4, 301, 301.6, 402.1, 403.1, 559.8, 560.5, 701.3,
701.4, 790.9, 791.3, 895.3, 943.4, 1046.5, 1137.4, 1193.3, 1242.6,
1292.7, 1345.1, 1405.4)), .Names = c("SA", "CT", "p"), row.names = c(NA,
-24L), class = "data.frame")
dh <- gsw::gsw_geo_strf_dyn_height(df$SA, df$CT, df$p)
print(data.frame(df$SA, df$CT, df$p, dh))
TEOS-10/GSW-R documentation built on April 14, 2025, 12:46 p.m.