As a deep time modeler, you are responsible for creating a new topography file (bnd_topo) which you will point to in the user_nl_cam namelist. You also specify trace gas constituents and the solar constant in user_nl_cam, if they differ from pre-industrial. The CESM User Guide describes the user namelist functionality in CESM1.
A note about topography
We recommend that you smooth your high resolution (10min) topography to reduce the liklihood of model instabilities. We provide a rudimentary smoothing program for this purpose: smth9_topo.ncl. The code allows you to set the number of smoothing passes (niter) you make over your initial topography, and produces a latitudinal cross section for comparison.
Atmospheric boundary conditions ~ bnd_topo
The atmosphere component model requires initial boundary conditions, including fractional landcover (LANDFRAC), topography (PHIS). surface roughness, as the standard deviation of topography (SGH), and of high resolution surface roughness (SGH30), and an aerosol scaling factor (LANDM_COSLAT). For deep time, SGH and SGH30 are estimated by multiplying the input surface topography by 0.3 and 0.1, respectively.
Initial atmospheric state ~ ncdata
We suggest two approaches for providing the initial atmospheric state ("cami") the model needs.
1) We provide a script that will produce a basic physical atmospheric state suitable for initialization that allows the atmospheric model to spin-up during the first few decades of integration. The tools compute a cosine weighted equator-pole latitudinal temperature distribution based on user-specified polar and equator temperature estimates. Default values are set at 28oC (equator) and 12oC (pole). Surface pressure and the surface geopotential are computed based on the topography (topobathy.nc), and the zonal velocity is based on the thermal wind balance. The initial condition can also be used to remove the sea ice fraction.
src: create_cesm1_cami-InitialConditions.ncl, create_cesm1_cami-ZERO-ICEFRAC-Only.ncl
2) Alternatively, we successfully initialize deep time simulations using the fully coupled model from a default cesm initial condition file (e.g., cami_0000-01-01_1.9x2.5_L30_c070703.nc). The model equilibrates to the deep time conditions within a few months to years.
- Src: run2create_cesm1_cami.csh
- Calls: create_cesm1_cami.ncl
- input: topobathy.nc
- output: cami_bnd_topo_RESOLN_MYRUN.nc (bnd_topo)
The c-shell script run2create_cesm1_cami.csh drives an NCL script called create_cesm1_cami.ncl which requires topobathy.nc as well as a cami skeleton file at your atm/land resolution. The skeleton file is simply used as a template for the required cami variables names and will be overwritten by the basic physical atmospheric state computed by create_cesm1_cami.ncl. You can use any cesm1 cami file from a previous CAM run as a template for this process.
Edit the environmental variables in run2create_cesm1_cami.csh to your local environment. Details on this script can be found in the comments of the script itself.
NOTE: Angular, highly variable topography can cause instability in the CAM atmosphere model that is difficult to detect but will cause the model to crash without errors. We recommend smoothing PHIS in (cami_bnd_topo) if your topography is angular or highly variable.
Other atmoshere forcing files:
- The present day absorption/emissivity forcing dataset was built with wide constraints and is therefore flexible and can be used for paleoclimate cases.
- Typically, present day or pre-industrial ozone mixing ratio boundary forcing files are used for paleoclimate cases. Choice of dataset will depend on your control experiment.
- Although it is necessary to input present day or pre-industrial aerosol boundary forcing files into the model, it is recommended that deep time paleoclimate researches use the namelist functionality to flag CAM to internally compute aerosol optical depths.
- [CESM1 DEFAULT]
- [CESM1 DEFAULT]
Namelist adjustments for physical forcing:
Group A: Solar constant and trace gases
Define appropriate values for the solar constant and for trace gas concentrations for your time period. Orbital parameters are set in the coupler.
|CO2VM2||CO2 volume mixing ratio|
|CH4VMR||CH4 volume mixing ratio|
|N20VMR||N2O volume mixing ratio|
|F11VMR1||CFC11 volume mixing ratio|
|F12VMR1||CFC12 volume mixing ration|
1 For deep time paleo experiments, F11VMR and F12VMR should be set to 0.
Group B: Aerosol Optical Depth
CAM radiation code requires a boundary forcing dataset for aerosol mass mixing ratios and aerosol optical properties. Because this is unknown for many paleoclimate cases, we must flag the code to use the CCM3.6 method for computing aerosol optical depths, which does not require spatial knowledge for aerosol mixing ratios. This method will compute a uniform optical depth across all grid points. Group B parameters must all be specified in the namelist.
A tauback value of 0.28 is equivalent to CCM3.6 default value of TAUVIS =0.14. The various scaling parameters must be set to zero to render the present day aerosol mixing ratio values in the code obsolete and allow the uniform optical depth to be used.
Default values for all namelist parameters can be found in the CAM3 User’s Guide.
|TAUBACK||uniform optical depth||0.28|
|SSLTSCL||sea salt scaling||0.0|