implemented ndem functionality into Darcy solver, not tested when ndem!=1