Using a C++ productivity framework for the shallow water equations/Modeling Glacial Isostatic Adjustment and Sea Level changes in CESM

Recently, members of the CGD Oceangraphy section and CISL participated in the NOAA/NCAR/NREL hackathon.  The hackathon aimed to create a version of the shallow water equations model written in AMReX.  AMReX is a C++ productivity framework that provides a block-structured adaptive mesh refinement capability.  AMReX has several useful features, including single-source support for CPU and GPU,  simple boundary exchange and I/O capabilities,  and interoperability with both legacy Fortran and C codebases.  We describe the performance characteristics of our shallow water model on both CPU and GPU platforms.  We also discuss how AMReX could be used to simultaneously enable code modernization while minimizing disruption to scientific development.

CESM is a complex and powerful model that includes prognostic components for most of the major parts of earth’s climate system, including the atmosphere, oceans, sea ice, land surfaces, vegetation, ocean biogeochemistry, and ice sheets. But there is one major part of the Earth system that is not prognosed: the actual earth. This short exchange talk gives a brief overview of early efforts to couple the Citcom SVE3 solid earth model with CESM, and the benefits of including prognostic glacial isostatic adjustment and sea level changes in our earth system model.