[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Small error in TA_ITER_1




There seems to be a small sign error in the calculation of
the derivative of alkalinity with respect to H+ in TA_ITER_1;
the term that comes from the alkalinity of the H+-ion has
the wrong sign:

In fn everything is o.k.:
	fn = k1*x*dic/b +
[snip]
     &	     sit/(1.0 + x/ksi) -  <------ 
     &	     x/c - [snip]

but in df = dfn/dx the term is wrong:
	df = ((k1*dic*b) - k1*x*dic*db)/b2 -
[snip]
     &	     sit/ksi/(1.0+x/ksi)**2. +  <------ 
     &	     1.0/c + [snip]          

But since [h+] is of the order of 10^-5 mol/m^3, this probably
makes no big difference.

Cheers, Christoph

**************************************************
Christoph Voelker
Dept. of Theoretical Oceanography
Institut fuer Meereskunde an der Universitaet Kiel
**************************************************