Instructions to run CESM1.5: "28"

 

CESM1.5 (mini-breck version)

This is the cesm1.5 discussed at the mini breck. The so called "28".

The control runs for the "28" simulations (pre-industrial control, 20th century, standalone simualtions) can be found at here,
All the other cesm1.5 development simulations can be found at: CESM1.5 simulations ("mini Breck" simulations)

Instructions to reproduce the preindustrial control for 28 are below.

Experimental tag

On the svn repository
https://svn-ccsm-models.cgd.ucar.edu/cesm1/exp/branch_tags/cesm1_5_beta03_cntlexp_tags/cesm1_5_beta03_cntlexp06

There is a local copy of the tag at:

/glade/p/work/hannay/cesm_tags/cesm1_5_beta03_cntlexp06

Compset and create_newcase command

In this tag, the B1850 compset uses:
- ATM: cam5.5 physics
- SEA ICE: prognostic cice5
- LAND ICE: cism1 land ice
- LAND: clm5.0 physics + clm5.0 bgc (cn and methane)
- RIVER MODEL: MOSART
- OCN: POP2 default+ BGC (CO2=diag)

An example of the create_newcase command is:
./create_newcase -case $yourcasename -compset B1850 -res f09_g16 -mach yellowstone -user_pes_setby allactive -project $project_number

pe_layout

The throughput/cost for the following pe_layout
/glade/p/cesmdata/cseg/runs/cesm1_5/b.e15.B1850G.f09_g16.pi_control.28/env_mach_pes.xml

Model Cost: 4526.52 pe-hrs/simulated_year
Model Throughput: 20.02 simulated_years/day

code modification

Copy in SourceMods/src/cam
/glade/p/cesmdata/cseg/runs/cesm1_5/b.e15.B1850G.f09_g16.pi_control.28/SourceMods/src.cam/parameters_tunable.F90

xml changes

# run-off mapping file
./xmlchange ROF2OCN_RMAPNAME=cpl/gridmaps/r05/map_r05_to_gx1v6_e1000r300_151109.nc

# Determines provenance of surface Chl for radiative penetration,computations.
./xmlchange OCN_CHL_TYPE=diagnostic

Change to user_nl_** files

# user_nl_cam

bnd_topo = '/glade/p/cesmdata/cseg/inputdata/atm/cam/topo/fv_0.9x1.25-gmted2010_modis-cam_fv_smooth-intermediate_ncube3000-no_anisoSGH_c151029.nc'
dust_emis_fact = 0.60D0
micro_mg_dcs = 140.D-6


# user_nl_clm

suplnitro = 'NONE'
use_flexiblecn=.true.
use_luna=.true.
use_fun=.true.
fsurdat = '/glade/p/cesmdata/cseg/inputdata/lnd/clm2/surfdata_map/surfdata_0.9x1.25_16pftsmodwshrub_simyr1850_c151210.nc'
paramfile = '/glade/p/cesmdata/cseg/inputdata/lnd/clm2/paramdata/clm5_params.c151217.nc'
leaf_mr_vcm = 0.015d00
use_dynroot = .false.
lower_boundary_condition = 2
use_bedrock = .true.
rooting_profile_method=1
soil_layerstruct = '20SL_8.5m'
limit_irrigation = .false.
use_init_interp = .true.
init_interp_fill_missing_with_natveg = .true.

# user_nl_cpl

flux_diurnal = .true.

# user_nl_cice
# (by default the initial condition for sea-ice is coming from a previous SOM run. You want to use the initial condition from the hybrid run instead)

ice_ic = 'b.e15.B1850G.f09_g16.pi_control.25.cice.r.0041-01-01-00000.nc'

# user_nl_pop

alk_riv_flux_input%filename = '/glade/p/cesm/bgcwg/forcing/BEC_gx1v6_forcing_no_ms/river_nutrients_GNEWS2000_gx1v6_no_ms_c150702.nc'
dfe_riv_flux_input%filename = '/glade/p/cesm/bgcwg/forcing/BEC_gx1v6_forcing_no_ms/river_nutrients_GNEWS2000_gx1v6_no_ms_c150702.nc'
dic_riv_flux_input%filename = '/glade/p/cesm/bgcwg/forcing/BEC_gx1v6_forcing_no_ms/river_nutrients_GNEWS2000_gx1v6_no_ms_c150702.nc'
din_riv_flux_input%filename = '/glade/p/cesm/bgcwg/forcing/BEC_gx1v6_forcing_no_ms/river_nutrients_GNEWS2000_gx1v6_no_ms_c150702.nc'
dip_riv_flux_input%filename = '/glade/p/cesm/bgcwg/forcing/BEC_gx1v6_forcing_no_ms/river_nutrients_GNEWS2000_gx1v6_no_ms_c150702.nc'
doc_riv_flux_input%filename = '/glade/p/cesm/bgcwg/forcing/BEC_gx1v6_forcing_no_ms/river_nutrients_GNEWS2000_gx1v6_no_ms_c150702.nc'
don_riv_flux_input%filename = '/glade/p/cesm/bgcwg/forcing/BEC_gx1v6_forcing_no_ms/river_nutrients_GNEWS2000_gx1v6_no_ms_c150702.nc'
dop_riv_flux_input%filename = '/glade/p/cesm/bgcwg/forcing/BEC_gx1v6_forcing_no_ms/river_nutrients_GNEWS2000_gx1v6_no_ms_c150702.nc'
dsi_riv_flux_input%filename = '/glade/p/cesm/bgcwg/forcing/BEC_gx1v6_forcing_no_ms/river_nutrients_GNEWS2000_gx1v6_no_ms_c150702.nc'
dust_flux_input%filename = '/glade/p/cesm/bgcwg/forcing/BEC_gx1v6_forcing_no_ms/dst79gnx_gx1v6_090416_no_ms_c150702.nc'
fesedflux_input%filename = '/glade/p/cesm/bgcwg/forcing/BEC_gx1v6_forcing_no_ms/fesedflux_gx1v6_etopo2v2_Nov2015_vents_no_ms_c151218.nc'
iron_flux_input%filename = '/glade/p/cesm/bgcwg/forcing/BEC_gx1v6_forcing_no_ms/solFe_scenario4_current_gx1v6_no_ms_c150702.nc'
nhy_flux_monthly_input%filename = '/glade/p/cesm/bgcwg/forcing/BEC_gx1v6_forcing_no_ms/ndep_ocn_1850_gx1v6_no_ms_c150702.nc'
nox_flux_monthly_input%filename = '/glade/p/cesm/bgcwg/forcing/BEC_gx1v6_forcing_no_ms/ndep_ocn_1850_gx1v6_no_ms_c150702.nc'

env_mach_specific.xml

Remove the following line from the file env_mach_specific.xml
<env name="MP_EAGER_LIMIT">0</env>

It will increase significantly your throughput.