算法[最好是fortran]来插入来自2D非结构化网格的数据以形成笛卡尔网格

Question

我正在使用2D动态非结构化网格对裂缝传播进行建模.随着裂缝随时间传播,元素相应地移动.对于给定的时间步长,我想插入我的非结构化网格的数据zi = f(xi,yi)(其中(xi,yi)是非结构化三角形的节点)以获得函数z的值.笛卡尔网格(其中x和y有规律地间隔).我的网格大小约为100x100个节点.有很多方法可以实现这种插值方案,但我对该主题知之甚少,无法确定哪种方案足够健壮和简单.我正在使用Fortran 90(不要问为什么......).有没有可用的开源算法？我不想重新发明轮子.

非常感谢!

Answer 1

您可能对地球系统建模框架 (ESMF) 感兴趣：

http://www.earthsystemmodeling.org/

它配备了 Fortran、C 和最近的 Python（功能有限）API 和优秀的文档。它通常用于地球物理模型及其耦合中的数据抽象对象，但它还提供了在结构化和非结构化网格之间重新网格化的例程，可用于离线（独立命令行实用程序）或在线（通过子例程调用）插值权重一代。

上次我研究它时，ESMF 提供了双线性、双三次和数量守恒重新网格化方法。