## b. Radiation Scheme

Here, the radiation scheme used in the numerical calculations and the characteristics of the scheme will be described. The key is to discretize the radiative flux equation in two ways and then join the two discretized equations, which allows radiation calculations for a gray atmosphere with an arbitrary optical depth profile.

## Finite-Difference Representation of Radiative Flux

Net upward radiative flux for the case of a gray atmosphere is given by the following equation:

where is optical depth at the top of the atmosphere, thus is .

πBis equal toσT.^{4}Here, the radiative flux equation is discretized in the same way as in Nakajima

et al. (1992), and thusHere,

Krepresents the upper-most layer, and_{MAX}K= 32 in a typical model calculation. The radiative flux is calculated by joining two sets of discretized equations. Specifically, the integrals for the contributions to net upward radiative flux from optically thin layers and optically thick layers are evaluated using the upper and lower options in the equation, respectively. The fraction of the usage of the upper and lower options in the equation changes according to the value of . If a significantly large value is taken for , the integrals are evaluated entirely with the upper option, which results in a radiative scheme identical to that used in Numaguti and Hayashi (1991)._{MAX}There exists no that makes the upper and lower options identical. If such hypothetically exists, the relationship

needs to hold. Rearranging this condition by letting leads to

Unfortunately, no solutions exist for this equation. More accurately, a discontinuity emerges in the evaluated fluxes at the altitude at which the transition between the two options occurs. However, the actual difference between the evaluated fluxes is almost negligibly small. In fact, no evident visible discontinuities exist in the vertical profiles of the radiative flux. (For more information, see Figure 5 of Characteristics of One-Dimensional Radiative-Convective Equilibrium Solutions.)

The above-mentioned complex procedure is necessary because problems arise if either the upper or lower option in the discretized radiative flux equation alone is used. The reasons that the problem arises will be described below, and the values of that are suitable for 3D calculations will be examined.

## For the Case of Large

C_{Δτ}

First, the extreme limits are considered in which is set to a significantly large value. This corresponds to using the upper option of the discretized equation for the entire layer, that is, using the radiative scheme of Numaguti and Hayashi (1991). As an example, results will be shown for = 10

^{6},P= 10_{n0}^{5}, and a vertical grid resolution of L1000. The relationship betweenTand OLR from this case is almost identical to that from = 10_{g}^{-1}(not shown). However, the vertical profile of radiative flux at high temperatures is not calculated properly. Figure 1 shows the vertical profile of net upward radiative flux for the case ofT= 550 K. Unlike in the case of = 10_{g}^{-1}, net upward radiative flux increases toward the surface in the layer below the vicinity of log_{10}σ= - 2. Net upward radiative flux decreases from the top to a certain level in the middle layer, and below the level flux increases again. This tendency is attributable to the following reason. If the atmosphere becomes optically thick, only the terms associated withj=kremain in the lower option of the discretized equation, which then leads tofor the evaluation of net upward radiative flux. Properly,

Fshould be represented as^{net}

(a)

(b)

Figure 1:The results for the case of = 10^{6}. For L1000,p= 10_{n0}^{5}, andT= 550 K. (a) The vertical profile of net upward radiative flux (W/m_{s}^{2}). (b) The vertical profile of radiative heating rate (K/sec).

In this example calculation, the values of

σTand net upward radiative flux are^{4}

For the case of L32, the above-described phenomenon is more amplified. Net upward radiative flux starts increasing near log

_{10}σ= - 3 from higher to lower altitudes, however, it starts decreasing again around log_{10}σ= -1 (Figure 2(a)). This profile results because the net upward radiative flux is determined byσT-^{4}_{k-1}σTbelow log^{4}_{k}_{10}σ= -3.0. The decrease in net upward radiative flux from higher to lower altitudes near the surface is attributable to small vertical grid intervals near the surface. The net upward radiative flux values around log_{10}σ= - 1 for the case of = 10^{6}are larger than those for the case of = 10^{-1}by 400 W/m^{2}. The differences in both upward and downward radiative fluxes in the middle layer between the case of = 10^{6}and the case of = 10^{-1}are approximately 200 W/m^{2}. On the other hand, the error in the radiative heating profile which is important for time evolution problems in 3D model is relatively small. (Figure 2(b)).(a)

(b)

Figure 2:The results for the case of = 10^{6}. For L32,p= 10_{n0}^{5}Pa, andT= 550 K. (a) The vertical profile of net upward radiative flux (W/m_{s}^{2}). (b) The vertical profile of radiative heating rate (K/sec).In the case of L32, a discontinuity arises in the middle layer of the vertical profile of radiative flux for ≥ 10.0. Though not shown, in the case of L32, the radiative flux can be calculated with few problems as long as ≤ 1.0.

## For the Case of Small

C_{Δτ}

The behaviors of the results from the extreme limit of a significantly small will be examined. The results under consideration correspond to those that are obtained by evaluating radiative flux integrals for the entire layer with the use of the lower option of the discretized equation. As an example, results will be shown for = 0.0,

P= 10_{n0}^{5}, and the vertical grid resolution of L32. When 1D equilibrium solutions for this case are obtained according to the procedure for the 1D convective-radiative equilibrium calculation, the final values of net upward radiative flux appear to be accurate (not shown). However, in the case ofT= 250 K, the radiation calculation for states in which an adiabatic structure is adopted from the bottom to the top of the atmosphere cannot be accurately performed. Discontinuities arise in the upward and downward radiative flux profiles as shown in Figures 3(a) and 3(b). Incidentally, such discontinuities are clearly present up to approximately_{g}T= 300 K. The reason is that the value of_{g}becomes zero in the optically thin layer, that is, in the upper, low-temperature layer. When this option is used, the values of

τand_{j-1/2}τbecome significantly small. For this reason, the calculation yields a numerically equal value for and within the range of numerical error, resulting in zero for the difference between the two terms. Properly, the values of both and become quite small, leading to_{j+1/2}In the layer above the altitude at which becomes zero, the downward radiative flux is determined by

is thus smaller than the correct value of downward radiative flux. Similarly, for the upward radiative flux, some of the additive terms drop away above the altitude at which becomes zero.

(a)

(b)

Figure 3:Radiative flux profiles obtained with = 0.0 for a state in which an adiabatic structure is adopted from the bottom to the top of the atmosphere. (a) Vertical profile of upward radiative flux (W/m^{2}). (b) Vertical profile of downward radiative flux (W/m^{2}). For the case of L32,p=10_{n0}^{5}Pa, andT= 250 K._{s}Furthermore, the above-mentioned problem occurs in a more extreme manner with increasing vertical grid resolution. In the case of L1000, the value of OLR can no longer be calculated correctly, and thus an accurate relationship between

Tand OLR cannot be obtained (not shown)._{g}Consequently, the upper option of the discretized equation cannot be applied for the entire atmospheric layer. However, it has been confirmed that radiative flux can be calculated without difficulties for ≥ 10

^{-7}(not shown).## Summary

If the upper option of the discretized radiative flux equation alone is used, a problem arises in the radiation calculations for an optically thin atmosphere. Furthermore, the lower option of the discretized radiative flux equation alone cannot be used for an optically thick atmosphere. However, a combination of the two options of the discretized equation allows radiation calculations of a gray atmosphere with an arbitrary profile of optical depth. In this regard, the value of

C_{Δτ}that prescribes the switchover between the discretized equations needs to take a value between 1.0 and 10^{-7}. Based on the results above, a value of 0.1 is adopted forC_{Δτ}in the present study.