Radiative transfer: Simple gray and multi-band radiation schemes.
Documentation for two_stream_gray_rad_mod
.
Summary
This module includes different simple configurations for solving the two stream radiative transfer equations. There are three semi-gray configurations available (gray in shortwave and longwave), and one intermediate multi-band configuration (gray in shortwave, multi-band in longwave). These configurations are selected by setting the rad_scheme
namelist option.
|
Short description |
Reference |
|
|
Semi-gray scheme with prescribed longwave and shortwave
optical depths.
|
||
|
Semi-gray scheme with longwave optical depth dependent on
water vapour content and
optical depth is prescribed.
|
||
|
Multi-band scheme with two longwave bands and one shortwave
band. One longwave band corresponds to an infrared window
region (
infrared wavelengths (
depths depend on water vapour content and
Note: the Geen scheme is currently only stable when run in a fixed SST
configuration. A version stable with a slab ocean is expected by 30/11/2020.
|
||
|
Semi-gray scheme for use in giant planet simulations. Longwave
and shortwave optical depths are prescribed. Does not require a
surface temperature as input, and allows specification of an
interior heat flux.
|
For each of these schemes diurnally-and-seasonally averaged, and time-dependent options for the incoming solar radiation are available.
In the next few sections we will describe the main options available to the user for solving the radiative transfer equations, and prescribing the insolation. Namelist options and their defaults, and available output diagnostics are summarised at the end of the documentation.
Frierson/Byrne schemes
Each of these schemes solve the same equations to compute the longwave and shortwave fluxes. The difference between them lies in the definition of the longwave optical depths.
where mixed_layer_mod
.
rad_scheme
set to 'FRIERSON'
, the longwave optical depth is specified as a function of pressurewhere
When rad_scheme
is set to 'BYRNE'
, the longwave optical depth is calculated as a function of specific humidity
In the equation above,
where
where
Schneider giant planet scheme
This scheme is suitable to use for giant planet experiments where there is no solid surface at the bottom of the model atmosphere. It is selected by setting rad_scheme
to 'SCHNEIDER'
.
surface_flux_mod
.where
and
is the Bond albedo. Here,
The only reflection of shortwave radiation occurs at the top of atmosphere where a fraction
There is no modelled upward flux of shortwave radiation in the atmosphere
where
where
Geen scheme
The Geen scheme provides an intermediate option between gray radiation and more complete descriptions of radiative transfer (e.g., the correlated-
Note: the Geen scheme is currently only stable when run in a fixed SST configuration. A version stable with a slab ocean is expected by 30/11/2020.
and in the window region:
The superscripts
As with the Frierson and Byrne schemes, the implementation of the shortwave band assumes that there is no scattering within the atmosphere, and shortwave radiation is only absorbed from the downward stream. All reflection occurs at the surface. The upward flux of shortwave
and in the non-window region
where the
where
Incoming solar radiation
For each of the radiative transfer schemes, diurnally-and-seasonally averaged, and time-dependent options for the incoming solar radiation are avaialable. The user selects whether diurnally-and-seasonally averaged or time-dependent solar forcing is used by setting the namelist option do_seasonal
.
do_seasonal = False
.For the Frierson, Byrne, and Geen schemes, setting do_seasonal = False
imposes a
where
For the Schneider giant planet scheme, setting do_seasonal = False
imposes an insolation profile that varies with
which corresponds to the latitudinal distribution of radiation received by a planet with no obliquity (perpetual equinox). As with the simple
do_seasonal = True
. In this case, astronomy_mod
is used to calculate the zenith angle and planet-star distance at each location and time. This information is then used by the two_stream_gray_rad
module to calculate the top of atmosphere insolation as a function of time.astronomy_mod
calculates the zenith angle and planet-star distance as a function of the following orbital and planetary parameters: obliquity, eccentricity, semi-major axis, longitude of perihelion (w.r.t NH autumn equinox), orbital period, and rotation rate. These options must be specified to astronomy_mod_nml
. An input parameter in two_stream_gray_rad
, equinox_day
determines the time of year when NH autumn equinox occurs.
two_stream_gray_rad
then calculates the insolation as
where rrsun
in the code.
If the user wishes, they may average incoming solar radiation over a period
There is also an option to run perpetually on one day by setting the namelist variable solday
. For example if Northern Hemisphere autumn equinox was set to occur on day solday=180
. This can be used in conjunction with an appropriate choice for
Namelist options
The namelist options for two_stream_gray_rad_nml are listed below.
Namelist option to choose scheme
- rad_scheme:
String choosing the radiation scheme. Options are
'FRIERSON'
,'BYRNE'
,'GEEN'
,'SCHNEIDER'
. Default option is'FRIERSON'
.
Namelist options for Frierson scheme longwave optical depth
- ir_tau_eq:
Surface longwave optical depth at equator. Default
.- ir_tau_pole:
Surface longwave optical depth at pole. Default
.- odp:
Frierson optical depth scaling parameter
. Default .- linear_tau:
. Determines partitioning between linear term and term in Frierson longwave optical depth. Default .- wv_exponent:
Pressure exponent,
in definition of optical depth. Default .
Namelist options for Byrne scheme longwave optical depth
- bog_a:
Absorption coefficient
in Byrne longwave optical depth. Default .- bog_b:
Absorption coefficient
in Byrne longwave optical depth. Default 1997.9.- bog_mu:
Scaling parameter
in Byrne longwave optical depth. Default 1.0.
Namelist options for Frierson/Byrne scheme shortwave optical depth
- atm_abs:
Shortwave optical depth at the equator
. Default .- sw_diff:
Amplitude of latitudinal optical depth variation
. Default .- solar_exponent:
Pressure exponent
. Default value is .
- window:
Window fraction
for longwave radiative transfer. Default value .- ir_tau_co2_win:
Absorption coefficient
in longwave optical depth. Default value .- ir_tau_wv_win1:
Absorption coefficient
in longwave optical depth. Default value .- ir_tau_wv_win2:
Absorption coefficient
in longwave optical depth. Default value .- ir_tau_co2:
Absorption coefficient
in longwave optical depth. Default value .- ir_tau_wv1:
Absorption coefficient
in longwave optical depth. Default value .- ir_tau_wv2:
Absorption coefficient
in longwave optical depth. Default value .
- diabatic_acce:
Multiplicative scaling factor for temperature tendency due to radiation in Schneider scheme. This can be used to speed up model spin-up. It will speed up the model spin-up if greater than 1. Default value is 1.
Namelist options for incoming solar radiation
- do_seasonal:
Sets whether insolation is diurnally-and-seasonally averaged (
FALSE
), or time dependent (TRUE
). DefaultFALSE
.- solar_constant:
The solar constant
. Default value: .
The following namelist options are used when rad_scheme
is set to 'FRIERSON'
, 'BYRNE'
or 'GEEN'
and do_seasonal=FALSE
:
- del_sol:
Parameter
determining insolation amplitude. Default value: .- del_sw:
Parameter
defining magnitude of modification to insolation profile. Default value: .
The following namelist options are used when do_seasonal=TRUE
:
- use_time_average_coszen:
TRUE
orFALSE
. IfTRUE
, average ( is zenith angle) over the perioddt_rad_avg
. For example, for the Earth’s diurnal period,use_time_average_coszen=TRUE
anddt_rad_avg=86400.
would achieve diurnally averaged insolation.- dt_rad_avg:
Averaging period (units: seconds) for time-dependent insolation
. Default=-1 sets averaging period to model timestep.- solday:
Day of year to run time-dependent insolation perpetually. If negative, the option to run perpetually on a specific day is not used. Default -10.
- equinox_day:
Fraction of year [0,1] where Northern Hemisphere autumn equinox occurs. Default = 0.75 (e.g. end of September for 360 day year).
Namelist options for setting carbon dioxide concentration
- do_read_co2:
TRUE
orFALSE
. IfTRUE
, reads time-varying concentration from an input file [needs to be 4D (3 spatial dimensions and time), but no spatial variation should be defined (the code only reads in maximum value at a given time)]. DefaultFALSE
.- carbon_conc:
Prescribed
(units: ppmv). Used ifdo_read_co2=FALSE
. Default value .- co2_file:
Name of
file to read (note, should be specified without .nc appendix). Default'co2'
.- co2_variable_name:
Name of
variable in file. Default'co2'
.
Important parameters not set in two_stream_gray_rad_nml
- pstd_mks:
This is used as the reference pressure
in the Frierson/Byrne/Schneider shortwave optical depth, and the Frierson/Schneider longwave optical depth. It is set inconstants_mod_nml
. Default value is . Note: this should be changed to for the giant planet configuration.- albedo_value:
The surface albedo
used by the Frierson/Byrne/Geen schemes at the lower boundary is set inmixed_layer_nml
. Default value is for a simple homogeneous slab ocean surface. The albedo can vary spatially if land or ice is introduced. For more details, see the documentation formixed_layer_mod
.- flux_heat_gp:
A prescribed heat flux through the lower boundary can be added for the giant planet case. This is set in
surface_flux_nml
. Default value is .
The following astronomical parameters are set in astronomy_mod_nml
. They are used if do_seasonal=True
.
- ecc:
Orbital eccentricity. Default
.- obliq:
Obliquity. Default
degrees.- per:
Longitude of perihelion (point in orbit when planet is closest to star) with respect to autumnal equinox in Northern Hemisphere. Default
degrees.
The following astronomical parameters are set in constant_nml
. They will be used to calculate the diurnal period if do_seasonal=True
.
- orbital_period:
Orbital period in seconds. Default is
. Only used ifcalendar
is set to'no_calendar'
inmain_nml
.- omega:
Planetary rotation rate in
. Default value is
Diagnostics
These are the diagnostics associated with the two_stream_gray_rad
module.
Name |
Description |
Units |
---|---|---|
olr |
Outgoing longwave radiation. |
|
swdn_sfc |
Absorbed shortwave at surface. |
|
swdn_toa |
Shortwave flux down at top of atmosphere. |
|
net_lw_surf |
Net upward longwave flux at surface. |
|
lwdn_sfc |
Longwave flux down at surface. |
|
lwup_sfc |
Longwave flux up at surface. |
|
tdt_rad |
Temperature tendency due to radiation. |
|
tdt_solar |
Temperature tendency due to solar radiation. |
|
flux_rad |
Total radiative flux (positive up). |
|
flux_lw |
Net longwave radiative flux (positive up). |
|
flux_sw |
Net shortwave radiative flux (positive up). |
|
coszen |
Cosine of zenith angle. |
none |
fracsun |
Daylight fraction of time interval. |
none |
co2 |
Carbon dioxide concentration. |
|
lw_dtrans |
Longwave (non-window) transmission. |
none |
lw_dtrans_win |
Longwave window transmission.
Note: only for
'GEEN' scheme. |
none |
sw_dtrans |
Shortwave transmission.
Note: only for
'GEEN' scheme. |
none |
Relevant modules and subroutines
Modules relevant to this one include:
- astronomy_mod:
Module that performs astronomical calcuations used for insolation.
- mixed_layer_mod:
Surface albedo is set here. This is also where the surface temperature is updated.
- surface_flux_mod:
An internal heat flux for giant planets can be set here.
- constants_mod:
Planetary rotation rate and orbital period are set here. These are used in the calculations made by
astronomy_mod
.
Other radiative transfer schemes are included in the following modules:
- rrtm_radiation:
Correlated-
scheme tuned for Earth-like applications.- socrates_interface_mod:
Interface for flexible Met-Office correlated-
scheme used for Earth-like and exoplanetary atmospheres.