gcmfaces icon indicating copy to clipboard operation
gcmfaces copied to clipboard

Published 20 hours ago • verified publisher icon MITgcm

Flux conserving, bin average of 2-D field

Open dcarrollsci opened this issue 5 years ago • 3 comments

Hi Gaël,

Does gcmfaces allow for flux conserving, bin averaging from LLC output to lon/lat grid? The application of this is bin averaging of ECCO-Darwin air-sea CO2 fluxes for various end users in the carbon community.

Thanks, Dustin

dcarrollsci avatar Sep 13 '19 19:09 dcarrollsci

I do not know of anyone who might have implemented this type of regridding but would definitely be a very useful addition. Are you thinking about schemes that deal with c-grid vector fields (U,V,W) or more common "conservative interpolation" schemes?

gaelforget avatar Sep 23 '19 22:09 gaelforget

It is for surface fluxes. Dustin and I are writing code for llc270 to 1-deg-grid bin-averaging. Happy to share when done.

menemenlis avatar Sep 24 '19 16:09 menemenlis

Would be great. Here is what I would suggest:

    1. take a look at sample_processing/example_interp.m which is mentioned in the gcmfaces manual under tutorial examples. interp=gcmfaces_interp_coeffs(lon(:),lat(:)); should work with any of our llc, cs, etc grids, and give you the four nearest neighbors. No need to redo this part -- build on it as much as possible.
    1. Fork, branch, and put your proposed implementation in the gcmfaces_devel/ folder.
    1. PR back to this repo once you have an implementation that you want me to review.
    1. Discuss and iterate on the PR if needed

Just like what we routinely do with MITgcm (see https://mitgcm.readthedocs.io/en/latest/contributing/contributing.html for detailed instructions if needed).

ok?

ps: could you please add, in this thread, the paper reference for the scheme you are looking at with Dustin.

gaelforget avatar Sep 25 '19 15:09 gaelforget