Run a coupled biogeochemistry (MARBL) & ocean (MOM6) model

1. MARBL is part of the CESM

2. Add MARBL to the compset

3. Add relevant information to configure_forcingss

4. See below for an example.

5. Adding runoff to a BGC example and getting river nutrients requires an offline mapping. Add the GLOFAS runoff mesh as well as the river nutrients (which is on the GLOFAS grid)

from CrocoDash.case import Case
case = Case(
    cesmroot = cesmroot,
    caseroot = caseroot,
    inputdir = inputdir,
    ocn_grid = grid,
    ocn_vgrid = vgrid,
    ocn_topo = topo,
    project = 'NCGD0011',
    override = True,
    machine = "derecho",
    compset = "CR1850MARBL_JRA"
)

case.configure_forcings(
    date_range = ["2000-01-01 00:00:00", "2000-02-01 00:00:00"],
    data_input_path = "", # Path to CESM OUTPUT data
    product_name = "CESM_OUTPUT",
    marbl_ic_filepath = "" # Path to MARBL Global IC File
)

# For River Nutrients, which should be used when runoff is also on, add
#    runoff_esmf_mesh_filepath = "", # For river nutrient mapping 
#    global_river_nutrients_filepath = "", # Path to global river nutrients mapping
 

Note

If you’re working on Derecho or Casper, you can use the following paths:

data_input_path = "/glade/campaign/collections/cmip/CMIP6/CESM-HR/FOSI_BGC/HR/g.e22.TL319_t13.G1850ECOIAF_JRA_HR.4p2z.001/ocn/proc/tseries/month_1"
product_name = "CESM_OUTPUT"
runoff_esmf_mesh_filepath = "/glade/campaign/cesm/cesmdata/cseg/inputdata/ocn/mom/croc/rof/glofas/dis24/GLOFAS_esmf_mesh_v4.nc"
global_river_nutrients_filepath = "/glade/campaign/cesm/cesmdata/cseg/inputdata/ocn/mom/croc/rof/river_nutrients/river_nutrients.GNEWS_GNM.glofas.20250916.64bit.nc"
marbl_ic_filepath = "/glade/campaign/collections/gdex/data/d651077/cesmdata/inputdata/ocn/mom/tx0.66v1/ecosys_jan_IC_omip_latlon_1x1_180W_c231221.nc"