Addressing numerical challenges in introducing a reactive transport code into a land surface model: A biogeochemical modeling proof-of-concept with CLM-PFLOTRAN 1.0: Modeling Archive

 Modeling Archive Citation

G. Tang, F. Yuan, G. Bisht, G. E. Hammond, P. C. Lichtner, J. Kumar, R. T. Mills, X. Xu, B. Andre, F. M. Hoffman, S. L. Painter, and P. E. Thornton. Addressing numerical challenges in introducing a reactive transport code into a land surface model: A biogeochemical modeling proof-of-concept with CLM-PFLOTRAN 1.0: Modeling Archive. 2016. Next Generation Ecosystem Experiments Arctic Data Collection, Carbon Dioxide Information Analysis Center, Oak Ridge National Laboratory, Oak Ridge, Tennessee, USA. Accessed at http://dx.doi.org/10.5440/1239799

Abstract

This Modeling Archive is in support of an NGEE Arctic discussion paper under review and available at doi:10.5194/gmd-9-927-2016.

The purpose is to document the simulations to allow verification, reproducibility, and follow-up studies. This dataset contains shell scripts to create the CLM-PFLOTRAN cases,  specific input files for PFLOTRAN and CLM, outputs, and python scripts to make the figures using the outputs  in the publication. Through these results, we demonstrate that CLM-PFLOTRAN can approximately reproduce CLM results in selected cases for the Arctic, temperate and tropic sites. In addition, the new framework facilitates mechanistic representations of soil biogeochemistry processes in the land surface model.

Related NGEE Arctic Publication Citation

This Modeling Archive is provided in support of the following paper. Please cite this paper in addition to the modeling archive for full attribution of the modeling endeavour.

Tang, G., Yuan, F., Bisht, G., Hammond, G. E., Lichtner, P. C., Kumar, J., Mills, R. T., Xu, X., Andre, B., Hoffman, F. M., Painter, S. L., and Thornton, P. E.: Addressing numerical challenges in introducing a reactive transport code into a land surface model: a biogeochemical modeling proof-of-concept with CLM–PFLOTRAN 1.0, Geosci. Model Dev., 9, 927-946, doi:10.5194/gmd-9-927-2016, 2016.


Modeling Archive Contents

 

Model:

PFLOTRAN-DEV

https://bitbucket.org/pflotran/pflotran-dev/

changeset: 86fb10a809ea

CLM-PFLOTRAN

clm-pflotran: https://www.bitbucket.org/clm_pflotran/clm-pflotran-ngee-sci (changeset 3837c12)

pflotran-clm: https://www.bitbucket.org/clm_pflotran/pflotran-clm-ngee-sci (changeset 6347d01)

Input Data:

https://www.bitbucket.org/clm_pflotran/clm-pflotran-data-trunk-testing (changeset 1fe9844)

Parameters:

CLM parameters:

brw1.nc

pit1.nc

cax1.nc

PFLOTRAN parameters

hanford-clm.dat

mass_transefer_decomp.h5

sgrid-1x1x10-clm2pf.meshmap

sgrid-1x1x10-clm2pf_surf.meshmap

sgrid-1x1x10-pf2clm.meshmap

US-Brw_I1850CLM45CN_ad_spinup.in

Initializations:

Scripts

brw

pit

cax

Output:

brw.tgz

pit.tgz

cax.tgz

Configurations:

PETSC (changeset: 821a7925fede8aa3b3b482fc9ccb2d087e2f80fa) on oic

./configure \

   PETSC_DIR=/projects/cesm/devtools/petsc \

   PETSC_ARCH=arch-linux2-gcc4.8.1-mpich3.0.4-opt \

   --download-cmake=yes \

   --download-parmetis=yes \

   --download-metis=yes \

   --with-c2html=no \

   --with-debugging=0 \

   COPTFLAGS=-O1 \

   FOPTFLAGS=-O1 \

   --download-hdf5=yes \

   --with-mpi-dir=/projects/cesm/devtools/mpich-3.0.4-gcc4.8.1

 

~                                                                  

Post Processing:

Figure 4: brw500yl.py (tar xzf brw.tgz)

Figure 5: pit300yl.py (tar xzf pit.tgz)

Figure 6: cax300yl.py (tar xzf cax.tgz)

Results:

http://www.geosci-model-dev.net/9/927/2016/