Numerical issues for coupling biological models with isopycnal mixing schemes