Yes, thanks! Well, putting NO3CONC = 0 and Kfraction = 0 (as is in the KEDS screenshot) makes things a bit simpler. I also think that, using screenshot conditions, exhal must be = 0.1. Interesting noting that juncpot goes to infinite as z2 tends to zero. (I think that often taking limits gives some insight in what is going on..) I'll ponder on this...
Do you happen to have something about activity coefficient, which may explain why my results of post #23 are always a bit off?