Modify the rain_threshold in CLM

Modify the rain_threshold in CLM#

In the below exercise, we will change the rain threshold for stress deciduous vegetation, which includes C3 grasses. The rain threshold is the amount of rain required to initiate leaf onset. Reaching the rain threshold is one of several requirements for stress deciduous vegetation leaf onset. If you are interested, you can find more information about the stress deciduous phenology representation in the CLM Technical Note.

The current value of rain_threshold is 20mm as specified in the code:

rain_threshold = 20._r8
Exercise: Modify the CLM rain_threshold parameter

Create a case called b1850_rain_threshold using the compset B1850 at f19_g17 resolution.

  • Change the rain_threshold from 20mm to 1mm.

  • Make a 1-month run.

Click here for hints

Tip to locate the ‘rain_threshold’ parameter

  • Find the subroutine containing rain_threshold. For instance, in the CESMROOT, use the command:

    grep –r rain_threshold *
Click here for the solution

Create a new case b1850_rain_threshold with the command:

cd /glade/u/home/$USER/code/my_cesm_code/cime/scripts
./create_newcase --case ~/cases/b1850_rain_threshold --compset B1850 --res f19_g17

Case setup:

cd ~/cases/b1850_rain_threshold
./case.setup

Locate the file where rain_threshold is defined and copy it into SourceMods/src.lnd:

cp /glade/u/home/$USER/code/my_cesm_code/components/clm/src/biogeochem/CNPhenologyMod.F90 SourceMods/src.clm

Edit the file SourceMods/src.clm/CNPhenologyMod.F90 and add the lines:

      ! specify rain threshold for leaf onset                                   
      rain_threshold = 20._r8

to

      ! specify rain threshold for leaf onset                                   
      rain_threshold = 1._r8

Change the run length:

./xmlchange STOP_N=1,STOP_OPTION=nmonths

If needed, change job queue and account number For instance:

./xmlchange JOB_QUEUE=tutorial,PROJECT=UESM0013 --force

Build and submit:

qcmd -- ./case.build
./case.submit