CESM initialization

Simulations

SE dycore

ID Casename Case Description Diagnostics Details/Comments
AA0 b.e11.B1850C5CN.f09_g16.005
Large-ensemble

timeseries (1-200)
timeseries (1-100)
w/ obs (300-319)

ocn timeseries (1-1000)
ocn timeseries (1-100)

sea-ice diags

Comments:
AA1 b.e12.B1850C5CN.ne30_g16.init.ch.001
Coupled run with CESM1.2.0
out of the box + hybrid from b.e12.B1850C5CN.ne30_g16.control.022 at yr 145
timeseries
w/ obs
Comments:
AA2 b.e12.B1850C5CN.ne30_g16.init.ch.002
Same as AA1
+ ocean coupling frequency set to 2 hours
timeseries
w/ obs
w/ AA1
Comments:
./xmlchange -file env_run.xml -id OCN_NCPL -val 12
./xmlchange -file env_run.xml -id ROF_NCPL -val 12

Gokhan hopes this configuration will produce a slight warming of the SSTs (by adding a diurnal cycle)
This doesn't have the expected effect. SSTS are slightly warmer but the TOA balance is off by almost 0.5 Wm/2.
AA3 b.e12.B1850C5CN.ne30_g16.init.ch.003

Same as AA1
+ rhminl = 0.8975D0->0.8935D0
+ rhminl = 0.8935D0->0.89D0 (at year 112)

timeseries
w/ obs
w/ Large ens
w/ AA1

ocn mean
ocn timeseries

sea-ice diags

Comments:
Retuning to balance TOA
AA4 b.e12.B1850C5CN.ne30_g16.init.ch.004
Same as AA3
+ the wind stress seen by the ocean is reduced by 30% south of 35S.
timeseries
w/ obs
w/ AA3
Comments:
the wind stress seen by the ocean is reduced by 30% south of 35S.

AA5 b.e12.B1850C5CN.ne30_g16.init.ch.005
Same as AA3
+ start from Levitus

timeseries
w/ obs
w/ Large ens
w/ AA3


Comments:
None

AA6 b.e12.B1850C5CN.ne30_g16.init.ch.006
Same as AA3
+ start from Large ensemble at year 400

timeseries
w/ obs
w/ Large ens
w/ AA3


Comments:
None

AA9 b.e12.B1850C5CN.ne30_g16.init.ch.009
Same as AA5
start from Levitus
+ rhminl = 0.8935D0->0.8895D0
+ rhminl = 0.8895D0->0.884D0 (at year 49)

timeseries
w/ obs
w/ Large ens
w/ AA5
w/ AA3

ocn mean
ocn timeseries

sea-ice diags


Comments:
None

AA12 b.e12.B1850C5CN.ne30_g16.init.ch.012
Same as AA9
start using spunup ocean at year 402
+ rhminl = 0.8935D0->0.885D0

timeseries
w/ obs
w/ Large ens
w/ AA3
w/ WW9

ocn mean
ocn timeseries


Comments:
None

AA13 b.e12.B1850C5CN.ne30_g16.init.ch.013
Same as AA9
start from core
+ rhminl = 0.8935D0->0.885D0

timeseries
w/ obs
w/ Large ens
w/ AA3
w/ AA9

ocn mean
ocn timeseries


Comments:
None

AA14 b.e12.B1850C5CN.ne30_g16.init.ch.014
Same as AA3
hybrid from AA3
+ using spunup ocean at year 402
+ rhminl = 0.8935D0->0.885D0

timeseries
w/ obs
w/ Large ens
w/ AA3

ocn mean
ocn timeseries


Comments:
None

AA15 b.e12.B1850C5CN.ne30_g16.init.ch.015
Same as AA9
+ the wind stress seen by the ocean is reduced by 30% south of 35S.

timeseries
w/ obs
w/ Large ens
w/ AA3
w/ AA9

ocn mean
ocn timeseries

 


Comments:
None

AA16 b.e12.B1850C5CN.ne30_g16.init.ch.016
Start from Levitus in SE
+ spunup land from AA3 at year 147
+ rhimnl = 0.884 (value from AA9)

+ extra output from ocean

timeseries
Sfc budget
w/ obs


Comments:
None

AA17 b.e12.B1850C5CN.ne30_g16.init.ch.017

Start from Levitus in FV
+ spunup land from large ensemble
+ rhminl = 0.8925 (value of large ensemble)

+ extra output from ocean

timeseries
Sfc budget
w/ obs


Comments:
I forgot to change the casename to "f09"

AA19 b.e12.B1850C5CN.ne30_g16.init.ch.019

Start from Levitus in SE
+ spunup land from AA3 at year 147
+ rhimnl = 0.884 (value from AA9)

+ attempt to change vertical mixing

+ extra output from ocean

timeseries
Sfc budget
w/ obs


Comments:


Here are the changes you need to make to POP to get it to run with constant mixing instead of KPP:

1) Add the following to user_nl_pop2

vmix_choice = 'const'
ltidal_mixing = .false.
lhoriz_varying_bckgrnd = .false.
overflows_on = .false.
overflows_interactive = .false.
const_vdc = 0.2
const_vvc = 0.2

2) Copy the following files from /glade/p/work/mlevy/cases/C.const_mixing.cesm1_2_0/SourceMods/src.pop2 to your SourceMods/src.pop2 directory

i) gx1v6_tavg_contents -- this is the list of all variables output by POP, and I removed the ones that are specific to KPP

ii) ocn.base.tavg.csh -- this file controls how many stream functions POP should output... by default there are 3, but I think the 3rd file only contains KPP-specific variables, so I changed it to 2 streams.

For the experiment you described (changing the kappa term in the diffusion equation), add the following two lines to user_nl_pop2:

const_vdc = [new value; old value is 0.25]
const_vvc = [new value; old value is 0.25]

I was able to run for 5 days with the default (0.25) values as well as dropping down to 0.2; that's all the testing I did to make sure these code-mods wouldn't break anything.

 

AA20 b.e12.B1850C5CN.ne30_g16.init.ch.020
Start from Levitus in SE
+ spunup land from AA3 at year 147
+ rhimnl = 0.884 (value from AA9)
+ dust_emis_fact = 0.35D0 (same as in FV)

+ extra output from ocean

timeseries
Sfc budget
w/ obs


Comments:
None

AA21 b.e12.B1850C5CN.ne30_g16.init.ch.021

Start from Levitus in SE
+ spunup land from AA3 at year 147
+ rhimnl = 0.884 (value from AA9)

+ mapping files (see details)


+ extra output from ocean

timeseries
Sfc budget
w/ obs


Comments:
atm2ocn_fmapname = '/glade/p/cesmdata/cseg/inputdata/cpl/gridmaps/ne30np4/map_ne30np4_TO_gx1v6_aave.130522.nc'
atm2ocn_smapname = '/glade/p/cesmdata/cseg/inputdata/cpl/gridmaps/ne30np4/map_ne30np4_TO_gx1v6_blin.130522.nc'
atm2ocn_vmapname = '/glade/p/cesmdata/cseg/inputdata/cpl/gridmaps/ne30np4/map_ne30np4_TO_gx1v6_patc.130522.nc'
ocn2atm_fmapname = '/glade/p/cesmdata/cseg/inputdata/cpl/gridmaps/gx1v6/map_gx1v6_TO_ne30np4_aave.130522.nc'
ocn2atm_smapname = '/glade/p/cesmdata/cseg/inputdata/cpl/gridmaps/gx1v6/map_gx1v6_TO_ne30np4_aave.130522.nc'

AA22 b.e12.B1850C5CN.ne30_g16.init.ch.022

Start from Levitus in SE
+ spunup land from AA3 at year 147
+ rhimnl = 0.884 (value from AA9)

+ changes from gohkan to modify TAUX seen by ocean

+ extra output from ocean

timeseries
Sfc budget
w/ obs


Comments:
doesn't work properly

AA23 b.e12.B1850C5CN.ne30_g16.init.ch.023

Start from Levitus in SE
+ spunup land from AA3 at year 147
+ rhimnl = 0.884 (value from AA9)

+ turn TMS off

+ extra output from ocean

timeseries
Sfc budget
w/ obs

w/ control


Comments:


AA24 b.e12.B1850C5CN.ne30_g16.init.ch.024

Start from Levitus in SE
+ spunup land from AA3 at year 147
+ rhimnl = 0.884 (value from AA9)

+ tms_orocnst = 1.0->2.0D0

+ extra output from ocean

timeseries
Sfc budget
w/ obs

vs /control


Comments:


AA25 b.e12.B1850C5CN.ne30_g16.init.ch.025

Start from Levitus in SE
+ spunup land from AA3 at year 147
+ rhimnl = 0.884 (value from AA9)

+ old SE topo (smooth):
USGS-gtopo30_0.9x1.25_smooth2000-100_ne30np4_c051027.nc


+ extra output from ocean

timeseries
Sfc budget
w/ obs

w/ control


Comments:


AA26 b.e12.B1850C5CN.ne30_g16.init.ch.026

Start from Levitus in SE
+ spunup land from AA3 at year 147
+ rhimnl = 0.884 (value from AA9)

+ hdsp(i) = 4.0_r8 * sgh(i)


+ extra output from ocean

timeseries
Sfc budget
w/ obs

w/ control


Comments:


AA27 b.e12.B1850C5CN.ne30_g16.init.ch.027

Start from Levitus in SE
+ spunup land from AA3 at year 147
+ rhminl and rpen like FV
+ cldfrc_rhminl = 0.8925D0
+ uwshcu_rpen = 10.0


+ extra output from ocean

timeseries
Sfc budget
w/ obs


Comments:


AA28 b.e12.B1850C5CN.ne30_g16.init.ch.028

Start from Levitus in SE
+ spunup land from AA3 at year 147
+ nugding


+ extra output from ocean

timeseries
Sfc budget
w/ obs

w/ control


Comments:


AA29 b.e12.B1850C5CN.ne30_g16.init.ch.029

Start from Levitus in SE
+ spunup land from AA3 at year 147
+ strong nudging


+ extra output from ocean

timeseries
Sfc budget
w/ obs

w/ control


Comments:


AA31 b.e12.B1850C5CN.poles_30x4_g16.init.ch.031

Start from Levitus in SE
+ spunup land from AA3 at year 147
+ refined poles


+ extra output from ocean

timeseries
Sfc budget
w/ obs


Comments:


AA32 b.e12.B1850C5CN.ne30_g16.init.ch.032

Start from Levitus in SE
+ spunup land from AA3 at year 147
+ nudging U and V only


+ extra output from ocean

timeseries
Sfc budget


Comments:


AA34 b.e12.B1850C5CN.poles_30x4_g16.init.ch.034

Same as AA16:
Start from Levitus in SE
+ spunup land from AA3 at year 147
+ rhimnl = 0.884 (value from AA9)
+ extra output from ocean

+ min(kappa_{iso,thic}) = 600 m^2/s

timeseries
w/ obs


Comments:


AA35 b.e12.B1850C5CN.poles_30x4_g16.init.ch.035

Same as AA16:
Start from Levitus in SE
+ spunup land from AA3 at year 147
+ rhimnl = 0.884 (value from AA9)
+ extra output from ocean

+ min(kappa_thic) = 600 m^2/s, min(kappa_iso) = 1200 m^2/s

timeseries
w/ obs


Comments:


AA36 b.e12.B1850C5CN.poles_30x4_g16.init.ch.036

Same as AA16:
Start from Levitus in SE
+ spunup land from AA3 at year 147
+ rhimnl = 0.884 (value from AA9)
+ extra output from ocean

+ min(kappa_{iso,thic}) = 600 m^2/s and KPP modification (stability-dependent surface layer depth)

timeseries
w/ obs


Comments: