Kitchen Sink - Oh My! It’s Everything!

This project is about setting up everything we can and seeing what happens!

Download this project by either:

Running the CrocoDash CLI command: crocodash template --machine derecho --notebook crocodash.projects.KitchenSink
Copying the file from your CrocoDash checkout: gallery/crocodash/projects/KitchenSink.ipynb

Grid Generation¶

Let’s make a small domain in a northerly area

from pathlib import Path

from CrocoDash.grid import Grid
from CrocoDash.vgrid import VGrid
from CrocoDash.topo import Topo
from CrocoDash.case import Case


grid = Grid.from_projection(
    crs="EPSG:3338",        # Alaska Albers (std parallels 55N/65N, CM 154W)
    x_min=-1_000_000,       # meters, Albers origin at (154W, 50N)
    x_max= 1_200_000,
    y_min=  400_000,
    y_max= 1_600_000,
    resolution_m=10_000,    # 10 km
    name="alaska_10km",
)


topo = Topo(
    grid = grid,
    min_depth = 9.5, # in meters
)

topo.set_from_dataset(    
    bathymetry_path = "<GEBCO_LOWRES>",
    longitude_coordinate_name="lon",
    latitude_coordinate_name="lat",
    vertical_coordinate_name="elevation"
)

topo.depth.plot()
vgrid  = VGrid.hyperbolic(
    nk = 75, # number of vertical levels
    depth = topo.max_depth,
    ratio=20.0 # target ratio of top to bottom layer thicknesses
)

# CESM case (experiment) name
casename = "ks_alaska_10km"

# CESM source root (Update this path accordingly!!!)
cesmroot ="<CESM>"

# Place where all your input files go 
inputdir = Path("<inputdir>") / casename
    
# CESM case directory
caseroot = Path("<casedir>") / casename


case = Case(
    cesmroot = cesmroot,
    caseroot = caseroot,
    inputdir = inputdir,
    ocn_grid = grid,
    ocn_vgrid = vgrid,
    ocn_topo = topo,
    project = '<PROJECT>',
    override = True,
    machine = "derecho",
    compset = "1850_DATM%JRA_SLND_CICE_MOM6%REGIONAL%MARBL-BIO_DROF%GLOFAS_SGLC_WW3" )

Forcings¶

DROF is the only new component here. All we need to do is provide a mesh file to get the runoff going. All CrocoDash does is generate a mapping file for the runoff data. The CESM coupler actually passes the data in (which we specify through the compset).

case.configure_forcings(
    date_range=["2020-01-01 00:00:00", "2020-01-09 00:00:00"],
    product_name="reference_ocean",
    function_name="get_reference_ocean_data",
    marbl_ic_filepath="/glade/campaign/collections/gdex/data/d651077/cesmdata/inputdata/ocn/mom/tx0.66v1/ecosys_jan_IC_omip_latlon_1x1_180W_c231221.nc",   
    cice_product_name="reference_ice",            
    cice_function_name="get_reference_ice_data", 
    n_halo_cells=1, 
    ww3_obc_product_name="reference_waves",              
    ww3_obc_function_name="get_reference_wave_spectra",
    rof_esmf_mesh_filepath="/glade/campaign/cesm/cesmdata/cseg/inputdata/ocn/mom/croc/rof/glofas/dis24/GLOFAS_esmf_mesh_v4.nc",                       
)

case.process_forcings()

Let’s see if this works! If you’re interested in this setup, you can swap in real data for each product as you go along! It’ll take some serious time to run!

With this many components active, I’d mess with the NTASKS/ROOTPE layout for each one to make sure everything looks reasonable, then a quick build & submit.

Feel free to iterate with process forcings! You can run it on the command line with crocodash process --caseroot YOURCASEROOT!