search

giovannellilab / geochem

Geochemistry pipeline for processing and visualising data from IRMS, IC and ICP-MS
https://github.com/giovannellilab/geochem
Other
0 stars 0 forks source link

Add Gibbs plot #38

Closed guillecg closed 2 months ago

guillecg commented 3 months ago

Check Tamburello et al. 2022

alebastianoni commented 3 months ago

ggplot(subset(data, Funding %in% c("BMS", "ERC","ERC/NASA Habitable Worlds","ERC/FEAMP")), aes(x=(na/(na+ca)), y=tds,label=SiteID))+ geom_point(aes(color=ExpID, size=T, alpha=0.99))+ geom_text(data=subset(data, SiteID %in% c("SW_CF")))+ scale_color_viridis_d()+ guides(fill=guide_legend(override.aes=list(shape=21)))+ geom_segment(aes(x = 0.05, y= 1e+2, xend = 0.9, yend = 10), colour = "lightgrey", alpha=0.3) + geom_segment(aes(x = 0.05, y= 1e+3, xend = 0.9, yend = 1e+5), colour = "lightgrey", alpha=0.3) + geom_segment(aes(x = 0.6, y= 1e+3, xend = 0.9, yend = 1e+2), colour = "lightgrey", alpha=0.3) + geom_segment(aes(x = 0.6, y= 1e+3, xend = 0.9, yend = 1e+4), colour = "lightgrey", alpha=0.3) +

annotate(geom="text", x=0.2, y=2e+2, label="Rock", color="black", size=3, alpha=0.6) + annotate(geom="text", x=0.6, y=1e+2, label="Precipitation", color="black", size=3, alpha=0.6) + annotate(geom="text", x=0.4, y=2e+4, label="Evaporation", color="black", size=3, alpha=0.6) + scale_y_log10()+ annotation_logticks(sides = "l") + labs(title="", shape="", fill="Temperature (°C)", x="Na/(Na+Ca)", y="TDS (ppm)")+ theme_glab(base_size=18)+ theme(legend.position="right")

guillecg commented 2 months ago

Solved in 0516f541100d25df3e765fb902caab5a4381cfd6