Search

Scattering by infinite hexagonal ice prisms is calculated using Maxwell's equations in the discrete dipole approximation for size parameters up to x=400. Birefringence is included in the calculations. Applicability of the geometric optics approximation is investigated. Excellent agreement between wave optics and geometric optics is observed for large size parameter in the outer part of the 22 degree halo feature. For smaller ice crystals halo broadening is predicted, and there is appreciable "spillover" of the halo into shadow scattering angles <22 degrees. Ways to retrieve ice crystal sizes are suggested based on the full width at half-maximum of the halo, the power at <22 deg, and the halo polarization.

DDSCAT 7.3 is an open-source Fortran-90 software package applying the discrete dipole approximation to calculate scattering and absorption of electromagnetic waves by targets with arbitrary geometries and complex refractive index. The targets may be isolated entities (e.g., dust particles), but may also be 1-d or 2-d periodic arrays of "target unit cells", allowing calculation of absorption, scattering, and electric fields around arrays of nanostructures. The theory of the DDA and its implementation in DDSCAT is presented in Draine (1988) and Draine & Flatau (1994), and its extension to periodic structures in Draine & Flatau (2008), and efficient near-field calculations in Flatau & Draine (2012). DDSCAT 7.3 includes support for MPI, OpenMP, and the Intel Math Kernel Library (MKL). DDSCAT supports calculations for a variety of target geometries. Target materials may be both inhomogeneous and anisotropic. It is straightforward for the user to "import" arbitrary target geometries into the code. DDSCAT automatically calculates total cross sections for absorption and scattering and selected elements of the Mueller scattering intensity matrix for user-specified scattering directions. DDSCAT 7.3 can efficiently calculate E and B throughout a user-specified volume containing the target. This User Guide explains how to use DDSCAT 7.3 to carry out electromagnetic scattering calculations, including use of DDPOSTPROCESS, a Fortran-90 code to perform calculations with E and B at user-selected locations near the target. A number of changes have been made since the last release, DDSCAT 7.2 ., Comment: 102 pages, 17 figures. arXiv admin note: substantial text overlap with arXiv:1202.3424, arXiv:1002.1505, arXiv:0809.0337, arXiv:astro-ph/0409262, arXiv:astro-ph/0309069, arXiv:astro-ph/0008151

On March 11, 2013 I talked with Warren Wiscombe about his contributions to scientific computer programming, atmospheric science and radiative transfer. Our conversation is divided into three parts related to light scattering, radiative transfer and his general thoughts about scientific programming. There are some reflections on how radiative transfer parameterizations gradually sneaked in to modern climate and atmospheric Global Circulation Models. Why some software programs such as light scattering code MIEV and DISORT are very successful and why some of them fizzle. We talked about the role of tools in modern science, open source movement, repeatability of scientific results, computer languages, computer programs as objects of arts, and even if programs can be revolutionary.

This manual describes the library of conjugate gradients codes CCGPAK, which solves system of complex linear system of equations. The library is written in FORTRAN90 and is highly portable. The codes are general and provide mechanism for matrix times vector multiplication which is separated from the conjugate gradient iterations itself. It is simple to switch between single and double precisions. All codes follow the same naming conventions.

DDSCAT 7.2 is a freely available open-source Fortran-90 software package applying the discrete dipole approximation (DDA) to calculate scattering and absorption of electromagnetic waves by targets with arbitrary geometries and complex refractive index. The targets may be isolated entities (e.g., dust particles), but may also be 1-d or 2-d periodic arrays of "target unit cells", which can be used to study absorption, scattering, and electric fields around arrays of nanostructures. The DDA approximates the target by an array of polarizable points. The theory of the DDA and its implementation in DDSCAT is presented in Draine (1988) and Draine & Flatau (1994), and its extension to periodic structures in Draine & Flatau (2008). Efficient near-field calculations are enabled following Flatau & Draine (2012). DDSCAT 7.2 allows accurate calculations of electromagnetic scattering from targets with size parameters 2*pi*aeff/lambda < 25 provided the refractive index m is not large compared to unity (|m-1| < 2). DDSCAT 7.2 includes support for MPI, OpenMP, and the Intel Math Kernel Library (MKL). DDSCAT 7.2 supports calculations for a variety of target geometries (e.g., ellipsoids, regular tetrahedra, rectangular solids, finite cylinders, hexagonal prisms, etc.). Target materials may be both inhomogeneous and anisotropic. It is straightforward for the user to import new target geometries into the code. DDSCAT 7.2 calculates total cross sections for absorption and scattering and selected elements of the Mueller scattering intensity matrix for specified orientation of the target relative to the incident wave, and for specified scattering directions. DDSCAT 7.2 calculates E throughout a user-specified rectangular volume containing the target. A Fortran-90 code READNF to read E and P from files created by DDSCAT 7.2 is included in the distribution., Comment: 95 pages. Expanded and updated to conform to latest release, DDSCAT 7.2.1. Complete software package available at http://code.google.com/p/ddscat/

DDSCAT 7.1 is an open-source Fortran-90 software package applying the discrete dipole approximation to calculate scattering and absorption of electromagnetic waves by targets with arbitrary geometries and complex refractive index. The targets may be isolated entities (e.g., dust particles), but may also be 1-d or 2-d periodic arrays of "target unit cells", allowing calculation of absorption, scattering, and electric fields around arrays of nanostructures. The theory of the DDA and its implementation in DDSCAT is presented in Draine (1988) and Draine & Flatau (1994), and its extension to periodic structures (and near-field calculations) in Draine & Flatau (2008). DDSCAT 7.1 includes support for MPI, OpenMP, and the Intel Math Kernel Library (MKL). DDSCAT supports calculations for a variety of target geometries. Target materials may be both inhomogeneous and anisotropic. It is straightforward for the user to "import" arbitrary target geometries into the code. DDSCAT automatically calculates total cross sections for absorption and scattering and selected elements of the Mueller scattering intensity matrix. This User Guide explains how to use DDSCAT 7.1 to carry out electromagnetic scattering calculations. DDfield, a Fortran-90 code to calculate E and B at user-selected locations near the target, is included in the distribution. A number of changes have been made since the last release, DDSCAT 7.0 ., Comment: 83 pages, 11 figures

DDSCAT 7.0 is an open-source Fortran-90 software package applying the discrete dipole approximation to calculate scattering and absorption of electromagnetic waves by targets with arbitrary geometries and complex refractive index. The targets may be isolated entities (e.g., dust particles), but may also be 1-d or 2-d periodic arrays of "target unit cells", allowing calculation of absorption, scattering, and electric fields around arrays of nanostructures. The theory of the DDA and its implementation in DDSCAT is presented in Draine (1988) and Draine & Flatau (1994), and its extension to periodic structures (and near-field calculations) in Draine & Flatau (2008). DDSCAT 7.0 includes support for MPI, OpenMP, and the Intel Math Kernel Library (MKL). DDSCAT supports calculations for a variety of target geometries. Target materials may be both inhomogeneous and anisotropic. It is straightforward for the user to "import" arbitrary target geometries into the code. DDSCAT automatically calculates total cross sections for absorption and scattering and selected elements of the Mueller scattering intensity matrix. This User Guide explains how to use DDSCAT 7.0 to carry out electromagnetic scattering calculations. DDfield, a Fortran-90 code DDfield to calculate E and B at user-selected locations near the target, is included in the distribution., Comment: 78 pages, 10 figures. v4 has some typos corrected, and includes more information for installing DDSCAT on non-Linux platforms. v5 corrects some typos in section 8.8, and updates the citation to Draine & Flatau (2008)

The discrete-dipole approximation (DDA) is a powerful method for calculating absorption and scattering by targets that have sizes smaller than or comparable to the wavelength of the incident radiation. The DDA can be extended to targets that are singly- or doubly-periodic. We generalize the scattering amplitude matrix and the 4 x 4 Mueller matrix to describe scattering by singly- and doubly-periodic targets, and show how these matrices can be calculated using the DDA. The accuracy of DDA calculations using the open-source code DDSCAT is demonstrated by comparison to exact results for infinite cylinders and infinite slabs. A method for using the DDA solution to obtain fields within and near the target is presented, with results shown for infinite slabs., Comment: 19 pages, 7 figures, submitted to J. Opt. Soc. Am. A

DDSCAT 6.1 is a software package which applies the discrete dipole approximation (DDA) to calculate scattering and absorption of electromagnetic waves by targets with arbitrary geometries and complex refractive index. DDSCAT 6.1 allows accurate calculations of electromagnetic scattering from targets with size parameters 2 pi a_eff/lambda < 15 provided the refractive index m is not large compared to unity (|m-1| < 2). DDSCAT 6.1 includes support for MPI and FFTW. We also make available a "plain" distribution of DDSCAT 6.1 that does not include support for MPI, FFTW, or netCDF, but is much simpler to install than the full distribution. The DDSCAT package is written in Fortran and is highly portable. The program supports calculations for a variety of target geometries (e.g., ellipsoids, regular tetrahedra, rectangular solids, finite cylinders, hexagonal prisms, etc.). Target materials may be both inhomogeneous and anisotropic. It is straightforward for the user to import arbitrary target geometries into the code, and relatively straightforward to add new target generation capability to the package. DDSCAT automatically calculates total cross sections for absorption and scattering and selected elements of the Mueller scattering intensity matrix for specified orientation of the target relative to the incident wave, and for specified scattering directions., Comment: 59 pages. Updated to describe new target option CYLCAP, to correct error in formula for linear polarization P of scattered light, and provide more explicit instructions for use of multiple dielectric functions

DDSCAT.6.0 is a freely available software package (http://www.astro.princeton.edu/~draine/DDSCAT.6.0.html) which applies the "discrete dipole approximation" (DDA) to calculate scattering and absorption of electromagnetic waves by targets with arbitrary geometries and complex refractive index. DDSCAT.6.0 allows accurate calculations of electromagnetic scattering from targets with ``size parameters'' 2*pi*a/lambda < 15 provided the refractive index m is not large compared to unity (|m-1| < 1). DDSCAT.6.0 includes the option of using the FFTW (Fastest Fourier Transform in the West) package. DDSCAT.6.0 also includes MPI support, permitting parallel calculations on multiprocessor systems. DDSCAT package is written in Fortran and is highly portable. The program supports calculations for a variety of target geometries (e.g., ellipsoids, regular tetrahedra, rectangular solids, finite cylinders, hexagonal prisms, etc.). Target materials may be both inhomogeneous and anisotropic. It is straightforward for the user to ``import'' arbitrary target geometries into the code, and relatively straightforward to add new target generation capability to the package. DDSCAT automatically calculates total cross sections for absorption and scattering and selected elements of the Mueller scattering intensity matrix for specified orientation of the target relative to the incident wave, and for specified scattering directions. This User Guide explains how to use DDSCAT.6.0 to carry out electromagnetic scattering calculations. CPU and memory requirements are described., Comment: 46 pages, 6 figures

DDSCAT.5a is a freely available software package which applies the "discrete dipole approximation" (DDA) to calculate scattering and absorption of electromagnetic waves by targets with arbitrary geometries and complex refractive index. The DDA approximates the target by an array of polarizable points. DDSCAT.5a requires that these polarizable points be located on a cubic lattice. DDSCAT.5a10 allows accurate calculations of electromagnetic scattering from targets with "size parameters" 2 pi a/lambda < 15 provided the refractive index m is not large compared to unity (|m-1| < 1). The DDSCAT package is written in Fortran and is highly portable. The program supports calculations for a variety of target geometries (e.g., ellipsoids, regular tetrahedra, rectangular solids, finite cylinders, hexagonal prisms, etc.). Target materials may be both inhomogeneous and anisotropic. It is straightforward for the user to import arbitrary target geometries into the code, and relatively straightforward to add new target generation capability to the package. DDSCAT automatically calculates total cross sections for absorption and scattering and selected elements of the Mueller scattering intensity matrix for specified orientation of the target relative to the incident wave, and for specified scattering directions. This User Guide explains how to use DDSCAT.5a10 to carry out EM scattering calculations. CPU and memory requirements are described., Comment: Added more information for non-Unix sites, and new target option PRISM3 (triangular prism). 42 pages, 5 figures

The discrete-dipole approximation (DDA) is a flexible technique for computing scattering and absorption by targets of arbitrary geometry. In this paper we perform systematic study of various non-stationary iterative (conjugate gradient) methods in search for the most efficient one in order to solve the system of equations arising in DDA. We document implementation of these methods in our public domain code DDSCAT.5a

We report that the submerged microbubbles are an efficient source of diffuse radiance and may contribute to a rapid transition to the diffuse asymptotic regime. In this asymptotic regime an average cosine is easily predictable and measurable., Comment: 4 pages, 3 Postscript figures, opex2.sty (enclosed), also available from the Optical Society of America htpp://epubs.osa.org/oearchive/pdf/11948.pdf

A Monte Carlo model was used to study the scattering error of an absorption meter with a divergent light beam and a limited acceptance angle of the receiver. Reflections at both ends of the tube were taken into account. Calculations of the effect of varying optical properties of water, as well as the receiver geometry, were performed. A weighting function showing the scattering error quantitatively as a function of angle was introduced. Some cases of the practical interests are discussed., Comment: 17 pages, 7 Postscript figures, also available from http://atol.ucsd.edu/~pflatau. In Press, Journal of Atmospheric and Oceanic Technology, 2000

We report on the influence of submerged bubble clouds on the remote sensing properties of water. We show that the optical effect of bubbles on radiative transfer and on the estimate of the ocean color is significant. We present a global map of the volume fraction of air in water derived from daily wind speed data. This map, together with the parameterization of the microphysical properties, shows the possible significance of bubble clouds on the albedo of incoming solar energy, Comment: 17 pages, 6 Postscript figures, see also http://atol.ucsd.edu/~pflatau publications for related papers. Q. J. Roy. Met. Soc. in press 2000

ACE-Asia's extensive measurements from land, ocean, air and space quantified aerosol-radiation interactions. Results from each platform type, plus satellite-suborbital combinations, include: 1. Time series of multiwavelength aerosol optical depth (ADD), Angstrom exponent (alpha), single-scattering albedo (SSA), and size distribution from AERONET radiometry at 13 stations. In China and Korea AOD and alpha were strongly anticorrelated (reflecting transient dust events); dust volume-size modes peaked near 8 microns diameter; and SSA(dust) greater than SSA(pollution). 2. Calculations and measurements of photosynthetically active radiation and aerosols in China yield 24-h average downward surface radiative forcing per AOD(500 nm) of -27 W/sq m (400-700 nm). 3. The Hawaii-Japan cruise sampled a gradient with AOD(500 nm) extremes of 0.1 and 1.1. Shipboard measurements showed that adding dust to pollution increased SSA(550 nm, 55% RH), typically from -0.91 to approx. 0.97. Downwelling 8-12 micron radiances showed aerosol effects, especially in the major April dust event, with longwave forcing estimated at -5 to 15 W/sq m. 4. Extinction profiles from airborne sunphotometry and total-direct-diffuse radiometry show wavelength dependence often varying strongly with height, reflecting layering of dust-dominated over pollution-dominated aerosols. Comparing sunphotometric extinction profiles to those from in situ measurements (number and composition vs size, or scattering and absorption) shows layer heights agree, but extinction sometimes differs. 5. Airborne solar spectral flux radiometry yields absorption spectra for layers. Combining with AOD spectra yields best-fit aerosol single scattering albedo spectra. 6. Visible, NIR and total solar fluxes combined with AOD give radiative forcing efficiencies at surface and aloft.

A primary, ACE-Asia objective was to quantify the interactions between aerosols and radiation in the Asia-Pacific region. Toward this end, radiometric and related aerosol measurements were made from ocean, land, air and space platforms. Models that predict aerosol fields guided the measurements and are helping integrate and interpret results. Companion overview's survey these measurement and modeling components. Here we illustrate how these components were combined to determine aerosol radiative. impacts and their relation to aerosol properties. Because clouds can obscure or change aerosol direct radiative effects, aircraft and ship sorties to measure these effects depended on predicting and finding cloud-free areas and times with interesting aerosols present. Pre-experiment satellite cloud climatologies, pre-flight aerosol and cloud forecasts, and in-flight guidance from satellite imagery all helped achieve this. Assessments of aerosol regional radiative impacts benefit from the spatiotemporal coverage of satellites, provided satellite-retrieved aerosol properties are accurate. Therefore, ACE-Asia included satellite retrieval tests, as part of many comparisons to judge the consistency (closure) among, diverse measurements. Early results include: (1) Solar spectrally resolved and broadband irradiances and optical depth measurements from the C-130 aircraft and at Kosan, Korea yielded aerosol radiative forcing efficiencies, permitting comparisons between efficiencies of ACE-Asia and INDOEX aerosols, and between dust and "pollution" aerosols. Detailed results will be presented in separate papers. (2) Based on measurements of wavelength dependent aerosol optical depth (AOD) and single scattering albedo the estimated 24-h a average aerosol radiative forcing efficiency at the surface for photosynthetically active radiation (400 - 700 nm) in Yulin, China is approx. 30 W sq m per AOD(500 nm). (3) The R/V Brown cruise from Honolulu to Sea of Japan sampled an aerosol optical depth gradient, with AOD(500 nm) extremes from 0.1 to 1.1. On the Pacific transit from Honolulu to Hachijo AOD(500 nm) averaged 0.2, including increases to 0.4 after several storms, suggesting the strong impact of wind-generated seasalt. The AOD maximum, found in the Sea of Japan, was influenced by dust and anthropogenic sources. (4) In Beijing, single scattering albedo retrieved from AERONET sun-sky radiometry yielded midvisible SSA=0.88 with strong wavelength dependence, suggesting a significant black carbon component. SSA retrieved during dust episodes was approx. 0.90 and variable but wavelength neutral reflecting the presence of urban haze with the dust. Downwind at Anmyon Island SSA was considerably higher, approx. 0.94, but wavelength neutral for dust episodes and spectrally dependent during non dust periods. (5) Satellite retrievals show major aerosol features moving from Asia over the Pacific; however, determining seasonal-average aerosol effects is hampered by sampling frequency and large-scale cloud systems that obscure key parts of aerosol patterns. Preliminary calculations using, satellite-retrieved AOD fields and initial ACE-Asia aerosol properties (including sulfates, soot, and dust) yield clear-sky aerosol radiative effects in the seasonal-average ACE-Asia plume exceeding those of manmade greenhouse gases. Quantifying all-sky direct aerosol radiative effects is complicated by the need to define the height of absorbing aerosols with respect to cloud decks.

The topics discussed include the following: an overview of radiation and mesoscale flows and lessons learned from the intercomparison of GCM radiative codes.

A diffusion-like description of radiative transfer in clouds and the free atmosphere is often used. The two stream model is probably the best known example of such a description. The main idea behind the approach is that only the first few moments of radiance are needed to describe the radiative field correctly. Integration smooths details of the angular distribution of specific intensity and it is assumed that the closure parameters of the theory (diffusivity factors) are only weakly dependent on the distribution. The diffusivity factors are investigated using the results obtained from both Stratocumulus and Cirrus phases of FIRE experiment. A new theoretical framework is described in which two (upwards and downwards) diffusivity factors are used and a detailed multistream model is used to provide further insight about both the diffusivity factors and their dependence on scattering properties of clouds.