How to build the coupled CM2Mc-LPJmL5.1 model

Stefan Petri and the POEM team
poem@pik-potsdam.de
Wed Apr 14 16:09:34 CEST 2021

This is supplementary material to Drüke, M., von Bloh, W., Petri, S., Sakschewski, B., Schaphoff, S., Forkel, M., Huiskamp, W., Feulner, G., and Thonicke, K.: CM2Mc-LPJmL v1.0: Biophysical coupling of a process-based dynamic vegetation model with managed land to a general circulation model, Geosci. Model Dev. Discuss. [preprint], https://doi.org/10.5194/gmd-2020-436, in review, 2021. A short note about the naming: CM2Mc is the project name of the coarse-grid atmosphere-ocean general circulation model developed by Eric Galbraith et al [1] [2]. This setup found its way into the ``official'' MOM5 repository [3] as test case CM2M_coarse_BLING [4]. For better readability, we mostly use the name ``CM2Mc'' as meaning MOM5 with CM2M_coarse_BLING configuration.

CM2Mc-LPJmL5.1 is an experiment setup, where LPJmL is plugged into CM2M_coarse_BLING. The land model component land_lad resp. land_atlantes and LPJmL are run at 0.5x0.5 deg resolution. Ocean and atmosphere are unchanged from CM2M_coarse_BLING.

Prepare the MOM5 source code

Prepare the LPJmL source code

Download and unpack the LPJmL code. We assume that it will be located in the directory POEM-GMD-CODE/LPJmL.
cd POEM-GMD-CODE/LPJmL
./configure.sh
That creates a file named Makefile.inc , which you could adapt to your computing environment, if necessary.

Configure and build the coupled model executable

We use FMS coupler that is distributed with MOM5 as the main program for the coupled model.
cd POEM-GMD-CODE/MOM5
For that, create an environment definition and a makefile template for your platform according to the instructions on https://mom-ocean.github.io/docs/quick-start-guide/ .

We choose pik-hlrs2015-ifort as our platform name. We define a Shell variable with that name, and create files

export PLATFORM=pik-hlrs2015-ifort
bin/environs.$PLATFORM
bin/mkmf.template.$PLATFORM
The build process in exp/lpj_compile.csh expects to find the to-be-used LPJmL version in a subdirectory src/land_atlantes/LPJmL. Thus, create a symbolic link src/land_atlantes/LPJmL which points to the LPJmL source code, which was unpacked and configure above:
cd src/land_atlantes
ln -s ../../../LPJmL LPJmL
cd ../..

Now compile an executable:

cd exp 
./POEM_compile.csh --type MOM_ATLANTES_LPJ_CM2
Eventually, this creates the executable model exec/$PLATFORM/MOM_ATLANTES_LPJ_CM2/fms_MOM_ATLANTES_LPJ_CM2.x .

Input Data

Input data for the MOM5 components

The CM2M_coarse_BLING test case [4] that accompanies the ``official'' MOM5 code is used as input data for the ocean, sea-ice, atmosphere and land_lad model components, with the following few modifications.

CM2M_coarse_BLING runs the land model component on 96x60 cells. For LPJmL we use 0.5deg resolution. Thus build a new land grid definition, and then new exchange grid specifications, including a new land-sea mask, using the tools provided with the MOM5 distribution, following the instructions in the MOM5 documentation.

For some reason, the orography height definition used by Galbraith et al [2] is not included in the ``official'' CM2M_coarse_BLING input data [4]. Following [5], the surface definition file named surf.res is extracted from the restart file fv_rst.res.nc provided by Eric Galbraith [6] [7].

The field_table [9] contains some cross-land mixing parameterisations which are not contained in the ``official'' CM2M_coarse_BLING test case [4].

Following the discussion in [5] we changed some namelist parameters to the values used for Galbraith et al. (2011) [2].

Input data for the LPJmL component

The above-created land-sea mask is used to generate a grid definition file for LPJmL. Soil data is taken from the Harmonized World Soil Database (HWSD) [10]. Lake fractions are taken from the Global Lakes and Wetlands Database (GLWD) [11]. The drainage map is taken from DDM30 [12]. Land use data is taken from [13].

We run the MOM5 components, with the above-stated input and configuration, without the LPJmL component, for 50 years to aproximate a steady state. Then we run it for another 30 years with daily output of the atmosphere-to-land fluxes. That is used to construct a climatology for a LPJmL spinup run. LPJmL is run with bare soil as starting condition for 5000 years, using a 30-year climatology as forcing.

References

[1]
CM2Mc Simulation Library https://earthsystemdynamics.org/models/cm2mc-simulation-library/ Last access 15.4.2021
[2]
Galbraith, Eric D., Eun Young Kwon, Anand Gnanadesikan, Keith B. Rodgers, Stephen M. Griffies, Daniele Bianchi, Jorge L. Sarmiento, John P. Dunne, Jennifer Simeon, Richard D. Slater, Andrew T. Wittenberg, and Isaac M. Held. "Climate Variability and Radiocarbon in the CM2Mc Earth System Model". Journal of Climate 24(16) (2011): 4230-4254. https://doi.org/10.1175/2011JCLI3919.1
[3]
MOM5 Modular Ocean Model https://mom-ocean.github.io Last access 15.4.2021
[4]
CM2M_coarse_BLING example setup input data. http://portal.sf.utas.edu.au/thredds/catalog/momtest/test_data/catalog.html?dataset=momtest/test_data/CM2M_coarse_BLING.input.tar.gz Last access 15.4.2021 via https://mom-ocean.github.io/downloads/
[5]
North Pacific deep water with CM2M_coarse_BLING test case.
Posting in a discussion thread on the mom-users mailing list. Last access 15.4.2021
[6]
North Pacific deep water with CM2M_coarse_BLING test case.
Another posting in a discussion thread on the mom-users mailing list. Last access 15.4.2021
[7]
Initial conditions used for Galbraith et al. (2011) [2].Link obtained via [6]
[8]
Boundary conditions used for Galbraith et al. (2011) [2].Link obtained via [6]
[9]
Runscript used for Galbraith et al. (2011) [2].Link obtained via [6]
[10]
FAO/IIASA/ISRIC/ISSCAS/JRC, 2012. Harmonized World Soil Database (version 1.21). FAO, Rome, Italy and IIASA, Laxenburg, Austria. Data set downloaded 26.4.2012
[11]
Bernhard Lehner, Petra Döll, Development and validation of a global database of lakes, reservoirs and wetlands, Journal of Hydrology 296(1–4), 2004, pp 1-22. https://doi.org/10.1016/j.jhydrol.2004.03.028 (https://www.sciencedirect.com/science/article/pii/S0022169404001404)
[12]
Döll, P., Lehner, B. (2002): Validation of a new global 30-min drainage direction map. Journal of Hydrology, 258(1-4), 214-231. https://doi.org/10.1016/S0022-1694(01)00565-0
[13]
Fader, M., Rost, S., Mueller, C., Bondeau, A., and Gerten, D.: Virtual water content of temperate cereals and maize: Present and potential future patterns, J. Hydrol., 384, 218–231, https://doi.org/10.1016/j.jhydrol.2009.12.011, 2010.