martens73/scarface

ScarFace - seacarb calculations with R Shiny user interface

scarface's Issues

seacarb::carb() and shiny app results differ

Hi

With flag=15, ALK=2200 µmol kg-1 and DIC=2000 µmol kg-1 the shiny app yields

Sample	flag	S	T	Patm	P	pH	CO2	fCO2	pCO2	fCO2pot	pCO2pot	fCO2insitu	pCO2insitu	HCO3	CO3	DIC	ALK	OmegaAragonite	OmegaCalcite
NA	15.00	35	25	1.00	0	7.87	17.62	620.53	622.52	620.53	622.52	620.53	622.52	1836.73	145.66	2000.00	2200.00	2.31	3.51

When I do it directly in R the result is

library(seacarb)
seacarb::carb(flag=15, var1=2200/1e6, var2=2000/1e6, S=35, T=25) # flag=15 --> var1 = Talk in mol kg-1; var2 = DIC in mol kg-1

flag	S	T	Patm	P	pH	CO2	fCO2	pCO2	fCO2pot	pCO2pot	fCO2insitu	pCO2insitu	HCO3	CO3	DIC	ALK	OmegaAragonite	OmegaCalcite
15	35	25	1	0	7.870532	1.739097e-05	612.5331	614.4945	612.5331	614.4945	612.5331	614.4945	0.001835282	0.0001473273	0.002	0.0022	2.338158	3.547316

I am aware that the shiny app returns different units than the carb() function, e.g. for ALK and DIC.

However, pCO2_app = 622.52 while pCO2_carb = 614.4945. Could you please explain these differences?

Thanks for any help and kind regards,
Chris

packageVersion("seacarb")
[1] ‘3.3.0’