Add Adusumilli et al. (2020) melt rates

[xylar]

This merge switches climatology maps of Antarctic melt rates to use Adusumilli et al. (2020) instead of Rignot et al. (2013). It also adds the Adusumilli et al. melt rates and fluxes to the time-series plots as a 3rd observational point (keeping Rignot et al. 2013 present-day and steady-state estimates).

To do:

• run test suite
• run for SORRM simulation
• update documentation

[xylar]

@cbegeman and @milenaveneziani, although I've marked you to review, there's not a big rush on this. It likely makes sense to wait until LCRC is back up next week and maybe until after the Hackathon to look this over.

[xylar]

Testing

I ran the test suite on Perlmutter. Results can be found here:
https://portal.nersc.gov/cfs/e3sm/xylar/analysis_testing/pm-cpu/add-adusumilli-melt-rates/

In particular, the Adusumilli et al. (2020) melt rates can be seen in both climatology maps:

and time series:

[cbegeman]

@xylar This all looks great! What testing would you like from me?

[xylar]

@cbegeman, good question. If you have time, you could run MPAS-Analysis using this branch on a decade of a recent SORRM simulation.

I plan to use this branch on the 3 runs for the paper (I just haven't got to it since LCRC came back up). Once that's done, I'll ask you to have a look at that, too.

[xylar]

@cbegeman, here's the function I would appreciate you sanity-checking if you have time.

[cbegeman]

@xylar I looked it over and I think it's correct. Thanks!

[xylar]

Here is the analysis from a 5-year run on SORRM that used prescribed melt rates from Adusumilli, which makes for a nice sanity check both on the analysis and on the prescribed melt fluxes:
https://web.lcrc.anl.gov/public/e3sm/diagnostic_output/ac.xasay-davis/analysis/20230516.GMPAS-IAF-DISMF.T62_SOwISC12to60E2r4.chrysalis/

Melt fields are nearly identical, with small differences due to interpolation artifacts:

Regionally averaged melt time series are similar but slightly smaller from MPAS-Ocean, e.g.:

[xylar]

@cbegeman, if the analysis I've done so far is sufficient for you to be happy, I don't think you need to do anything else in your review except give it a look. (A second pair of eyes on the regional preprocessing was the main thing I wanted.)

[xylar]

@milenaveneziani, do you have time to give this a look?

[milenaveneziani]

Looks great to me @xylar.

Sort of unrelated question: what does the 'SS' stand for in one of the Rignot's labels?

[milenaveneziani]

Approved based on visual inspection and developer testing.

[xylar]

Sort of unrelated question: what does the 'SS' stand for in one of the Rignot's labels?

Steady state. This would be the melt rates if the ice sheet were in equilibrium instead of losing mass. (So an approximation of what melt rates might have been in preindustrial times.)

[xylar]

Thanks @milenaveneziani!