Skip to content

AK K64G Demo Quick Start Guide 2b: C Cycle (adjusting Vcmax via flnr)

Fengming Yuan edited this page Sep 1, 2022 · 7 revisions

NGEE-Arctic: AK-K64G simulation example 2 - Parameterization

This example 2C runs a single-point, extracted from half-degree Global data, ELM v2 simulation, within which NGEE Kougarok Mile 64 study site located, Using user-tuned parameter, flnr, Fraction of Leaf Nitrogen in Rubisco, to adjust specific PFT's photosynthesis parameter Vcmax, namely VCMAX25TOP (top canopy Vcmax at 25 degree C ) in ELM.

(1) It will start from an initial restart file from full-run's spinup stage

(2) It will modify 'clm_params.nc', parameter file for ELM to read, in a temporary work directory.

(3) Re-run ELM, using newly-generated 'clm_params.nc'.

(4) Check specified PFT VCMAX25TOP (or GPP/NEP/NEE) output in jupyter notebook.

REPEAT (2)-(4), for purpose to obtain a reasonbly well GPP simulation for site, or Vcmax if observed.

I. Parametering ELM Demo Example B

NOTE, assuming that

(1) you've already done with AK-K64G Demo Quick Start Guide 2. one year simulation of transient run.

AND,

(2)you should have TWO (2) containers, one for running ELM model (but stopped), another for visualizing outputs in jupyter notebook which is active in your internet brower. Your laptop layout should be in some ways like in Parametering ELM Demo Example A

STEP 1. Modifying ELM PFT's leaf rubisco N fraction, flnr, to adjust VCMAX25TOP

  • Click Jupyter Notebook for visualizing's Left Panel Folder brower so that it back to folder /work/, like following:

Screen Shot 2022-08-30 at 2 04 53 PM

  • Click open folder work, and one file, clm_params.txt, like following:

Screen Shot 2022-08-30 at 2 10 26 PM

NOTE that clm_params.txt is an ASCII file converted from clm_params.nc which actually read by ELM, and editable in juypter notebook text editor.

  • Scrolling down the file clm_params.txt to Line 375-378 (shown in shade), where are the header information of parameter flnr.
    Screen Shot 2022-08-30 at 2 35 19 PM

  • Further Scrolling down to Line 1104-1129 (shown in shade), where are the PFT names and order for current ELM model. Screen Shot 2022-08-30 at 2 39 07 PM

NOTE that Line 1116 and 1117 are two primary PFTs possibly exit in high-latitude region, including Arctic, in default ELM model. REMEMBER their orders, 12th and 13th respectively. NEXT we will try to edit the 13th (c3_arctic_grass) value of parameter flnr.

  • Further Scrolling down to Line 1485-1487 (shown in shaded are the 12th and 13th values for what we're interested in this Demo).

Screen Shot 2022-08-30 at 2 53 10 PM

  • For an example, modify default value of 0.1365 of the 13th (c3_arctic_grass), to 0.2365.

Screen Shot 2022-08-30 at 2 56 21 PM

NOTE that Jupyter Notebook text editor will automatically save this modified clm_params.txt file. For sure, you may click open menu File and chose save text to save your change(s). Screen Shot 2022-08-30 at 3 01 23 PM

AND, this modified clm_params.txt will be converted back to overwrite clm_params.nc when running OLMT example script in STEP 2 next.



STEP 2. Re-Run the model of 1 year transient simulation using OLMT

  • In one of Terminal (for running ELM docker container), Type the following command, OR, copy/paste it:
docker run -t -i --hostname=docker --user=modeluser -v elmdata:/inputdata/ -v elmoutput:/output/ -v elmwork:/tmp/ fasstsimulation/elm-builds:elm_v2-for-ngee /bin/bash /tools/OLMT/examples/site_fullrun_docker_transient_oneyear.sh

OR, by pressing UP-ARROW key one time (to get the same command as used before), as following:

Screen Shot 2022-08-30 at 3 13 23 PM

STEP 3. Visualizing new ELM output variable, VCMAX25TOP, and/or others.

  • AFTER STEP 2 (re-run) finished, Click-open UPPER Tab Plot_ELM_variable.ipynb file again. And move starting cell to [1].

Screen Shot 2022-08-30 at 3 21 30 PM

  • THEN, for an example, if click that restart the kernel & run all cells button, image , it will run script in all-steps at once, and will plot **GPP** by default, for individual PFT(s) if exist(s).

  • Scrolling up to move your mouse pointer to cell [4], like


![Screen Shot 2022-08-30 at 11 25 58 AM](https://user-images.githubusercontent.com/1952960/187477648-d77c02fe-60ee-4cf3-bcaf-7f6dcbfa7cd5.png)
  • By clicking the dropdown arrow, it will show whole list of variables from ELM outputs, and Choose one, such as VCMAX25TOP, Canopy Top Vcmax at 25 degree C:

Screen Shot 2022-08-30 at 11 39 38 AM


  • AND, then Clicking one-step advance button, image in the upper bar, until end of script,

OR, Clicking open menu Run, and clicking Run Selected Cell and All Below:

Screen Shot 2022-08-30 at 3 42 08 PM

THEN, it will plot like following:

Screen Shot 2022-08-30 at 3 37 59 PM

YOU may notice what have changed, comparing to default run shown in the previous Demo.



REPEAT STEP 1 - 3, to obtain reasonably well VCMAX for this 13th (c3_arctic_grass) PFT. THEN changing other concerned PFT(s).

FINALLY, save modified clm_params.nc into your file systems, like following (back to Docker App) Screen Shot 2022-08-30 at 11 05 45 PM