b.e30_alpha08b.B1850C_LTso.ne30_t232_wgx3.323

Open
#307 6 comments 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

Assessment

Difficulty
5/5
Estimated time
Over a week
Newbie friendliness
25/100
Issue type
Feature
Clarity
Mostly clear
Activity status
Quiet
Tech stack
fortran
Domain
backend

Research direction

Start with the b.e30_alpha08b.B1850C_LTso.ne30_t232_wgx3.323 case directory and compare issue #304, then review the user_nl_cam, user_nl_cice, user_nl_mom, and user_nl_cpl settings and required INPUT files. Run the case from the specified CESM sandbox and check the linked diagnostics and archive output. Done means the case runs with the requested startup, grid, runoff, volcano, and diagnostic configuration.

Written by the indexing model from the issue text.

Description

Purpose:
Repeat 320

  • update ocean configuration (tx2_3v3),
  • frozen runoff mapping file is a merge of two mapping files: 100 km spreading for the SH and 250 km spreading for the NH.
  • Start from a G case at year 8.
  • Add volcano bugfix as in 319.

Description:
Same as 320 (#304) except:

  • change run type
./xmlchange RUN_TYPE=startup
  • change runoff
./xmlchange ROF2OCN_ICE_RMAPNAME=/glade/work/gmarques/cesm/tx2_3/runoff_mapping/map_r05_to_tx2_3_nnsm_e100r100sh_e250r250nh_merged_260307.nc
./xmlchange ROF2OCN_LIQ_RMAPNAME=/glade/work/gmarques/cesm/tx2_3/runoff_mapping/map_r05_to_tx2_3_nnsm_e100r100_260306.nc
  • change grid
./xmlchange ICE_DOMAIN_MESH=/glade/work/gmarques/cesm/tx2_3/mesh/ESMF_mesh_tx2_3v3_260305_cdf5.nc
./xmlchange OCN_DOMAIN_MESH=/glade/work/gmarques/cesm/tx2_3/mesh/ESMF_mesh_tx2_3v3_260305_cdf5.nc
./xmlchange MASK_MESH=/glade/work/gmarques/cesm/tx2_3/mesh/ESMF_mesh_tx2_3v3_260305_cdf5.nc
  • add sourcemods for MAM5 volcanos:
/glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/b.e30_alpha08b.B1850C_LTso.ne30_t232_wgx3.321/SourceMods/src.cam/
aerosol_optics_cam.F90  
phys_prop.F90
  • add to user_nl_cam corrected rad_climate and extra variables to look at MAM5 volcanos
add to fincl1: 'CME', 'CMEIOUT', 'AODSSdn','CME', 'CMEIOUT', 'AODSSdn', 'EXTINCTdn','VOLC_MMR1','VOLC_MMR2','VOLC_MMR3','VOLC_MMR5','VOLC_SAD',
'VOLC_MASS5','BURDENSO4dn'

rad_climate = 'A:Q:H2O', 'N:O2:O2', 'A:CO2:CO2', 'N:ozone:O3', 'A:N2O:N2O', 'A:CH4:CH4', 'N:CFC11STAR:CFC11', 'A:CFC12:CFC12',          'M:mam4_mode1:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/mam4_mode1_rrtmg_aeronetdust_sig1.6_dgnh.48_c140304.nc',           'M:mam4_mode2:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/mam4_mode2_rrtmg_aitkendust_c141106.nc',          'M:mam4_mode3:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/mam4_mode3_rrtmg_aeronetdust_c141106.nc',         'M:mam4_mode4:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/mam4_mode4_rrtmg_c130628.nc',           'N:VOLC_MMR1:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/volc_camRRTMG_byradius_sigma1.6_mode1_c210211.nc',           'N:VOLC_MMR2:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/volc_camRRTMG_byradius_sigma1.6_mode2_c210211.nc',           'N:VOLC_MMR3:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/volc_camRRTMG_byradius_sigma1.2_mode3_c210211.nc',          'N:VOLC_MMR5:/glade/campaign/cesm/cesmdata/inputdata/atm/cam/physprops/volc_camRRTMG_byradius_sigma1.2_mode5_c260228.nc'
  • user_nl_cice
grid_file = "/glade/work/gmarques/cesm/tx2_3/seaice_files/cice_grid_tx2_3v3_260305_cdf5.nc"
grid_format = 'nc'
ice_ic      = "/glade/derecho/scratch/gmarques/archive/g.e30_a08g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2026.005/rest/0008-01-01-00000/g.e30_a08g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2026.005.cice.r.0008-01-01-00000.nc"
kmt_file = "/glade/work/gmarques/cesm/tx2_3/seaice_files/cice_grid_tx2_3v3_260305_cdf5.nc"

tr_fsd = .true.
restart_fsd = .true. 
nfsd = 12
wave_spec_type = 'alt'
histfreq = "m", "x", "x", "x", "x"
histfreq_n = 1, 1, 0, 0, 0

f_fresh = 'mxxxx'
f_fresh_ai = 'mxxxx'
f_frzmlt = 'mdxxx'
f_afsd = "mdxxx"
f_afsdn = "mdxxx"
f_dafsd_latg = "mdxxx"
f_dafsd_latm = "mdxxx"
f_dafsd_newi = "mdxxx"
f_dafsd_wave = "mdxxx"
f_dafsd_weld = "mdxxx"
f_wave_sig_ht = "mdxxx"
f_diam_ww = "mdxxx"
f_hice_ww = "mdxxx"
f_fsdrad = "mdxxx"
f_fsdperim = "mdxxx"
f_uvel = "mdxxx"
f_vvel = "mdxxx"
r_snw = 1.0
  • user_nl_mom
INPUTDIR = "INPUT"
GRID_FILE = "ocean_hgrid_250930.nc"
TOPO_FILE = "ocean_topo_tx2_3v3_260305.nc"
TOPO_EDITS_FILE = ""
TIDAL_ENERGY_FILE = "tidal_energy_tx2_3v3_conserve_260306.nc"
CHL_FILE = "seawifs-clim-1997-2010-tx2_3v3_260305.nc"
CHANNEL_LIST_FILE = "channels_tx2_3v3_260304.txt"
GEOTHERMAL_FILE = "geothermal_davies2013_tx2_3v3_260305.nc"
MLE%
CR_FILE = "mle_cr_labsea_mask_tx2_3v3_260306.nc"
%MLE

MEKE_GEOMETRIC_ALPHA = 0.07
SMAGORINSKY_KH = False
MAXTRUNC = 1.0E6

! Leith+E Settings
TAPER_LEITHY = True
LEITH_BI_CONST = 85.0
LEITHY_CK = 1.0
LEITHY_DEPTH = 1000.0
LEITHY_WIDTH = 500.0

! Background KH
USE_KH_BG_2D = True
KH_BG_2D_FILENAME = "KH_BG_2D_tx2_3v3_260306.nc"
KH_BG_2D_VARNAME = "KH"

! Resolution scaling:
!! GM & Redi the same
!! Don't apply scaling to Laplacian viscosity, which is the 2D background
RESOLN_SCALED_KH = False
RESOLN_SCALED_KHTH = True  ! Already using, here for clarity
RESOLN_SCALED_KHTR = True  ! Already using, here for clarity
KH_RES_SCALE_COEF = 0.5    ! Tighter around the equator
KH_RES_FN_POWER = 6        ! Already using, here for clarity
VISC_RES_SCALE_COEF = 0.5  ! Used for KHTR
VISC_RES_FN_POWER = 6      ! Used for KHTR
INTERPOLATE_RES_FN = True  ! Only computes resolution fn once, then interpolates

KHTR_MIN = 50.

THICKNESS_CONFIG = thickness_file
VELOCITY_CONFIG = file
TS_CONFIG = file

INIT_LAYERS_FROM_Z_FILE = False

THICKNESS_FILE = "/glade/derecho/scratch/gmarques/archive/g.e30_a08g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2026.005/rest/0008-01-01-00000/g.e30_a08g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2026.005.mom6.r.0008-01-01-00000.nc"

TS_FILE = "/glade/derecho/scratch/gmarques/archive/g.e30_a08g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2026.005/rest/0008-01-01-00000/g.e30_a08g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2026.005.mom6.r.0008-01-01-00000.nc"
VELOCITY_FILE = "/glade/derecho/scratch/gmarques/archive/g.e30_a08g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2026.005/rest/0008-01-01-00000/g.e30_a08g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2026.005.mom6.r.0008-01-01-00000.nc"

TEMP_IC_VAR = "Temp"
SALT_IC_VAR = "Salt"
  • user_nl_cpl
histaux_l2x1yrg = .true.
histaux_atm2med_file1_enabled = .true.
histaux_atm2med_file2_enabled = .true.
histaux_atm2med_file3_enabled = .true.
histaux_atm2med_file4_enabled = .true.
histaux_atm2med_file5_enabled = .true.
histaux_atm2med_file5_history_n = 1
histaux_atm2med_file5_history_option = 'ndays'
histaux_atm2med_file5_ntperfile = 1
  • create INPUT with all required files
mkdir /glade/derecho/scratch/hannay/b.e30_alpha08b.B1850C_LTso.ne30_t232_wgx3.321/run/INPUT 
cd /glade/derecho/scratch/hannay/b.e30_alpha08b.B1850C_LTso.ne30_t232_wgx3.321/run/INPUT 
cp -L /glade/derecho/scratch/gmarques/g.e30_a08g.GW_JRA.TL319_t232_wgx3_hycom1_N75.2026.005/run/INPUT/* .

Case directory:

Sandbox:

Diagnostics:

Output:

  • Raw output on derecho (if still available): /glade/derecho/scratch/hannay/archive/b.e30_alpha08b.B1850C_LTso.ne30_t232_wgx3.323/
  • CAM Climos: /glade/campaign/cgd/amp/amwg/climo/b.e30_alpha08b.B1850C_LTso.ne30_t232_wgx3.323

Contacts:
@cecilehannay @dlawrenncar @adamrher @PeterHjortLauritzen @JulioTBacmeister @swrneale @slevis-lmwg @wwieder @olyson @gustavo-marques @iangrooms @fobryan3 @klindsay28 @dabail10 @duvivier @tilmes @mbramberger @gunterl

Extra details:

Dominant language
No language data
Stars
11
Forks
3
No merged PRs in 30d

Contributor guide

No contributing guide indexed for this repository

First steps

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. Open a pull request that references the issue number.

More from NCAR/cesm_dev

All issues in NCAR/cesm_dev

Similar issues

More Backend & API Design issues

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.