Salinity data needs to be taken from the dat files and added to the data frame first,
tdata$salinity <- as.numeric(txdata$Salinity)
in the .dat import loop, same place temperature is added, then replace color = ~tempC with color = ~salinity in the plotly setup (it takes columns from the data frame specified at the start)
One problem I'm running into is that you can only use one colorscale per plot, still trying to hack my way around that, until then not sure how to display both data, maybe as circle lines, or offset slightly in one direction to make a two sided curtain..
I made the bathy work on ITP 116 which is floating around the north pole and just crossing the Lomonosov ridge, you can see the buoy got stuck for a while in spot where there's an almost 2km vertical cliff / cutout in the ridge. Takes a bit to load but worth it:
https://rpubs.com/macid/ITP116_TempC_Bathy (hover over a bathy point to see depth in meters, not to z-axis scale..)