Hi Romain, are you sure, that you use the right x-values in the method boundaryTypes(..)? Try using the onUpperBoundary_(..) as several other test cases use them.
You do not need to change anything in neumann(..), because in bounadryTypes(..) you specified the injectors to be of Dirichlet type. Add the right value for dirichlet(..) and it should work. Salut Christoph _______________________________________________ Dumux mailing list [email protected] https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
