FET010 = readtable("FET010_Tris_9.4.25.csv");
[FET010.pH_free_verify, FET010.pH_tot_verify] = phcalc_jp(FET010.Vrse,0,FET010.Temperature,FET010.Salinity,-0.92716570766144,7.406790E-04,[0, 0, 0 ,0 , 0, 0]')
FET010.diff = FET010.pHtot - FET010.pH_tot_verify

FET011 = readtable("FET011_Tris_9.4.25.csv");
[FET011.pH_free_verify, FET011.pH_tot_verify] = phcalc_jp(FET011.Vrse,0,FET011.Temperature,FET011.Salinity,-0.99686991849884,3.260250E-04,[0, 0, 0 ,0 , 0, 0]')
FET011.diff = FET011.pHtot - FET011.pH_tot_verify