Thanks gdidio
Actually we figured that out, by applying a in/out mass transfer rate, in addition to the in/out flow transfer rate that goes with the 3rd kd bc. This is necessary, unless the C that you specified in the river bc is also the C in the groundwater, which is usually not the case by several orders of magnitude. Fortunately we knew which mass had to be loaded from the river to the aquifer, so we could calibrate the in/out mass transfer rates.
Cheers !