4: Turn on the ecosystem

You can also explore setting up a similar case but using the GW1850MARBL_JRA component set. Note how this differs from the previous GW_JRA component set we used in Exercise 1.

Exercise: Run a control case

Create a case called g.GW1850MARBL_JRA.TL319_t233_wg37.001 using the compset GW1850MARBL_JRA at TL319_t233_wg37 resolution.

Set the run length to 1 year.

Note that the BGC support is not fully incorporated yet, which means you cannot run it “out of the box” but have to make the following modifications:

RIV_FLUX_FILE = "riv_nut.gnews_gnm.rJRA025_to_tx2_3v3_nnsm_e100r100_260306.20260424.nc"
MARBL_FESEDFLUX_FILE = "fesedflux_2024algo_tx2_3v3.c260608.nc"
MARBL_FESEDFLUXRED_FILE = "fesedfluxRed_2024algo_tx2_3v3.c260608.nc"
MARBL_FEVENTFLUX_FILE = "feventflux_2026algo_5gmol_tx2_3v3.c260612.nc"

Refer to Exercise 2 for how to make these modifications, and then build and run the model.

Click here for hints

How do I modify the runtime parameters?

In user_nl_mom add the following lines:

RIV_FLUX_FILE = "riv_nut.gnews_gnm.rJRA025_to_tx2_3v3_nnsm_e100r100_260306.20260424.nc"
MARBL_FESEDFLUX_FILE = "fesedflux_2024algo_tx2_3v3.c260608.nc"
MARBL_FESEDFLUXRED_FILE = "fesedfluxRed_2024algo_tx2_3v3.c260608.nc"
MARBL_FEVENTFLUX_FILE = "feventflux_2026algo_5gmol_tx2_3v3.c260612.nc"

How do I compile and run?

You can compile with the command:

qcmd -- ./case.build

You can run with the command:

./case.submit

How do I check the lenght of the run?

Use xmlquery to search for the variables that control the run length

Click here for the solution

Create a new case g.GW1850MARBL_JRA.TL319_t233_wg37.001 with the command:

cd /glade/u/home/$USER/code/my_cesm_code/cime/scripts/
./create_newcase --case /glade/u/home/$USER/cases/g.GW1850MARBL_JRA.TL319_t233_wg37.002 --compset GW1850MARBL_JRA --res TL319_t233_wg37 --run-unsupported

Case setup:

cd /glade/u/home/$USER/cases/g.GW1850MARBL_JRA.TL319_t233_wg37.001 
./case.setup

Change the run length:

./xmlchange STOP_N=1,STOP_OPTION=nyears

Edit the file user_nl_mom and add the lines:

RIV_FLUX_FILE = "riv_nut.gnews_gnm.rJRA025_to_tx2_3v3_nnsm_e100r100_260306.20260424.nc"
MARBL_FESEDFLUX_FILE = "fesedflux_2024algo_tx2_3v3.c260608.nc"
MARBL_FESEDFLUXRED_FILE = "fesedfluxRed_2024algo_tx2_3v3.c260608.nc"
MARBL_FEVENTFLUX_FILE = "feventflux_2026algo_5gmol_tx2_3v3.c260612.nc"

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

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

Build and submit:

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

When the run is completed, look into the archive directory for: g.GW1850MARBL_JRA.TL319_t233_wg37.001.

cd /glade/derecho/scratch/$USER/archive/g.GW1850MARBL_JRA.TL319_t233_wg37.001/ocn/hist
ls