-
Notifications
You must be signed in to change notification settings - Fork 3
AK K64G Demo Quick Start Guide 2a: Modifying the vegetation parameter Specific Leaf Area
This example 2B runs a single-point, extracted from half-degree Global data, ELM v2 simulation, within which NGEE Kougarok Mile 64 study site located, Using user-modified plant physiological parameter, slatop
, Specific Leaf Area (SLA) at top canopy (unit: m2/gC).
(1) It will start from an initial restart file from after-final-spinup in regular run of a site, but with option --site_name=kougarok
(2) It will modify clm_params.nc
, and Customize ELM outputs in user_nl_elm
, in a temporary work directory.
(3) Re-run ELM, using newly-generated clm_params.nc
and user_nl_elm
.
(4) Check a couple of ELM output variables, e.g. TLAI
, LEAFC
, etc., in jupyter notebook.
NOTE, assuming that you've already done with AK-K64G Demo Quick Start Guide 2. one year simulation of transient run.
AND, 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 following:
TWO Terminal windows:
![Screen Shot 2022-08-30 at 12 29 14 PM](https://user-images.githubusercontent.com/1952960/187508089-79e8d70d-03ee-4078-94ed-8a30b8250226.png)
ELM visualizing Jupyter notebook in internet brower:
![Screen Shot 2022-10-10 at 2 41 05 PM](https://user-images.githubusercontent.com/1952960/194932772-eca2aa88-1a77-49e8-98f0-b20374c053cf.png)
STEP 1. Modifying ELM parameter, slatop
, Specific Leaf Area at top canopy.
- If NOT YET, Click Jupyter Notebook for visualizing's Left Panel Folder brower, and Click open folder
work/kougarok
, and one of files for ELM parameters,clm_params.txt
, like following:
![Screen Shot 2022-09-05 at 4 25 25 PM](https://user-images.githubusercontent.com/1952960/188509912-8cf5506a-da5f-4b80-b47b-754202ae08cf.png)
NOTE that, in this Demo, we will aslo show how to customize user_nl_elm
, for output variable LEAFN
, leaf biomass N, which originally not meant to output in AK-K64G Demo Quick Start Guide 2. one year simulation of transient run.
- Scrolling down the file
clm_params.txt
toLine 170-172
(shown in shaded), where are the header information of parameterslatop
.
- Further Scrolling down to
Line 1327-1329
(shown in shaded are the 12th and 13th values for what we're interested in this Demo).
- For an example, modify default value of
0.03
of the13th (c3_arctic_grass)
, to0.0165
(m2/gC leaf), i.e. thicker leaf per biomass.
NOTE: this value here is totally assumed, implying USER should obtain its real value from field/lab works or literature studies.
![Screen Shot 2022-09-05 at 4 31 02 PM](https://user-images.githubusercontent.com/1952960/188510067-69a55e60-897a-4bfb-ae41-920e30751165.png)
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-10-10 at 3 03 06 PM](https://user-images.githubusercontent.com/1952960/194935902-b7c50c97-d6cf-41c6-acd5-5e432c2cdcbd.png)
AND, this modified clm_params.txt
will be converted back to overwrite clm_params.nc
when running OLMT example script in STEP 3 of this Demo.
STEP 2. (OPTIONALLY) Add ELM output of LEAFN
and other vegetation N varibles into user_nl_elm
, which will control ELM's actual output variable list.
- Click-open file
user_nl_elm
in folder /work/kougarok, which had already customized but NO vegetation N variables output originally, as shown:
![Screen Shot 2022-10-10 at 2 53 02 PM](https://user-images.githubusercontent.com/1952960/194934395-b4af5fc5-8aa1-4d88-a056-113c94dbb374.png)
- Add
TOTVEGN
,LEAFN
, etc. (separated by comma) into the list. NOTE again the file will be automactically saved, but you may manully save it for sure.
![Screen Shot 2022-10-10 at 2 50 17 PM](https://user-images.githubusercontent.com/1952960/194934416-1a9a872f-49cc-4179-9056-45a6519232fb.png)
STEP 3. Re-Run the model of 1 year transient simulation
using OLMT
- In the Terminal for running ELM docker container, Type the following command, OR, copy/paste it. NOTE the command line is appended one more option
--case_prefix=MySLA
so that this demo won't erase AK-K64G Demo Quick Start Guide 2. one year simulation of transient run:
docker run -t -i --hostname=docker --user=modeluser -v elmdata:/inputdata/ -v elmoutput:/output/ -v elmwork:/tmp/ fasstsimulation/elm-builds:elm_v2-for-ngee_multiarch /bin/bash /tools/OLMT/examples/site_fullrun_docker_transient_nyears.sh --transient_years=1 --case_prefix=MySLA
STEP 4. Visualizing newly-simulated ELM output variable, TLAI
, LEAFC
, LEAFN
, and/or others.
- AFTER STEP 3 (re-run) finished, Click-open UPPER Tab
Plot_ELM_variable_v2.ipynb
file again.
![Screen Shot 2022-10-10 at 3 29 29 PM](https://user-images.githubusercontent.com/1952960/194939526-13611c16-99ac-4d9a-987b-293d98f4bb3c.png)
-
THEN, click that
restart the kernel & run all cells
button,, it will run script in all-steps at once, and will plot
**GPP**
by default. NOTE this rerun script is Needed to refresh simulations just finished. -
Scrolling up to move your mouse pointer to cell
[3]
, and clicking the dropdown arrow, choose case nameMySLA_...
:
![Screen Shot 2022-10-10 at 3 15 10 PM](https://user-images.githubusercontent.com/1952960/194937743-1063c20a-bfd8-4c00-bba9-548af20b82e0.png)
- Click the
one-step advance
button (OR CTRL+Enter) and move cursor down to cell[5]
, choseIndividual PFTs
:
![Screen Shot 2022-10-10 at 3 15 26 PM](https://user-images.githubusercontent.com/1952960/194938071-39a46ef5-dd98-4d4f-80c8-35a7606238f6.png)
-
- Click the
one-step advance
button (OR CTRL+Enter) and move cursor down further to Choose Variable Name dropdow, and clicking the dropdown arrow, it will show whole list of variables from ELM outputs, and Choose one, such asTLAI
, total leaf area index:
- Click the
![Screen Shot 2022-10-10 at 3 13 53 PM](https://user-images.githubusercontent.com/1952960/194938480-3ac40dc6-4951-4bed-bbd6-675bef76d5bc.png)
- AND, then Clicking open menu
Run
, and clickingRun Selected Cell and All Below
:
![Screen Shot 2022-09-05 at 5 11 55 PM](https://user-images.githubusercontent.com/1952960/188512893-dfc727c1-a4bc-41de-a0d2-187449631c34.png)
THEN, it will plot like following:
![Screen Shot 2022-10-10 at 3 17 30 PM](https://user-images.githubusercontent.com/1952960/194938698-8afac4a6-94e1-4bc6-b34a-63a0677579e6.png)
FOR LEAFC
, we can do similarly to show as following:
![Screen Shot 2022-10-10 at 3 26 34 PM](https://user-images.githubusercontent.com/1952960/194938949-9390f93e-c2cb-47cb-9f74-7908c49b6400.png)
REPEAT STEP 1 - 4, to obtain reasonably well LAI and leaf biomass 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)