Purpose:
This simulation is testing the CESM with a fully dynamically evolving Greenland ice sheet with the removing ice cap option turned off.
Description:
This is a hybrid run branching from year 73 of #320. Prior to running the model, the Greenland ice sheet (GrIS) was initialized following the workflow:
Run an I compset using data atmosphere from run 334 between years 60 and 88.
Spin-up the GrIS for 10,000 years with CISM standalone using coupler history files from the I compset for SMB and temperature.
Use the last time slice of the CISM spinup as CISM input file for this run.
This issue is similar to #338 except that the configuration option in CISM remove_ice_caps = .false.
Case directory:
- Locally (if still available):
/glade/campaign/cesm/development/liwg/leguy/CESM3/cases/bg.e30_beta08.B1850C_MTso.ne30_t232_wgx3.334.gris2way_init_CLMicase_GLCwithIceCap
Sandbox:
- Locally (if still available): /glade/campaign/cesm/development/liwg/leguy/CESM3/code/cesm3_0_beta08
Output:
- Raw output on derecho (if still available): /glade/derecho/scratch/gunterl/archive/bg.e30_beta08.B1850C_MTso.ne30_t232_wgx3.334.gris2way_init_CLMicase_GLCwithIceCap
Contacts:
@gunterl
Extra details:
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 = .true.
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.
cisminputfile='/glade/campaign/cesm/development/liwg/leguy/CESM3/CESM3_input/GrIS/input_for_BGrun334/gris_init_334.cesm.nc'
user_nl_clm
reseed_dead_plants = .true.
use_excess_ice_streams=.false.
use_init_interp = .true.
user_nl_cam
mfilt = 0, 5, 20, 40, 12, 120, 1, 1
nhtfrq = 0, -24, -24, -3, 0, -2, 0, -8760
ndens = 2, 2, 2, 2, 2, 1, 2, 1
interpolate_output = .true., .true., .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'
fincl3 = 'PRECT', 'PRECC', 'FLUT', 'U850', 'U200', 'V850', 'V200', 'OMEGA500', 'TS', 'SST', 'PSL', 'Z500'
fincl4 = 'PRECC','PRECL'
phys_grid_ctem_nfreq=-6
phys_grid_ctem_zm_nbas=120
phys_grid_ctem_za_nlat=90
seasalt_emis_scale= 1.55D0
user_nl_cice
grid_file = "/glade/work/gmarques/cesm/tx2_3/seaice_files/cice_grid_tx2_3v3_260305_cdf5.nc"
grid_format = 'nc'
kmt_file = "/glade/work/gmarques/cesm/tx2_3/seaice_files/cice_grid_tx2_3v3_260305_cdf5.nc"
r_snw = 1.0
restart_fsd = .true.
histfreq = "m", "x", "x", "x", "x"
user_nl_mom
INPUTDIR = "/glade/derecho/scratch/hannay/b.e30_alpha08o.B1850C_MTso.ne30_t232_wgx3.330/run/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.
CORRECT_BBL_BOUNDS = True
CHANNEL_DRAG_MAX_BBL_THICK = 5.0
user_nl_cpl
histaux_l2x1yrg = .true.
Purpose:
This simulation is testing the CESM with a fully dynamically evolving Greenland ice sheet with the removing ice cap option turned off.
Description:
This is a hybrid run branching from year 73 of #320. Prior to running the model, the Greenland ice sheet (GrIS) was initialized following the workflow:
Run an I compset using data atmosphere from run 334 between years 60 and 88.
Spin-up the GrIS for 10,000 years with CISM standalone using coupler history files from the I compset for SMB and temperature.
Use the last time slice of the CISM spinup as CISM input file for this run.
This issue is similar to #338 except that the configuration option in CISM remove_ice_caps = .false.
Case directory:
/glade/campaign/cesm/development/liwg/leguy/CESM3/cases/bg.e30_beta08.B1850C_MTso.ne30_t232_wgx3.334.gris2way_init_CLMicase_GLCwithIceCap
Sandbox:
Output:
Contacts:
@gunterl
Extra details:
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 = .true.
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.
cisminputfile='/glade/campaign/cesm/development/liwg/leguy/CESM3/CESM3_input/GrIS/input_for_BGrun334/gris_init_334.cesm.nc'
user_nl_clm
reseed_dead_plants = .true.
use_excess_ice_streams=.false.
use_init_interp = .true.
user_nl_cam
mfilt = 0, 5, 20, 40, 12, 120, 1, 1
nhtfrq = 0, -24, -24, -3, 0, -2, 0, -8760
ndens = 2, 2, 2, 2, 2, 1, 2, 1
interpolate_output = .true., .true., .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'
fincl3 = 'PRECT', 'PRECC', 'FLUT', 'U850', 'U200', 'V850', 'V200', 'OMEGA500', 'TS', 'SST', 'PSL', 'Z500'
fincl4 = 'PRECC','PRECL'
phys_grid_ctem_nfreq=-6
phys_grid_ctem_zm_nbas=120
phys_grid_ctem_za_nlat=90
seasalt_emis_scale= 1.55D0
user_nl_cice
grid_file = "/glade/work/gmarques/cesm/tx2_3/seaice_files/cice_grid_tx2_3v3_260305_cdf5.nc"
grid_format = 'nc'
kmt_file = "/glade/work/gmarques/cesm/tx2_3/seaice_files/cice_grid_tx2_3v3_260305_cdf5.nc"
r_snw = 1.0
restart_fsd = .true.
histfreq = "m", "x", "x", "x", "x"
user_nl_mom
INPUTDIR = "/glade/derecho/scratch/hannay/b.e30_alpha08o.B1850C_MTso.ne30_t232_wgx3.330/run/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.
CORRECT_BBL_BOUNDS = True
CHANNEL_DRAG_MAX_BBL_THICK = 5.0
user_nl_cpl
histaux_l2x1yrg = .true.