How do I implement a boundary condition \partial T /\partial y = 0 at facesTop (for instance)?
It seems to me that doing var.faceGrad.constrain(((0,),(0,)), where=mesh.facesTop) (the initial example in the manual) I would be like doing \partial T /\partial x + \partial T /\partial y = 0... Thanks a lot Francisco Vega _______________________________________________ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]