The incoming solar radiation at the top of the atmosphere is computed from the formulas given in Smithsonian Meteorological Tables (List 1958), with an attenuation by clouds as in Reed (1977). The formula adopted here is the same as described in Rosati and Miyakoda (1988):
The parameter , the incident radiation under clear skies, is composed of direct and diffuse components of short wave radiation:
The direct component of short wave radiation is
where the atmospheric transmission coefficient , is the zenith angle of the sun, and is the incident radiation at the top of the atmosphere. The diffuse component of net short wave radiation is
where is the absorption due to water vapor and ozone. The solar noon altitude in degrees (), , and are computed from present day astronomy (see Appendix B).
The ocean surface albedo is taken from Payne (1972). In order to compute monthly means of , results from the formula are averaged for each day of the month, using analyzed monthly means of cloudiness for .
The net long wave radiation is based on an empirical formula which includes the effect of air-sea temperature difference (e.g., Rosati and Miyakoda 1988)
The dimensionless cloud coefficient is allowed to vary with latitude in order to account for different cloud types (Budyko 1974, Oberhuber 1988). This bulk parameterization is based on a formula originally given by Brunt (1932), later modified to include the term proportional to which attempts to include the effects of the temperature difference between the surface layer and the lower layer of air (Budyko 1974). Although several authors have adopted the proportionality term as above (e.g., Simpson and Paulson 1979, Rosati and Miyakoda 1988), others (Bunker 1976, Hsiung 1986, Oberhuber 1988, Isemer et al. 1989) have used instead.