Current version of carbonate.F90 consider hydrogen ions concentration (AHPLUS) as local variable and therefore this variable can't be saved.