Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

reading remapping weights generated by fregrid/ESMF into data_override #1399

Open
raphaeldussin opened this issue Oct 25, 2023 · 3 comments
Open
Labels
enhancement Issue/PR for a modification that increases performance, improves syntax, or adds functionality.

Comments

@raphaeldussin
Copy link

Is your feature request related to a problem? Please describe.

The horizontal remapping (horiz_interp) called from data_override is likely to become more and more
costly as horizontal resolution of forcing sets increases and has to be done for every job. Reading pre-computed
remapping weights could save time AND allow to use weights generated by other methods (e.g. ESMF's patch,...)

Describe the solution you'd like

In the data_table the interpolation field (usually bilinear/bicubic) would be able to receive a filename from
which the weights would be read and passed to data_override

Describe alternatives you've considered

Remapping to the model grid results in duplication of large files and need to be done for each model resolution so
is not terribly practical.

@raphaeldussin raphaeldussin added the enhancement Issue/PR for a modification that increases performance, improves syntax, or adds functionality. label Oct 25, 2023
@uramirez8707
Copy link
Contributor

In the data_table the interpolation field (usually bilinear/bicubic) would be able to receive a filename from
which the weights would be read and passed to data_override

Are you referring to the weights used by horiz_interp to do the interpolation?

Do you have an example of what this weight file will look like?

@raphaeldussin
Copy link
Author

correct!

I do have some old examples of comparing weights between fregrid and esmf here:
/work/Raphael.Dussin/sandbox/FREEDOM-NCtools/FREEDOM-NCtools/FRE2ESMF/examples

they both are variations over storing a sparse matrix. There might be some choices to make on what format to support (or both).

@theresa-morrison
Copy link

This feature would be helpful for some of the regional domains being developed. In particular, our Arctic configurations have a problem because the regional ocean domain is not re-entrant but the forcing from the data atmosphere should be interpolated across the zonal boundary.

Currently to solve this problem, all the forcing data must be interpolated onto the ocean grid offline, which for high resolution or long runs is a large volume of data to produce an manage. What Raf has suggested here would allow the correct weights for the atmosphere interpolation to be generated offline and used by FMS solving any issues at the seam in the forcing data.

I am really interested in having this feature in order to develop accurate regional Arctic simulations.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
enhancement Issue/PR for a modification that increases performance, improves syntax, or adds functionality.
Projects
None yet
Development

No branches or pull requests

3 participants