Thank you. I implemented as you suggested and it works very good now. This approach seems also applicable for the mesh that are not aligned in the interface.
Best, Zhekai On Mon, Oct 24, 2016 at 9:00 AM, Daniel Wheeler <[email protected]> wrote: > On Mon, Oct 24, 2016 at 12:20 AM, Zhekai Deng > <[email protected]> wrote: > > > > This is not a issue if mesh system is 1D, and boundary cell is just 1. Or > > maybe on the rectangular mesh system. However, does any one has any idea > to > > how to transfer it properly ? > > Hi Zhekai; I haven't read through everything but I would guess that > you need to interpolate from one mesh to another. To do that you can > use one of Scipy's inerpolation functions such as > > https://docs.scipy.org/doc/scipy-0.14.0/reference/ > generated/scipy.interpolate.griddata.html > > We also built that capability into FiPy, but it's probably better to > use Scipy's interpolation at this point. > > https://github.com/usnistgov/fipy/blob/develop/fipy/ > variables/cellVariable.py#L167 > > If the mesh points are aligned, it is still a good idea to use the > interpolation approach as it removes the issue of worrying about mesh > ordering. > > -- > Daniel Wheeler > _______________________________________________ > fipy mailing list > [email protected] > http://www.ctcms.nist.gov/fipy > [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] >
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
