CaCO3flux_l = (caco3flux_u+CaCO3_uptake(i,j,k)*drF(k)
& *drF(k) - CaCO3flux_l)*recip_drF(k)
CaCO3flux_u = CaCO3flux_l
_RL CaCO3flux_l