Technical note: The libRadtran software package for radiative transfer calculations - description and examples of use

The libRadtran software package is a suite of tools for radiative transfer calculations in the Earth's atmosphere. Its main tool is the uvspec program. It may be used to compute radiances, irradiances and actinic fluxes in the solar and terrestrial part of the spectrum. The design of uvspec allows simple problems to be easily solved using defaults and included data, hence making it suitable for educational purposes. At the same time the flexibility in how and what input may be specified makes it a powerful and versatile tool for research tasks. The uvspec tool and additional tools included with libRadtran are described and realistic examples of their use are given. The libRadtran software package is available from http://www.libradtran.org

Recursive formulas for the motivic Milnor basis

We give recursive formulas for the generating elements in the Milnor basis of the mod 2 motivic Steenrod algebra.Comment: 10 page

Radiation transport in cloudy and aerosol loaded atmospheres

Thesis (Ph.D.) University of Alaska Fairbanks, 1992The equation for radiation transport in vertical inhomogeneous absorbing, scattering, and emitting atmospheres is derived from first principles. It is cast in a form amenable to solution, and solved using the discrete ordinate method. Based on the discrete ordinate solution a new computationally efficient and stable two-stream algorithm which accounts for spherical geometry is developed. The absorption and scattering properties of atmospheric molecules and particulate matter is discussed. The absorption cross sections of the principal absorbers in the atmosphere, H\sb2O, CO\sb2 and O\sb3, vary erratically and rapidly with wavelength. To account for this variation, the correlated-k distribution method is employed to simplify the integration over wavelength necessary for calculation of warming/cooling rates. The radiation model, utilizing appropriate absorption and scattering cross sections, is compared with ultraviolet radiation measurements. The comparison suggests that further experiments are required. Ultraviolet (UV) and photosynthetically active radiation (PAR) is computed for high and low latitudes for clear and cloudy skies under different ozone concentrations. An ozone depletion increases UV-B radiation detrimental to life. Water clouds diminish UV-B, UV-A and PAR for low surface albedos and increase them for high albedos. The relative amount of harmful UV-B increases on overcast days. The daily radiation doses exhibit small monthly variations at low latitudes but vary by a factor of 3 at high latitudes. Photodissociation and warming/cooling rates are calculated for clear skies, aerosol loaded atmospheres, and atmospheres with cirrus and water clouds. After major volcanic explosions aerosols change O\sb3 and NO\sb2 photodissociation rates by 20%. Both aged aerosols and cirrus clouds have little effect on photodissociation rates. Water clouds increase $(\sim$100%) photodissociation rates that are sensitive to visible radiation above the cloud. Solar warming rates vary by 50% in the stratosphere due to changing surface albedo. Water clouds have a similar effect. The net effect of cirrus clouds is to warm the troposphere and the stratosphere. Only extreme volcanic aerosol loadings affect the terrestrial warming rate, causing warming below the aerosol layer and cooling above it. Aerosols give increased solar warming above the aerosol layer and cooling below it

Effects of atmospheric sphericity on stratospheric chemistry and dynamics over Antarctica

Atmospheric sphericity is an important factor that must be considered in order to evaluate an accurate ozone loss rate in the polar stratosphere. The built-in plane-parallel radiative transfer scheme of a nudging chemical transport model (CTM) and an atmospheric general circulation model (AGCM) with coupled chemistry is modified by a pseudospherical approximation. The plane-parallel atmosphere radiative transfer version (PPA version) is compared with the pseudospherical atmosphere radiative transfer version (SA version) for both the nudging CTM and AGCM. The nudging CTM can isolate the chemical effects for a given dynamical field, while the interaction among the chemical, radiative, and dynamical processes can be studied with the AGCM. The present analysis focuses on Antarctica during an ozone hole period. In the ozone loss period over Antarctica, ozone starts to decrease earlier and minimum value of total ozone becomes lower in the SA versions of both the nudging CTM and the AGCM than in the corresponding PPA versions. The ozone mixing ratio decreases earlier in the SA version because of an earlier increase of ClO concentration initiated by the upward actinic flux at solar zenith angles greater than 90°. Dynamics plays an important role as well as the chemical processes. During the ozone recovery period, the ozone distribution becomes almost the same in the SA and PPA versions of the nudging CTM, while in the AGCM the ozone amount in the SA version remains at lower values compared to those of the PPA version. In the AGCM, a decrease of ozone over Antarctica enhances the latitudinal gradient of temperature and thus strengthens the polar vortex in the SA version. A resultant delay of the polar vortex breakup causes the delay of the ozone recovery. For the AGCM, ensemble runs are performed. The ensemble experiment exhibits large ozone variances after the middle of December, when the ozone recovery is dynamically controlled. Most ensemble members of the AGCM show a delay of the polar vortex breakup in the SA version, while a few members show opposite results. In the latter members, the polar vortex breakup is strongly affected by the enhanced EP flux from the troposphere around 100 hPa, which causes the variances in the ozone recovery period. Most members, however, do not show large statistical variances; that justifies the conclusions from the ensemble means

Optimizing separator pressure in a multistage crude oil production plant

This thesis contains a steady state simulation model of a multistage crude oil production plant. The unstable reservoir fluid is separated into crude oil and gas in a three stage gas oil separation plant. Peng Robinson Equation of State and successive substitution is used to model the separators. The liquid output for the third stage separators is stable crude oil and the crude oil is sold on the international market. The gas liberated in the separators is compressed by centrifugal compressors and injected into the reservoir. Several methods for modelling centrifugal compressors are discussed. The discussion concludes that Dimensionless Parameters is the best suited method for this case and this method is therefore used. The plant studied has three separation stages at decreasing pressure and the standard control method is to use fixed separator pressure. This is not optimal and brute force optimization is used to find the optimal pressure set-points. Adaptive pressure set-points will improve product separation and therefore increase the profit. The profit is calculated as the income from crude oil sale minus the energy cost. Using this new method, the optimum separator pressures for Snorre A are determined. As a result, the profit is increased with 0.17, 0.11 and 0.07 present for three different inlet conditions. For one year this small increase in profit gives an extra two million dollars in income. The adaptive pressure set-points give an extension of the plant s maximum capacity for a reservoir fluid with a high gas oil ratio