ai2cm / fv3gfs-fortran

FV3GFS Fortran for internal development at AI2
Other
6 stars 9 forks source link

Add option to pressure level coarsen outputs with limited extrapolation #366

Closed spencerkclark closed 1 year ago

spencerkclark commented 1 year ago

Summary

Here the extrapolation is using a limited nearest neighbor approach. The extrapolation is limited in the sense that we have only shifted the threshold for masking to the coarse grid's vertical layer midpoint pressures instead of the coarse grid's vertical layer bottom interface pressures. Using the latter results in zero extrapolation; using the former results in allowing values filled with nearest neighbor extrapolation to be included in averages if extrapolation was not required too far. This is mostly relevant over ocean where our default no-extrapolation approach requires masking roughly 50% of the grid cells in the lowest model level; this limited extrapolation approach enables retaining 100% of the grid cells over most ocean regions while maintaining the mask over land regions with variable topography.

2023-06-01-C3072-to-C48-unmasked-area-comparison

The above figure shows the fraction of the area of the C3072 grid cells left unmasked in the lowest model level when coarsening to C48 resolution for three different approaches. The left and middle panels represent our default and this new limited extrapolation approach. Unlimited extrapolation, the right panel, obviously results in no masking ever occurring. This plot was produced in this notebook.

Links

The PR implementing the equivalent functionality in fv3net is ai2cm/fv3net#2237. This notebook checks that the online approach implemented here approximately produces the same results as the offline approach implemented in https://github.com/ai2cm/fv3net/pull/2237.