Using ModularVariable in a mesh generated by gmsh

Hi everyone, I've just started using FiPy. I'm particularely interested in the phase field problems. I tried the examples of solidification and grain impingement given in the documentation. What I intend to do is to use a non-rectangular domain with triangular mesh for the problem. After

Re: How to declare diffusion function

Kevin - You are correct that there's no way to project a FaceVariable back to a CellVariable. Since this is only for the initial condition, I would just explicitly set it, e.g., phi.value = (1 + x) / (gamma * psi.grad[0]) As for the boundary condition, I think FiPy's default no-flux