Skip to content

b.e30_alpha08b.B1850C_LTso.ne30_t232_wgx3.322 #306

@cecilehannay

Description

@cecilehannay

Purpose:
Repeat case 321 but start from a spun-up ocean instead of a G-case initial condition.

Case 321 shows a strong cooling at the beginning of the run and a negative RESTOM. We want to confirm whether this behavior is caused by starting from a G case.

Image

Description:
Same as 321 (#305) except:

  • start from 287 at year 125
./xmlchange RUN_REFCASE=b.e30_alpha07g.B1850C_LTso.ne30_t232_wgx3.287
./xmlchange RUN_REFDATE=0125-01-01
./xmlchange RUN_TYPE=hybrid
./xmlchange GET_REFCASE=true
./xmlchange RUN_REFDIR=cesm2_init
  • 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"

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.

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.322/
  • CAM Climos: /glade/campaign/cgd/amp/amwg/climo/b.e30_alpha08b.B1850C_LTso.ne30_t232_wgx3.322

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


Extra details:

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions