Hi, I have an additional question about this section of the docs (on conservative interpolation).

The example in the docs says:

Note that this is not a perfectly precise calculation. However, it’s sufficient to illustrate the basic principles of the calculation.

and below:

Finally, we attempt to plot a crude meridional overturning streamfunction for a single timestep.

What makes this calculation imprecise/crude, and what is missing to make it physically rigorous?
I see that the zonal integration assumes fixed zonal spacing of gridcells (dx), what else is missing?

Maybe this information could be mentioned in this section of the documentation to make it directly applicable to "real-world" analyses.

Please let me know if this is not the right place to ask this, I didn't want to open a new issue since it's not really an "issue" with xgcm but more of a science question.

Hi, I have an additional question about this section of the docs (on conservative interpolation).

The example in the docs says:

Note that this is not a perfectly precise calculation. However, it’s sufficient to illustrate the basic principles of the calculation.

and below:

Finally, we attempt to plot a crude meridional overturning streamfunction for a single timestep.

What makes this calculation imprecise/crude, and what is missing to make it physically rigorous? I see that the zonal integration assumes fixed zonal spacing of gridcells (dx), what else is missing?

Maybe this information could be mentioned in this section of the documentation to make it directly applicable to "real-world" analyses.

Please let me know if this is not the right place to ask this, I didn't want to open a new issue since it's not really an "issue" with xgcm but more of a science question.

Hi @linus-0 could you open a separate issue for this so we can keep issues tightly scoped? Thank you.