Skip to content

Radiation

frierson_atmospheric_heating(ds, albedo=0)

Returns the atmospheric radiative heating rate from the surface and top of atmosphere energy fluxes. A negative value indicates that the atmosphere is cooling.

This takes into account any radiation that is absorbed by the atmosphere on its way down from space to the surface, as specified through tau_equator and the amount reflected at the surface through the albedo.

In Isca, there is no absorption of the shortwave radiation as it moves back up through the atmosphere to space after being reflected at the surface.

Parameters:

Name Type Description Default
ds Dataset

Dataset for particular experiment, must contain:

  • swdn_toa - Incident shortwave radiation at the top of atmosphere. This is saved by Isca if the variable swdn_toa in the two_stream module is specified in the diagnostic table.
  • swdn_sfc - Net shortwave radiation absorbed at the surface i.e. incident - reflected. This is the negative of net upward shortwave radiation at the surface. This is saved by Isca if the variable swdn_sfc in the two_stream module is specified in the diagnostic table.
  • lwup_sfc - Upward longwave flux at the surface. This is saved by Isca if the variable lwdn_sfc in the two_stream module is specified in the diagnostic table.
  • lwdn_sfc - Downward longwave flux at the surface. This is saved by Isca if the variable lwdn_sfc in the two_stream module is specified in the diagnostic table.
  • olr - Outgoing longwave radiation at the top of atmosphere. This is saved by Isca if the variable lwdn_sfc in the two_stream module is specified in the diagnostic table.
 required
albedo float

Fraction of incident shortwave radiation reflected by the surface. It is specified through the option albedo_value in the mixed_layer_nml namelist.

 0

Returns: Atmospheric radiative heating rate in \(W/m^2\).

Source code in isca_tools/utils/radiation.py 
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
def frierson_atmospheric_heating(ds: Dataset, albedo: float = 0) -> xr.DataArray:
    """
    Returns the atmospheric radiative heating rate from the surface and top of atmosphere energy fluxes. A negative
    value indicates that the atmosphere is cooling.

    This takes into account any radiation that is absorbed by the atmosphere on its way down from space to the surface,
    as specified through `tau_equator` and the amount reflected at the surface through the `albedo`.

    In *Isca*, there is no absorption of the shortwave radiation as it moves back up through the atmosphere to space
    after being reflected at the surface.

    Args:
        ds: Dataset for particular experiment, must contain:

            * `swdn_toa` - Incident shortwave radiation at the top of atmosphere.
                This is saved by *Isca* if the variable `swdn_toa` in the `two_stream` module is specified in the
                diagnostic table.
            * `swdn_sfc` - Net shortwave radiation absorbed at the surface i.e. incident - reflected.
                This is the negative of net upward shortwave radiation at the surface.
                This is saved by *Isca* if the variable `swdn_sfc` in the `two_stream` module is specified in the
                diagnostic table.
            * `lwup_sfc` - Upward longwave flux at the surface.
                This is saved by *Isca* if the variable `lwdn_sfc` in the `two_stream` module is specified in the
                diagnostic table.
            * `lwdn_sfc` - Downward longwave flux at the surface.
                This is saved by *Isca* if the variable `lwdn_sfc` in the `two_stream` module is specified in the
                diagnostic table.
            * `olr` - Outgoing longwave radiation at the top of atmosphere.
                This is saved by *Isca* if the variable `lwdn_sfc` in the `two_stream` module is specified in the
                diagnostic table.
        albedo: Fraction of incident shortwave radiation reflected by the surface.
            It is specified through the option `albedo_value` in the `mixed_layer_nml` namelist.
    Returns:
        Atmospheric radiative heating rate in $W/m^2$.

    """
    # #This is the full method of doing it, but we can do it simpler without the sw absorption stuff
    # swup_net_sfc = -ds.swdn_sfc
    # # need to account for SW absorbed by atmosphere and reflected at surface
    # swup_net_toa = -frierson_net_toa_sw_dwn(ds.swdn_toa, ds.ps, albedo, tau_equator, tau_lat_var, pressure_exponent,
    #                                         ref_pressure)
    # flux_surf = swup_net_sfc + ds.lwup_sfc - ds.lwdn_sfc
    # flux_toa = swup_net_toa + ds.olr
    # return ds.swdn_toa - ds.swdn_sfc / (1 - albedo) + ds.lwup_sfc - ds.lwdn_sfc - ds.olr

    swup_toa = albedo/(1-albedo) * ds.swdn_sfc
    flux_toa = swup_toa - ds.swdn_toa + ds.olr
    flux_surf = -ds.swdn_sfc + ds.lwup_sfc - ds.lwdn_sfc
    return flux_surf - flux_toa

frierson_net_toa_sw_dwn(insolation, surface_pressure, albedo=0, tau_equator=0, tau_lat_var=0, pressure_exponent=4, ref_pressure=101325)

Function to calculate the net downward shortwave radiation at the top of atmosphere for the Frierson and Byrne radiation schemes.

This takes into account any radiation that is absorbed by the atmosphere on its way down from space to the surface, as specified through tau_equator and the amount reflected at the surface through the albedo.

In Isca, there is no absorption of the shortwave radiation as it moves back up through the atmosphere to space after being reflected at the surface.

All default values are the default values in Isca if the option is not specified in the relavent namelist.

Parameters:

Name Type Description Default
insolation DataArray

Incident shortwave radiation at the top of atmosphere with dimensions of latitude, longitude and time. This is saved by Isca if the variable swdn_toa in the two_stream module is specified in the diagnostic table.

 required
surface_pressure DataArray

Surface pressure in Pa with dimensions of latitude, longitude and time, \(p_s\). This is saved by Isca if the variable ps in the dynamics module is specified in the diagnostic table.

 required
albedo float

Fraction of incident shortwave radiation reflected by the surface. It is specified through the option albedo_value in the mixed_layer_nml namelist.

 0
tau_equator float

Surface optical depth at the equator, \(\tau_{0e}^*\). It is specified through the option atm_abs in the two_stream_gray_rad_nml namelist.

 0
tau_lat_var float

Variation of optical depth with latitude, \(\Delta \tau^*\). It is specified through the option sw_diff in the two_stream_gray_rad_nml namelist.

 0
pressure_exponent float

Determines the variation of optical depth with pressure, \(\kappa^*\). It is specified through the option solar_exponent in the two_stream_gray_rad_nml namelist.

 4
ref_pressure float

Reference pressure used by Isca in Pa. It is specified through the option pstd_mks in the constants_nml namelist.

 101325

Returns:

Type Description
DataArray

Net downward shortwave radiation at the top of atmosphere with dimensions of latitude, longitude and time.
Source code in isca_tools/utils/radiation.py 
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
def frierson_net_toa_sw_dwn(insolation: xr.DataArray, surface_pressure: xr.DataArray, albedo: float = 0,
                            tau_equator: float = 0, tau_lat_var: float = 0, pressure_exponent: float = 4,
                            ref_pressure: float = 101325) -> xr.DataArray:
    """
    Function to calculate the net downward shortwave radiation at the top of atmosphere for the *Frierson*
    and *Byrne* radiation schemes.

    This takes into account any radiation that is absorbed by the atmosphere on its way down from space to the surface,
    as specified through `tau_equator` and the amount reflected at the surface through the `albedo`.

    In *Isca*, there is no absorption of the shortwave radiation as it moves back up through the atmosphere to space
    after being reflected at the surface.

    All default values are the default values in *Isca* if the option is not specified in the relavent `namelist`.

    Args:
        insolation: Incident shortwave radiation at the top of atmosphere
            with dimensions of latitude, longitude and time.
            This is saved by *Isca* if the variable `swdn_toa` in the `two_stream` module is specified in the
            diagnostic table.
        surface_pressure: Surface pressure in *Pa* with dimensions of latitude, longitude and time, $p_s$.
            This is saved by *Isca* if the variable `ps` in the `dynamics` module is specified in the diagnostic table.
        albedo: Fraction of incident shortwave radiation reflected by the surface.
            It is specified through the option `albedo_value` in the `mixed_layer_nml` namelist.
        tau_equator: Surface optical depth at the equator, $\\tau_{0e}^*$.
            It is specified through the option `atm_abs` in the `two_stream_gray_rad_nml` namelist.
        tau_lat_var: Variation of optical depth with latitude, $\Delta \\tau^*$.
            It is specified through the option `sw_diff` in the `two_stream_gray_rad_nml` namelist.
        pressure_exponent: Determines the variation of optical depth with pressure, $\\kappa^*$.
            It is specified through the option `solar_exponent` in the `two_stream_gray_rad_nml` namelist.
        ref_pressure: Reference pressure used by Isca in *Pa*.
            It is specified through the option `pstd_mks` in the `constants_nml` namelist.

    Returns:
        Net downward shortwave radiation at the top of atmosphere with dimensions of latitude, longitude and time.
    """
    tau = frierson_sw_optical_depth(surface_pressure, tau_equator, tau_lat_var, pressure_exponent, ref_pressure)
    return insolation*(1-albedo*np.exp(-tau))

frierson_sw_optical_depth(surface_pressure, tau_equator=0, tau_lat_var=0, pressure_exponent=4, ref_pressure=101325)

Function to calculate shortwave surface optical depth, \(\tau_s\), as a function of latitude, \(\phi\), as performed in Isca for the Frierson and Byrne radiation schemes:

\[ \tau_s(\phi) = (1-\Delta \tau^* \sin^2 \phi)\tau_{0e}^*\big(\frac{p_s}{p_0}\big)^{\kappa^*} \]

All default values are the default values in Isca if the option is not specified in the relavent namelist.

Parameters:

Name Type Description Default
surface_pressure DataArray

Surface pressure in Pa with dimensions of latitude, longitude and time, \(p_s\). This is saved by Isca if the variable ps in the dynamics module is specified in the diagnostic table.

 required
tau_equator float

Surface optical depth at the equator, \(\tau_{0e}^*\). It is specified through the option atm_abs in the two_stream_gray_rad_nml namelist.

 0
tau_lat_var float

Variation of optical depth with latitude, \(\Delta \tau^*\). It is specified through the option sw_diff in the two_stream_gray_rad_nml namelist.

 0
pressure_exponent float

Determines the variation of optical depth with pressure, \(\kappa^*\). It is specified through the option solar_exponent in the two_stream_gray_rad_nml namelist.

 4
ref_pressure float

Reference pressure used by Isca in Pa. It is specified through the option pstd_mks in the constants_nml namelist.

 101325

Returns:

Type Description
DataArray

Shortwave surface optical depth with dimensions of latitude, longitude and time.
Source code in isca_tools/utils/radiation.py 
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
def frierson_sw_optical_depth(surface_pressure: xr.DataArray, tau_equator: float = 0, tau_lat_var: float = 0,
                              pressure_exponent: float = 4, ref_pressure: float = 101325) -> xr.DataArray:
    """
    Function to calculate shortwave surface optical depth, $\\tau_s$, as a function of latitude, $\\phi$, as performed in
    [Isca](https://execlim.github.io/Isca/modules/two_stream_gray_rad.html#frierson-byrne-schemes) for the *Frierson*
    and *Byrne* radiation schemes:

    $$
    \\tau_s(\\phi) = (1-\Delta \\tau^* \sin^2 \\phi)\\tau_{0e}^*\\big(\\frac{p_s}{p_0}\\big)^{\\kappa^*}
    $$

    All default values are the default values in *Isca* if the option is not specified in the relavent `namelist`.

    Args:
        surface_pressure: Surface pressure in *Pa* with dimensions of latitude, longitude and time, $p_s$.
            This is saved by *Isca* if the variable `ps` in the `dynamics` module is specified in the diagnostic table.
        tau_equator: Surface optical depth at the equator, $\\tau_{0e}^*$.
            It is specified through the option `atm_abs` in the `two_stream_gray_rad_nml` namelist.
        tau_lat_var: Variation of optical depth with latitude, $\Delta \\tau^*$.
            It is specified through the option `sw_diff` in the `two_stream_gray_rad_nml` namelist.
        pressure_exponent: Determines the variation of optical depth with pressure, $\\kappa^*$.
            It is specified through the option `solar_exponent` in the `two_stream_gray_rad_nml` namelist.
        ref_pressure: Reference pressure used by Isca in *Pa*.
            It is specified through the option `pstd_mks` in the `constants_nml` namelist.

    Returns:
        Shortwave surface optical depth with dimensions of latitude, longitude and time.
    """
    tau_surface = (1-tau_lat_var * np.sin(np.deg2rad(surface_pressure.lat))**2) * tau_equator
    return tau_surface * (surface_pressure/ref_pressure)**pressure_exponent

get_heat_capacity(c_p, density, layer_depth)

Given heat capacity un units of \(JK^{-1}kg^{-1}\), this returns heat capacity in units of \(JK^{-1}m^{-2}\).

Parameters:

Name Type Description Default
c_p float

Specific heat at constant pressure.
Units: \(JK^{-1}kg^{-1}\)

 required
density float

Density of substance (usually air or water).
Units: \(kgm^{-3}\)

 required
layer_depth float

Depth of layer.
Units: \(m\)

 required

Returns:

Type Description
float

Heat capacity in units of \(JK^{-1}m^{-2}\).
Source code in isca_tools/utils/radiation.py 
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
def get_heat_capacity(c_p: float, density: float, layer_depth: float) -> float:
    """
    Given heat capacity un units of $JK^{-1}kg^{-1}$, this returns heat capacity in units of $JK^{-1}m^{-2}$.

    Args:
        c_p: Specific heat at constant pressure.</br>
            Units: $JK^{-1}kg^{-1}$
        density: Density of substance (usually air or water).</br>
            Units: $kgm^{-3}$
        layer_depth: Depth of layer.</br>
            Units: $m$

    Returns:
        Heat capacity in units of $JK^{-1}m^{-2}$.
    """
    return c_p * density * layer_depth

opd_lw_gray(lat, pressure=None, kappa=1, tau_eq=6, tau_pole=1.5, pressure_ref=10 ** 5, frac_linear=0.1, k_exponent=4)

Returns the longwave optical depth used in the Frierson Isca rad_scheme.

If pressure not provided, will return surface value. Otherwise, will return value at give pressure.

Parameters:

Name Type Description Default
lat ndarray

float [n_lat]
Latitude in degrees.

 required
pressure Optional[float]

Pressure in Pa.

 None
kappa float

Frierson optical depth scaling parameter.
opd in two_stream_gray_rad_nml namelist.

 1
tau_eq float

Surface longwave optical depth at equator.
ir_tau_eq in two_stream_gray_rad_nml namelist.

 6
tau_pole float

Surface longwave optical depth at pole.
ir_tau_pole in two_stream_gray_rad_nml namelist.

 1.5
pressure_ref float

Reference pressure in Pa.

 10 ** 5
frac_linear float

Determines partitioning between linear term and \(p^k\) term.
linear_tau in two_stream_gray_rad_nml namelist.

 0.1
k_exponent float

Pressure exponent.
wv_exponent in two_stream_gray_rad_nml namelist.

 4

Returns:

Name Type Description
opd ndarray

float [n_lat]
Longwave optical depth
Source code in isca_tools/utils/radiation.py 
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
def opd_lw_gray(lat: np.ndarray, pressure: Optional[float] = None,
                kappa: float = 1, tau_eq: float = 6, tau_pole: float = 1.5,
                pressure_ref: float = 10**5, frac_linear: float=0.1, k_exponent: float=4) -> np.ndarray:
    """
    Returns the longwave optical depth used in the
    [Frierson](https://execlim.github.io/Isca/modules/two_stream_gray_rad.html#frierson-byrne-schemes)
    Isca `rad_scheme`.

    If `pressure` not provided, will return surface value. Otherwise, will return value at give `pressure`.

    Args:
        lat: `float [n_lat]`</br>
            Latitude in degrees.
        pressure: Pressure in Pa.
        kappa: Frierson optical depth scaling parameter.</br>
            `opd` in `two_stream_gray_rad_nml` namelist.
        tau_eq: Surface longwave optical depth at equator.</br>
            `ir_tau_eq` in `two_stream_gray_rad_nml` namelist.
        tau_pole: Surface longwave optical depth at pole.</br>
            `ir_tau_pole` in `two_stream_gray_rad_nml` namelist.
        pressure_ref: Reference pressure in Pa.
        frac_linear: Determines partitioning between linear term and $p^k$ term.</br>
            `linear_tau` in `two_stream_gray_rad_nml` namelist.
        k_exponent: Pressure exponent.</br>
            `wv_exponent` in `two_stream_gray_rad_nml` namelist.

    Returns:
        opd: `float [n_lat]`</br>
            Longwave optical depth
    """
    opd_surf = kappa * (tau_eq + (tau_pole - tau_eq) * np.sin(np.deg2rad(lat)) ** 2)
    if pressure is None:
        return opd_surf
    else:
        pressure_factor = frac_linear * (pressure/pressure_ref) + (1-frac_linear) * (pressure/pressure_ref)**k_exponent
        return opd_surf * pressure_factor