b.e23_alpha17f.BLT1850.ne30_t232.091

[gustavo-marques]

Purpose:

Evaluate new out-of-the-box baseline with cesm2_3_alpha17f (or cesm2_3_beta17).

---

Description:

New baseline with out-of-the-box settings from cesm2_3_alpha17f (or cesm2_3_beta17). The vertical coordinate in MOM6 is z* with 65 layers. The ICs for MOM and CICE are from corresponding G case NCAR/omwg_dev#48; restart file 0016-01-01.

---

Case directory:

• Locally (if still available):
  /glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/b.e23_alpha17f.BLT1850.ne30_t232.091

---

Sandbox:

• Locally (if still available): /glade/work/gmarques/cesm.sandboxes/cesm2_3_alpha17f

• On github: https://github.com/ESCOMP/cesm/tree/cesm2_3_alpha17f

---

Diagnostics:

• Diags (if available):
  https://webext.cgd.ucar.edu/BLT1850/b.e23_alpha17f.BLT1850.ne30_t232.091/

---

Output:

• Raw output on derecho (if still available): /glade/derecho/scratch/gmarques/archive/b.e23_alpha17f.BLT1850.ne30_t232.091/

---

Contacts:
@gustavo-marques, @cecilehannay, @PeterHjortLauritzen, @adamrher, @brianpm, @JulioTBacmeister, @swrneale, @fobryan3, @bstephens82, @klindsay28, @dabail10, @wwieder, @slevis-lmwg

---

Extra details:

In user_nl_cice:

ice_ic = '/glade/p/cesmdata/cseg/inputdata/cesm2_init/g.e23_a17f.GJRAv4.TL319_t232_zstar_N65.2024.056/0016-01-01-00000/g.e23_a17f.GJRAv4.TL319_t232_zstar_N65.2024.056.cice.r.0016-01-01-00000.nc'
dumpfreq = 'x','x','x','x','x'
f_fresh = 'mxxxx'
f_fresh_ai = 'mxxxx'

In user_nl_mom:

MIN_SALINITY = 1.e-06
INIT_LAYERS_FROM_Z_FILE = False
THICKNESS_CONFIG = "thickness_file"
THICKNESS_FILE = "/glade/p/cesmdata/cseg/inputdata/cesm2_init/g.e23_a17f.GJRAv4.TL319_t232_zstar_N65.2024.056/0016-01-01-00000/g.e23_a17f.GJRAv4.TL319_t232_zstar_N65.2024.056.mom6.r.0016-01-01-00000.nc"
TS_CONFIG = "file"
TS_FILE = "/glade/p/cesmdata/cseg/inputdata/cesm2_init/g.e23_a17f.GJRAv4.TL319_t232_zstar_N65.2024.056/0016-01-01-00000/g.e23_a17f.GJRAv4.TL319_t232_zstar_N65.2024.056.mom6.r.0016-01-01-00000.nc"
TEMP_IC_VAR = "Temp"
SALT_IC_VAR = "Salt"
VELOCITY_CONFIG = "file"
VELOCITY_FILE = "/glade/p/cesmdata/cseg/inputdata/cesm2_init/g.e23_a17f.GJRAv4.TL319_t232_zstar_N65.2024.056/0016-01-01-00000/g.e23_a17f.GJRAv4.TL319_t232_zstar_N65.2024.056.mom6.r.0016-01-01-00000.nc"

[adamrher]

@gustavo-marques since the atm/hist/ output is on the unstructured mesh, would you mind executing this regridding script to map all the .cam.h0a. files to a lat-lon grid? Execute it in the atm/hist/ directory and it will create a sub-dir atm/hist/regrid/ with the same filenames but all mapped to f09:

#!/bin/bash

case=b.e23_alpha17f.BLT1850.ne30_t232.091
wgts=/glade/work/aherring/CESM2/TSTEPS/tseries/ne30pg3_TO_f09-bilin.nc

mkdir regrid

for file in ${case}.cam.h0a.*; do
  echo ${file}
  ncremap -m ${wgts} ${file} regrid/${file}
done

exit

[gustavo-marques]

@adamrher: sure! I am running it now, but remapping all the files will take some time.

[wwieder]

Can you do this for all the land files too?  That would make our post processing much easier. Now that I ask this, I wonder if we need to be a bit more nuanced in how we create the wgts file, to handle land/ocean masks more intelligently? @olyson how have you been remapping files for other cases?

Also, @olyson, maybe @slevis-lmwg and I can take a look at this case with you to see if the land component is set up as expected?

[olyson]

When required, for individual history files, I have been using mkprocdata, which uses conservative regridding and is pretty fast really. Currently, I've been using this command to create the mapping file required:

qcmd -- ./create_ESMF_map.sh --maptype aave --filesrc /glade/p/cesmdata/cseg/inputdata/atm/cam/coords/ne30pg3_scrip_170604.nc --filedst /glade/p/cesm/amwg/amwg_diagnostics/map_files/0.9x1.25.nc --namesrc ne30pg3 --namedst fv0.9x1.25

I think ncremap --sgs_frc=landfrac -m /glade/work/aherring/CESM2/TSTEPS/tseries/ne30pg3_TO_f09-bilin.nc ... will run on the land history files as well, not sure how long it takes.
It would be a good idea to compare the two methods. I'm not sure how well the "bilin" works on our fields that are not global.
Within the land diagnostics package, it computes the climatological, means, and trends files on the SE grid, then regrids them using the ESMF regridding capability within NCL. Then plots are produced.

[adamrher]

@cecilehannay can we try to run ADF diagnostics on the remapped output in the /atm/hist/regrid/ directory?

[olyson]

Land diagnostics are here:

https://webext.cgd.ucar.edu/BLT1850/b.e23_alpha17f.BLT1850.ne30_t232.091/lnd/b.e23_alpha17f.BLT1850.ne30_t232.091_15_24-b.e23_alpha16g.BLT1850.ne30_t232.082b_15_24/setsIndex.html

[gustavo-marques]

@wwieder and @olyson, I am happy to also regrid/remap the land output. What script should I run?

[wwieder]

Thanks @gustavo-marques the LMWG diagnostics package can work on the raw ne30 grid and then regrids for contour plots.

[megandevlan]

@gustavo-marques Would it be possible to extend this run to 100 years? It's a really nice clean comparison between 091, 092, 093, and 098 for assessing the impacts of each change relative to this control baseline! I think it's 34 years now, but for ENSO longer would be better.

[gustavo-marques]

Hi @megandevlan! This case has already been run for 100 years.

[megandevlan]

Ugh, got my cases mixed up - thank you for checking this @gustavo-marques!!