## c. Examination on the Vertical Resolution

In this section, we examine the vertical grid resolution that is necessary for radiation calculations. Through comparisons of 1D calculation results obtained from the use of various vertical grid resolutions, the minimum vertical grid resolution required for radiation calculations will be identified.

## Positions of the Vertical Layers

Nakajima et al. (1992) conducts calculations using 400 to 700 vertical layers in order to acquire 1D radiative equilibrium solutions. However, it is practically infeasible to conduct 3D time evolutions with that number of vertical grid resolutions. In order to conduct 3D calculations, it is necessary to identify the minimum number of vertical layers that allows radiation calculations with sufficient accuracy. To accomplish this, we will compare calculations from four scenarios: 16, 21, 32, and 1000 vertical layer resolutions. The models based on these four resolutions will be referred to as L16, L21, L32, and L1000, respectively. Table 1 shows the positions of the vertical layers for the case of L32.In L16, the 15 lowest layers listed in Table 1 are used , and the 16th layer is set to 0. In L21,

σ= 10^{-1}, 10^{-2}, 10^{-3}, 10^{-4}, and 10^{-5}are added to the vertical layers of L16. In L1000, the entire atmospheric layer (σ= 1 to 10^{-6}) is divided into equal intervals of lnp.Incidentally, the entire atmospheric layer is divided into equal intervals of ln p in Nakajima

et al. (1992). In their study, the number of the vertical layers is set to

The actual numbers of vertical layers are 560 for

c= 4_{pv}R,c= 3.5_{pn}R,p= 10_{n0}^{5}Pa,T= 500 K and 702 for_{g}c= 4_{pv}R,c= 4.5_{pn}R,p= 10_{n0}^{8}Pa,T= 500 K._{g}

Table 1:Positions of the vertical grids.k: the vertical grid number, numbered in order 0, 1, 2, ... from the surface up.σ:the value of theσ-coordinate.z: the global mean values of the geopotential heights (m) obtained from the case ofS= 1380 W/m^{2}.

Tg-OLR RelationshipHere, it will be determined if an accurate relationship between

Tand OLR can be obtained with various vertical grid resolutions._{g}Figure 1 illustrates the relationships between surface temperature and OLR that were acquired for various vertical resolutions for the case of

p= 10_{n0}^{5}Pa.

The dependency of OLR on Tg in the case of L32 is identical to that in the case of L1000. However, when the resolution is reduced to L21, the

T-OLR curve becomes undulating. Since this case does not allow an accurate evaluation of the tropopause height, radiative flux cannot be accurately calculated. When the vertical resolution is further reduced to that of L16, no range exists in which the_{g}T-OLR curve is characterized by a flat line. Such a result is obtained because the top of the troposphere is situated at an altitude higher than the top of the model domain at values of_{g}Texceeding 350 K. Although the highest layer corresponds to_{g}σ= 0.05 in the case of L16, the altitude of the tropopause becomes even higher than this layer whenTexceeds 350 K._{g}

Figure 1:The relationship between surface temperature (K) and OLR (W/m2) forC= 0.1. For the case of_{Δ τ}p= 10_{n0}_{5}Pa. Black: the results from L1000. Red: the results from L32. Green: the results from L21. Blue: the results from L16.For the purpose of comparing these values between the cases of L32 and L1000, Tables 2 and 3 show the values of the OLR from the flat portion and near the peak of the

T-OLR curve, respectively. The difference in the pairs of values between L32 and L1000 is as large as approximately 5 W/m_{g}^{2}. Thus, it can be claimed that radiation calculations can be conducted with L32 with an accuracy of approximately 5 W/m^{2}.

Table 2:The values of OLR from the flat portion of the surface temperature-OLR curve . The values in the rows labeled L1000 and L32 are the values of OLR (W/m^{2}) obtained from setting the number of vertical layers to 1000 and 32, respectively.

Table 3:The values of OLR from near the peak of the surface temperature-OLR curve. The values in the rows labeled L1000 and L32 are the values of OLR (W/m_{2}) obtained from setting the number of vertical layers to 1000 and 32, respectively.

The discussions above reveals that L32 allows calculations of OLR with sufficient accuracy for the case of p

_{n0}= 10^{5}Pa. The subsequent examination addresses whether 1D equilibrium solutions can be accurately obtained when partial pressure of the non-condensable constituents at the surface is varied. Figure 2 shows theT-OLR relationship evaluated with L32 for various partial pressures of the non-condensable constituents at the surface, p_{g}_{n0}= 10^{4}, 10^{5}, 10^{6},10^{7}Pa. The figure suggests that the relationship obtained from L32 is quite similar to that obtained from L1000 for surface partial pressures up to pn0 = 105 Pa. However, small undulations become evident in the curve in the case ofp= 10_{n0}^{6}Pa. In the case ofp= 10_{n0}^{7}Pa, significant undulations are present on the peak portion of the curve. For all values ofp, the values of OLR are larger in L32 than in L1000. In particular, for the case of_{n0}p= 10_{n0}^{7}Pa, the value of OLR even exceeds the Komabayashi-Ingersoll limit (385 W/m^{2}). It can be presumed that the value of OLR can be accurately calculated by L32 up to approximatelyp= 10_{n0}^{6}Pa.

Figure 2:The relationship between surface temperature (K) and OLR (W/m^{2}) for cases with various amounts of the non-condensable constituents are given. Blue:p= 10_{n0}^{4}Pa. Aqua:p= 10_{n0}^{5}Pa. Green:p= 10_{n0}^{6}Pa. Red:p= 10_{n0}^{7}Pa. The vertical grid resolution is that of L32.## Vertical Structures

It is confirmed that the value of OLR can be calculated accurately for

p< 10_{n0}^{6}Pa by L32. Here, we examine whether the vertical profiles of physical variables from L32 are also similar to those from L1000. Figures 3 and 4 show the vertical profiles of temperature and net upward radiative flux, respectively. The vertical profiles in both figures are similar to those from L1000, suggesting that L32 is also able to calculate the vertical structures with sufficient accuracy. In conclusion, we find that the use of L32 is sufficient for the purpose of calculating the atmospheric conditions with the present scheme forp= 10_{n0}^{5}Pa andT≤ 560 K._{g}

(a) (b)

Figure 3:Vertical profiles of temperature (K) of 1D radiative solutions. (a) Results from L32. (b) Results from L1000. ForpPa. The displayed results are for_{n0}=10^{5}T= 250, 300, 350, 400, 450, 500, and 550 K._{g}

(a) (b)

Figure 4:Vertical profiles of net upward radiative flux (W/m^{2}) based on 1D radiative solutions. (a) Results from L32. (b) Results from L1000. ForC= 0.1 and_{Δ τ}p= 10_{n0}^{5}Pa. The displayed results are forT= 250, 300, 350, 400, 450, 500, and 550 K._{g}## Summary

For radiation calculations, the use of 32 vertical layers is sufficient if

p< 10_{n0}^{6}Pa. However, this result represents merely the minimum number of vertical layers that is required in radiation calculation to obtain 1D equilibrium solutions. Additional investigations will be necessary to determine if the above result is applicable for atmospheric dynamics calculations. This topic will be addressed in future research.