trefread/0000755006307300054160000000000013047376462012027 5ustar islascgdamptrefread/namelist_definition/0000755006307300054160000000000013047371545016050 5ustar islascgdamptrefread/namelist_definition/namelist_definition.xml0000644006307300054160000075424513047371545022637 0ustar islascgdamp Set to true to read in relaxation temperature profile from file Path for relaxation temperature profile file Full pathname of time-variant boundary dataset for aerosol masses. Default: set by build-namelist. Add CAM3 prescribed aerosols to the physics buffer. Default: FALSE Dynamics/physics transpose method for nonlocal load-balance. 0: use mpi_alltoallv. 1: use point-to-point MPI-1 two-sided implementation. 2: use point-to-point MPI-2 one-sided implementation if supported, otherwise use MPI-1 implementation. 3: use Co-Array Fortran implementation if supported, otherwise use MPI-1 implementation. 11-13: use mod_comm, choosing any of several methods internal to mod_comm. The method within mod_comm (denoted mod_method) has possible values 0,1,2 and is set according to mod_method = phys_alltoall - modmin_alltoall, where modmin_alltoall is 11. -1: use option 1 when each process communicates with less than half of the other processes, otherwise use option 0 (approx.). Default: -1 Select target number of chunks per thread. Must be positive. Default: 1 Physics grid decomposition options. -1: each chunk is a dynamics block. 0: chunk definitions and assignments do not require interprocess comm. 1: chunk definitions and assignments do not require internode comm. 2: optimal diurnal, seasonal, and latitude load-balanced chunk definition and assignments. 3: chunk definitions and assignments only require communication with one other process. 4: concatenated blocks, no load balancing, no interprocess communication. Default: 2 Physics grid decomposition options. 0: assign columns to chunks as single columns, wrap mapped across chunks 1: use (day/night; north/south) twin algorithm to determine load-balanced pairs of columns and assign columns to chunks in pairs, wrap mapped Default: 0 for unstructured grid dycores, 1 for lat/lon grid dycores Output constituent tendencies due to convection. Set to 'none', 'q_only' or 'all'. Default: 'q_only', but 'none' for WACCM configurations. Turns on TEM circulation diagnostics history output. Only valid for FV dycore. Default: .true. for WACCM with interactive chemistry, otherwise .false. Turn on verbose output identifying columns that fail energy/water conservation checks. Default: FALSE Number of layers from the top of the model over which to do dry convective adjustment. Must be less than plev (the number of vertical levels). Default: 3 Number of dynamics timesteps per physics timestep. If zero, a best-estimate will be automatically calculated. Default: 0 Number of tracer advection timesteps per physics timestep. Nsplit is partitioned into nspltrac and nsplit/nspltrac, with the latter being the number of dynamics timesteps per tracer timestep, possibly rounded upward; after initialization, the code quantity nsplit is redefined to be the number of dynamics timesteps per tracer timestep. Default: 0 Number of vertical re-mapping timesteps per physics timestep. Nspltrac is partitioned into nspltvrm and nspltrac/nspltvrm, with the latter being the number of tracer timesteps per re-mapping timestep, possibly rounded upward; after initialization, the code quantity nspltrac is redefined to be the number of tracer timesteps per re-mapping timestep. Default: 0 Order (mode) of X interpolation (1,..,6). East-West transport scheme. Default: 4 Order (mode) of Y interpolation (1,..,6). North-South transport scheme. Default: 4 Scheme to be used for vertical mapping. Default: 4 Flag indicating whether the dynamics uses internal algorithm for energy conservation. Default: .false. Enables optional filter for intermediate c-grid winds, (courtesy of Bill Putman). Default: 0 1 for FFT filter always, 0 for combined algebraic/FFT filter. The value 0 is used for CAM3, otherwise it is using the value 1. Default: set by build-namelist Chooses type of divergence damping and velocity diffusion. div24del2flag = 2 for ldiv2 (default), = 4 for ldiv4, = 42 for ldiv4 + ldel2 where ldiv2: 2nd-order divergence damping everywhere and increasing in top layers ldiv4: 4th-order divergence damping ldel2: 2nd-order velocity-component damping targetted to top layers, with coefficient del2coef Default: set by build-namelist Chooses level of velocity diffusion. Default: 3.0e5 Set to 1 to force the 2D transpose computation when a 1D decomposition is used. This is intended for debugging purposes only. Default: 0 Geopotential method (routines geopk, geopk16, or geopk_d). 0 for transpose method, 1 for method using semi-global z communication with optional 16-byte arithmetic, 2 for method using local z communication; method 0, method 1 with 16-byte arithmetic and method 2 are all bit-for-bit across decompositions; method 0 scales better than method 1 with npr_z, and method 1 is superior to method 0 for small npr_z. The optimum speed is attained using either method 1 with 8-byte arithmetic (standard for geopk16) or method 2 when utilizing the optimal value for the associated parameter geopkblocks; for the last two subcycles of a timestep, method 0 is automatically used; see geopk.F90 and cd_core.F90. Default: 0 Geopotential method 2 pipeline parameter (routine geopk_d). geopk_d implements a pipeline algorithm by dividing the information that must be moved between processes into blocks. geopkblocks specifies the number of blocks to use. The larger the number of blocks, the greater the opportunity for overlapping communication with computation and for decreasing instantaneous bandwidth requirements. The smaller the number of blocks, the fewer MPI messages sent, decreasing MPI total latency. See geopk_d within geopk.F90. Default: 1 Mod_comm irregular underlying communication method for dyn_run/misc. 0 for original mp_sendirr/mp_recvirr 1 for mp_swapirr and a point-to-point implementation of communication pattern 2 for mp_swapirr and a collective (MPI_Alltoallv) implementation of communication pattern Default: 0 True for mod_comm irregular communication handshaking for dyn_run/misc Default: .true. True for mod_comm irregular communication blocking send for dyn_run/misc, false for nonblocking send Default: .true. Maximum number of outstanding nonblocking MPI requests to allow when using mp_swapirr and point-to-point communications for dyn_run/misc. Setting this less than the maximum can improve robustness for large process count runs. If set to less than zero, then do not limit the number of outstanding send/receive requests. Default: -1 (so no limit) Mod_comm irregular underlying communication method for cd_core/geopk 0 for original mp_sendirr/mp_recvirr 1 for mp_swapirr and a point-to-point implementation of communication pattern 2 for mp_swapirr and a collective (MPI_Alltoallv) implementation of communication pattern Default: 0 True for mod_comm irregular communication handshaking for cd_core/geopk Default: .true. True for geopk_d and mod_comm irregular communication blocking send for cd_core/geopk; false for nonblocking send. Default: .true. Maximum number of outstanding nonblocking MPI requests to allow when using mp_swapirr and point-to-point communications for cd_core/geopk. Setting this less than the maximum can improve robustness for large process count runs. If set to less than zero, then do not limit the number of outstanding send/receive requests. Default: -1 (so no limit) Mod_comm irregular underlying communication method for gather 0 for original mp_sendirr/mp_recvirr 1 for mp_swapirr and a point-to-point implementation of communication pattern 2 for mp_swapirr and a collective (MPI_Alltoallv) implementation of communication pattern Default: 0 True for mod_comm irregular communication handshaking for gather Default: .true. True for mod_comm irregular communication blocking send for gather, false for nonblocking send Default: .true. Maximum number of outstanding nonblocking MPI requests to allow when using mp_swapirr and point-to-point communications for gather. Setting this less than the maximum can improve robustness for large process count runs. If set to less than zero, then do not limit the number of outstanding send/receive requests. Default: -1 (so no limit) Mod_comm irregular underlying communication method for scatter 0 for original mp_sendirr/mp_recvirr 1 for mp_swapirr and a point-to-point implementation of communication pattern 2 for mp_swapirr and a collective (MPI_Alltoallv) implementation of communication pattern Default: 0 True for mod_comm irregular communication handshaking for scatter Default: .true. True for mod_comm irregular communication blocking send for scatter, false for nonblocking send Default: .true. Maximum number of outstanding nonblocking MPI requests to allow when using mp_swapirr and point-to-point communications for scatter. Setting this less than the maximum can improve robustness for large process count runs. If set to less than zero, then do not limit the number of outstanding send/receive requests. Default: -1 (so no limit) Mod_comm irregular underlying communication method for multiple tracers 0 for original mp_sendtrirr/mp_recvtrirr 1 for mp_swaptrirr and point-to-point communications 2 for mp_swaptrirr and all-to-all communications Default: 0 True for mod_comm irregular communication handshaking for multiple tracers Default: .true. True for mod_comm irregular communication blocking send for multiple tracers, false for nonblocking send Default: .true. Maximum number of outstanding nonblocking MPI requests to allow when using mp_swaptrirr and point-to-point communications for multiple tracers. Setting this less than the maximum can improve robustness for large process count runs. If set to less than zero, then do not limit the number of outstanding send/receive requests. Default: -1 (so no limit) One or two simultaneous mod_comm irregular communications (excl. tracers) Default: 1 Max number of tracers for simultaneous mod_comm irregular communications Default: 1 For mod_comm gather/scatters, 0 for temporary contiguous buffers; 1 for mpi derived types. Default: 0 For geopk (geopktrans=1) messages, 0 for temporary contiguous buffers; 1 for mpi derived types. Default: 0 For mod_comm transposes, 0 for temporary contiguous buffers; 1 for mpi derived types. Default: 0 A four element integer array which specifies the YZ and XY decompositions. The first two elements are the number of Y subdomains and number of Z subdomains in the YZ decomposition. The second two elements are the number of X subdomains and the number of Y subdomains in the XY decomposition. Note that both the X and Y subdomains must contain at least 3 grid points. For example, a grid with 96 latitudes can contain no more than 32 Y subdomains. There is no restriction on the number of grid points (levels) in a Z subdomain, but note that the threading parallelism in the FV dycore is over levels, so for parallel efficiency it is best to have at least the number of levels in each Z subdomain as there are threads available. There are a couple of rough rules of thumb to follow when setting the 2D decompositions. The first is that the number of Y subdomains in the YZ decomposition should be the same as the number of Y subdomains in the XY decomposition (npr_yz(1) == npr_yz(4)). The second is that the total number of YZ subdomains (npr_yz(1)*npr_yz(2)) should equal the total number of XY subdomains (npr_yz(3)*npr_yz(4)). Default: ntask,1,1,ntask where ntask is the number of MPI tasks. This is a 1D decomposition in latitude. Overlapping of trac2d and cd_core subcycles. Default: 0 Size of tracer domain decomposition for trac2d. Default: 1 TRUE => the offline meteorology winds are defined on the model grid cell walls. Default: FALSE Name of file that contains the offline meteorology data. Default: none Name of directory that contains the offline meteorology data. Default: none Name of file that contains names of the offline meteorology data files. Default: none TRUE => the offline meteorology file will be removed from local disk when no longer needed. Default: FALSE (km) top of relaxation region of winds for offline waccm Default: (km) bottom of relaxation region of winds for offline waccm Default: Relaxation time (hours) applied to specified meteorology. - positive values less then time step size gives 100% nudging - negative values gives 0.0% nudging (infinite relaxation time) Default: 0.0 switch to turn on/off mass fixer for offline driver Default: true srf heat flux field name in met data file Default: water vapor flux field name in met data file Default: multiplication factor for srf heat flux Default: 1.0 multiplication factor for water vapor flux Default: 1.0 multiplication factor for snow hieght Default: 1.0 if false then do not allow surface models feedbacks influence climate Default: true if true nudge meteorology surface fields TAUX, TAUY, SHFLX, QFLX rather than force Default: true del^2 horizontal diffusion coefficient. This is used above the Nth order diffusion. Default: set by build-namelist Order (N) of horizontal diffusion operator used below the sponge layers. N must be a positive multiple of 2. Default: 4 The order N horizontal diffusion operator will be used in and below the layer specified by this variable. Default: 4 Nth order horizontal diffusion coefficient. Default: set by build-namelist Number of days (from timestep 0) to run divergence damper. Use only if spectral model becomes dynamicallly unstable during initialization. Suggested value: 2. (Value must be >= 0.) Default: 0. Time filter coefficient. Default: 0.06 Number of levels over which to apply Courant limiter, starting at top of model. Default: 5 Number of dynamics timesteps per physics timestep. If zero, a best-estimate will be automatically calculated. Default: 1 del^2 horizontal diffusion coefficient. This is used above the 4th order diffusion. Default: set by build-namelist del^4 horizontal diffusion coefficient used below the sponge layers. Default: set by build-namelist Number of days (from timestep 0) to run divergence damper. Use only if spectral model becomes dynamicallly unstable during initialization. Suggested value: 2. (Value must be >= 0.) Default: 0. Time filter coefficient. Default: 0.06 Number of levels over which to apply Courant limiter, starting at top of model. Default: 5 Spectral dynamics gather option. Default: 0 Spectral dynamics transpose option. Default: 0 Flag indicating whether to assign latitudes to equidistribute columns or latitudes. This only matters when using a reduced grid. Default: TRUE Number of processes assigned to dynamics (SE, EUL and SLD dycores). Default: Total number of processes assigned to job. Stride for dynamics processes (EUL and SLD dycores). E.g., if stride=2, assign every second process to the dynamics. Default: 1 Whether or not to enable gravity waves produced by orography. Default: set by build-namelist. Whether or not to enable gravity waves produced by frontogenesis. Default: set by build-namelist. Whether or not to enable inertial gravity waves produced by frontogenesis. Default: set by build-namelist. Whether or not to enable gravity waves produced by deep convection. Default: set by build-namelist. Whether or not to enable gravity waves produced by shallow convection. Default: .false. Gravity wave spectrum dimension (wave numbers are from -pgwv to pgwv). Default: set by build-namelist. Width of speed bins (delta c) for gravity wave spectrum (reference wave speeds are from -pgwv*dc to pgwv*dc). Default: set by build-namelist. Dimension for long wavelength gravity wave spectrum (wave numbers are from -pgwv_long to pgwv_long). Default: set by build-namelist. Width of speed bins (delta c) for long wavelength gravity wave spectrum (reference wave speeds are from -pgwv_long*dc_long to pgwv_long*dc_long). Default: set by build-namelist. Force the stress due to gravity waves to be zero at the top of the model. In the low-top model, this helps to conserve momentum and produce a QBO. Default: set by build-namelist. Apply limiters to tau before applying the efficiency factor, rather than afterward. Default: .false. Apply limiter on maximum wind tendency from stress divergence in gravity wave drag scheme. Default: .true. Efficiency associated with convective gravity waves from the Beres scheme (deep convection). Default: set by build-namelist. Efficiency associated with convective gravity waves from the Beres scheme (shallow convection). Default: set by build-namelist. Efficiency associated with gravity waves from frontogenesis. Default: set by build-namelist. Efficiency associated with inertial gravity waves from frontogenesis. Default: set by build-namelist. Efficiency associated with orographic gravity waves. Default: set by build-namelist. Whether or not to enable gravity waves produced by meso-Beta Ridges. Default: FALSE (set by build-namelist) Number of meso-Beta ridges (per gridbox) to invoke. Default: 10 (set by build-namelist) Efficiency scaling factor associated with anisotropic OGW. Default: set by build-namelist. Max efficiency associated with anisotropic OGW. Default: 1.0 Drag coefficient for obstacles in low-level flow. Default: 1.0 Whether or not to allow trapping for meso-Beta Ridges. Default: FALSE (set by build-namelist) Whether or not to enable gravity waves produced by meso-gamma Ridges. Default: FALSE (set by build-namelist) Number of meso-gamma ridges (per gridbox) to invoke. Default: -1 (set by build-namelist) Efficiency scaling factor associated with anisotropic OGW. Default: set by build-namelist. Max efficiency associated with anisotropic OGW. Default: 1.0 Drag coefficient for obstacles in low-level flow. Default: 1.0 Whether or not to allow trapping for meso-gamma Ridges. Default: FALSE (set by build-namelist) Full pathname of boundary dataset for meso-gamma ridges. Default: set by build-namelist. Critical Froude number squared (used only for orographic waves). Default: set by build-namelist. Factor to multiply tau by, for orographic waves in the southern hemisphere. Default: 1._r8 Inverse Prandtl number used in gravity wave diffusion Default: 0.5 Scaling factor for heating depth in gravity waves from convection. If less than 1.0 this acts as an effective reduction of the gravity wave phase speeds needed to drive the QBO. Default: 1.0 Scale SGH by land fraction in gravity wave drag Default: TRUE Frontogenesis function critical threshold. Default: set by build-namelist. Full pathname of Beres lookup table data file for gravity waves sourced from deep convection. Default: set by build-namelist. Full pathname of Beres lookup table data file for gravity waves sourced from shallow convection. Default: set by build-namelist. Background source strength (used for waves from frontogenesis). Default: set by build-namelist. Background source strength (used for inertial waves from frontogenesis). Default: set by build-namelist. Whether or not to use tapering at the poles to reduce the effects of mid-scale gravity waves from frontogenesis. Default: set by build-namelist. Full pathname of time-variant boundary dataset for greenhouse gas surface values. Default: set by build-namelist. CH4 volume mixing ratio. This is used as the time invariant surface value of CH4 if no time varying values are specified. Default: set by build-namelist. CO2 volume mixing ratio. This is used as the time invariant surface value of CO2 if no time varying values are specified. Default: set by build-namelist. User override for the prescribed CO2 volume mixing ratio used by the radiation calculation. Note however that the prescribed value of CO2 which is sent to the surface models is still the one that is set using either the co2vmr or the scenario_ghg variables. Default: not used CFC11 volume mixing ratio adjusted to reflect contributions from many GHG species. This is used as the time invariant surface value of F11 if no time varying values are specified. Default: set by build-namelist. CFC12 volume mixing ratio. This is used as the time invariant surface value of CFC12 if no time varying values are specified. Default: set by build-namelist. N2O volume mixing ratio. This is used as the time invariant surface value of N2O if no time varying values are specified. Default: 0.311e-6 Data start year. Use in conjunction with ghg_yearstart_model. Default: 0 Model start year. Use in conjunction with ghg_yearstart_data. Default: 0 Amount of co2 ramping per year (percent). Only used if scenario_ghg = 'RAMP_CO2_ONLY' Default: 1.0 CO2 cap if > 0, floor otherwise. Specified as multiple or fraction of inital value; e.g., setting to 4.0 will cap at 4x initial CO2 setting. Only used if scenario_ghg = 'RAMP_CO2_ONLY' Default: boundless if ramp_co2_annual_rate > 0, zero otherwise. Date on which ramping of co2 begins. The date is encoded as an integer in the form YYYYMMDD. Only used if scenario_ghg = 'RAMP_CO2_ONLY' Default: 0 If scenario_ghg is set to "RAMPED" then the greenhouse gas surface values are interpolated between the annual average values read from the file specified by bndtvghg. In that case, the value of this variable (> 0) fixes the year of the lower bounding value (i.e., the value for calendar day 1.0) used in the interpolation. For example, if rampyear_ghg = 1950, then the GHG surface values will be the result of interpolating between the values for 1950 and 1951 from the dataset. Default: 0 Controls treatment of prescribed co2, ch4, n2o, cfc11, cfc12 volume mixing ratios. May be set to 'FIXED', 'RAMPED', 'RAMP_CO2_ONLY', or 'CHEM_LBC_FILE'. FIXED => volume mixing ratios are fixed and have either default or namelist input values. RAMPED => volume mixing ratios are time interpolated from the dataset specified by bndtvghg. RAMP_CO2_ONLY => only co2 mixing ratios are ramped at a rate determined by the variables ramp_co2_annual_rate, ramp_co2_cap, and ramp_co2_start_ymd. CHEM_LBC_FILE => volume mixing ratios are set from the chemistry lower boundary conditions dataset specified by flbc_file. Default: FIXED Full pathname of time-variant boundary dataset for greenhouse gas production/loss rates. Only used by the simple prognostic GHG chemistry scheme that is enabled via the argument "-prog_species GHG" to configure. Default: set by build-namelist. This variable should not be set by the user. It is set by build-namelist when the user specifies the argument "-prog_species GHG" to configure which turns on a simple prognostic chemistry scheme for CH4, N2O, CFC11 and CFC12. Default: set by build-namelist Flag to set rad_climate variable so that the prognostic CO2 controlled by the co2_cycle module is radiatively passive. Default: FALSE If TRUE turn on CO2 code. Default: FALSE If TRUE read co2 flux from fuel. Default: FALSE If TRUE read co2 flux from ocn. Default: FALSE Filepath for dataset containing CO2 flux from ocn. Default: none Filepath for dataset containing CO2 flux from fossil fuel. Default: none Sets the averaging flag for all variables on a particular history file series. Valid values are: A ==> Average B ==> GMT 00:00:00 average I ==> Instantaneous M ==> Minimum X ==> Maximum L ==> Local-time S ==> Standard deviation The default is to use the averaging flags for each variable that are set in the code via calls to subroutine addfld. Defaults: set in code via the addfld and add_default subroutine calls. If true don't put any of the variables on the history tapes by default. Only output the variables that the user explicitly lists in the fincl# namelist items. Default: FALSE List of fields to exclude from the 1st history file (by default the name of this file contains the string "h0"). Default: none List of fields to exclude from the 2nd history file (by default the name of this file contains the string "h1"). Default: none List of fields to exclude from the 3rd history file (by default the name of this file contains the string "h2"). Default: none List of fields to exclude from the 4th history file (by default the name of this file contains the string "h3"). Default: none List of fields to exclude from the 5th history file (by default the name of this file contains the string "h4"). Default: none List of fields to exclude from the 6th history file (by default the name of this file contains the string "h5"). Default: none List of fields to exclude from the 7th history file (by default the name of this file contains the string "h6"). Default: none List of fields to exclude from the 8th history file (by default the name of this file contains the string "h7"). Default: none List of fields to exclude from the 9th history file (by default the name of this file contains the string "h8"). Default: none List of fields to exclude from the 10th history file (by default the name of this file contains the string "h9"). Default: none List of fields to include on the first history file (by default the name of this file contains the string "h0"). The added fields must be in Master Field List. The averaging flag for the output field can be specified by appending a ":" and a valid averaging flag to the field name. Valid flags are: A ==> Average B ==> GMT 00:00:00 average I ==> Instantaneous M ==> Minimum X ==> Maximum L ==> Local-time S ==> Standard deviation Default: set in code via the addfld and add_default subroutine calls. Same as fincl1, but for the 2nd history file (by default the name of this file contains the string "h1"). Default: none. Same as fincl1, but for the 3rd history file (by default the name of this file contains the string "h2"). Default: none. Same as fincl1, but for the 4th history file (by default the name of this file contains the string "h3"). Default: none. Same as fincl1, but for the 5th history file (by default the name of this file contains the string "h4"). Default: none. Same as fincl1, but for the 6th history file (by default the name of this file contains the string "h5"). Default: none. Same as fincl1, but for the 7th history file (by default the name of this file contains the string "h6"). Default: none. Same as fincl1, but for the 8th history file (by default the name of this file contains the string "h7"). Default: none. Same as fincl1, but for the 9th history file (by default the name of this file contains the string "h8"). Default: none. Same as fincl1, but for the 10th history file (by default the name of this file contains the string "h9"). Default: none. if .true. then output CLUBBs history statistics if .true. then output CLUBBs radiative history statistics Same as fincl1, but for CLUBB statistics on zt grid. Default: none. Same as fincl1, but for CLUBB statistics on zm grid. Default: none. Same as fincl1, but for CLUBB statistics on radiation zt grid. Default: none. Same as fincl1, but for CLUBB statistics on radiation zm grid. Default: none. Same as fincl1, but for CLUBB statistics on surface. Default: none. Name of the IOP case so case specific adjustments can be made in CLUBB. Default: none. Collect all column data into a single field and output in ncol format, much faster than default when you have a lot of columns. List of columns or contiguous columns at which the fincl1 fields will be output. Individual columns are specified as a string using a longitude degree (greater or equal to 0.) followed by a single character (e)ast/(w)est identifer, an underscore '_' , and a latitude degree followed by a single character (n)orth/(s)outh identifier. For example, '10e_20n' would pick the model column closest to 10 degrees east longitude by 20 degrees north latitude. A group of contiguous columns can be specified using bounding latitudes and longitudes separated by a colon. For example, '10e:20e_15n:20n' would select the model columns which fall with in the longitude range from 10 east to 20 east and the latitude range from 15 north to 20 north. Default: none Same as fincl1lonlat, but for 2nd history file. Same as fincl1lonlat, but for 3rd history file. Same as fincl1lonlat, but for 4th history file. Same as fincl1lonlat, but for 5th history file. Same as fincl1lonlat, but for 6th history file. Same as fincl1lonlat, but for 7th history file. Same as fincl1lonlat, but for 8th history file. Same as fincl1lonlat, but for 9th history file. Same as fincl1lonlat, but for 10th history file. Specific fields which will be written using the non-default precision on the 1st history file. Default: none Specific fields which will be written using the non-default precision on the 2nd history file. Default: none Specific fields which will be written using the non-default precision on the 3rd history file. Default: none Specific fields which will be written using the non-default precision on the 4th history file. Default: none Specific fields which will be written using the non-default precision on the 5th history file. Default: none Specific fields which will be written using the non-default precision on the 6th history file. Default: none Specific fields which will be written using the non-default precision on the 7th history file. Default: none Specific fields which will be written using the non-default precision on the 8th history file. Default: none Specific fields which will be written using the non-default precision on the 9th history file. Default: none Specific fields which will be written using the non-default precision on the 10th history file. Default: none Array of history filename specifiers. The filenames of up to six history output files can be controlled via this variable. Filename specifiers give generic formats for the filenames with specific date and time components, file series number (0-5), and caseid, filled in when the files are created. The following strings are expanded when the filename is created: %c=caseid; %t=file series number (0-5); %y=year (normally 4 digits, more digits if needed); %m=month; %d=day; %s=seconds into current day; %%=% symbol. Note that the caseid may be set using the namelist variable case_name. For example, for a simulation with caseid="test" and current date and time of 0000-12-31 0:00UT, a filename specifier of "%c.cam2.h%t.%y-%m.nc" would expand into "test.cam2.h0.0000-12.nc" for the first history file. The filename specifier "%c.cam2.h%t.%y-%m-%d-%s.nc" would expand to "test.cam2.h1.0000-12-31-00000.nc" for the second history file. Spaces are not allowed in filename specifiers. Although the character "/" is allowed in the specifier, it will be interpreted as a directory name and the corresponding directories will have to be created in the model execution directory (directory given to configure with -cam_exedir option) before model execution. The first element is for the primary history file which is output by default as a monthly history file. Entries 2 through 6 are user specified auxilliary output files. Defaults: "%c.cam2.h0.%y-%m.nc", "%c.cam2.h1.%y-%m-%d-%s.nc", ..., "%c.cam2.h5.%y-%m-%d-%s.nc" Full pathname of the satellite track data used by the satellite track history output feature. Default: none Satellite track history filename specifier. See hfilename_spec Default: "%c.cam2.sat.%y-%m-%d-%s.nc" List of history fields to output along the satellite track specified by sathist_track_infile Default: none Sets the number of observation columns written the satellite track history file series. Default: 100000 Sets the number of columns closest to the observation that should be output. Setting this to a number greater than 1 allows for spatial interpolation in the post processing. Default: 1 Sets the number of timesteps closest to the observation that should be output. Setting this to a number greater than 1 allows for temporal interpolation in the post processing. Default: 1 Frequency that initial files will be output: 6-hourly, daily, monthly, yearly, or never. Valid values: 'NONE', '6-HOURLY', 'DAILY', 'MONTHLY', 'YEARLY', 'CAMIOP', 'ENDOFRUN'. Default: 'YEARLY' If false then include only REQUIRED fields on IC file. If true then include required AND optional fields on IC file. Default: FALSE Array containing the maximum number of time samples written to a history file. The first value applies to the primary history file, the second through tenth to the auxillary history files. Default: 1,30,30,30,30,30,30,30,30,30 Array containing the starting time of day for local time history averaging. Used in conjuction with lcltod_stop. If lcltod_stop is less than lcltod_start, then the time range wraps around 24 hours. The start time is included in the interval. Time is in seconds and defaults to 39600 (11:00 AM). The first value applies to the primary hist. file, the second to the first aux. hist. file, etc. Default: none Array containing the stopping time of day for local time history averaging. Used in conjuction with lcltod_start. If lcltod_stop is less than lcltod_start, then the time range wraps around 24 hours. The stop time is not included in the interval. Time is in seconds and defaults to 0 (midnight). The first value applies to the primary hist. file, the second to the first aux. hist. file, etc. Default: none Array specifying the precision of real data written to each history file series. Valid values are 1 or 2. '1' implies output real values are 8-byte and '2' implies output real values are 4-byte. Default: 2,2,2,2,2,2,2,2,2,2 Array of write frequencies for each history file series. If nhtfrq(1) = 0, the file will be a monthly average. Only the first file series may be a monthly average. If nhtfrq(i) > 0, frequency is specified as number of timesteps. If nhtfrq(i) < 0, frequency is specified as number of hours. Default: 0,-24,-24,-24,-24,-24,-24,-24,-24,-24 If interpolate_output(k) = .true., then the k'th history file will be interpolated to a lat/lon grid before output. Default: .false. Size of latitude dimension of grid for interpolated output. If interpolate_nlat and interpolate_nlon are zero, reasonable values will be chosen by the dycore based on the run resolution. Default: 0 Size of longitude dimension of grid for interpolated output. If interpolate_nlat and interpolate_nlon are zero, reasonable values will be chosen by the dycore based on the run resolution. Default: 0 Selects interpolation method for output on lat/lon grid. 0: Use SE's native high-order method. 1: Use a bilinear method. Default: 1 (bilinear) Selects output grid type for lat/lon interpolated output. 1: Equally spaced, including poles (FV scalars output grid). 2: Gauss grid (CAM Eulerian). 3: Equally spaced, no poles (FV staggered velocity). Default: 1 Full pathname of initial atmospheric state dataset (NetCDF format). Default: set by build-namelist. Perturb the initial conditions for temperature randomly by up to the given amount. Only applied for initial simulations. Default: 0.0 Full pathname of master restart file from which to branch. Setting is Required for branch run. Default: none If TRUE, try to initialize data for all consituents by reading from the initial conditions dataset. If variable not found then data will be initialized using internally-specified default values. If FALSE then don't try reading constituent data from the IC file; just use the internally-specified defaults. Default: TRUE If true, the COSP cloud simulator is run. Setting this namelist variable happens automatically if you compile with COSP. COSP will not run unless this is set to .true. in the namelist! Turn on the desired simulators using lXXX_sim namelist vars If no specific simulators are specified, all of the simulators are run on all columns and all output is saved. (useful for testing). COSP is available with CAM4 and CAM5 physics. This default logical is set in cospsimulator_intr.F90. Default: FALSE If true, COSP cloud simulators are run to produce all output required for the COSP plots in the AMWG diagnostics package. sets cosp_ncolumns=10 and cosp_nradsteps=3 (appropriate for COSP statistics derived from seasonal averages), and runs MISR, ISCCP, MODIS, CloudSat radar and CALIPSO lidar simulators (cosp_lmisr_sim=.true.,cosp_lisccp_sim=.true., cosp_lmodis_sim=.true.,cosp_lradar_sim=.true.,cosp_llidar_sim=.true.). This default logical is set in cospsimulator_intr.F90. Default: TRUE If true, the COSP cloud simulators are run to produce select output for the AMWG diagnostics package. sets cosp_ncolumns=10 and cosp_nradsteps=3 (appropriate for COSP statistics derived from seasonal averages), and runs MISR, ISCCP, MODIS, and CALIPSO lidar simulators (cosp_lmisr_sim=.true.,cosp_lisccp_sim=.true., cosp_lmodis_sim=.true.,cosp_llidar_sim=.true.). This default logical is set in cospsimulator_intr.F90. Default: FALSE If true, the passive COSP cloud simulators are run to produce select output for the AMWG diagnostics package. sets cosp_ncolumns=10 and cosp_nradsteps=3 (appropriate for COSP statistics derived from seasonal averages), and runs MISR, ISCCP, and MODIS simulators (cosp_lmisr_sim=.true.,cosp_lisccp_sim=.true.,cosp_lmodis_sim=.true.). This default logical is set in cospsimulator_intr.F90. Default: FALSE If true, the active COSP cloud simulators are run to produce select output for the AMWG diagnostics package. sets cosp_ncolumns=10 and cosp_nradsteps=3 (appropriate for COSP statistics derived from seasonal averages), and runs CloudSat radar and CALIPSO lidar simulators (cosp_lradar_sim=.true.,cosp_llidar_sim=.true.). This default logical is set in cospsimulator_intr.F90. Default: FALSE If true, the ISCCP cloud simulator is run to produce select output for the AMWG diagnostics package. sets cosp_ncolumns=10 and cosp_nradsteps=3 (appropriate for COSP statistics derived from seasonal averages), and runs ISCCP simulator (cosp_lmisr_sim=.false.,cosp_lisccp_sim=.true., cosp_lmodis_sim=.false.,cosp_lradar_sim=.false.,cosp_llidar_sim=.false.). This default logical is set in cospsimulator_intr.F90. 1236: Default: FALSE If true, run all simulators using the default values cosp_ncolumns=50 and cosp_nradsteps=1. This option is mainly intended for testing, but it also must be used in order to output the input fields needed to run the simulator in an offline mode (via setting cosp_histfile_aux=.true.). Default: FALSE If true, COSP radar simulator will be run and all non-subcolumn output will be saved. Default: FALSE If true, COSP lidar simulator will be run and all non-subcolumn output will be saved Default: FALSE If true, COSP ISCCP simulator will be run and all non-subcolumn output will be saved. ISCCP simulator is run on only daylight columns. Default: FALSE If true, MISR simulator will be run and all non-subcolumn output will be saved. MISR simulator is run on only daylight columns. Default: FALSE If true, MODIS simulator will be run and all non-subcolumn output will be saved. Default: FALSE If true, the COSP cloud simulator is run for CFMIP 3-hourly experiments. This default logical is set in cospsimulator_intr.F90 Default: FALSE If true, the COSP cloud simulator is run for CFMIP daily experiments. This default logical is set in cospsimulator_intr.F90 Default: FALSE If true, the COSP cloud simulator is run for CFMIP off-line monthly experiments. This default logical is set in cospsimulator_intr.F90 Default: FALSE If true, the COSP cloud simulator is run for CFMIP monthly experiments. This default logical is set in cospsimulator_intr.F90 Default: FALSE Number of subcolumns in SCOPS This default logical is set in cospsimulator_intr.F90 Default: 50 Turns on sampling along a-train orbit for radar and lidar simulators. This default logical is set in cospsimulator_intr.F90 Default: FALSE Path for the Atrain orbit data file provided by CFMIP. There is no default for this, but sample_atrain = .true. will not work if this namelist variable is undefined. currently /project/cms/jenkay/SCAM/repository/cloudsat_orbit_08921_14250.nc Default: no default set now, need to specify in namelist This specifies the CAM history tape where COSP diagnostics will be written. Ignored/not used if any of the cosp_cfmip_* namelist variables are invoked. This default is set in cospsimulator_intr.F90 Default: 1 If true, additional output is added to make it possible to run COSP off-line. This default is set in cospsimulator_intr.F90 Default: FALSE This specifies the CAM history tape where extra COSP diagnostics will be written. This default is set in cospsimulator_intr.F90 Default: 2 This specifies the frequency at which is COSP is called, every cosp_nradsteps radiation timestep. This default is set in cospsimulator_intr.F90 Default: 1 Turns on sub-column output from COSP. If both the isccp/misr simulators and the lidar/radar simulators are run, lfrac_out is from the isccp/misr simulators columns. This default logical is set in cospsimulator_intr.F90 Default: FALSE Number of macrophysics/microphysics substeps. Default: 1 Threshold for autoconversion of cold ice in RK microphysics scheme. Default: set by build-namelist Threshold for autoconversion of warm ice in RK microphysics scheme. Default: set by build-namelist Tunable constant for evaporation of precip in RK microphysics scheme. Default: set by build-namelist Critical radius at which autoconversion become efficient in RK microphysics scheme. Default: set by build-namelist Relative humidity threshold for stratospheric cloud water condensation in RK microphysics poleward of 50 degrees. Default: set by build-namelist Switch to control whether Park macrophysics should prognose cloud ice (cldice). Default: .true. Switch to control whether Park macrophysics should prognose cloud liquid (cldliq). Default: .true. Switch to control whether Park macrophysics should perform detrainment into the stratiform cloud scheme. Default: .true. Version number for MG microphysics Default: 1 Sub-version number for MG microphysics Default: 0 Autoconversion size threshold Default: set by build-namelist Switch to control whether MG microphysics should prognose cloud ice (cldice). Default: .true. Switch to control whether MG microphysics should prognose cloud liquid (cldliq). Default: .true. Number of substeps over MG microphysics. Default: 1 Type of precipitation fraction. Default: max_overlap Efficiency factor for berg Default: 1 Do Seifert and Behang (2001) autoconversion and accretion physics when set to true. Default: .false. Switch to control whether MG microphysics performs a uniform calculation or not (useful for sub-columns) Default: .false. prescribed aerosol bulk sulfur scale factor Default: 2 Switch to turn on heterogeneous freezing code. Default: .false. Add diagnostic output for heterogeneous freezing code. Default: .false. Switch to turn on treatment of pre-existing ice in the ice nucleation code. Default: .false. Add diagnostics for pre-existing ice option in ice nucleation code to history output. Default: .false. Subgrid scaling factor for relative humidity in ice nucleation code. Default: set by build-namelist Switch to determine whether ice nucleation happens using the incloud (true) or the gridbox average (false) relative humidity. When true, it is assumed that the incloud relative humidity for nucleation is 1. Default: set by build-namelist Fraction of Aitken mode sulfate particles assumed to nucleate ice in the polar stratospheric. Provides an increase in homogeneous freezing over the Liu&Penner method. Temporary solution to adjust ice surface area density and dehydration in the polar stratosphere where there doesn't seem to be enough nucleation. A value of zero means Liu&Penner is used. Default: set by build-namelist Characteristic adjustment time scale for Hack shallow scheme. Default: set by build-namelist Rain water autoconversion coefficient for Hack shallow scheme. Default: set by build-namelist Penetrative entrainment efficiency in UW shallow scheme. Default: set by build-namelist Switch for Vavrus "freeze dry" adjustment in cloud fraction. Set to FALSE to turn the adjustment off. Default: set by build-namelist Switch for ice cloud fraction calculation. Default: set by build-namelist Minimum rh for low stable clouds. Default: set by build-namelist Adjustment to rhminl for land without snow cover. Default: 0.10 Minimum rh for high stable clouds. Default: set by build-namelist parameter for shallow convection cloud fraction. Default: set by build-namelist parameter for shallow convection cloud fraction. Default: set by build-namelist parameter for deep convection cloud fraction. Default: set by build-namelist parameter for deep convection cloud fraction. Default: set by build-namelist top pressure bound for mid level cloud. Default: set by build-namelist Bottom height (Pa) for mid-level liquid stratus fraction. Default: 700.e2 Scheme for ice cloud fraction: 1=wang & sassen, 2=schiller (iciwc), 3=wood & field, 4=Wilson (based on smith), 5=modified slingo (ssat & empyt cloud) Default: 5 Critical RH for ice clouds (Wilson & Ballard scheme). Default: 0.93 Minimum rh for ice cloud fraction > 0. Default: set by build-namelist rhi at which ice cloud fraction = 1. Default: set by build-namelist Use cloud fraction to determine whether to do growth of ice clouds below RHice of 1 down to RHice = rhmini. Default: set by build-namelist Convective momentum transport parameter (upward) Default: set by build-namelist Convective momentum transport parameter (downward) Default: set by build-namelist Autoconversion coefficient over land in ZM deep convection scheme. Default: set by build-namelist Autoconversion coefficient over ocean in ZM deep convection scheme. Default: set by build-namelist Tunable evaporation efficiency for land in ZM deep convection scheme. Default: set by build-namelist Tunable evaporation efficiency in ZM deep convection scheme. Default: set by build-namelist Include organization parameterization in ZM Default: .false. The number of negative buoyancy regions that are allowed before the convection top and CAPE calculations are completed. Default: 5 Factor applied to the ice fall velocity computed from Stokes terminal velocity. Default: set by build-namelist Type of water vapor saturation vapor pressure scheme employed. 'GoffGratch' for Goff and Gratch (1946); 'MurphyKoop' for Murphy & Koop (2005) Default: GoffGratch Control use of sub-columns within macro/micro physics; 'false' for no subcolumns. Default: 'false' Type of sub-column generator scheme employed. 'SIHLS' Sub-columns generated with Latin Hypercube sampling of the CLUBB PDF; Type of sub-column generator scheme employed. 'CloudObj' Create sub-columns where most water is assigned to cloud sub-columns; Type of sub-column generator scheme employed. 'tstcp' testing; Type of sub-column generator scheme employed. 'vamp' Variation Across Microphysics Profiles simple deterministic scheme; 'off' for none. Default: 'off' Turns off averaging and assigns first subcolumn back to grid. Needed for BFB comparisons 'true' for no averaging. Default: '.false.' Turns on/off filtering during averaing in tstcp 'true' to use filtering. Default: '.false.' Turns on/off use of weights during averaging in tstcp 'true' to use weights. Default: '.false.' Perturbs the temperatures in state after copying for testing purposes 'true' to perturb temperatures. Default: '.false.' Tests the restart capabilities of weights with a more adequate test 'true' to set the weights to a slightly more complicated pattern for restart testing Default: '.false.' Turns on/off use of weights during averaging in tstcp 'true' to use weights. Default: '.true.' Number of subcolumns/samples to use in this simulation. Must be less than psubcols. Default: 4 Type of condensate to assume in VAMP Generator 1 Uniform Condensate 2 Variable Condensate Uniform Number 3 Variable Condensate Variable Number Default: 3 Type of overlap to assume in VAMP Generator 1 Maximum Default: 1 Number of subcolumns in VAMP Generator Default: 10 Type of deep convection scheme employed. 'ZM' for Zhang-McFarlane; 'off' for none; or 'UNICON' which doesn't distinquish shallow and deep. Default: 'ZM' unless using 'UNICON' Type of microphysics scheme employed. 'RK' for Rasch and Kristjansson (1998); 'MG' for Morrison and Gettelman (2008), Gettelman et al (2010) two moment scheme for CAM5 and CAM6 SPCAM has two different microphysics schemes: SPCAM_m2005 (Morrison et al 2005), SPCAM_sam1mom (Khairoutinov 2003) Default: set by build-namelist (depends on value set in configure). Type of macrophysics scheme employed. 'park' for Park (1998); 'RK' for Rasch and Kristjansson (1998); 'CLUBB_SGS' clubb. Default: set by build-namelist Switch for CLUBB_SGS Default: FALSE Type of shallow convection scheme employed. 'Hack' for Hack shallow convection; 'UW' for original McCaa UW pbl scheme, modified by Sungsu Park; 'CLUBB_SGS' for CLUBB_SGS; or UNICON which doesn't distinquish shallow and deep. Default: set by build-namelist (depends on eddy_scheme). Logical switch to turn on the beljaars scheme Default: set by build-namelist Logical switch to turn on turbulent mountain stress calculation in vertical diffusion routine. Default: set by build-namelist Turbulent mountain stress parameter used when turbulent mountain stress calculation is turned on. See do_tms. Default: 1.0 for CAM, set by build-namelist for WACCM, T31 Factor determining z_0 from orographic standard deviation [ no unit ] Used when turbulent mountain stress calc is turned on. See do_tms. set by build-namelist for WACCM, T31 Maximum master length scale designed to address issues in diag_TKE outside the boundary layer. In order not to disturb turbulence characteristics in the lower troposphere, this should be set at least larger than a few km. However, this does not significantly improve the values outside of the boundary layer. Smaller values make some improvement, but it is also noisy. Better results are seen using eddy_leng_max or kv_freetrop_scale. Default: 40.e3 (m) Maximum dissipation length scale designed to address issues with diag_TKE outside the boundary layer, where the default value generates large diffusivities. A value of 30 m is consistent with the length scales used in the HB scheme; however, this will also reduce value in the boundary layer. Default: 40.e3 (m) Bottom pressure level at which namelist values for eddy_leng_max and eddy_lbulk_max are applied. Default values are used at lower levels (i.e. the boundary layer). Default: 100.e3 (hPa) Moist entrainment enhancement parameter. Default: set by build-namelist Pressure (Pa) that defined the upper atmosphere for adjustment of eddy diffusivities from diag_TKE using kv_top_scale. Default: 0. Scaling factor that is applied (multiplied) to the eddy diffusivities in the upper atmosphere (see kv_top_pressure). Default: 1.0 Scaling factor that is applied (multiplied) to the eddy diffusivities in the free troposphere (boundary layer to kv_top_pressure) Default: 1.0 Perform mass conservation check on eddy diffusion operation. Default: FALSE Logical switch to turn on implicit turbulent surface stress calculation in diffusion solver routine. Default: set by build-namelist Produce output for the offline unicon driver. Default: .false. History file number for offline unicon driver output. Default: 2 (i.e., h1 history file) Apply cloud top radiative cooling parameterization Default: .false. Include effects of precip evaporation on turbulent moments Default: .false. Explicit diffusion on temperature and moisture when CLUBB is on Default: .false. CLUBB do explicit diffusion with a stability correction Default: .false. CLUBB timestep. Default: set by build-namelist Rain evaporation efficiency factor. Default: set by build-namelist Switch for CLUBB_ADV Default: FALSE Low Skewness in C11 Skw. Function High Skewness in C11 Skw. Function Coef. applied to log(avg dz/thresh) Low Skw.: gamma coef. Skw. Fnct. Momentum coefficient of Kh_zm Thermo of Kh_zm Plume widths for theta_l and rt C2 coef. for the rtp2_dp1 term C2 coef. for the thlp2_dp1 term C2 coef. for the rtpthlp_dp1 term Coef. #1 in C8 Skewness Equation Low Skewness in C7 Skw. Function High Skewness in C7 Skw. Function Factor to decrease sensitivity in the denominator of Skw calculation Intensity of stability correction applied to C1 and C6 Uses PDF to compute perturbed values for l_avg_Lscale code Include the effects of ice latent heating in turbulence terms Apply liquid supersaturation adjustment code Default: false The name of the active CARMA microphysics model or none when CARMA is not active. Default: none A fraction that scales how tight the convergence criteria are to determine that the substepping has resulted in a valid solution. Smaller values will force more substepping. CARMA particles. Default: 0.1 When non-zero, the largest change in temperature (K) allowed per substep. Default: 0. Flag indicating that the CARMA model is an aerosol model, and should be called in tphysac. Default: TRUE Flag indicating that CARMA is a cloud ice model and should be called in tphysbc. Default: FALSE Flag indicating that CARMA is a cloud liquid model and should be called in tphysbc. Default: FALSE Flag indicating that CARMA should do clear sky calculations for particles that are not part of a cloud in addition to doing a separate calculation for incloud particles. Only valid when carma_do_incloud is true. Default: FALSE Flag indicating whether the coagulation process is enabled for CARMA particles. Default: FALSE Flag indicating that CARMA is responsible for detrain condensate from convection into the model. Default: FALSE Flag indicating that the dry deposition process is enabled for CARMA particles. Default: FALSE Flag indicating that the emission of particles is enabled for CARMA. Default: FALSE Flag indicating that sedimentation should be calculated using an explicit technique where the substepping is used to keep the CFL condition from being violated rather than the default PPM scheme. Default: FALSE Flag indicating CARMA coefficients should only be initialized once from a fixed temperature profile rather than recomputed for each column. This improves performance, but reduces accuracy. By default the temperature profile used is calculated as the average of the initial condition file, but a predefined profile can be provided. Default: FALSE Flag indicating used in cunjunction with carma_do_fixedinit to indicate that only the coagulation coefficients should only be initialized from a fixed temperature profile and all other coeeficients will be recalculated. Coagulation is the slowest initialization, so this improves performance while still retaining accuracy for most processes. Default: FALSE Flag indicating that the condensational growth process is enabled for CARMA particles. Default: FALSE Flag indicating that CARMA sulfate mass mixing ratio will be used in radiation calculation. Default: FALSE Flag indicating that CARMA sulfate surface area density will be used in heterogeneous chemistry rate calculation. Default: FALSE Flag indicating that CARMA should treat cloud particles as incloud rather than gridbox average calculations. Default: FALSE Flag indicating that carma should generate optical properties files for the CAM radiation code. Default: FALSE Flag indicating that particle heating will be used for the condensational growth process. Default: FALSE Flag indicating that particle heating will affect the atmospheric temperature. Default: FALSE Flag indicating that substepping will be used for the condensational growth process. Default: FALSE Flag indicating that changes in heating will be calculated as a result CARMA processes and will affect the CAM heating tendency. Default: FALSE Flag indicating that the wet deposition process is enabled for CARMA particles. Default: FALSE Flag indicating that the effect of Brownian diffusion will be calculated for CARMA particles. NOTE: This needs to be used in conjunction with CARMA sedimentation. Default: FALSE Flag indicating that the sedimentation process is enabled for CARMA particles. Default: FALSE Flag indicating whether CARMA is enabled. If CARMA has been included in the build (configure -carma with something other than none), then this will cause all of the CARMA constituents and field names to be registered, but no other CARMA process will be preformed. This overrides the individual CARMA process flags. Default: FALSE Specifies the maximum number of retry attempts to be used when condensational growth requires substepping, but the original estimate for the amount of substepping was insufficient. Default: 8 Specifies the maximum number of substeps that could be used for the first guess when condensational growth requires substepping. Default: 1 Specifies the name of the reference temperature file that will be used (and created if necessary) for initialization of CARMA to a fixed temperature profile. Default: carma_reft.nc Accommodation coefficient for coagulation. Default: 1.0 Accommodation coefficient for growth with ice. Default: 0.93 Accommodation coefficient for growth with liquid. Default: 1.0 Accommodation coefficient for temperature. Default: 1.0 Critical relative humidity for liquid cloud formation, used for sub-grid scale in-cloud saturation. Default: 1.0 Global mass of dust emission for the event. Default: 0. (kg) Global mass of dust emission for the event. Default: 0. (kg) Starting date for emissions in the form of (yyyyddd) where yyyy is a year and ddd is a day of year. Default: 1 (yyyyddd) Starting time for the emission event in GMT. Default: 0. (s Z) Stopping date for emissions in the form of (yyyyddd) where yyyy is a year and ddd is a day of year. Default: 1 (yyyyddd) Stoping time for the emission event in GMT. Default: 0. (s) Minimum latitude of the area for emssions from the event. Default: -90. (degrees north) Maximum latitude of the area for emssions from the event. Default: 90. (degrees north) Minimum longitude of the area for emssions from the event. Default: 0. (degrees east) Maximum longitude of the area for emssions from the event. Default: 360. (degrees east) Are the soot particles treated as fractals? Default: FALSE Flag indicating that meteor smoke emission will be scaled by a global relative flux based upon the carma_escale_file. Default: FALSE The total meteor smoke emission rate in kt/year. The flux will be scaled to total that value. Default: 16.0 Specifies the name of the file containing the meteor smoke emission (ablation) profile. Default: set by build-namelist. Specifies the name of the file containing the global realtive flux specification. Default: set by build-namelist. Specifies the day of year when tracers will start being emitted for the tracer test. Default: 1 The emission rate of inert tracers used in the test. A positive value indicates that the rate is a column mass (kg/m2/s) and a negative value indicate that it is a mass mixing ratio (kg/kg/s). Default: 1e-9 Flag indicating that h2so4 vapor pressures should be calculated as if they were over sulfates that have been totally neutralized. Default: FALSE Specifies the method to use to get the prescribed sulfate aerosols for use with nucleation of cirrus clouds. This can be different than the sulfate aerosols that are used with the climate. Default: fixed Specifies the name of the file containing ice refrative indicies as a function of wavelength used for the particle heating calculation. Default: set by build-namelist. Specifies the name of the file containing soil erosion factors. This is used by the dust model. Default: set by build-namelist. Flag indicating that a calculated Weibull K should be used. Default: FALSE Specifies the name of the sea salt emission parameterization. Default: Gong ======= Full pathname of time-variant ozone mixing ratio boundary dataset. Default: set by build-namelist. Add CAM3 prescribed ozone to the physics buffer. Default: FALSE Flag for yearly cycling of ozone data. If set to FALSE, a multi-year dataset is assumed, otherwise a single-year dataset is assumed, and ozone will be cycled over the 12 monthly averages in the file. Default: TRUE String identifying a hardware counter to the papi library. Default: PAPI_TOT_CYC String identifying a hardware counter to the papi library. Default: PAPI_FP_OPS String identifying a hardware counter to the papi library. Default: PAPI_FP_INS String identifying a hardware counter to the papi library. Default: PAPI_NO_CTR Flag indicating whether the mpi_barrier in t_barrierf should be called. Default: FALSE Maximum number of levels of timer nesting . Default: 99999 Maximum detail level to profile. Default: 1 Flag indicating whether timers are disabled. Default: FALSE Collect and print out global performance statistics (for this component communicator). Default: FALSE Maximum number of processes writing out timing data (for this component communicator). Default: -1 Separation between process ids for processes that are writing out timing data (for this component communicator). Default: 1 Flag indicating whether the PAPI namelist should be read and HW performance counters used in profiling. Default: FALSE Flag indicating whether the performance timer output should be written to a single file (per component communicator) or to a separate file for each process. Default: TRUE Initialization of GPTL timing library. Default: GPTLmpiwtime Swap communication protocol option (reduced set): 3, 5: nonblocking send 2, 3, 4, 5: nonblocking receive 4, 5: ready send Default: 4 Swap communication maximum request count: <=0: do not limit number of outstanding send/receive requests >0: do not allow more than swap_comm_maxreq outstanding nonblocking send requests or nonblocking receive requests Default: 128 fc_gather flow control option: < 0 : use MPI_Gather >= 0: use point-to-point with handshaking messages and preposting receive requests up to max(min(1,fc_gather_flow_cntl),max_gather_block_size) ahead. Default value is defined by private parameter max_gather_block_size, which is currently set to 64. Default: 64 Allocate all buffers as global. This is a performance optimization on machines for which allocation/deallocation of physpkg scope buffers on every timestep was slow (Cray-X1). Default: TRUE Name of the CAM physics package. N.B. this variable may not be set by the user. It is set by build-namelist via information in the configure cache file to be consistent with how CAM was built. Default: set by build-namelist Name of the CAM chemistry package. N.B. this variable may not be set by the user. It is set by build-namelist via information in the configure cache file to be consistent with how CAM was built. Default: set by build-namelist Runtime options of upper thermosphere WACCM-X. 'ionosphere' for full ionopshere and neutral thermosphere, 'neutral' for just neutral thermosphere, and off for no WACCM-X. Default: 'off' Troposphere cloud physics will be done only below the top defined by this pressure (Pa). Default: set by build-namelist MAM affects climate only below the top defined by this pressure (Pa). Default: 0 for non-MAM cases, otherwise set by build-namelist Molecular diffusion will be done only if the lowest pressure is below this limit (Pa). Default: 0.1 The level closest to this pressure (Pa) is the bottom of the region where molecular diffusion is done. Default: 50. Use this variable to specify the latitude (in degrees) of a column to debug. The closest column in the physics grid will be used. Default: none Use this variable to specify the longitude (in degrees) of a column to debug. The closest column in the physics grid will be used. Default: none If set to .true., turns on extra validation of physics_state objects in physics_update. Used mainly to track down which package is the source of invalid data in state. Default: .false. Switch to turn on adjustment of the surface fluxes to reduce instabilities in the surface layer. Set to 1 to turn on the adjustments. Default: 0 if eddy_scheme is 'HB', 1 otherwise. Type of eddy scheme employed by the vertical diffusion package. 'HB' for Holtslag and Boville; 'diag_TKE' for diagnostic tke version of Grenier and Bretherton; 'HBR' for Rasch modified version of 'HB'. Default: 'HB' Switch to use new convective scavenging for modal aerosols. This scheme replaces the call to ZM's convtran for the the modal aerosol number and mass mixing ratio constituents. Default: .false. Produce output for the AMWG diagnostic package. Default: .true. Produce output for the AMWG variability diagnostics. Default: .false. Switch for diagnostic output of the aerosol tendencies Default: .false. Switch for diagnostic output of the aerosol optics Default: .false. Switch for diagnostic output of eddy variables Default: .false. Switch for cam4 T/Q budget diagnostic output Default: .false. History tape number T/Q budget output is written to. Default: 1 Switch for diagnostic output used primarily for WACCM runs. Default: .true. if WACCM physics is on, .false. otherwise. Switch for diagnostic output used primarily for WACCM-X runs. Default: .true. if WACCM-X is on, .false. otherwise. Switch for diagnostics specific to the current chemistry package or configuration. Default: .true. Switch for diagnostics specific to the current CARMA model. Default: .true. Switch for diagnostics specific to CLUBB. Default: .true. True when model is configured to use an offline driver. Default: Set by build-namelist. Type of radiation scheme employed. Default: set by build-namelist Convective water used in radiation? 0 ==> No 1 ==> Yes - Arithmetic average. 2 ==> Yes - Average in emissivity. Default: set by build-namelist Lower limit of cumulus cloud fraction. Default: set by build-namelist Full pathname of absorption/emission dataset. Used only by camrt scheme. It consists of terms used for determining the absorptivity and emissivity of water vapor in the longwave parameterization of radiation. Default: set by build-namelist. Frequency of absorptivity/emissivity calculations in time steps (if positive) or model hours (if negative). To avoid having the abs/ems values saved on the restart output, make sure that the interval of the abs/ems calculation evenly divides the restart interval. Default: -12 Frequency of long-wave radiation calculation in timesteps (if positive) or model hours (if negative). Default: -1 Frequency of short-wave radiation calculation in timesteps (if positive) or model hours (if negative). Default: -1 Specifies length of time in timesteps (positive) or hours (negative) SW/LW radiation will be run for every timestep from the start of an initial run. Default: 0 Return fluxes per band in addition to the total fluxes. Default: FALSE If true, then average the zenith angle over the radiation timestep rather than using instantaneous values. Default: FALSE Definitions for the aerosol modes that may be used in the rad_climate and rad_diag_* variables. Default: set by build-namelist A list of the radiatively active species, i.e., species that affect the climate simulation via the radiative heating rate calculation. Default: set by build-namelist A list of species to be used in the first diagnostic radiative heating rate calculation. These species are not the ones affecting the climate simulation. This is a hook for performing radiative forcing calculations. Default: none Analogous to rad_diag_1, but for the 2nd diagnostic calculation. Default: none Analogous to rad_diag_1, but for the 3rd diagnostic calculation. Default: none Analogous to rad_diag_1, but for the 4th diagnostic calculation. Default: none Analogous to rad_diag_1, but for the 5th diagnostic calculation. Default: none Analogous to rad_diag_1, but for the 6th diagnostic calculation. Default: none Analogous to rad_diag_1, but for the 7th diagnostic calculation. Default: none Analogous to rad_diag_1, but for the 8th diagnostic calculation. Default: none Analogous to rad_diag_1, but for the 9th diagnostic calculation. Default: none Analogous to rad_diag_1, but for the 10th diagnostic calculation. Default: none output data needed for off-line radiation calculations Default: FALSE History tape number radiation driver output data is written to. Default: 0 Averaging flag for adiation driver output data. Default: 'A' Switch to turn on Fixed Dynamical Heating in the offline radiation tool (PORT). Default: false Full pathname of dataset for water refractive indices used in modal aerosol optics Default: none Dry deposition surface values interpolated to model grid, required for unstructured atmospheric grids with modal chemistry. Default: none filepath and name for ice optics data for rrtmg Default: none filepath and name for ice optics data for rrtmg Default: none filepath and name for ice optics data for rrtmg Default: none filepath and name for ice optics data for rrtmg Default: none filepath and name for liquid cloud (gamma distributed) optics data for rrtmg Default: none Variable to specify the vertical index at which the Rayleigh friction term is centered (the peak value). Default: 2 Rayleigh friction parameter to determine the width of the profile. If set to 0 then a width is chosen by the algorithm (see rayleigh_friction.F90). Default: 0. Rayleigh friction parameter to determine the approximate value of the decay time (days) at model top. If 0.0 then no Rayleigh friction is applied. Default: 0. Full pathname of IOP dataset. Default: set by build-namelist. Column radiation mode. Default: FALSE If using diurnal averaging or not. Default: FALSE Use the specified surface properties. Default: FALSE IOP name for CLUBB running in single column mode Default: "" Use relaxation. Default: FALSE Total solar irradiance (W/m2). Default: set by build-namelist Full pathname of dataset for file that contains the solar photon enerspectra or TSI data as a time series Default: none Full pathname of dataset for file that contains the solar EUV data as a time series Default: none Type of time interpolation for data in solar_data_file. Can be set to "FIXED" or "SERIAL". Default: SERIAL If solar_data_type is "FIXED" then solar_data_ymd is the date the solar data is fixed to. If solar_data_type is "SERIAL" the solar_data_ymd is the start date of the time series of solar data. Format: YYYYMMDD Default: none Seconds of the day corresponding to solar_data_ymd Default: current model time of day Use spectral scaling in the radiation heating Default: false User can specify names for test tracers to be read from the initial file. The number of names specified should be given as the value of the -nadv_tt option to configure. Default: '' This variable should not be set by the user. If configure has been invoked with the '-nadv_tt N' option then build-namelist will set this variable to the value N. If test_tracer_names have been specified then N should be the number of names supplied. If test_tracer_names have not been specified, then the tracer_suite module generates the tracer names and supplies the initial values. Default: set by configure If true age of air tracers are included. This variable should not be set by the user. It will be set by build-namelist to be consistent with the '-age_of_air_trcs' argument specified to configure. Default: set by configure If true age of air tracers are read from the initial conditions file. If this is not specified then they are not read from IC file. Default: TRUE The length (in seconds) of the atm time step, i.e., the driver calls the atm component once every dtime seconds. This is also the coupling interval between the dynamics and physics packages. This variable is not actually used in the atm model, but rather is used by build-namelist to set the value of atm_cpl_dt. So it will have an effect only when running CAM using standalone scripts. The CESM scripts have their own method for setting atm_cpl_dt. Default: is resolution and dycore dependent and is set by build-namelist. Full pathname of time-invariant boundary dataset for topography fields. Default: set by build-namelist. Setting use_topo_file=.false. allows the user to specify that PHIS, SGH, SGH30, and LANDM_COSLAT are all zero without having to supply a topo file full of zeros. Default: TRUE Full pathname of boundary dataset for tropopause climatology. Default: set by build-namelist. Flag to tell build-namelist to use time-dependent external forcing files for the aircraft emissions. Default: FALSE Flag to set rad_climate variable so that the chemical tracers are radiatively passive. Default: FALSE Wet deposition method used MOZ --> mozart scheme is used NEU --> J Neu's scheme is used OFF --> wet deposition is turned off Default: MOZ List of gas-phase species that undergo wet deposition via the wet deposition scheme. Default: NONE Turns on accumulation to coarse mode exchange appropriate for the stratosphere. This also changes the default mode definitions (widths and edges) via default aerosol property input files. Default: FALSE Turns on prognostic modal sulfate aerosols in the stratosphere. Default: FALSE List of aerosol species that undergo wet deposition. Default: set by build-namelist. In-cloud solubility factor used in SO4 wet removal Default: set by build-namelist. Below-cloud solubility factor used in SO4 wet removal Default: set by build-namelist. In-cloud solubility factor used in NH4 wet removal Default: set by build-namelist. Below-cloud solubility factor used in NH4 wet removal Default: set by build-namelist. In-cloud solubility factor used in NH4NO3 wet removal Default: set by build-namelist. Below-cloud solubility factor used in NH4NO3 wet removal Default: set by build-namelist. In-cloud solubility factor used in CB2 wet removal Default: set by build-namelist. Below-cloud solubility factor used in CB2 wet removal Default: set by build-namelist. In-cloud solubility factor used in OC2 wet removal Default: set by build-namelist. Below-cloud solubility factor used in OC2 wet removal Default: set by build-namelist. In-cloud solubility factor used in wet removal of BULK dust Default: set by build-namelist. Below-cloud solubility factor used in wet removal of BULK dust Default: set by build-namelist. In-cloud solubility factor used in wet removal of BULK sea salt Default: set by build-namelist. Below-cloud solubility factor used in wet removal of BULK sea salt Default: set by build-namelist. List of aerosol species that undergo sediment (dry deposition). Default: set by build-namelist. Tuning for below cloud scavenging of interstitial modal aerosols. Default: set by build-namelist. Tuning for in-cloud scavenging of interstitial modal aerosols. Default: set by build-namelist. Tuning for in-cloud scavenging of cloud-borne modal aerosols. Default: set by build-namelist. Tuning for seasalt_emis Default: set by build-namelist. In-cloud solubility factor used in BULK aerosol wet removal Default: set by build-namelist. Below-cloud solubility factor used in BULK aerosol wet removal Default: set by build-namelist. Scavenging coefficient used in BULK aerosol wet removal Default: set by build-namelist. Scavenging coefficient used in the wet removal of SO4 Default: set by build-namelist. Scavenging coefficient used in the wet removal of NH4 Default: set by build-namelist. group="aerosol_nl" valid_values="" > Scavenging coefficient used in the wet removal of NH4NO3 Default: set by build-namelist. Scavenging coefficient used in the wet removal of CB2 Default: set by build-namelist. Scavenging coefficient used in the wet removal of OC2 Default: set by build-namelist. Scavenging coefficient used in the wet removal of DST01 Default: set by build-namelist. Scavenging coefficient used in the wet removal of DST02 Default: set by build-namelist. Scavenging coefficient used in the wet removal of DST03 Default: set by build-namelist. Scavenging coefficient used in the wet removal of DST04 Default: set by build-namelist. Scavenging coefficient used in the wet removal of SSLT01 Default: set by build-namelist. Scavenging coefficient used in the wet removal of SSLT02 Default: set by build-namelist. Scavenging coefficient used in the wet removal of SSLT03 Default: set by build-namelist. Scavenging coefficient used in the wet removal of SSLT04 Default: set by build-namelist. Full pathname of boundary dataset for airplane emissions. Default: set by build-namelist. Full pathname of dataset containing soil moisture fraction information used in 'xactive_atm' method of calculating dry deposition of chemical tracers. Default: set by build-namelist. Full pathname of dataset which contains the prescribed deposition velocities used in the 'table' method of calculating dry deposition of chemical tracers. Default: set by build-namelist. Full pathname of dataset which contains land vegitation information used in 'xactive_atm' method of calculating dry deposition of chemical tracers. Default: set by build-namelist. Full pathname of dataset which contains season information used in 'xactive_atm' method of calculating dry deposition of chemical tracers. Default: set by build-namelist. Tuning parameter for dust emissions. Default: set by build-namelist. Full pathname of dataset for coefficient data used in WACCM to calculate ion drag for high solar fluxes from the Scherliess low latitude electric potential model. Default: set by build-namelist. Full pathname of dataset for coefficient data used in WACCM to calculate ion drag for low solar fluxes from the Scherliess low latitude electric potential model. Default: set by build-namelist. Full pathname of dataset for coefficient data used in WACCM to calculate ion drag from the Weimer96 high latitude electric potential model. Default: set by build-namelist. Full pathname of dataset for the neutral species absorption cross sections for EUV photo reactions producing electrons. Default: set by build-namelist. Type of time interpolation of emission datasets specified. Can be set to 'CYCLICAL', 'SERIAL', 'INTERP_MISSING_MONTHS', or 'FIXED'. by srf_emis_specifier. Default: 'CYCLICAL' The cycle year of the surface emissions data if srf_emis_type is 'CYCLICAL'. Format: YYYY Default: 0 The date at which the surface emissions are fixed if srf_emis_type is 'FIXED'. Format: YYYYMMDD Default: 0 The time of day (seconds) corresponding to srf_emis_fixed_ymd at which the surface emissions are fixed if srf_emis_type is 'FIXED'. Default: 0 seconds Full pathname of dataset for EUVAC solar EUV model (0.05-121nm). Default: set by build-namelist. The cycle year of the external forcings (3D emissions) data if ext_frc_type is 'CYCLICAL'. Format: YYYY Default: 0 Default: current model date The date at which the external forcings are fixed if ext_frc_type is 'FIXED'. Format: YYYYMMDD Default: 0 The time of day (seconds) corresponding to ext_frc_fixed_ymd at which the external forcings are fixed if ext_frc_type is 'FIXED'. Default: 0 seconds List of full pathnames of elevated emission (or external chemical forcings) datasets. The chemistry package reads in elevated emission data from a set of netcdf files in units of "molecules/cm3/s". Each tracer species emissions is read from its own file as directed by the namelist variable ext_frc_specifier. The ext_frc_specifier variable tells the model which species have elevated emissions and the file path for the corresponding species. That is, the ext_frc_specifier variable is set something like: ext_frc_specifier = 'SO2 -> /path/vrt.emis.so2.nc', 'SO4 -> /path/vrt.emis.so4.nc', etc... Each emission file can have more than one source. When the emission are read in the sources are summed to give a total emission field for the corresponding species. The emission can be read in as time series of data, cycle over a given year, or be fixed to a given date. The vertical coordinate in these emissions files should be 'altitude' (km) so that the vertical redistribution to the model layers is done using a mass conserving method. If the vertical coordinate is altitude then data needs to be ordered from the surface to the top (increasing altitude). Default: set by build-namelist. Type of time interpolation for fixed lower boundary data. Can be set to 'CYCLICAL', 'SERIAL', 'INTERP_MISSING_MONTHS', or 'FIXED'. Default: 'CYCLICAL' The cycle year of the fixed lower boundary data if flbc_type is 'CYCLICAL'. Format: YYYY Default: 0 The date at which the fixed lower boundary data is fixed if flbc_type is 'FIXED'.. Format: YYYYMMDD Default: 0 The time of day (seconds) corresponding to flbc_fixed_ymd at which the fixed lower boundary data is fixed if flbc_type is 'FIXED'. Default: 0 seconds Full pathname of dataset for fixed lower boundary conditions. Default: set by build-namelist. List of species that are fixed at the lower boundary. Default: set by build-namelist. Type of time interpolation for fixed lower boundary data. Default: 'CYCLICAL' File name of dataset for NOy upper boundary conditions. Default: set by build-namelist. Full pathname of the directory that contains the NOy upper boundary conditions files specified in noy_ubc_filelist. Default: set by build-namelist. Filename of file that contains a sequence of filenames for prescribed aerosols. The filenames in this file are relative to the directory specied by prescribed_aero_datapath. Default: set by build-namelist. Type of time interpolation for NOy upper boundary conditions. Default: 'SERIAL' The cycle year of the NOy upper boundary data if flbc_type is 'CYCLICAL'. Format: YYYY Default: 0 The date at which the NOy upper boundary data is fixed if flbc_type is 'FIXED'.. Format: YYYYMMDD Default: 0 The time of day (seconds) corresponding to noy_ubc_fixed_ymd at which the NOy upper boundary data is fixed if flbc_type is 'FIXED'. Default: 0 seconds Full pathname of dataset for chemical tracers constrained in the stratosphere Default: set by build-namelist. List of species that are constrained in the stratosphere. Default: set by build-namelist. Full pathname of dataset for land mask applied to the lighting NOx production Default: set by build-namelist. Multiplication factor applied to the lighting NOx production Default: 1.0. Multiplication factor applied to the upper boundary NO mass mixing ratio. Default: 1.0 Full pathname of dataset for the neutral species absorption cross sections. Default: set by build-namelist. Full pathname of dataset for fast-tuv photolysis cross sections Default: set by build-namelist. Full pathname of dataset of O2 cross sections for fast-tuv photolysis Default: set by build-namelist. Full pathname of dataset of O2 and 03 column densities above the model for look-up-table photolysis Default: set by build-namelist. Full pathname of the aircraft input file list Default: none Type of time interpolation for data in aircraft aerosol files. Default: 'CYCLICAL_LIST' Full pathname of the directory that contains the files specified in gcr_ionization_filelist. Default: set by build-namelist. Filename of dataset for ionization rates by galactic cosmic rays. Default: set by build-namelist. Filename of file that contains a sequence of filenames for ionization rates by galactic cosmic rays. The filenames in this file are relative to the directory specied by gcr_ionization_datapath. Default: set by build-namelist. Names of variables containing ionization rates (/cm3/sec) in the cosmic rays datasets. Default: none Type of time interpolation for data in gcr_ionization files. Can be set to 'CYCLICAL', 'SERIAL', 'INTERP_MISSING_MONTHS', or 'FIXED'. Default: 'SERIAL' The cycle year of the prescribed green house gas data if gcr_ionization_type is 'CYCLICAL'. Format: YYYY Default: 0 The date at which the prescribed green house gas data is fixed if gcr_ionization_type is 'FIXED'. Format: YYYYMMDD Default: 0 The time of day (seconds) corresponding to gcr_ionization_fixed_ymd at which the prescribed green house gas data is fixed if gcr_ionization_type is 'FIXED'. Default: 0 seconds Full pathname of the directory that contains the files specified in prescribed_aero_filelist. Default: set by build-namelist. Switch used to indicate which type of aerosols are prescribed -- bulk or modal. This is used to set the default prescribed_aero_specifier and aerodep_flx_specifier namelist variables. Default: set by build-namelist Filename of dataset for prescribed aerosols. Default: set by build-namelist. Filename of file that contains a sequence of filenames for prescribed aerosols. The filenames in this file are relative to the directory specied by prescribed_aero_datapath. Default: set by build-namelist. Remove the file containing prescribed aerosol concentrations from local disk when no longer needed. Default: FALSE A list of variable names of the concentration fields in the prescribed aerosol datasets and corresponding names used in the physics buffer seperated by colons. For example: prescribed_aero_specifier = 'pbuf_name1:ncdf_fld_name1','pbuf_name2:ncdf_fld_name2', ... If there is no colon seperater then the specified name is used as both the pbuf_name and ncdf_fld_name, Default: none Type of time interpolation for data in prescribed_aero files. Can be set to 'CYCLICAL', 'SERIAL', 'INTERP_MISSING_MONTHS', or 'FIXED'. Default: 'SERIAL' The cycle year of the prescribed aerosol data if prescribed_aero_type is 'CYCLICAL'. Format: YYYY Default: 0 The date at which the prescribed aerosol data is fixed if prescribed_aero_type is 'FIXED'. Format: YYYYMMDD Default: 0 The time of day (seconds) corresponding to prescribed_aero_fixed_ymd at which the prescribed aerosol data is fixed if prescribed_aero_type is 'FIXED'. Default: 0 seconds Full pathname of the directory that contains the files specified in aerodep_flx_filelist. Default: set by build-namelist. Filename of dataset for prescribed aerosols. Default: set by build-namelist. Filename of file that contains a sequence of filenames for prescribed aerosols. The filenames in this file are relative to the directory specied by aerodep_flx_datapath. Default: set by build-namelist. Remove the file containing prescribed aerosol deposition fluxes from local disk when no longer needed. Default: FALSE Names of variables containing aerosol data in the prescribed aerosol datasets. Default: none Type of time interpolation for data in aerodep_flx files. Can be set to 'CYCLICAL', 'SERIAL', 'INTERP_MISSING_MONTHS', or 'FIXED'. Default: 'SERIAL' The cycle year of the prescribed aerosol flux data if aerodep_flx_type is 'CYCLICAL'. Format: YYYY Default: 0 The date at which the prescribed aerosol flux data is fixed if aerodep_flx_type is 'FIXED'. Format: YYYYMMDD Default: 0 The time of day (seconds) corresponding to ssaerodep_flx_fixed_ymd at which the prescribed aerosol flux data is fixed if saerodep_flx_type is 'FIXED'. Default: 0 seconds Full pathname of the directory that contains the files specified in prescribed_ghg_filelist. Default: set by build-namelist. Filename of dataset for prescribed GHGs. Default: set by build-namelist. Filename of file that contains a sequence of filenames for prescribed GHGs. The filenames in this file are relative to the directory specied by prescribed_ghg_datapath. Default: set by build-namelist. Remove the file containing prescribed green house gas concentrations from local disk when no longer needed. Default: FALSE Names of variables containing GHG data in the prescribed GHG datasets. Default: none Type of time interpolation for data in prescribed_ghg files. Can be set to 'CYCLICAL', 'SERIAL', 'INTERP_MISSING_MONTHS', or 'FIXED'. Default: 'SERIAL' The cycle year of the prescribed green house gas data if prescribed_ghg_type is 'CYCLICAL'. Format: YYYY Default: 0 The date at which the prescribed green house gas data is fixed if prescribed_ghg_type is 'FIXED'. Format: YYYYMMDD Default: 0 The time of day (seconds) corresponding to prescribed_ghg_fixed_ymd at which the prescribed green house gas data is fixed if prescribed_ghg_type is 'FIXED'. Default: 0 seconds Full pathname of the directory that contains the files specified in prescribed_ozone_filelist. Default: set by build-namelist. Filename of dataset for prescribed ozone. Default: set by build-namelist. Filename of file that contains a sequence of filenames for prescribed ozone. The filenames in this file are relative to the directory specied by prescribed_ozone_datapath. Default: set by build-namelist. Name of variable containing ozone data in the prescribed ozone datasets. Default: 'ozone' Remove the file containing prescribed ozone concentrations from local disk when no longer needed. Default: FALSE Type of time interpolation for data in prescribed_ozone files. Can be set to 'CYCLICAL', 'SERIAL', 'INTERP_MISSING_MONTHS', or 'FIXED'. Default: 'SERIAL' The cycle year of the prescribed ozone data if prescribed_ozone_type is 'CYCLICAL'. Format: YYYY Default: 0 The date at which the prescribed ozone data is fixed if prescribed_ozone_type is 'FIXED'. Format: YYYYMMDD Default: 0 The time of day (seconds) corresponding to prescribed_ozone_fixed_ymd at which the prescribed ozone data is fixed if prescribed_ozone_type is 'FIXED'. Default: 0 seconds Full pathname of the directory that contains the files specified in prescribed_volcaero_filelist. Default: set by build-namelist. Filename of dataset for prescribed volcaero. Default: set by build-namelist. Filename of file that contains a sequence of filenames for prescribed volcaero. The filenames in this file are relative to the directory specied by prescribed_volcaero_datapath. Default: set by build-namelist. Name of variable containing volcaero data in the prescribed volcaero datasets. Default: 'MMRVOLC' Remove the file containing prescribed volcanic aerosol concentrations from local disk when no longer needed. Default: FALSE Type of time interpolation for data in prescribed_volcaero files. Can be set to 'CYCLICAL', 'SERIAL', 'INTERP_MISSING_MONTHS', or 'FIXED'. Default: 'SERIAL' The cycle year of the prescribed volcanic aerosol data if prescribed_volcaero_type is 'CYCLICAL'. Format: YYYY Default: 0 The date at which the prescribed volcanic aerosol data is fixed if prescribed_volcaero_type is 'FIXED'. Format: YYYYMMDD Default: 0 The time of day (seconds) corresponding to prescribed_volcaero_fixed_ymd at which the prescribed volcanic aerosol data is fixed if prescribed_volcaero_type is 'FIXED'. Default: 0 seconds Full pathname of the directory that contains the files specified in prescribed_strataero_filelist. Default: set by build-namelist. Filename of dataset for prescribed volcaero. Default: set by build-namelist. Filename of file that contains a sequence of filenames for prescribed volcaero. The filenames in this file are relative to the directory specied by prescribed_strataero_datapath. Default: set by build-namelist. Name of variable containing volcaero data in the prescribed volcaero datasets. Default: 'MMRVOLC' Remove the file containing prescribed volcanic aerosol concentrations from local disk when no longer needed. Default: FALSE Type of time interpolation for data in prescribed_strataero files. Can be set to 'CYCLICAL', 'SERIAL', 'INTERP_MISSING_MONTHS', or 'FIXED'. Default: 'SERIAL' The cycle year of the prescribed volcanic aerosol data if prescribed_strataero_type is 'CYCLICAL'. Format: YYYY Default: 0 The date at which the prescribed volcanic aerosol data is fixed if prescribed_strataero_type is 'FIXED'. Format: YYYYMMDD Default: 0 The time of day (seconds) corresponding to prescribed_strataero_fixed_ymd at which the prescribed volcanic aerosol data is fixed if prescribed_strataero_type is 'FIXED'. Default: 0 seconds Switch to turn on climate feed backs due to prescribed stratospheric aerosols via the rad_climate namelist variable. Default: false Full pathname of dataset for radiative source function used in look up table photloysis Default: set by build-namelist. Full pathname of dataset for the coefficients of the NOEM nitric oxide model used to calculate its upper boundary concentration. Default: set by build-namelist. Full pathname of boundary dataset for soil erodibility factors. Default: set by build-namelist. Full pathname of time-variant boundary dataset for the time-dependent proxies for solar and geomagnetic activity( F10.7, F10.7a, Kp, Ap ). Default: set by build-namelist. List of full pathnames of surface emission datasets. The chemistry package reads in emission data from a set of netcdf files in units of "molecules/cm2/s". Each tracer species emissions is read from its own file as directed by the namelist variable srf_emis_specifier. The srf_emis_specifier variable tells the model which species have emissions and the file path for the corresponding species. That is, the srf_emis_specifier variable is set something like: srf_emis_specifier = 'CH4 -> /path/emis.ch4.nc', 'CO -> /path/emis.co.nc', etc... Each emission file can have more than one source. When the emission are read in the sources are summed to give a total emission field for the corresponding species. The emission can be read in as time series of data, cycle over a given year, or be fixed to a given date. Default: set by build-namelist. Full pathname of dataset containing tropopheric sulfate aerosols Default: set by build-namelist. Full pathname of the directory that contains the files specified in sulf_filelist. Default: set by build-namelist. Filename of file that contains a sequence of filenames for prescribed sulfate. The filenames in this file are relative to the directory specied by sulf_datapath. Default: set by build-namelist. Name of variable containing sulfate data in the prescribed sulfate datasets. Default: 'SULFATE' Remove the file containing prescribed sulfate concentrations from local disk when no longer needed. Default: FALSE Type of time interpolation for data in prescribed sulfate files. Can be set to 'CYCLICAL', 'SERIAL', 'INTERP_MISSING_MONTHS', or 'FIXED'. Default: 'SERIAL' The cycle year of the prescribed sulfate data if sulf_type is 'CYCLICAL'. Format: YYYY Default: 0 The date at which the prescribed sulfate data is fixed if sulf_type is 'FIXED'. Format: YYYYMMDD Default: 0 The time of day (seconds) corresponding to sulf_fixed_ymd at which the prescribed sulfate data is fixed if sulf_type is 'FIXED'. Default: 0 seconds Full pathname of dataset for TGCM upper boundary Default: set by build-namelist. Type of time interpolation for data in TGCM upper boundary file. Can be set to 'CYCLICAL', 'SERIAL', 'INTERP_MISSING_MONTHS', or 'FIXED'. Default: 'SERIAL' The cycle year of the TGCM upper boundary data if tgcm_ubc_type is 'CYCLICAL'. Format: YYYY Default: 0 The date at which the TGCM upper boundary data is fixed if tgcm_ubc_type is 'FIXED'. Format: YYYY Default: 0 The time of day (seconds) corresponding to tgcm_ubc_fixed_ymd at which the TGCM upper boundary data is fixed if tgcm_ubc_type is 'FIXED'. Default: 0 seconds Perturbation applied to the upper boundary temperature. Default: 0.0 Frequency in time steps at which the chemical equations are solved. Default: 1 Filename of dataset for linoz cholirine loading. Default: none. Type of time interpolation type for data in chlorine_loading_file Default: 'SERIAL' The time of day (seconds) corresponding to chlorine_loading_fixed_ymd at which the chlorine loading data is fixed if chlorine_loading_type is 'FIXED'. Default: 0 seconds The date at which the chlorine loading data is fixed if chlorine_loading_type is 'FIXED'. Format: YYYYMMDD Default: 0 Full pathname of the directory that contains the files specified in linoz_data_filelist. Default: none. Filename of dataset for LINOZ data. Default: none. Filename of file that contains a sequence of filenames of the linoz data. The filenames in this file are relative to the directory specied by linoz_data_datapath. Default: set by build-namelist. Type of time interpolation for data in linoz_data files. Can be set to 'CYCLICAL', 'SERIAL', 'INTERP_MISSING_MONTHS', or 'FIXED'. Default: 'SERIAL' Remove the file containing LINOZ data from local disk when no longer needed. Default: FALSE The cycle year of the LINOZ data if linoz_data_type is 'CYCLICAL'. Format: YYYY Default: 0 The date at which the LINOZ data is fixed if linoz_data_type is 'FIXED'. Format: YYYYMMDD Default: 0 The time of day (seconds) corresponding to linoz_data_fixed_ymd at which the LINOZ data is fixed if linoz_data_type is 'FIXED'. Default: 0 seconds Full pathname of the directory that contains the files specified in tracer_cnst_filelist. Default: set by build-namelist. Filename of dataset for the prescribed chemical constituents. Default: set by build-namelist. Filename of file that contains a sequence of filenames for the prescribed chemical constituents. The filenames in this file are relative to the directory specied by tracer_cnst_datapath. Default: set by build-namelist. Remove the file containing prescribed chemical constituents from local disk when no longer needed. Default: FALSE List of prescribed chemical constituents. Default: set by build-namelist. Type of time interpolation for data in tracer_cnst files. Default: 'SERIAL' The cycle year of the prescribed chemical constituents data if tracer_cnst_type is 'CYCLICAL'. Format: YYYY Default: 0 The date at which the chemical constituents data is fixed if tracer_cnst_type is 'FIXED'. Format: YYYYMMDD Default: 0 The time of day (seconds) corresponding to tracer_cnst_fixed_ymd at which the chemical constituents data is fixed if tracer_cnst_type is 'FIXED'. Default: 0 seconds Full pathname of the directory that contains the files specified in tracer_srcs_filelist. Default: set by build-namelist. Filename of dataset for the prescribed chemical sources. Default: set by build-namelist. Filename of file that contains a sequence of datasets for the prescribed chemical sources. The filenames in this file are relative to the directory specied by tracer_srcs_datapath. Default: set by build-namelist. Remove the file containing prescribed chemical sources from local disk when no longer needed. Default: FALSE List of prescribed chemical sources Default: set by build-namelist. Type of time interpolation for data in tracer_srcs files. Default: 'SERIAL' The cycle year of the prescribed chemical sources data if tracer_srcs_type is 'CYCLICAL'. Format: YYYY Default: 0 The date at which the chemical sources data is fixed if tracer_srcs_type is 'FIXED'. Format: YYYYMMDD Default: 0 The time of day (seconds) corresponding to tracer_srcs_fixed_ymd at which the chemical sources data is fixed if tracer_srcs_type is 'FIXED'. Default: 0 seconds If TRUE then use the FTUV method to calculate the photolysis reactions rates, otherwise use the look up table method. Default: FALSE Full pathname of dataset for Chebyshev polynomial Coeff data used for photolysis cross sections. Default: set by build-namelist. Full pathname of cross section dataset for long wavelengh photolysis Default: set by build-namelist. Full pathname of cross section dataset for short wavelengh photolysis Default: set by build-namelist. List of species that undergo dry deposition. Default: set by build-namelist. Dry deposition method used. This specifies the method used to calculate dry deposition velocities of gas-phase chemical species. The available methods are: 'table' - prescribed method in CAM 'xactive_atm' - interactive method in CAM 'xactive_lnd' - interactive method in CLM Default: set by build-namelist Give the user the ability to specify rate families (or groupings) diagnostics based on reaction tag names. These group names can be added to history fincl variables. Example: rate_sums = 'OX_P = NO_HO2 + CH3O2_NO + 2*jo2_b ... ', 'OX_L = NO2_O_M + HO2_O3 + CLO_O ...', fincl1 = 'OX_P','OX_L', ... Default: none File containing fire emissions factors. Default: none Fire emissions specifier. Default: none If true fire emissions are input into atmosphere as elevated forcings. Otherwise they are treated as surface emissions. Default: TRUE File containing MEGAN emissions factors. Default: set by build-namelist. MEGAN specifier. Default: set by build-namelist. MEGAN mapped isoprene emissions facters switch If true then use mapped MEGAN emissions facters for isoprene. Default: .false. List of fluxes needed by the CARMA model, from CLM to CAM. Default: set by build-namelist. Filename of the prescribed waccm forcing data used with waccm_sc chemistry. This contains prescribed constituents for non-LTE calculations and heating rates for wavelengths less than 200 nm. Default: set by build-namelist. Full pathname of the directory that contains the files specified in waccm_forcing_filelist. Default: set by build-namelist. A file that contains a sequence of filenames for prescribed waccm forcing data. The filenames in this file are relative to the directory specied by waccm_forcing_datapath. Default: set by build-namelist. Remove the file containing prescribed waccm forcing data from local disk when no longer needed. Default: FALSE Names of variables containing concentrations and heating rate in the prescribed waccm forcing datasets. Default: none Type of time interpolation for data in waccm_forcing files. Can be set to 'CYCLICAL', 'SERIAL', 'INTERP_MISSING_MONTHS', or 'FIXED'. Default: 'CYCLICAL' The cycle year of the prescribed waccm forcing data if waccm_forcing_type is 'CYCLICAL'. Format: YYYY Default: 0 The date at which the prescribed waccm forcing data is fixed if waccm_forcing_type is 'FIXED'. Format: YYYYMMDD Default: 0 The time of day (seconds) corresponding to waccm_forcing_fixed_ymd at which the prescribed waccm forcing data is fixed if waccm_forcing_type is 'FIXED'. Default: 0 seconds Full pathname of time-variant boundary dataset for H2O production/loss rates. Default: set by build-namelist. Determines which constituents are used from NLTE calculations. TRUE implies use prognostic constituents. FALSE implies use constituents from dataset specified by waccm_forcing_file. Default: TRUE for full WACCM code; FALSE for WACCM_GHG. TRUE implies assume cyclic qbo data. Default: FALSE Filepath for qbo forcing dataset. Default: Set by build-namelist. TRUE implies qbo package is active. Default: FALSE If set this year is used for setting geomagnetic coordinates through out the simulation. If not set the model simulation year is used. Default: none International Geomagnetic Reference Field (IGRF) coefficients. Default: None. Filepath input dataset for ionization due to energetic particle precipitation. Default: None. Variable name in netCDF file epp_all_filepath which contains ion pairs production rates. Default: None. Filepath input dataset for ionization due to solar proton events. Default: None. Variable name in netCDF file epp_spe_filepath which contains ion pairs production rates. Default: None. Filepath input dataset for ionization due to medium energy electrons. Default: None. Variable name in netCDF file epp_mee_filepath which contains ion pairs production rates. Default: None. Filepath input dataset for ionization due to galactic cosmic rays. Default: None. Variable name in netCDF file epp_gcr_filepath which contains ion pairs production rates. Default: None. Set for refined exodus meshes (variable viscosity). Viscosity in namelist specified for regions with a resolution equivilant to a uniform grid with se_ne = se_fine_ne. Default: -1 (not used) CAM physics forcing option: 0: tendencies 1: adjustments 2: hybrid Default: Set by build-namelist. Hyperviscosity operator is the Laplacian^hypervis_order. The only supported value in CAM is "2". Default: Set by build-namelist. Scalar viscosity with variable coefficient. Use variable hyperviscosity based on element area limited by se_max_hypervis_courant. Use tensor hyperviscosity. Citation: Guba, O., Taylor, M. A., Ullrich, P. A., Overfelt, J. R., and Levy, M. N.: The spectral element method (SEM) on variable-resolution grids: evaluating grid sensitivity and resolution-aware numerical viscosity, Geosci. Model Dev., 7, 2803-2816, doi:10.5194/gmd-7-2803-2014, 2014. Default: 0 (i.e., not used) Number of hyperviscosity subcycles per dynamics timestep. Default: Set by build-namelist (probably 2). Number of hyperviscosity subcycles done in tracer advection code. The only supported value in CAM is 1. Default: Set by build-namelist. Limiter used for horizontal tracer advection: 0: None 4: Sign-preserving limiter. 8: Monotone limiter. Default: Set by build-namelist (probably 8). Upper bound for Courant number, used to limit se_hypervis_power. Default: 1.0e99 (i.e., not used) Filename of exodus file to read grid from (generated by CUBIT or SQuadGen). Default: "" Number of elements along a cube edge. Must match value of grid. Set this to zero to use a refined mesh. Default: Set by build-namelist. Number of PEs to be used by SE dycore. Default: Number of PEs used by CAM. Number of dynamics steps per physics timestep. Default: Set by build-namelist. Dynamics hyperviscosity [m^4/s]. Default: Set by build-namelist. Weighting of divergence component when applying hyperviscosity. If < 0, uses nu. Default: value of nu Hyperviscosity applied to dp3d (layer thickness when running vertically lagrangian dynamics) [m^4/s]. Default: Set by build-namelist. Hyperviscosity applied during tracer advection [m^4/s]. If < 0, uses nu. Default: Set by build-namelist. Second-order viscosity applied only near the model top [m^2/s]. Default: Set by build-namelist. Tracer advection is done every qsplit dynamics timesteps. Default: Set by build-namelist. TRUE specified use of a refined grid (mesh) for this run. Default: FALSE Vertically lagrangian code updates every rsplit tracer timesteps. If rsplit=0, vertically lagrangian code is off. Default: Set by build-namelist. Frequency with which diagnostic output is written to log (output every statefreq dynamics timesteps). Default: Set by build-namelist. Time stepping method for SE dycore se_tstep_type=0 pure Leapfrog except for very first timestep se_tstep_type=1 RK2 followed by qsplit-1 Leapfrog steps; second-order accurate in time (CESM1.2.0 setting) se_tstep_type=2 RK2-SSP 3 stage (as used by tracers) se_tstep_type=3 classic Runga-Kutta (RK) 3 stage se_tstep_type=4 Kinnmark&Gray Runga-Kutta (RK) 4 stage se_tstep_type=5 Kinnmark&Gray Runga-Kutta (RK) 5 stage (3rd order accurate in time) Default: 5 CAM-SE vertical remap algorithm 0: default value, Zerroukat monotonic splines 1: PPM vertical remap with mirroring at the boundaries (solid wall bc's, high-order throughout) 2: PPM vertical remap without mirroring at the boundaries (no bc's enforced, first-order at two cells bordering top and bottom boundaries) Default: Set by build-namelist. Default: 4 Default: 1 Default: -1 Default: "netcdf" Default: .false. Default: 0 Default: -1 Default: 0 pio_rearranger = 1 for box rearranger = 2 for subset rearranger Default: box rearranger used by pio1, subset is default for pio2. Full pathname of docn restart file. Default: set by build-namelist. Full pathname of docn restart file. Default: set by build-namelist. Full pathname of docn restart file. Default: set by build-namelist. Full pathname of time-variant sea-surface temperature and sea-ice concentration boundary dataset. Default: set by build-namelist. Full pathname of Default: set by build-namelist. Full pathname of grid file for time-variant sea-surface temperature and sea-ice concentration boundary dataset. Default: set by build-namelist. Flag for yearly cycling of SST data. If set to FALSE, a multi-year dataset is assumed, otherwise a single-year dataset is assumed, and SSTs will be cycled over the first 12 values in the file. Not used if running with CCSM. Default: TRUE The first year of the multi-year SST dataset which is read by CICE for the prescribed ice fraction. This needs to be set for AMIP simulations. Default: 0 The last year of the multi-year SST dataset which is read by CICE for the prescribed ice fraction. This needs to be set for AMIP simulations. Default: 0 Full pathname of master restart file from which to branch. Setting is Required for branch run. Default: none Name of file that the atmosphere component log messages will be written to. By default all log messages are written to stdout. Default: "" Absolute pathname of directory that the file specified by atm_logfile will be written to. Default: "." Name of file that the driver component log messages will be written to. By default all log messages are written to stdout. Default: "" Absolute pathname of directory that the file specified by cpl_logfile will be written to. Default: "." Name of file that the land component log messages will be written to. By default all log messages are written to stdout. Default: "" Absolute pathname of directory that the file specified by lnd_logfile will be written to. Default: "." Name of file that the runoff component log messages will be written to. By default all log messages are written to stdout. Default: "" Absolute pathname of directory that the file specified by rof_logfile will be written to. Default: "." Default: FALSE Default: FALSE Default: FALSE Default: FALSE Default: Stride used in selecting the processes in the atm communicator group. Default: 1 Root process of the atm communicator group. Default: 0 Number of atm tasks. Default: total number of tasks assigned to job. Number of threads in each atm task. Default: 1 Stride used in selecting the processes in the lnd communicator group. Default: 1 Root process of the lnd communicator group. Default: 0 Number of lnd tasks. Default: total number of tasks assigned to job. Number of threads in each lnd task. Default: 1 Stride used in selecting the processes in the ice communicator group. Default: 1 Root process of the ice communicator group. Default: 0 Number of ice tasks. Default: total number of tasks assigned to job. Number of threads in each ice task. Default: 1 Stride used in selecting the processes in the ocn communicator group. Default: 1 Root process of the ocn communicator group. Default: 0 Number of ocn tasks. Default: total number of tasks assigned to job. Number of threads in each ocn task. Default: 1 Stride used in selecting the processes in the rof communicator group. Default: 1 Root process of the rof communicator group. Default: 0 Number of rof tasks. Default: total number of tasks assigned to job. Number of threads in each rof task. Default: 1 Stride used in selecting the processes in the cpl communicator group. Default: 1 Root process of the cpl communicator group. Default: 0 Number of cpl tasks. Default: total number of tasks assigned to job. Number of threads in each cpl task. Default: 1 If true, run model in "aqua planet" mode. Only one of atm_adiabatic, atm_ideal_phys, or aqua_planet can be true. Default: FALSE Set the sst to a particular analytic solution. **Not currently functional** The type of analytic solution is currently hardcoded with a parameter in ocn_comp.F90. Default: 1 If TRUE, do not run model physics, only run the dynamical core. Only one of atm_adiabatic, atm_ideal_phys, or aqua_planet can be true. Default: FALSE If FALSE then CAM will set the deposition fluxes to zero before sending them to the coupler. A side effect of setting the variable chem_rad_passive to TRUE is that this variable will be set to FALSE (the deposition fluxes must be set to zero in order for the chemistry not to impact the climate). Default: TRUE If TRUE run the idealized Held&Suarez physics package. Only one of atm_adiabatic, atm_ideal_phys, or aqua_planet can be true. Default: FALSE bit for bit flag Default: FALSE If TRUE, use the pre-existing case name for a branch run. Default: FALSE annual budget level Default: 1 daily budget level Default: 0 instantaneous budget level Default: 0 long term budget level written at end of year Default: 1 long term budget level written at end of run Default: 0 monthly budget level Default: 1 Case title. Default: none. Case identifier. The value of case_name is used in the default filenames of both the history and restart files (see the hfilename_spec namelist option). The "%c" string in the hfilename_spec templates are expanded using the value of case_name when history filenames are created. Default: set by build-namelist. Use netcdf 64 bit offset, large file support. Default: FALSE T => do heat/water budget diagnostics Default: FALSE T => enable run time setting of thread count for each component Default: FALSE T => no diurnal cycle in ocn albedos. Default: FALSE Selects E,P,R adjustment technique. Default: 'off' Current machine. Default: 'unknown' Debug flag. Default: 1 Postfix for output log files. Default: '.log' Model version. Default: 'unknown' Coupler sequence option - CESM1_MOD_TIGHT is identical to old ocean_tight_coupling. Default: 'CESM1_MOD_TIGHT' Earth's eccentricity of orbit. (unitless: typically 0. to 0.1). Setting is Required if orb_iyear not set. Not used when running as part of CCSM. Default: none Mode to specify how orbital parameters are to be set. Not used when running as part of CCSM. Default: fixed_year Year (AD) used to compute earth's orbital parameters. If not set, then use the values from the orb_eccen, orb_mvelp, and orb_obliq namelist parameters. If only orb_iyear is set, orbital parameters will be computed automatically (based on Berger, 1977). If one of orb_eccen, orb_mvelp, or orb_obliq is set, all three must be set. If all four of the above are set by the user, orb_iyear takes precedence. Setting is Required unless orb_eccen, orb_mvelp, and orb_obliq are set. Not used when running as part of CCSM. Default: 1990. Earth's moving vernal equinox at perihelion (degrees: 0. to 360.0). Setting is Required if orb_iyear not set. Not used when running as part of CCSM. Default: none Earth's orbital angle of obliquity (degrees: -90. to +90., typically 22. to 26.). Setting is Required if orb_iyear not set. Not used when running as part of CCSM. Default: none Root output files Default: './' Set to TRUE to specify that the run will use a perpetual calendar, i.e., a diurnal cycle will be present for the fixed calendar day specified by perpetual_ymd. Default: FALSE Perpetual date encoded in an integer as (year*1000 + month*100 + day). If aqua_planet = .true. then perpetual_ymd is ignored and the perpetual date is set to 321. Default: none. Restart filename. Default: none List of namelist variables that may be overridden on a restart run. Default: none Restart pointer filename. Default: 'rpointer.drv' are ocean and ice grids same lat/lon/size Default: TRUE Latitude value of single column. Default: none. Longitude value of single column. Default: none. Set to TRUE to turn on single column mode. Default: FALSE Run type. 'startup' is an initial run. 'continue' is a restart run. 'branch' is a restart run in which properties of the output history files may be changed. Default: 'startup' Current user. Default: 'unknown' Invoke vector mapping option Default: 'cart3d' Flag to indicate whether to use the double-double distributed sum algorithm rather than the (almost) infinite precision reproducible distributed sum algorithm. Default: FALSE Flag to indicate whether a distributed sum that violates the difference tolerance specified by reprosum_diffmax should be recomputed using a floating point-based (but nonscalable) reproducible algorithm. Default: FALSE Relative difference between repro and nonrepro algorithms that will generate a warning. This will also force a recompute using a nonscalable algorithm if reprosum_recompute is true. If less than zero, then the difference will not be evaluated (and the nonrepro algorithm will not be computed). Default: -1.0 Coupling interval between the atmosphere and other system components. This is how frequently information can be communicated between the atmosphere and the surface models. Default: set by build-namelist. Default: Calendar type "NO_LEAP" for consistent 365-days per year or "GREGORIAN" to include leap-years. Note that if "GREGORIAN" is selected leap-years will be used in the time manager, but the calculation of the earth's orbit still assumes 365 day years. Valid values are "NO_LEAP" or "GREGORIAN". Default: "NO_LEAP" Write restart at end of run. Default: TRUE Default: set by build-namelist. Default: Default: set by build-namelist. Default: Default: set by build-namelist. Default: Default: set by build-namelist. Default: Reference time-of-day expressed as seconds past midnight. Used in conjuction with ref_ymd to set the reference time. Default: set to start_tod. Reference date encoded in an integer as (year*1000 + month*100 + day). Used in conjuction with ref_tod to set the reference time which is used to define a time coordinate for the output history files. The convention for the unit string of a time coordinate is of the form "time-unit since reference-time", for example, "days since 1990-01-01 00:00:00". The reference-time part of this string is specified by the ref_ymd and ref_tod variables. Default: set to start_ymd. Set the restart interval as a number of elapsed time units which are specified by restart_option. Default: 1 Set the interval between writing restart files using one of the options 'nsteps', 'ndays', 'nmonths', or 'nyears', in conjuction with stop_n to set the number of time units. A convenience option allows specifying that restart files be written at the end of each month or at the end of each year by using the options 'monthly' or 'yearly' respectively. It is also possible to request that no restart files be written via the option 'none', or that restart files be written only at the end of the run via the option 'end'. Default: 'monthly' Start time-of-day expressed as seconds past midnight. Used in conjuction with start_ymd to set the start time. Default: 0. Start date encoded in an integer as (year*1000 + month*100 + day). Used in conjuction with start_tod to set the start date of the simulation. Default: set by build-namelist. Set the length of run as a number of elapsed time units which are specified by stop_option. Default: 1 Set the length of run as an elapsed time using one of the options 'nsteps', 'ndays', 'nmonths', or 'nyears', in conjuction with stop_n to set the number of elapsed time units. Alternatively, set the final simulation time in absolute terms by using the option 'date' in conjuction with stop_ymd, and stop_tod to specify a date and time of day at which the simulation should stop. Default: 'ndays' Stop time-of-day expressed as seconds past midnight. Used in conjuction with stop_ymd to set the stop time. Default: 0. Stop date encoded in an integer as (year*1000 + month*100 + day). Used in conjuction with stop_tod to set the stop date of the simulation. Default: none. Clumps per processor. Default: Number of threads per process; or 1 if no OMP. Default: 355. Type of CO2 feedback. Default: 'constant' for BGC mode on if to create crop as separate landunits Default: FALSE Full pathname of grid/landfrac data file (on atm grid). Default: set by build-namelist. Full pathname of atmosphere topography data file (on atm grid). Default: set by build-namelist. Full pathname of initial conditions file. Default: set by build-namelist. Full pathname of land topography file. Default: set by build-namelist. Full pathname of PFT data file. Default: set by build-namelist. Full pathname of time varying PFT data file. Default: set by build-namelist. Full pathname of input datafile for rtm. Default: set by build-namelist. Full pathname of surface data file. Default: set by build-namelist. SNICAR (SNow, ICe, and Aerosol Radiative model) optical data file name SNICAR (SNow, ICe, and Aerosol Radiative model) snow aging data file name Per file averaging flag. Default: ' ',' ',' ',' ',' ',' ' TRUE implies do grid averaging. One value for each output history file. Default: All TRUE. TRUE indicates no default history fields. Default: FALSE Fields to exclude from history file 1. Default: none. Fields to exclude from history file 2. Default: none. Fields to exclude from history file 3. Default: none. Fields to exclude from history file 4. Default: none. Fields to exclude from history file 5. Default: none. Fields to exclude from history file 6. Default: none. Fields to add to history file 1. Default: none. Fields to add to history file 2. Default: none. Fields to add to history file 3. Default: none. Fields to add to history file 4. Default: none. Fields to add to history file 5. Default: none. Fields to add to history file 6. Default: none. Per file maximum number of time samples. Default: 30,30,30,30,30,30 Per file history output precision. Default: 2,2,2,2,2,2 Per file history write frequency (0=monthly) Default: 0,-24,-24,-24,-24,-24 Per file type1d Default: ' ',' ',' ',' ',' ',' ' Full pathname of master restart file for branch run. Default: none. number of segments per clump for decomposition Default: 20 on to output NetCDF files in large-file format Default: FALSE FALSE implies don't write any restart files. Default: TRUE If greater than one, average rtm over rtm_nsteps time steps This varible is only used internally by build-namelist to determine correct defaults for datasets such as surface datasets, initial conditions, dynamic PFT, aerosol-deposition, Nitrogen deposition rates etc. Default: set by build-namelist. Urban air conditioning/heating and wasteheat Default: 'off' TRUE if want urban traffic flux Default: FALSE true if want output written Default: FALSE River runoff Default: NULL Length of siderial day [seconds]. Default: set to shr_const value Radius of Earth [m]. Default: set to shr_const value Acceleration of gravity [m/s**2]. Default: set to shr_const value Molecular weight of dry air [g/mol] Default: set to shr_const value Heat capacity of dry air at constant pressure [J/kg/K]. Default: set to shr_const value Molecular weight of water [g/mol]. Default: set to shr_const value Heat capacity of water vapor at constant pressure [J/kg/K]. Default: set to shr_const value Freezing point of water [K]. Default: set to shr_const value Filepath for dataset for offline unit driver. Default: none List of filepaths for dataset for offline unit driver. Default: none trefread/NCL/0000755006307300054160000000000013047371145012434 5ustar islascgdamptrefread/NCL/make_tref_input.ncl0000644006307300054160000001252613047371145016314 0ustar islascgdamp; **************************************************************** ; * Script to generate a netcdf file containing a user specified * ; * relaxation temperature profile for input to the CAM dry * ; * dynamical core. * ; * Input: * ; * filein = File containig input tref (nlon,nlat,npre) * ; * expected that pressures are ordered bottom to top * ; * res = horizontal resolution * ; * nlev = number of vertical levels * ; * Output: * ; * fileout = location for output file * ; **************************************************************** ; * Isla Simpson (6th Feb 2017) * ; **************************************************************** ;***specify input relaxation tempreature profile file ;***expected to contain tref(lon,lat,pre) filein="/scratch/cluster/islas/trefread/inprofiles/trefprof002.nc" ;***create output file setfileoption("nc","Format","LargeFile") fileout = "/scratch/cluster/islas/trefread/outprofiles/tref002.nc" system("/bin/rm -f " fileout) ; remove if exists fileout = addfile( fileout, "c") ; open output file ;*** vertical level options nlev=30 ; nlev=40 ; nlev=60 ; nlev=80 load "./levels.ncl" ;*** horizontal resolution options res = "64x128" ;T42 ;res = "128x256" ;T85 ;res = "256x512" ;T255 ;res = "512x1024" ;T511 (?) load "./gaus_horizontal_grid.ncl" ;*** p0 P0 = 100000.0d0 P0!0="ncl_scalar" P0@long_name = "reference pressure" P0@units = "Pa" ;*** Read in input tref profile on arbitrary horizontal and vertical grid infile=addfile(filein,"r") lonin=infile->lon latin=infile->lat prein=infile->pre(::-1) trefin=infile->tref(::-1,:,:) ;*** horizontal interpolation onto gaussian grid trefnew = linint2(lonin,latin,trefin,False,lon,lat,0) ;*** vertical interpolation onto model pressure levels (assuming ps=1000) trefnew2=new( (/1,dimsizes(lev),dimsizes(lat),dimsizes(lon)/),"double") ;trefnew2=int2p(prein,trefnew,lev,0) trefnew2(0,:,:,:)=int2p_n(prein,trefnew,lev,-2,0) trefnew2!0 = "time" trefnew2!1 = "lev" trefnew2!2 = "lat" trefnew2!3 = "lon" ;***Time coordinate ntime=1 time = (/0.d0/) time!0 = "time" time@long_name = "time" time@units = "days since 1979-01-01 00:00:00" time@calendar = "noleap" time@bounds = "time_bnds" date = (/19790101/) date!0="time" date@long_name = "current date (YYYYMMDD)" datesec = (/0/) datesec!0="time" datesec@long_name = "current seconds of current date" ;***output interpolated tref setfileoption(fout,"DefineMode",True) ;***global attributes fAtt = True ; assign file attributes fAtt@title = "User specified relaxation temperature profile" fAtt@source_file = "ncl script" fAtt@Conventions = "None" fAtt@creation_date = systemfunc ("date") fileattdef( fout, fAtt ) ; copy file attributes ;***pre-defined coordinate variables dimNames = (/"time","lev","lat","lon"/) dimSizes= (/-1, dimsizes(lev), dimsizes(lat),dimsizes(lon)/) dimUnlim = (/True, False, False, False/) filedimdef(fout,dimNames,dimSizes,dimUnlim) ;***pre-define output variables filevardef(fout, "time", typeof(time), getvardims(time)) filevardef(fout, "date" ,typeof(date) ,getvardims(date)) filevardef(fout, "datesec" ,typeof(datesec) ,getvardims(datesec)) filevardef(fout, "lev" ,typeof(lev) ,getvardims(lev) ) filevardef(fout, "hyam" ,typeof(hyam) ,getvardims(hyam)) filevardef(fout, "hybm" ,typeof(hybm) ,getvardims(hybm)) filevardef(fout, "P0" ,typeof(P0) ,getvardims(P0)) filevardef(fout, "lat" ,typeof(lat) ,getvardims(lat)) filevardef(fout, "lon" ,typeof(lon) ,getvardims(lon)) filevardef(fout, "gw" ,typeof(gw) ,getvardims(gw)) filevardef(fout, "tref" ,typeof(trefnew2) ,getvardims(trefnew2)) ; copy attributes filevarattdef(fout,"time" ,time) ; copy time attributes filevarattdef(fout,"date" ,date) ; copy date attributes filevarattdef(fout,"datesec",datesec) ; copy datesec attributes filevarattdef(fout,"lev" ,lev) ; copy lev attributes filevarattdef(fout,"hyam" ,hyam) ; copy hyam attributes filevarattdef(fout,"hybm" ,hybm) ; copy hybm attributes filevarattdef(fout,"P0" ,P0) ; copy P0 attributes filevarattdef(fout,"lat" ,lat) ; copy lat attributes filevarattdef(fout,"lon" ,lon) ; copy lon attributes filevarattdef(fout,"gw" ,gw) ; copy gw attributes filevarattdef(fout,"tref" ,trefnew2) ; copy tref attributes setfileoption(fout,"DefineMode",False) fout->time = (/time/) fout->date = (/date/) fout->datesec = (/datesec/) fout->lev = (/lev/) fout->hyam = (/hyam/) fout->hybm = (/hybm/) fout->P0 = (/P0/) fout->lat = (/lat/) fout->lon = (/lon/) fout->gw = (/gw/) fout->tref = (/trefnew2/) trefread/NCL/levels.ncl0000755006307300054160000001645213047371145014437 0ustar islascgdamp;-------------------------------------------------------------------------- ; compute levels ;-------------------------------------------------------------------------- ; -----------note model top cannot be zero. Set to a small value ; -----------20 sigma levels. Evenly spaced in level interfaces if(nlev .eq. 20) then hyai=(/1e-5, 0.0, 0.0, 0.0, \ 0.0, 0.0, 0.0, 0.0, \ 0.0, 0.0, 0.0, 0.0, \ 0.0, 0.0, 0.0, 0.0, \ 0.0, 0.0, 0.0, 0.0, \ 0.0/) hybi=(/0.00000, 0.0500000, 0.100000, 0.150000, \ 0.200000, 0.250000, 0.300000, 0.350000, \ 0.400000, 0.450000, 0.500000, 0.550000, \ 0.600000, 0.650000, 0.700000, 0.750000, \ 0.800000, 0.850000, 0.900000, 0.950000, \ 1.00000/) end if ; -----------30 sigma levels. Evenly spaced in level interfaces if(nlev .eq. 30 ) then hyai = (/ 1e-5, 0.0, \ 0.0, 0.0, 0.0, 0.0, \ 0.0, 0.0, 0.0, 0.0, \ 0.0, 0.0, 0.0, 0.0, \ 0.0, 0.0, 0.0, 0.0, 0.0, \ 0.0, 0.0, 0.0, 0.0, \ 0.0, 0.0, 0.0, 0.0, \ 0.0, 0.0, 0.0, 0.0/) hybi = (/ 0.00000, 0.0333333, 0.0666667, 0.100000, \ 0.133333, 0.166667, 0.200000, 0.233333, \ 0.266667, 0.300000, 0.333333, 0.366667, \ 0.400000, 0.433333, 0.466667, 0.500000, \ 0.533333, 0.566667, 0.600000, 0.633333, \ 0.666667, 0.700000, 0.733333, 0.766667, \ 0.800000, 0.833333, 0.866667, 0.900000, \ 0.933333, 0.966667, 1/) end if ; -----------40 sigma levels. Evenly spaced in level interfaces if(nlev .eq. 40 ) then hyai = (/ 1e-5, 0.0, 0.0, 0.0, \ 0.0, 0.0, 0.0, 0.0, \ 0.0, 0.0, 0.0, 0.0, \ 0.0, 0.0, 0.0, 0.0, \ 0.0, 0.0, 0.0, 0.0, \ 0.0, 0.0, 0.0, 0.0, \ 0.0, 0.0, 0.0, 0.0, \ 0.0, 0.0, 0.0, 0.0, \ 0.0, 0.0, 0.0, 0.0, \ 0.0, 0.0, 0.0, 0.0, 0.0/) hybi = (/ 0.0, 0.0250000, 0.0500000, 0.0750000, \ 0.100000, 0.125000, 0.150000, 0.175000, \ 0.200000, 0.225000, 0.250000, 0.275000, \ 0.300000, 0.325000, 0.350000, 0.375000, \ 0.400000, 0.425000, 0.450000, 0.475000, \ 0.500000, 0.525000, 0.550000, 0.575000, \ 0.600000, 0.625000, 0.650000, 0.675000, \ 0.700000, 0.725000, 0.750000, 0.775000, \ 0.800000, 0.825000, 0.850000, 0.875000, \ 0.900000, 0.925000, 0.950000, 0.975000, 1.0/) end if ; -----------60 sigma levels. Evenly spaced in level interfaces if (nlev .eq. 60) then hyai = (/1e-5, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, \ 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, \ 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, \ 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, \ 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, \ 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00,\ 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 /) hybi = (/0.00000, 0.0166667, 0.0333334, 0.0500001, 0.0666668, \ 0.0833335, 0.100000, 0.116667, 0.133334, 0.150000, 0.166667, \ 0.183334, 0.200000, 0.216667, 0.233334, 0.250000, 0.266667, \ 0.283334, 0.300001, 0.316667, 0.333334, 0.350001, 0.366667, \ 0.383334, 0.400001, 0.416667, 0.433334, 0.450001, 0.466668, \ 0.483334, 0.500001, 0.516668, 0.533334, 0.550001, 0.566668, \ 0.583334, 0.600001, 0.616668, 0.633335, 0.650001, 0.666668, \ 0.683335, 0.700001, 0.716668, 0.733335, 0.750001, 0.766668, \ 0.783335, 0.800002, 0.816668, 0.833335, 0.850002, 0.866668, \ 0.883335, 0.900002, 0.916668, 0.933335, 0.950002, 0.966669, \ 0.983335, 1.00000 /) end if ; -----------80 sigma levels. Evenly spaced in level interfaces if (nlev .eq. 80) then hyai = (/1e-5, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, \ 0.00, 0.00, 0.00, \ 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, \ 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, \ 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, \ 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, \ 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, \ 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, \ 0.00, 0.00, 0.00, 0.00/) hybi = (/0.00000, 0.0125000, 0.0250000, 0.0375000, 0.0500000, 0.0625000, \ 0.0750000, 0.0875000, 0.100000, 0.112500, 0.125000, \ 0.137500, 0.150000, 0.162500, 0.175000, 0.187500, 0.200000, 0.212500, \ 0.225000, 0.237500, 0.250000, 0.262500, 0.275000, 0.287500, 0.300000, \ 0.312500, 0.325000, 0.337500, 0.350000, 0.362500, 0.375000, 0.387500, \ 0.400000, 0.412500, 0.425000, 0.437500, 0.450000, 0.462500, 0.475000, \ 0.487500, 0.500000, 0.512500, 0.525000, 0.537500, 0.550000, 0.562500, \ 0.575000, 0.587500, 0.600000, 0.612500, 0.625000, 0.637500, 0.650000, \ 0.662500, 0.675000, 0.687500, 0.700000, 0.712500, 0.725000, 0.737500, \ 0.750000, 0.762500, 0.775000, 0.787500, 0.800000, 0.812500, 0.825000, \ 0.837500, 0.850000, 0.862500, 0.875000, 0.887500, 0.900000, 0.912500, \ 0.925000, 0.937500, 0.950000, 0.962500, 0.975000, 0.987500, 1.00000/) end if ;-------------------------------------------------------------------------- ; derived variables and values that have to be set in CAM files ; everything below is derived from previous values ;-------------------------------------------------------------------------- nilev = nlev + 1 hyai!0="ilev" hyai@long_name = "hybrid A coefficient at layer interfaces" hybi!0="ilev" hybi@long_name = "hybrid B coefficient at layer interfaces" hyam=new(nlev,"double") do plev = 0, nlev-1 hyam(plev) = (hyai(plev)+hyai(plev+1))/2. end do hyam!0="lev" hyam@long_name = "hybrid A coefficient at layer midpoints" ; ------------------------------------------------------------------------ ; define mid-levels as the half way point between interfaces ; ------------------------------------------------------------------------ hybm=new(nlev,"double") do plev = 0, nlev-1 hybm(plev) = (hybi(plev)+hybi(plev+1))/2. end do hybm!0="lev" hybm@long_name = "hybrid B coefficient at layer midpoints" lev = (hyam+hybm)*1000. lev!0 = "lev" lev@long_name = "hybrid level at midpoints (1000*(A+B))" lev@units = "level" lev@positive = "down" lev@standard_name = "atmosphere_hybrid_sigma_pressure_coordinate" lev@formula_terms = "a: hyam b: hybm p0: P0 ps: PS" ilev = (hyai+hybi)*1000. ilev!0="ilev" ilev@long_name = "hybrid level at interfaces (1000*(A+B))" ilev@units = "level" ilev@positive = "down" ilev@standard_name = "atmosphere_hybrid_sigma_pressure_coordinate" ilev@formula_terms = "a: hyai b: hybi p0: P0 ps: PS" trefread/NCL/gaus_horizontal_grid.ncl0000755006307300054160000000207413047371145017355 0ustar islascgdamp ;=================================================================== ; Define horizontal grid ;=================================================================== if (res .eq. "8x16") then nlat = 8 nlon = 16 end if if (res .eq. "64x128") then nlat = 64 nlon = 128 end if if (res .eq. "128x256") then nlat = 128 nlon = 256 end if if (res .eq. "256x512") then nlat = 256 nlon = 512 end if if (res .eq. "512x1024") then nlat = 512 nlon = 1024 end if gaus_info = gaus(nlat/2) ;fspan(-90.0d0+90.0d0/nlat,90.0d0-90.d0/nlat,(/nlat/)) lat = gaus_info(:,0) lat!0="lat" lat@long_name = "latitude" lat@units = "degrees_north" gw = gaus_info(:,1) gw!0 = "lat" gw@long_name = "gauss weights" lon = fspan(0.d0,360.-360./nlon,(/nlon/)) lon!0 = "lon" lon@long_name = "longitude" lon@units = "degrees_east" trefread/README.txt0000644006307300054160000000772013047376462013533 0ustar islascgdampThe following is an example of modifying the CAM source code to (a) control the relaxation temperature profile used in the dry dynamical core via the namelist (b) read in a user defined relaxation temperature profile (Tref) from a netcdf file for the eulerian (pseudo-spectral) dynamical core.. The following describes how to set this up using a CESM release that is located in the directory $CESM within a run directory at the location $RUNDIR. The enclosed directory ./SourceMods/src.cam contains the source code modifications required to achieve this along with the original version of each file, denoted by filename-ORIG so that users can compare these to see the modifications that have to be made. The following files/modules require modification $CESM/components/cam/bld/namelist_files/namelist_definition.xml $CESM/components/cam/src/control/runtime_opts.F90 $CESM/components/cam/src/physics/simple/held_suarez.F90 $CESM/components/cam/src/physics/simple/held_suarez_cam.F90 For the CAM source code modifications i.e., the .F90 files. The best practise is to copy these files into the directory $RUNDIR/SourceMods/src.cam and make modifications there. The following modifications are applied within a particular CAM version. Users should translate them into the appropriate versions of these modules within their current CAM distribution. ----Step 1. Create the user defined (Tref) in required netcdf file format. The enclosed NCL scripts in ./NCL/ can be used to output a user defined relaxation temperature profile in the format required for input to the model. ----Step 2. Set up namelist parameters to control Tref reading. User defined namelist parameters need to be defined in the file $CESM/components/cam/bld/namelist_files/namelist_definition.xml The enclosed file ./namelist_defintion/namelist_definition.xml is a modified version of namelist_defition.xml that contains two additional namelist parameters treffromfile: type=logical. This controls whether or not the relaxation temperature profile is to be read in from a user specified file. .False. --> use the typical H-S relaxation temperature profil. .True. --> read in from file. treffile: type=char. Path for the file containing the user specified relaxation temperature profile for treffromfile=.True. These namelist parameters have been defined to be part of the group "trefread_nl" so these must be specified in user_nl_cam in the following format &trefread_nl treffile = 'myfilename.nc' treffromfile =.True. / ----Step 3. SourceMods to read in this new namelist. The modified version of held_suarez.F90 contains a new subroutine named trefread_readnl. This reads in the namelist parameters from user_nl_cam as defined above. Users who wish to add in additional namelist parameters for other reasons can adapt this to their needs. This subroutine has to be called from runtime_opts.F90 and so the enclosed version of runtime_opts.F90 has been modified to do this. ----Step 4. Modify held_suarez routines to contain the option to read in tref from a user specified file. The reading of the relaxation temperature profile is performed within the subroutine update_tref_eul This is called from within held_suarez_1994_init and results in an array tref that contains the user specified relaxation temperature profile. The subroutine held_suarez_1994 has been modified to contain the option to use this user specified tref if treffromfile=.True.. The additional parameter lchnk, which identifies the location within the state array where the relaxation is currently being performed, has to be passed from the subroutine held_suarez_tend within the module held_suarez_cam.F90 to held_suarez_1994. The enclosed held_suarez_cam.F90 module contains this modification. Many of the above modification are enclosed within the following comments. ! IRS mods are here ! END IRS but users can also see modifications by differenceing the new and original files e.g. xxdiff held_suarez.F90 held_suarez.F90-ORIG trefread/SourceMods/0000755006307300054160000000000013047373764014114 5ustar islascgdamptrefread/SourceMods/src.cam/0000755006307300054160000000000013047377074015440 5ustar islascgdamptrefread/SourceMods/src.cam/held_suarez.F90-ORIG0000644006307300054160000001371113047376377020773 0ustar islascgdampmodule held_suarez use shr_kind_mod, only: r8 => shr_kind_r8 implicit none private save public :: held_suarez_1994_init public :: held_suarez_1994 !! !! Forcing parameters !! real(r8), parameter :: efoldf = 1._r8 ! efolding time for wind dissipation real(r8), parameter :: efolda = 40._r8 ! efolding time for T dissipation real(r8), parameter :: efolds = 4._r8 ! efolding time for T dissipation real(r8), parameter :: sigmab = 0.7_r8 ! threshold sigma level real(r8), parameter :: t00 = 200._r8 ! minimum reference temperature real(r8), parameter :: kf = 1._r8/(86400._r8*efoldf) ! 1./efolding_time for wind dissipation real(r8), parameter :: onemsig = 1._r8 - sigmab ! 1. - sigma_reference real(r8), parameter :: ka = 1._r8/(86400._r8 * efolda) ! 1./efolding_time for temperature diss. real(r8), parameter :: ks = 1._r8/(86400._r8 * efolds) !! !! Model constants, reset in init call !! real(r8) :: cappa = 2.0_r8 / 7.0_r8 ! R/Cp real(r8) :: cpair = 1004.0_r8 ! specific heat of dry air (J/K/kg) real(r8) :: psurf_ref = 0.0_r8 ! Surface pressure ! pref_mid_norm are layer midpoints normalized by surface pressure ('eta' coordinate) real(r8), allocatable :: pref_mid_norm(:) integer :: pver ! Num vertical levels !======================================================================= contains !======================================================================= subroutine held_suarez_1994_init(cappa_in, cpair_in, psurf_ref_in, pref_mid_norm_in) !! Dummy arguments real(r8), intent(in) :: cappa_in real(r8), intent(in) :: cpair_in real(r8), intent(in) :: psurf_ref_in real(r8), intent(in) :: pref_mid_norm_in(:) pver = size(pref_mid_norm_in) allocate(pref_mid_norm(pver)) cappa = cappa_in cpair = cpair_in psurf_ref = psurf_ref_in pref_mid_norm = pref_mid_norm_in end subroutine held_suarez_1994_init subroutine held_suarez_1994(pcols, ncol, clat, pmid, & u, v, t, du, dv, s) !----------------------------------------------------------------------- ! ! Purpose: Implement idealized Held-Suarez forcings ! Held, I. M., and M. J. Suarez, 1994: 'A proposal for the ! intercomparison of the dynamical cores of atmospheric general ! circulation models.' ! Bulletin of the Amer. Meteor. Soc., vol. 75, pp. 1825-1830. ! !----------------------------------------------------------------------- ! ! Input arguments ! integer, intent(in) :: pcols ! Size of column dimension integer, intent(in) :: ncol ! Num active columns real(r8), intent(in) :: clat(pcols) ! latitudes(radians) for columns real(r8), intent(in) :: pmid(pcols,pver) ! mid-point pressure real(r8), intent(in) :: u(pcols,pver) ! Zonal wind (m/s) real(r8), intent(in) :: v(pcols,pver) ! Meridional wind (m/s) real(r8), intent(in) :: t(pcols,pver) ! Temperature (K) ! ! Output arguments ! real(r8), intent(out) :: du(pcols,pver) ! Zonal wind tend real(r8), intent(out) :: dv(pcols,pver) ! Meridional wind tend real(r8), intent(out) :: s(pcols,pver) ! Heating rate ! !---------------------------Local workspace----------------------------- ! integer :: i, k ! Longitude, level indices real(r8) :: kv ! 1./efolding_time (normalized) for wind real(r8) :: kt ! 1./efolding_time for temperature diss. real(r8) :: trefa ! "radiative equilibrium" T real(r8) :: trefc ! used in calc of "radiative equilibrium" T real(r8) :: cossq(ncol) ! coslat**2 real(r8) :: cossqsq(ncol) ! coslat**4 real(r8) :: sinsq(ncol) ! sinlat**2 real(r8) :: coslat(ncol) ! cosine(latitude) ! !----------------------------------------------------------------------- ! do i = 1, ncol coslat (i) = cos(clat(i)) sinsq (i) = sin(clat(i))*sin(clat(i)) cossq (i) = coslat(i)*coslat(i) cossqsq(i) = cossq (i)*cossq (i) end do ! !----------------------------------------------------------------------- ! ! Held/Suarez IDEALIZED physics algorithm: ! ! Held, I. M., and M. J. Suarez, 1994: A proposal for the ! intercomparison of the dynamical cores of atmospheric general ! circulation models. ! Bulletin of the Amer. Meteor. Soc., vol. 75, pp. 1825-1830. ! !----------------------------------------------------------------------- ! ! Compute idealized radiative heating rates (as dry static energy) ! ! do k = 1, pver if (pref_mid_norm(k) > sigmab) then do i = 1, ncol kt = ka + (ks - ka)*cossqsq(i)*(pref_mid_norm(k) - sigmab)/onemsig trefc = 315._r8 - (60._r8 * sinsq(i)) trefa = (trefc - 10._r8*cossq(i)*log((pmid(i,k)/psurf_ref)))*(pmid(i,k)/psurf_ref)**cappa trefa = max(t00,trefa) s(i,k) = (trefa - t(i,k))*kt*cpair end do else do i = 1, ncol trefc = 315._r8 - 60._r8*sinsq(i) trefa = (trefc - 10._r8*cossq(i)*log((pmid(i,k)/psurf_ref)))*(pmid(i,k)/psurf_ref)**cappa trefa = max(t00,trefa) s(i,k) = (trefa - t(i,k))*ka*cpair end do end if end do ! ! Add diffusion near the surface for the wind fields ! do k = 1, pver do i = 1, pcols du(i,k) = 0._r8 dv(i,k) = 0._r8 end do end do ! do k = 1, pver if (pref_mid_norm(k) > sigmab) then kv = kf*(pref_mid_norm(k) - sigmab)/onemsig do i = 1, ncol du(i,k) = -kv*u(i,k) dv(i,k) = -kv*v(i,k) end do end if end do end subroutine held_suarez_1994 end module held_suarez trefread/SourceMods/src.cam/runtime_opts.F90-ORIG0000644006307300054160000001666613047377062021222 0ustar islascgdampmodule runtime_opts !----------------------------------------------------------------------- ! ! Provide driver level routine for making calls to the namelist readers ! for the infrastructure and the dycore and physics parameterizations. ! !----------------------------------------------------------------------- use shr_kind_mod, only: r8=>shr_kind_r8 implicit none private save public :: read_namelist !======================================================================= contains !======================================================================= subroutine read_namelist(nlfilename, single_column, scmlat, scmlon) use cam_initfiles, only: cam_initfiles_readnl use constituents, only: cnst_readnl use phys_grid, only: phys_grid_readnl use chem_surfvals, only: chem_surfvals_readnl use check_energy, only: check_energy_readnl use radiation, only: radiation_readnl use carma_flags_mod, only: carma_readnl use co2_cycle, only: co2_cycle_readnl use scamMod, only: scam_readnl use spmd_utils, only: spmd_utils_readnl use cam_history, only: history_readnl use physconst, only: physconst_readnl use physics_buffer, only: pbuf_readnl use phys_control, only: phys_ctl_readnl use wv_saturation, only: wv_sat_readnl use ref_pres, only: ref_pres_readnl use cam3_aero_data, only: cam3_aero_data_readnl use cam3_ozone_data, only: cam3_ozone_data_readnl use dadadj_cam, only: dadadj_readnl use macrop_driver, only: macrop_driver_readnl use microp_driver, only: microp_driver_readnl use microp_aero, only: microp_aero_readnl use subcol, only: subcol_readnl use cloud_fraction, only: cldfrc_readnl use cldfrc2m, only: cldfrc2m_readnl use rk_stratiform, only: rk_stratiform_readnl use unicon_cam, only: unicon_cam_readnl use zm_conv_intr, only: zm_conv_readnl use hk_conv, only: hkconv_readnl use uwshcu, only: uwshcu_readnl use pkg_cld_sediment, only: cld_sediment_readnl use gw_drag, only: gw_drag_readnl use qbo, only: qbo_readnl use iondrag, only: iondrag_readnl use phys_debug_util, only: phys_debug_readnl use conv_water, only: conv_water_readnl use rad_constituents, only: rad_cnst_readnl use radiation_data, only: rad_data_readnl use modal_aer_opt, only: modal_aer_opt_readnl use clubb_intr, only: clubb_readnl use chemistry, only: chem_readnl use prescribed_volcaero, only: prescribed_volcaero_readnl use prescribed_strataero,only: prescribed_strataero_readnl use aerodep_flx, only: aerodep_flx_readnl use solar_data, only: solar_data_readnl use solar_euv_data, only: solar_euv_data_readnl use tropopause, only: tropopause_readnl use aoa_tracers, only: aoa_tracers_readnl use prescribed_ozone, only: prescribed_ozone_readnl use prescribed_aero, only: prescribed_aero_readnl use prescribed_ghg, only: prescribed_ghg_readnl use aircraft_emit, only: aircraft_emit_readnl use cospsimulator_intr, only: cospsimulator_intr_readnl use vertical_diffusion, only: vd_readnl use rayleigh_friction, only: rayleigh_friction_readnl use cam_diagnostics, only: diag_readnl use radheat, only: radheat_readnl #if ( defined OFFLINE_DYN ) use metdata, only: metdata_readnl #endif use offline_driver, only: offline_driver_readnl use rate_diags, only: rate_diags_readnl use tracers, only: tracers_readnl use dyn_comp, only: dyn_readnl !---------------------------Arguments----------------------------------- character(len=*), intent(in) :: nlfilename logical, intent(in) :: single_column real(r8), intent(in) :: scmlat real(r8), intent(in) :: scmlon !---------------------------Local variables----------------------------- character(len=*), parameter :: subname = "read_namelist" !----------------------------------------------------------------------- ! Call subroutines for modules to read their own namelist. ! In some cases namelist default values may depend on settings from ! other modules, so there may be an order dependence in the following ! calls. ! ***N.B.*** In particular, physconst_readnl should be called before ! the other readnl methods in case that method is used to set ! physical constants, some of which are set at runtime ! by the physconst_readnl method. ! Modules that read their own namelist are responsible for making sure ! all processes receive the values. call spmd_utils_readnl(nlfilename) call phys_grid_readnl(nlfilename) call physconst_readnl(nlfilename) !++bee 13 Oct 2015, need to fix the pbuf_global_allocate functionality, then ! can uncomment the pbuf_readnl line ! call pbuf_readnl(nlfilename) call cam_initfiles_readnl(nlfilename) call cnst_readnl(nlfilename) call history_readnl(nlfilename) call chem_surfvals_readnl(nlfilename) call phys_ctl_readnl(nlfilename) call wv_sat_readnl(nlfilename) call ref_pres_readnl(nlfilename) call cam3_aero_data_readnl(nlfilename) call cam3_ozone_data_readnl(nlfilename) call dadadj_readnl(nlfilename) call macrop_driver_readnl(nlfilename) call microp_driver_readnl(nlfilename) call microp_aero_readnl(nlfilename) call clubb_readnl(nlfilename) call subcol_readnl(nlfilename) call cldfrc_readnl(nlfilename) call cldfrc2m_readnl(nlfilename) call unicon_cam_readnl(nlfilename) call zm_conv_readnl(nlfilename) call rk_stratiform_readnl(nlfilename) call hkconv_readnl(nlfilename) call uwshcu_readnl(nlfilename) call cld_sediment_readnl(nlfilename) call gw_drag_readnl(nlfilename) call qbo_readnl(nlfilename) call iondrag_readnl(nlfilename) call phys_debug_readnl(nlfilename) call conv_water_readnl(nlfilename) call radiation_readnl(nlfilename) call rad_cnst_readnl(nlfilename) call rad_data_readnl(nlfilename) call modal_aer_opt_readnl(nlfilename) call chem_readnl(nlfilename) call prescribed_volcaero_readnl(nlfilename) call prescribed_strataero_readnl(nlfilename) call solar_data_readnl(nlfilename) call solar_euv_data_readnl(nlfilename) call carma_readnl(nlfilename) call tropopause_readnl(nlfilename) call aoa_tracers_readnl(nlfilename) call tracers_readnl(nlfilename) call aerodep_flx_readnl(nlfilename) call prescribed_ozone_readnl(nlfilename) call prescribed_aero_readnl(nlfilename) call prescribed_ghg_readnl(nlfilename) call co2_cycle_readnl(nlfilename) call aircraft_emit_readnl(nlfilename) call cospsimulator_intr_readnl(nlfilename) call diag_readnl(nlfilename) call check_energy_readnl(nlfilename) call radheat_readnl(nlfilename) call vd_readnl(nlfilename) call rayleigh_friction_readnl(nlfilename) #if ( defined OFFLINE_DYN ) call metdata_readnl(nlfilename) #endif call offline_driver_readnl(nlfilename) call rate_diags_readnl(nlfilename) call scam_readnl(nlfilename, single_column, scmlat, scmlon) call dyn_readnl(nlfilename) end subroutine read_namelist !======================================================================= end module runtime_opts trefread/SourceMods/src.cam/held_suarez.F900000644006307300054160000003364513047377052020234 0ustar islascgdampmodule held_suarez ! --------------------------------------------------------------------- ! Contains modification to read in relaxtation temperature profile from ! file. ! Isla Simpson 10th Feb 2017 ! --------------------------------------------------------------------- !IRS !use shr_kind_mod, only: r8 => shr_kind_r8 use shr_kind_mod,only:r8=>SHR_KIND_R8,cs=>SHR_KIND_CS,cl=>SHR_KIND_CL use cam_abortutils ,only:endrun use spmd_utils ,only:masterproc use cam_logfile ,only:iulog use spmd_utils, only: mpicom,mstrid=>masterprocid, mpi_integer, mpi_real8, & mpi_logical, mpi_character use dyn_grid ,only: get_horiz_grid_dim_d use ppgrid ,only: pcols,begchunk,endchunk ! END IRS implicit none private save public :: held_suarez_1994_init public :: held_suarez_1994 !IRS public :: trefread_readnl !END IRS !! !! Forcing parameters !! real(r8), parameter :: efoldf = 1._r8 ! efolding time for wind dissipation real(r8), parameter :: efolda = 40._r8 ! efolding time for T dissipation real(r8), parameter :: efolds = 4._r8 ! efolding time for T dissipation real(r8), parameter :: sigmab = 0.7_r8 ! threshold sigma level real(r8), parameter :: t00 = 200._r8 ! minimum reference temperature real(r8), parameter :: kf = 1._r8/(86400._r8*efoldf) ! 1./efolding_time for wind dissipation real(r8), parameter :: onemsig = 1._r8 - sigmab ! 1. - sigma_reference real(r8), parameter :: ka = 1._r8/(86400._r8 * efolda) ! 1./efolding_time for temperature diss. real(r8), parameter :: ks = 1._r8/(86400._r8 * efolds) !! !! Model constants, reset in init call !! real(r8) :: cappa = 2.0_r8 / 7.0_r8 ! R/Cp real(r8) :: cpair = 1004.0_r8 ! specific heat of dry air (J/K/kg) real(r8) :: psurf_ref = 0.0_r8 ! Surface pressure ! pref_mid_norm are layer midpoints normalized by surface pressure ('eta' coordinate) real(r8), allocatable :: pref_mid_norm(:) integer :: pver ! Num vertical levels ! IRS character(len=cl) :: treffile logical:: treffromfile real(r8),allocatable::tref(:,:,:) ! relaxation temperature read from file ! END IRS !======================================================================= contains !======================================================================= subroutine held_suarez_1994_init(cappa_in, cpair_in, psurf_ref_in, pref_mid_norm_in) !! Dummy arguments use error_messages,only: alloc_err real(r8), intent(in) :: cappa_in real(r8), intent(in) :: cpair_in real(r8), intent(in) :: psurf_ref_in real(r8), intent(in) :: pref_mid_norm_in(:) integer :: istat pver = size(pref_mid_norm_in) allocate(pref_mid_norm(pver)) cappa = cappa_in cpair = cpair_in psurf_ref = psurf_ref_in pref_mid_norm = pref_mid_norm_in allocate(tref(pcols,pver,begchunk:endchunk),stat=istat) call alloc_err(istat,'UPDATE_TREF_EUL','tref',pcols*pver*(endchunk-begchunk)) tref(:pcols,:pver,begchunk:endchunk)=0._r8 !IRS if (treffromfile) then call update_tref_eul(treffile) endif !END IRS ! allocate(tref(pcols,pver,begchunk:endchunk,1),stat=istat) ! call ! alloc_err(istat,'UPDATE_TREF_EUL','tref',pcols*pver*(endchunk-begchunk)) ! tref(:pcols,:pver,begchunk:endchunk,1)=0._r8 end subroutine held_suarez_1994_init subroutine held_suarez_1994(pcols, ncol, clat, pmid, & u, v, t, du, dv, s, lchnk) !----------------------------------------------------------------------- ! ! Purpose: Implement idealized Held-Suarez forcings ! Held, I. M., and M. J. Suarez, 1994: 'A proposal for the ! intercomparison of the dynamical cores of atmospheric general ! circulation models.' ! Bulletin of the Amer. Meteor. Soc., vol. 75, pp. 1825-1830. ! !----------------------------------------------------------------------- ! ! Input arguments ! integer, intent(in) :: pcols ! Size of column dimension integer, intent(in) :: ncol ! Num active columns real(r8), intent(in) :: clat(pcols) ! latitudes(radians) for columns real(r8), intent(in) :: pmid(pcols,pver) ! mid-point pressure real(r8), intent(in) :: u(pcols,pver) ! Zonal wind (m/s) real(r8), intent(in) :: v(pcols,pver) ! Meridional wind (m/s) real(r8), intent(in) :: t(pcols,pver) ! Temperature (K) ! ! Output arguments ! real(r8), intent(out) :: du(pcols,pver) ! Zonal wind tend real(r8), intent(out) :: dv(pcols,pver) ! Meridional wind tend real(r8), intent(out) :: s(pcols,pver) ! Heating rate !IRS integer, intent(in) :: lchnk ! chuhk identifier ! END IRS ! !---------------------------Local workspace----------------------------- ! integer :: i, k ! Longitude, level indices real(r8) :: kv ! 1./efolding_time (normalized) for wind real(r8) :: kt ! 1./efolding_time for temperature diss. real(r8) :: trefa ! "radiative equilibrium" T real(r8) :: trefc ! used in calc of "radiative equilibrium" T real(r8) :: cossq(ncol) ! coslat**2 real(r8) :: cossqsq(ncol) ! coslat**4 real(r8) :: sinsq(ncol) ! sinlat**2 real(r8) :: coslat(ncol) ! cosine(latitude) !IRS ! params for reading in Tref from file ! character(len=*) treffile ! logical treffromfile ! debugging integer :: icol,ichunk,ip !END IRS ! !----------------------------------------------------------------------- ! do i = 1, ncol coslat (i) = cos(clat(i)) sinsq (i) = sin(clat(i))*sin(clat(i)) cossq (i) = coslat(i)*coslat(i) cossqsq(i) = cossq (i)*cossq (i) end do ! !----------------------------------------------------------------------- ! ! Held/Suarez IDEALIZED physics algorithm: ! ! Held, I. M., and M. J. Suarez, 1994: A proposal for the ! intercomparison of the dynamical cores of atmospheric general ! circulation models. ! Bulletin of the Amer. Meteor. Soc., vol. 75, pp. 1825-1830. ! !----------------------------------------------------------------------- ! ! Compute idealized radiative heating rates (as dry static energy) ! ! ! IRS ! read in relaxation temperature profile from file if (treffromfile) then do k = 1, pver if (pref_mid_norm(k) > sigmab) then do i = 1, ncol kt = ka + (ks - ka)*cossqsq(i)*(pref_mid_norm(k) - sigmab)/onemsig trefa=tref(i,k,lchnk) s(i,k) = (trefa - t(i,k))*kt*cpair end do else do i = 1, ncol trefa=tref(i,k,lchnk) s(i,k) = (trefa - t(i,k))*ka*cpair end do end if end do ! end if else ! END IRS do k = 1, pver if (pref_mid_norm(k) > sigmab) then do i = 1, ncol kt = ka + (ks - ka)*cossqsq(i)*(pref_mid_norm(k) - sigmab)/onemsig trefc = 315._r8 - (60._r8 * sinsq(i)) trefa = (trefc - 10._r8*cossq(i)*log((pmid(i,k)/psurf_ref)))*(pmid(i,k)/psurf_ref)**cappa trefa = max(t00,trefa) s(i,k) = (trefa - t(i,k))*kt*cpair end do else do i = 1, ncol trefc = 315._r8 - 60._r8*sinsq(i) trefa = (trefc - 10._r8*cossq(i)*log((pmid(i,k)/psurf_ref)))*(pmid(i,k)/psurf_ref)**cappa trefa = max(t00,trefa) s(i,k) = (trefa - t(i,k))*ka*cpair end do end if end do end if !treffromfile ! ! Add diffusion near the surface for the wind fields ! do k = 1, pver do i = 1, pcols du(i,k) = 0._r8 dv(i,k) = 0._r8 end do end do ! do k = 1, pver if (pref_mid_norm(k) > sigmab) then kv = kf*(pref_mid_norm(k) - sigmab)/onemsig do i = 1, ncol du(i,k) = -kv*u(i,k) dv(i,k) = -kv*v(i,k) end do end if end do end subroutine held_suarez_1994 ! IRS - read in ideal physics namelist parameters subroutine trefread_readnl(nlfile) use shr_kind_mod,only:r8=>SHR_KIND_R8,cs=>SHR_KIND_CS,cl=>SHR_KIND_CL use namelist_utils ,only:find_group_name use units ,only:getunit,freeunit character(len=*),intent(in)::nlfile integer :: unitn, ierr character(len=*), parameter :: sub = 'trefread_readnl' namelist /trefread_nl/ treffromfile,treffile !Set default namelist values treffromfile=.False. treffile='randomfile' !Read in namelist values if (masterproc) then unitn = getunit() open(unitn,file=trim(nlfile),status='old') call find_group_name(unitn,'trefread_nl',status=ierr) if (ierr.eq.0) then read(unitn,trefread_nl,iostat=ierr) if (ierr.ne.0) then call endrun('trefread_readnl:: ERROR reading namelist') endif endif close(unitn) call freeunit(unitn) endif call mpi_bcast(treffromfile , 1, mpi_logical, mstrid, mpicom, ierr) if (ierr /= 0) call endrun(sub//": FATAL: mpi_bcast: held_suarez_1994") call mpi_bcast(treffile , len(treffile), mpi_character, mstrid, mpicom, ierr) if (ierr /= 0) call endrun(sub//": FATAL: mpi_bcast: held_suarez_1994") return end subroutine ! END IRS ! IRS subroutine update_tref_eul(treffilein) ! read in user specified relaxation temperature profile from file. use netcdf use dyn_grid ,only: get_horiz_grid_dim_d use error_messages,only: alloc_err use phys_grid ,only: scatter_field_to_chunk character(len=*), intent(in)::treffilein real(r8),allocatable::lonin(:),latin(:) real(r8),allocatable::trefin(:,:,:),trefintrans(:,:,:) integer istat,ncid,varid integer nlon,nlat,nlonin,nlatin,plevin integer ilev,ilat,ilon call get_horiz_grid_dim_d(nlon,nlat) if (masterproc) then allocate(lonin(nlon),stat=istat) call alloc_err(istat,'UPDATE_TREF_EUL','lonin',nlon) allocate(latin(nlat),stat=istat) call alloc_err(istat,'UPDATE_TREF_EUL','latin',nlat) allocate(trefin(nlon,nlat,pver)) call alloc_err(istat,'UPDATE_TREF_EUL','trefin',nlon*nlat*pver) allocate(trefintrans(nlon,pver,nlat)) call alloc_err(istat,'UPDATE_TREF_EUL','treintrans',nlon*nlat*pver) ! open tref file istat=nf90_open(trim(treffilein),NF90_NOWRITE,ncid) if(istat.ne.NF90_NOERR) then write(iulog,*)'NF90_OPEN: failed for file ',trim(treffile) write(iulog,*) nf90_strerror(istat) call endrun ('UPDATE_TREF_EUL') endif ! read in dimensions istat=nf90_inq_dimid(ncid,'lon',varid) if(istat.ne.NF90_NOERR) then write(iulog,*) nf90_strerror(istat) call endrun ('UPDATE_TREF_EUL') endif istat=nf90_inquire_dimension(ncid,varid,len=nlonin) if(istat.ne.NF90_NOERR) then write(iulog,*) nf90_strerror(istat) call endrun ('UPDATE_TREF_EUL') endif istat=nf90_inq_dimid(ncid,'lat',varid) if(istat.ne.NF90_NOERR) then write(iulog,*) nf90_strerror(istat) call endrun ('UPDATE_TREF_EUL') endif istat=nf90_inquire_dimension(ncid,varid,len=nlatin) if(istat.ne.NF90_NOERR) then write(iulog,*) nf90_strerror(istat) call endrun ('UPDATE_TREF_EUL') endif istat=nf90_inq_dimid(ncid,'lev',varid) if(istat.ne.NF90_NOERR) then write(iulog,*) nf90_strerror(istat) call endrun ('UPDATE_TREF_EUL') endif istat=nf90_inquire_dimension(ncid,varid,len=plevin) if(istat.ne.NF90_NOERR) then write(iulog,*) nf90_strerror(istat) call endrun ('UPDATE_TREF_EUL') endif istat=nf90_inq_varid(ncid,'lon',varid) if(istat.ne.NF90_NOERR) then write(iulog,*) nf90_strerror(istat) call endrun ('UPDATE_TREF_EUL') endif istat=nf90_get_var(ncid,varid,lonin) if(istat.ne.NF90_NOERR) then write(iulog,*) nf90_strerror(istat) call endrun ('UPDATE_TREF_EUL') endif istat=nf90_inq_varid(ncid,'lat',varid) if(istat.ne.NF90_NOERR) then write(iulog,*) nf90_strerror(istat) call endrun ('UPDATE_TREF_EUL') endif istat=nf90_get_var(ncid,varid,latin) if(istat.ne.NF90_NOERR) then write(iulog,*) nf90_strerror(istat) call endrun ('UPDATE_TREF_EUL') endif ! check dimensionality of input file if ((nlonin.ne.nlon).or.(nlatin.ne.nlat).or.(plevin.ne.pver)) then write(iulog,*) 'ERROR: update_tref_eul: nlon=',nlon,' nlonin=',nlonin write(iulog,*) 'ERROR: update_tref_eul: nlat=',nlat,' nlatin=',nlatin write(iulog,*) 'ERROR: update_tref_eul: plevin=',plevin,' pver=',pver call endrun('update_tref_eul: incorrect input tref dimensions') endif ! Read in tref file, transpose lat/lev indices and scatter istat=nf90_inq_varid(ncid,'tref',varid) if (istat.ne.NF90_NOERR) then write(iulog,*) nf90_strerror(istat) call endrun ('UPDATE_TREF_EUL') endif istat=nf90_get_var(ncid,varid,trefin) if (istat.ne.NF90_NOERR) then write(iulog,*) nf90_strerror(istat) call endrun ('UPDATE_TREF_EUL') endif do ilat=1,nlatin do ilev=1,plevin do ilon=1,nlonin trefintrans(ilon,ilev,ilat)=trefin(ilon,ilat,ilev) end do end do end do endif ! end (masterproc) call scatter_field_to_chunk(1,pver,1,nlon,trefintrans,tref(1,1,begchunk)) ! close tref file and deallocate arrays if (masterproc) then istat=nf90_close(ncid) if (istat.ne.NF90_NOERR) then write(iulog,*)'NF90_CLOSE: failed for file ',trim(treffile) write(iulog,*) nf90_strerror(istat) call endrun ('UPDATE_TREF_EUL') endif deallocate(lonin) deallocate(latin) deallocate(trefin) deallocate(trefintrans) endif ! (masterproc) end subroutine end module held_suarez trefread/SourceMods/src.cam/held_suarez_cam.F90-ORIG0000644006307300054160000001115713047376402021602 0ustar islascgdamp#define MODHS 1 #undef MODHS module held_suarez_cam use shr_kind_mod, only: r8 => shr_kind_r8 use ppgrid, only: pcols, pver implicit none private save public :: held_suarez_init, held_suarez_tend real(r8), parameter :: efoldf = 1._r8 ! efolding time for wind dissipation real(r8), parameter :: efolda = 40._r8 ! efolding time for T dissipation real(r8), parameter :: efolds = 4._r8 ! efolding time for T dissipation real(r8), parameter :: sigmab = 0.7_r8 ! threshold sigma level real(r8), parameter :: t00 = 200._r8 ! minimum reference temperature real(r8), parameter :: kf = 1._r8/(86400._r8*efoldf) ! 1./efolding_time for wind dissipation real(r8), parameter :: onemsig = 1._r8 - sigmab ! 1. - sigma_reference real(r8), parameter :: ka = 1._r8/(86400._r8 * efolda) ! 1./efolding_time for temperature diss. real(r8), parameter :: ks = 1._r8/(86400._r8 * efolds) !======================================================================= contains !======================================================================= subroutine held_suarez_init(pbuf2d) use physics_buffer, only: physics_buffer_desc use cam_history, only: addfld, add_default use physconst, only: cappa, cpair use ref_pres, only: pref_mid_norm, psurf_ref use held_suarez, only: held_suarez_1994_init type(physics_buffer_desc), pointer :: pbuf2d(:,:) ! Set model constant values call held_suarez_1994_init(cappa, cpair, psurf_ref, pref_mid_norm) ! This field is added by radiation when full physics is used call addfld('QRS', (/ 'lev' /), 'A', 'K/s', & 'Temperature tendency associated with the relaxation toward the equilibrium temperature profile') call add_default('QRS', 1, ' ') end subroutine held_suarez_init subroutine held_suarez_tend(state, ptend, ztodt) !----------------------------------------------------------------------- ! ! Purpose: ! algorithm 1: Held/Suarez IDEALIZED physics ! algorithm 2: Held/Suarez IDEALIZED physics (Williamson modified stratosphere ! algorithm 3: Held/Suarez IDEALIZED physics (Lin/Williamson modified strato/meso-sphere ! ! Author: J. Olson ! !----------------------------------------------------------------------- use physconst, only: cpairv use phys_grid, only: get_rlat_all_p use physics_types, only: physics_state, physics_ptend use physics_types, only: physics_ptend_init use cam_abortutils, only: endrun use cam_history, only: outfld use time_manager, only: get_nstep use held_suarez, only: held_suarez_1994 ! ! Input arguments ! type(physics_state), intent(inout) :: state real(r8), intent(in) :: ztodt ! Two times model timestep (2 delta-t) ! ! Output argument ! type(physics_ptend), intent(out) :: ptend ! Package tendencies ! !---------------------------Local workspace----------------------------- integer :: lchnk ! chunk identifier integer :: ncol ! number of atmospheric columns integer :: nstep ! timestep number real(r8) :: clat(pcols) ! latitudes(radians) for columns real(r8) :: pmid(pcols,pver) ! mid-point pressure integer :: i, k ! Longitude, level indices ! !----------------------------------------------------------------------- ! lchnk = state%lchnk ncol = state%ncol nstep = get_nstep() call get_rlat_all_p(lchnk, ncol, clat) do k = 1, pver do i = 1, ncol pmid(i,k) = state%pmid(i,k) end do end do ! initialize individual parameterization tendencies call physics_ptend_init(ptend, state%psetcols, 'held_suarez', ls=.true., lu=.true., lv=.true.) call held_suarez_1994(pcols, ncol, clat, pmid, & state%u, state%v, state%t, ptend%u, ptend%v, ptend%s) ! Note, we assume that there are no subcolumns in simple physics pmid(:ncol,:) = ptend%s(:ncol, :)/cpairv(:ncol,:,lchnk) if (pcols > ncol) then pmid(ncol+1:,:) = 0.0_r8 end if call outfld('QRS', pmid, pcols, lchnk) end subroutine held_suarez_tend end module held_suarez_cam trefread/SourceMods/src.cam/runtime_opts.F900000644006307300054160000001707113047377074020456 0ustar islascgdampmodule runtime_opts !----------------------------------------------------------------------- ! ! Provide driver level routine for making calls to the namelist readers ! for the infrastructure and the dycore and physics parameterizations. ! !----------------------------------------------------------------------- use shr_kind_mod, only: r8=>shr_kind_r8 implicit none private save public :: read_namelist !======================================================================= contains !======================================================================= subroutine read_namelist(nlfilename, single_column, scmlat, scmlon) use cam_initfiles, only: cam_initfiles_readnl use constituents, only: cnst_readnl use phys_grid, only: phys_grid_readnl use chem_surfvals, only: chem_surfvals_readnl use check_energy, only: check_energy_readnl use radiation, only: radiation_readnl use carma_flags_mod, only: carma_readnl use co2_cycle, only: co2_cycle_readnl use scamMod, only: scam_readnl use spmd_utils, only: spmd_utils_readnl use cam_history, only: history_readnl use physconst, only: physconst_readnl use physics_buffer, only: pbuf_readnl use phys_control, only: phys_ctl_readnl use wv_saturation, only: wv_sat_readnl use ref_pres, only: ref_pres_readnl use cam3_aero_data, only: cam3_aero_data_readnl use cam3_ozone_data, only: cam3_ozone_data_readnl use dadadj_cam, only: dadadj_readnl use macrop_driver, only: macrop_driver_readnl use microp_driver, only: microp_driver_readnl use microp_aero, only: microp_aero_readnl use subcol, only: subcol_readnl use cloud_fraction, only: cldfrc_readnl use cldfrc2m, only: cldfrc2m_readnl use rk_stratiform, only: rk_stratiform_readnl use unicon_cam, only: unicon_cam_readnl use zm_conv_intr, only: zm_conv_readnl use hk_conv, only: hkconv_readnl use uwshcu, only: uwshcu_readnl use pkg_cld_sediment, only: cld_sediment_readnl use gw_drag, only: gw_drag_readnl use qbo, only: qbo_readnl use iondrag, only: iondrag_readnl use phys_debug_util, only: phys_debug_readnl use conv_water, only: conv_water_readnl use rad_constituents, only: rad_cnst_readnl use radiation_data, only: rad_data_readnl use modal_aer_opt, only: modal_aer_opt_readnl use clubb_intr, only: clubb_readnl use chemistry, only: chem_readnl use prescribed_volcaero, only: prescribed_volcaero_readnl use prescribed_strataero,only: prescribed_strataero_readnl use aerodep_flx, only: aerodep_flx_readnl use solar_data, only: solar_data_readnl use solar_euv_data, only: solar_euv_data_readnl use tropopause, only: tropopause_readnl use aoa_tracers, only: aoa_tracers_readnl use prescribed_ozone, only: prescribed_ozone_readnl use prescribed_aero, only: prescribed_aero_readnl use prescribed_ghg, only: prescribed_ghg_readnl use aircraft_emit, only: aircraft_emit_readnl use cospsimulator_intr, only: cospsimulator_intr_readnl use vertical_diffusion, only: vd_readnl use rayleigh_friction, only: rayleigh_friction_readnl use cam_diagnostics, only: diag_readnl use radheat, only: radheat_readnl #if ( defined OFFLINE_DYN ) use metdata, only: metdata_readnl #endif use offline_driver, only: offline_driver_readnl use rate_diags, only: rate_diags_readnl use tracers, only: tracers_readnl use dyn_comp, only: dyn_readnl !IRS use held_suarez, only: trefread_readnl !END IRS !---------------------------Arguments----------------------------------- character(len=*), intent(in) :: nlfilename logical, intent(in) :: single_column real(r8), intent(in) :: scmlat real(r8), intent(in) :: scmlon !---------------------------Local variables----------------------------- character(len=*), parameter :: subname = "read_namelist" !----------------------------------------------------------------------- ! Call subroutines for modules to read their own namelist. ! In some cases namelist default values may depend on settings from ! other modules, so there may be an order dependence in the following ! calls. ! ***N.B.*** In particular, physconst_readnl should be called before ! the other readnl methods in case that method is used to set ! physical constants, some of which are set at runtime ! by the physconst_readnl method. ! Modules that read their own namelist are responsible for making sure ! all processes receive the values. call spmd_utils_readnl(nlfilename) call phys_grid_readnl(nlfilename) call physconst_readnl(nlfilename) !++bee 13 Oct 2015, need to fix the pbuf_global_allocate functionality, then ! can uncomment the pbuf_readnl line ! call pbuf_readnl(nlfilename) call cam_initfiles_readnl(nlfilename) call cnst_readnl(nlfilename) call history_readnl(nlfilename) call chem_surfvals_readnl(nlfilename) call phys_ctl_readnl(nlfilename) call wv_sat_readnl(nlfilename) call ref_pres_readnl(nlfilename) call cam3_aero_data_readnl(nlfilename) call cam3_ozone_data_readnl(nlfilename) call dadadj_readnl(nlfilename) call macrop_driver_readnl(nlfilename) call microp_driver_readnl(nlfilename) call microp_aero_readnl(nlfilename) call clubb_readnl(nlfilename) call subcol_readnl(nlfilename) call cldfrc_readnl(nlfilename) call cldfrc2m_readnl(nlfilename) call unicon_cam_readnl(nlfilename) call zm_conv_readnl(nlfilename) call rk_stratiform_readnl(nlfilename) call hkconv_readnl(nlfilename) call uwshcu_readnl(nlfilename) call cld_sediment_readnl(nlfilename) call gw_drag_readnl(nlfilename) call qbo_readnl(nlfilename) call iondrag_readnl(nlfilename) call phys_debug_readnl(nlfilename) call conv_water_readnl(nlfilename) call radiation_readnl(nlfilename) call rad_cnst_readnl(nlfilename) call rad_data_readnl(nlfilename) call modal_aer_opt_readnl(nlfilename) call chem_readnl(nlfilename) call prescribed_volcaero_readnl(nlfilename) call prescribed_strataero_readnl(nlfilename) call solar_data_readnl(nlfilename) call solar_euv_data_readnl(nlfilename) call carma_readnl(nlfilename) call tropopause_readnl(nlfilename) call aoa_tracers_readnl(nlfilename) call tracers_readnl(nlfilename) call aerodep_flx_readnl(nlfilename) call prescribed_ozone_readnl(nlfilename) call prescribed_aero_readnl(nlfilename) call prescribed_ghg_readnl(nlfilename) call co2_cycle_readnl(nlfilename) call aircraft_emit_readnl(nlfilename) call cospsimulator_intr_readnl(nlfilename) call diag_readnl(nlfilename) call check_energy_readnl(nlfilename) call radheat_readnl(nlfilename) call vd_readnl(nlfilename) call rayleigh_friction_readnl(nlfilename) #if ( defined OFFLINE_DYN ) call metdata_readnl(nlfilename) #endif call offline_driver_readnl(nlfilename) call rate_diags_readnl(nlfilename) call scam_readnl(nlfilename, single_column, scmlat, scmlon) call dyn_readnl(nlfilename) ! IRS call trefread_readnl(nlfilename) ! END IRS end subroutine read_namelist !======================================================================= end module runtime_opts trefread/SourceMods/src.cam/held_suarez_cam.F900000644006307300054160000001116613047376676021061 0ustar islascgdamp#define MODHS 1 #undef MODHS module held_suarez_cam use shr_kind_mod, only: r8 => shr_kind_r8 use ppgrid, only: pcols, pver implicit none private save public :: held_suarez_init, held_suarez_tend real(r8), parameter :: efoldf = 1._r8 ! efolding time for wind dissipation real(r8), parameter :: efolda = 40._r8 ! efolding time for T dissipation real(r8), parameter :: efolds = 4._r8 ! efolding time for T dissipation real(r8), parameter :: sigmab = 0.7_r8 ! threshold sigma level real(r8), parameter :: t00 = 200._r8 ! minimum reference temperature real(r8), parameter :: kf = 1._r8/(86400._r8*efoldf) ! 1./efolding_time for wind dissipation real(r8), parameter :: onemsig = 1._r8 - sigmab ! 1. - sigma_reference real(r8), parameter :: ka = 1._r8/(86400._r8 * efolda) ! 1./efolding_time for temperature diss. real(r8), parameter :: ks = 1._r8/(86400._r8 * efolds) !======================================================================= contains !======================================================================= subroutine held_suarez_init(pbuf2d) use physics_buffer, only: physics_buffer_desc use cam_history, only: addfld, add_default use physconst, only: cappa, cpair use ref_pres, only: pref_mid_norm, psurf_ref use held_suarez, only: held_suarez_1994_init type(physics_buffer_desc), pointer :: pbuf2d(:,:) ! Set model constant values call held_suarez_1994_init(cappa, cpair, psurf_ref, pref_mid_norm) ! This field is added by radiation when full physics is used call addfld('QRS', (/ 'lev' /), 'A', 'K/s', & 'Temperature tendency associated with the relaxation toward the equilibrium temperature profile') call add_default('QRS', 1, ' ') end subroutine held_suarez_init subroutine held_suarez_tend(state, ptend, ztodt) !----------------------------------------------------------------------- ! ! Purpose: ! algorithm 1: Held/Suarez IDEALIZED physics ! algorithm 2: Held/Suarez IDEALIZED physics (Williamson modified stratosphere ! algorithm 3: Held/Suarez IDEALIZED physics (Lin/Williamson modified strato/meso-sphere ! ! Author: J. Olson ! !----------------------------------------------------------------------- use physconst, only: cpairv use phys_grid, only: get_rlat_all_p use physics_types, only: physics_state, physics_ptend use physics_types, only: physics_ptend_init use cam_abortutils, only: endrun use cam_history, only: outfld use time_manager, only: get_nstep use held_suarez, only: held_suarez_1994 ! ! Input arguments ! type(physics_state), intent(inout) :: state real(r8), intent(in) :: ztodt ! Two times model timestep (2 delta-t) ! ! Output argument ! type(physics_ptend), intent(out) :: ptend ! Package tendencies ! !---------------------------Local workspace----------------------------- integer :: lchnk ! chunk identifier integer :: ncol ! number of atmospheric columns integer :: nstep ! timestep number real(r8) :: clat(pcols) ! latitudes(radians) for columns real(r8) :: pmid(pcols,pver) ! mid-point pressure integer :: i, k ! Longitude, level indices ! !----------------------------------------------------------------------- ! lchnk = state%lchnk ncol = state%ncol nstep = get_nstep() call get_rlat_all_p(lchnk, ncol, clat) do k = 1, pver do i = 1, ncol pmid(i,k) = state%pmid(i,k) end do end do ! initialize individual parameterization tendencies call physics_ptend_init(ptend, state%psetcols, 'held_suarez', ls=.true., lu=.true., lv=.true.) call held_suarez_1994(pcols, ncol, clat, pmid, & state%u, state%v, state%t, ptend%u, ptend%v, ptend%s, lchnk) ! Note, we assume that there are no subcolumns in simple physics pmid(:ncol,:) = ptend%s(:ncol, :)/cpairv(:ncol,:,lchnk) if (pcols > ncol) then pmid(ncol+1:,:) = 0.0_r8 end if call outfld('QRS', pmid, pcols, lchnk) end subroutine held_suarez_tend end module held_suarez_cam