b.e30_alpha09b.B1850C_MTso_Gris.ne30_t233_wgx3.367

Open
#372 1 comment 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, shell
Domain
build-system, hpc

Research direction

Start with ccs_config/modelgrid_aliases_nuopc.xml and the supplied user_nl_mom settings, then review the listed SourceMods and create_newcase commands. Build the case using the ne30pg3_t233_wg37_gris4 alias, initialize from 360a at year 139, and inspect the linked diagnostics to assess the 3000 m taper experiment.

Written by the indexing model from the issue text.

Description

Purpose:
Increase the taper depth to determine whether we can prevent Labrador Sea freezing while avoiding the reversal of the AMOC circulation within the Nordic Seas.

Same configuration as 360a, but with the taper depth increased to 3000 m.

Start from 360a at year 139.


Description:

Same as 360a but

  • taper depth 1500. -> 3000.
MLE_TAPER_FN_DEPTH = 1500. -> 3000. 
  • Start from 360a at year 139.
./xmlchange RUN_REFCASE=bg.e30_alpha09b.B1850C_MTso.ne30_t233_wgx3.360a
./xmlchange RUN_REFDATE=0139-01-01
./xmlchange RUN_TYPE=hybrid

**Compset
Compset: 1850C_CAM70%MT_CLM60%BGC-CROP_CICE_MOM6_MOSART_CISM2%GRIS-EVOLVE_WW3
res: ne30pg3_t233_wg37_gris4


Case directory:


Sandbox:

  • add this block to ccs_config/modelgrid_aliases_nuopc.xml (next to the
    other ne30pg3_t233 aliases):
  <model_grid alias="ne30pg3_t233_wg37_gris4">
    <grid name="atm">ne30np4.pg3</grid>
    <grid name="lnd">ne30np4.pg3</grid>
    <grid name="ocnice">tx2_3v3</grid>
    <grid name="glc">gris4</grid>
    <grid name="wav">wgx3v7</grid>
    <mask>tx2_3v3</mask>
  </model_grid>

Diagnostics:


Output:

  • Raw output on derecho (if still available): /glade/derecho/scratch/$USER/archive/b.e30_alpha09b.B1850C_MTso_Gris.ne30_t233_wgx3.367/
  • CAM Climos: /glade/campaign/cgd/amp/amwg/climo/b.e30_alpha09b.B1850C_MTso_Gris.ne30_t233_wgx3.367

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


Extra details:

  • create_newcase
./create_newcase --compset 1850C_CAM70%MT_CLM60%BGC-CROP_CICE_MOM6_MOSART_CISM2%GRIS-EVOLVE_WW3 --res ne30pg3_t233_wg37_gris4 \
--case /glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/b.e30_alpha09b.B1850C_MTso_Gris.ne30_t233_wgx3.367  \
--run-unsupported  
  • pe_layout
/bin/cp /glade/campaign/cesm/cesmdata/cseg/runs/cesm2_0/bg.e30_alpha09b.B1850C_MTso.ne30_t233_wgx3.365/env_mach_pes.xml .

xmlchange

./xmlchange CAM_CONFIG_OPTS=' -pcols 9' --append
./xmlchange CASE_GIT_REPOSITORY=git@github.com:NCAR/cesm_dev.git
./xmlchange RUN_REFCASE=bg.e30_alpha09b.B1850C_MTso.ne30_t233_wgx3.360a
./xmlchange RUN_REFDATE=0139-01-01
./xmlchange RUN_TYPE=hybrid
./xmlchange GET_REFCASE=true
./xmlchange RUN_REFDIR=cesm2_init
./xmlchange OCN_DIAG_MODE=spinup
./xmlchange JOB_PRIORITY=special --force
./xmlchange PROJECT=CESM0023,RESUBMIT=15,STOP_N=4,STOP_OPTION=nyears
./xmlchange PROJECT=CESM0023
./xmlchange JOB_WALLCLOCK_TIME=12:00:00  --subgroup case.run
  • user_nl_cam
mfilt               =       0,      0,     20,      40,      12,      120,      1,   1
nhtfrq              =       0,      0,    -24,      -3,       0,       -2,      0,  -8760
ndens               =       2,       2,      2,       2,       2,        1,      2,   1
interpolate_output  =  .true.,  .false., .true.,  .true., .false.,   .true.,     .true.
interpolate_nlat    =     192,     192,    192,     192,      192,     192,      192
interpolate_nlon    =     288,     288,    288,     288,      288,     288,      288

empty_htapes = .true.

fincl1 = 'ACTNI', 'ACTNL', 'ACTREI', 'ACTREL', 'AODDUST', 'AODDUSTdn','AODVIS', 'AODVISdn','BURDENBC',
'BURDENDUST', 'BURDENPOM', 'BURDENSEASALT',
'BURDENSO4', 'BURDENSOA', 'CAPE', 'CDNUMC',  'CLDHGH', 'CLDICE', 'CLDLIQ', 'CLDLOW',
'CLDMED', 'CLDTOT', 'CLOUD', 
'FCTI', 'FCTL', 'FLDS', 'FLNS', 'FLNSC', 'FLNT', 'FLNTC', 'FLUT',
'FREQZM', 'FSDS', 'FSDSC', 'FSNS', 'FSNSC', 'FSNT', 'FSNTC', 'FSNTOA', 'ICEFRAC', 'LANDFRAC', 'LHFLX', 'LWCF',
'OCNFRAC', 'OMEGA', 'OMEGA500', 'PBLH', 'PHIS', 'PINT', 'PMID', 'PRECC', 'PRECL', 'PRECSC', 'PRECSL', 'PRECT', 'PS', 'PSL',  'Q',
'QFLX', 'RELHUM', 'SHFLX', 'SOLIN', 'SST',
'SWCF',
'T', 'TAUX', 'TAUY', 'TFIX', 'TGCLDIWP', 'TGCLDLWP', 'TMQ', 'TREFHT', 'TS', 'U', 'U10',
'UBOT', 
'V', 'VBOT', 'Z3', 
'AODVISdn', 'CDNUMC',  'NUMICE', 'NUMLIQ','OMEGA500',
 'AQSO4_H2O2','AQSO4_O3', 
'ncl_a1SFWET',
'ncl_a2SFWET', 'ncl_a3SFWET', 'pom_a1SFWET', 'pom_a4SFWET', 'so4_a1SFWET', 'so4_a2SFWET', 'so4_a3SFWET', 'soa_a1SFWET',
'soa_a2SFWET', 'bc_c1SFWET', 'bc_c4SFWET', 'dst_c1SFWET', 'dst_c2SFWET', 'dst_c3SFWET', 'ncl_c1SFWET', 'ncl_c2SFWET',
'ncl_c3SFWET', 'pom_c1SFWET', 'pom_c4SFWET', 'so4_c1SFWET', 'so4_c2SFWET', 'so4_c3SFWET', 'soa_c1SFWET', 'soa_c2SFWET',
'bc_a1DDF', 'bc_a4DDF', 'dst_a1DDF', 'dst_a2DDF', 'dst_a3DDF', 'ncl_a1DDF', 'ncl_a2DDF', 'ncl_a3DDF',
'pom_a1DDF', 'pom_a4DDF', 'so4_a1DDF', 'so4_a2DDF', 'so4_a3DDF', 'soa_a1DDF', 'soa_a2DDF',
'so4_a1_CLXF', 'so4_a2_CLXF', 'SFbc_a4', 'SFpom_a4', 'SFso4_a1', 'SFso4_a2',
'so4_a1_sfgaex1', 'so4_a2_sfgaex1', 'so4_a3_sfgaex1', 'soa_a1_sfgaex1', 'soa_a2_sfgaex1',
'SFdst_a1','SFdst_a2', 'SFdst_a3', 'SFncl_a1', 'SFncl_a2', 'SFncl_a3',
'num_a2_sfnnuc1', 'SFSO2', 'OCN_FLUX_DMS', 'SAD_SULFC', 'SAD_TROP', 'SAD_AERO',
'AODVISstdn', 'REFF_AERO', 'so4_a1', 'so4_a2', 'SAD_AERO', 'PS', 'TROP_P', 'FSNTOAC',
'AODSSdn', 'AODSSdn', 'BURDENSO4dn', 'AODVISdn', 'AODVISstdn',
'AODBCdn','AODSOAdn','AODSSdn','AODPOMdn','AODSO4dn', 'AODDUSTdn','AODVISdn', 'AODVISstdn','BURDENSO4dn', 'SFDMS','BURDENBCdn', 'BURDENDUSTdn', 
'BURDENPOMdn', 'BURDENSEASALTdn', 'BURDENSO4dn', 'BURDENSOAdn', 'AODBCdn','AODSOAdn','AODSSdn','AODPOMdn','AODSO4dn', 'AODDUSTdn','AODVISdn', 'AODVISstdn','BURDENSO4dn', 'SFDMS','BURDENBCdn', 'BURDENDUSTdn', 'BURDENPOMdn', 'BURDENSEASALTdn', 'BURDENSO4dn', 'BURDENSOAdn', 
'CCN3', 'CFC11', 'CFC12', 'CH4','CO2', 'EXTINCTdn', 'N2O', 
'EFIX'

fincl3 = 'PRECT', 'PRECC', 'FLUT', 'U850', 'U200', 'V850', 'V200', 'OMEGA500', 'TS', 'SST', 'PSL', 'Z500'

fincl4 =  'PRECC','PRECL'


seasalt_emis_scale		= 1.6D0
tracer_cnst_datapath           = '/glade/campaign/acom/acom-climate/tilmes/CAM7/b.e30_alpha08o.B1850C_MTt4s.ne30_t232_wgx3.330_branch.fix/atm/proc/TracerCnstCyclical1850'

tracer_cnst_file               = 'tracer_cnst_MTt4s_halons_3Dmonthly_1850climo330_c20260423.nc'

prescribed_ozone_datapath      = '/glade/campaign/acom/acom-climate/tilmes/CAM7/b.e30_alpha08o.B1850C_MTt4s.ne30_t232_wgx3.330_branch.fix/atm/proc/ozone_strataero_day_5_zm'

prescribed_ozone_file          = 'ozone_strataero_b.e30_alpha08o.B1850C_MTt4s.ne30_t232_wgx3.330_branch.fix_1850_5_30_climo_zm_5day_c260424.nc'

prescribed_strataero_datapath  = '/glade/campaign/acom/acom-climate/tilmes/CAM7/b.e30_alpha08o.B1850C_MTt4s.ne30_t232_wgx3.330_branch.fix/atm/proc/ozone_strataero_day_5_zm'

 prescribed_strataero_file              = 'ozone_strataero_b.e30_alpha08o.B1850C_MTt4s.ne30_t232_wgx3.330_branch.fix_1850_5_30_climo_zm_5day_c260424.nc'

stream_ndep_data_filename              = '/glade/campaign/acom/acom-climate/tilmes/CAM7/b.e30_alpha08o.B1850C_MTt4s.ne30_t232_wgx3.330_branch.fix/atm/proc/NDepositionCyclical1850/fndep_clm_MTt4s_3Dmonthly_1850climo330_c20260424.nc'

clubb_gamma_coef                =  0.298     
clubb_gamma_coefb               =  0.298   

clubb_penta_solve_method = 2
clubb_tridiag_solve_method = 2
  • user_nl_clm
use_init_interp = .true.
reseed_dead_plants = .true.
use_excess_ice_streams=.false. 
  • 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
histaux_rof2med_file1_enabled = .true.
histaux_rof2med_file1_ntperfile = 8
  • user_nl_cice
r_snw = 0.75
restart_fsd = .true.
histfreq = "m", "x", "x", "x", "x"
  • user_nl_cism
dew = 4000.
dns = 4000.
ewn = 421
nsn = 721
upn = 11

idiag = 136
jdiag = 294

dycore = 2
flow_law = 2
evolution = 3
temperature = 1
temp_init = 3
marine_margin = 1
calving_domain = 1
calving_init = 1
limit_marine_cliffs = .false.
basal_mass_balance = 1

gthf = 1
smb_input = 0
smb_input_function = 0
artm_input_function = 0
bmlt_float = 0
dm_dt_diag = 1

compute_blocks = 0
which_ho_powerlaw_c = 2
which_ho_babc = 9
which_ho_bwat = 0
which_ho_effecpress = 0
which_ho_efvs = 2
which_ho_thermal_timestep = 0
which_ho_resid = 4
which_ho_sparse = 3
which_ho_nonlinear = 1
which_ho_gradient = 0
which_ho_gradient_margin = 1
which_ho_assemble_beta = 1
which_ho_assemble_bfric = 1
which_ho_assemble_taud = 1
which_ho_assemble_lateral = 1
which_ho_approx = 4
diva_slope_correction = .true.
which_ho_precond = 3
which_ho_calving_front = 0
which_ho_ground = 2
which_ho_flotation_function = 2
glissade_maxiter = 100
block_inception = .false.
remove_ice_caps = .false.

flow_factor_ground = 1.
flow_factor_float = 1.
ice_limit = 1.
thck_gradient_ramp = 100.
marine_limit = -200
max_slope = 0.5
rhoi = 917.
rhoo = 1026.
grav = 9.80616
shci = 2117.27
lhci = 3.337e5
trpt = 273.16
pmp_offset = 10.
beta_grounded_min = 100.
beta_powerlaw_umax = 5000.
pseudo_plastic_q = 0.25
pseudo_plastic_u0 = 100.
p_ocean_penetration = 0.5
cliff_timescale = 0.0
taumax_cliff = 1.e6
powerlaw_m = 3.0
powerlaw_c_const = 20000.
powerlaw_c_max = 100000.
powerlaw_c_min = 100.
  • user_nl_mom
DIAG_COORD_DEF_RHO2 = "FILE:ocean_rho2_190917.nc,interfaces=rho2"
MLE_TAPER_FN_DEPTH = 3000.
MLE_TAPER_FN_POWER = 4
BBL_EFFIC = 0.2 
EOF

-SourceMods

SourceMods/src.cam/ 
dimensions_mod.F90  
dyn_comp.F90  
fvm_mapping.F90  
global_norms_mod.F90  
prim_advance_mod.F90  
prim_driver_mod.F90  README  
viscosity_mod.F90
SourceMods/src.mom/ 
diag_table.unresolved  
MOM_mixed_layer_restrat.F90
Dominant language
No language data
Stars
11
Forks
3
PR merge metrics
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 Build System issues

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.