Converting ROMS data to finite element grid

[dwperkey]

Hello everyone. I am new to modeling so I apologize if this question is one that is generally known.

For my research I am in need of bringing bathy and hydro information from a ROMS model of the York River Estuary into the SMS supported Particle Tracking Model (PTM). PTM utilizes a finite element grid system, so I need a method to convert the ROMS curvilinear grid to finite element. In my search through the forum I noticed some people have mentioned utilizing ADCIRC data in their ROMS models (which I assume would require a conversion of finite element to curvilinear), but I was curious if commonly used methods/tools exist to go the opposite direction and convert the ROMS grid into finite element.

Any tips, hints, or resources would be greatly appreciated. This is my first venture into the modeling world so I'm not familiar with what might be a common practice within the community.

Thanks,

Dave

[kate]

This sort of problem is generally known as regridding or remapping. One remapping toolbox comes with ESMF which also has a Python interface. I'm sure there are Matlab tools for such things as well - we're always having to remap to get boundary conditions, for example.

[ngarg]

i routinely do something similar, where i remap data from a structured (stencil) type grid to unstructured (finite element) mesh. I use scrip for that purpose, also as Kate mentioned you can use EMSF (ESMpy). Also, you can use ncl (NCAR) which has some functions for working with both rectilinear, and unstructured grid.

NCL under the hood uses EMSF regridgen for the same purpose. There are some tools available for doing similar things for FVCOM, you can get started by looking at how they do it. I think they use MATLAB.

https://github.com/pwcazenave/fvcom-toolbox

and there is one more library, based in python from people at Stanford, where they use SUNTANS model, which is also FEM based.

https://github.com/mrayson/soda

[dwperkey]

Thanks for the tips. I'll look into both.