next up previous contents
Next: Running DDSCAT.5a Using the Up: User Guide for Previous: Moving the Executable

 

The Parameter File ddscat.par

The directory DDA should now contain a sample file ddscat.par which provides parameters to the program ddscat. As provided (see AppendixA), the file ddscat.par is set up to calculate scattering by a 8tex2html_wrap_inline32116tex2html_wrap_inline32114 rectangular array of 192 dipoles, with an effective radius tex2html_wrap_inline3237, at a wavelength of tex2html_wrap_inline3239 (for a ``size parameter'' tex2html_wrap_inline3241).

The dielectric function of the target material is provided in the file diel.tab, which is a sample file in which the refractive index is set to m=1.33+0.01i at all wavelengths; the name of this file is provided to ddscat by the parameter file ddscat.par.

The sample parameter file as supplied calls for the new GPFA FFT routine (GPFAFT) of Temperton (1992) to be employed and the PBCGST iterative method to be used for solving the system of linear equations. (See section §12 and §13 for discussion of choice of FFT algorithm and choice of equation-solving algorithm.)

The sample parameter file specifies (via option LATTDR) that the ``Lattice Dispersion Relation'' of Draine and Goodman (1993) be employed to determine the dipole polarizabilities. See §10 for discussion of other options.

The sample ddscat.par file specifies that the calculations be done for a single wavelength (tex2html_wrap_inline3239) and a single effective radius (tex2html_wrap_inline3247). Note that in DDSCAT.5a the ``effective radius'' tex2html_wrap_inline3249 is the radius of a sphere of equal volume - i.e., a sphere of volume tex2html_wrap_inline3251 , where d is the lattice spacing and N is the number of occupied (i.e., non-vacuum) lattice sites in the target. Thus the effective radius tex2html_wrap_inline3257 .

The incident radiation is always assumed to propagate along the x axis in the ``Lab Frame''. The sample ddscat.par file specifies incident polarization state tex2html_wrap_inline3261 to be along the y axis (and consequently polarization state tex2html_wrap_inline3265 will automatically be taken to be along the z axis). IORTH=2 in ddscat.par calls for calculations to be carried out for both incident polarization states (tex2html_wrap_inline3261 and tex2html_wrap_inline3265 - see 20).

The target is assumed to have two vectors tex2html_wrap_inline3221 and tex2html_wrap_inline3223 embedded in it; tex2html_wrap_inline3223 is perpendicular to tex2html_wrap_inline3221. In the case of the 8tex2html_wrap_inline32116tex2html_wrap_inline32114 rectangular array of the sample calculation, the vector tex2html_wrap_inline3221 is along the ``long'' axis of the target, and the vector tex2html_wrap_inline3223 is along the ``intermediate'' axis. The target orientation in the Lab Frame is set by three angles: tex2html_wrap_inline3289, tex2html_wrap_inline3291, and tex2html_wrap_inline3293, defined and discussed below in §16. Briefly, the polar angles tex2html_wrap_inline3291 and tex2html_wrap_inline3293 specify the direction of tex2html_wrap_inline3221 in the Lab Frame. The target is assumed to be rotated around tex2html_wrap_inline3221 by an angle tex2html_wrap_inline3289. The sample ddscat.par file specifies tex2html_wrap_inline3305 and tex2html_wrap_inline3307 (see lines in ddscat.par specifying variables BETA and PHI), and calls for three values of the angle tex2html_wrap_inline3291 (see line in ddscat.par specifying variable THETA). DDSCAT.5a chooses tex2html_wrap_inline3291 values uniformly spaced in tex2html_wrap_inline3313; thus, asking for three values of tex2html_wrap_inline3291 between 0 and tex2html_wrap_inline3317 yields tex2html_wrap_inline3319, tex2html_wrap_inline3321, and tex2html_wrap_inline3317.

Appendix A provides a detailed description of the file ddscat.par.gif


next up previous contents
Next: Running DDSCAT.5a Using the Up: User Guide for Previous: Moving the Executable

Bruce Draine
Thu Aug 10 09:34:16 EDT 2000