T-Matrix related workspace methods. More...

#include <cfloat>
#include <cmath>
#include <stdexcept>
#include "check_input.h"
#include "logic.h"
#include "math_funcs.h"
#include "messages.h"
#include "refraction.h"
#include "special_interp.h"
#include "tmatrix.h"

Functions
void	diameter_maxFromDiameter_volume_equ (Numeric &diameter_max, Numeric &diameter_aspect_area_max, const String &shape, const Numeric &diameter_volume_equ, const Numeric &aspect_ratio, const Verbosity &)
	WORKSPACE METHOD: diameter_maxFromDiameter_volume_equ. More...

void	diameter_volume_equFromDiameter_max (Numeric &diameter_volume_equ, Numeric &volume, const String &shape, const Numeric &diameter_max, const Numeric &aspect_ratio, const Verbosity &)
	WORKSPACE METHOD: diameter_volume_equFromDiameter_max. More...

void	scat_data_singleTmatrix (SingleScatteringData &scat_data_single, ScatteringMetaData &scat_meta_single, const GriddedField3 &complex_refr_index, const String &shape, const Numeric &diameter_volume_equ, const Numeric &aspect_ratio, const Numeric &mass, const String &ptype, const Vector &data_f_grid, const Vector &data_t_grid, const Vector &data_za_grid, const Vector &data_aa_grid, const Numeric &precision, const String &cri_source, const Index &ndgs, const Index &robust, const Index &quiet, const Verbosity &verbosity)
	WORKSPACE METHOD: scat_data_singleTmatrix. More...

void	TMatrixTest (const Verbosity &verbosity)
	WORKSPACE METHOD: TMatrixTest. More...

Variables
const Numeric	PI

Detailed Description

T-Matrix related workspace methods.

Author: Oliver Lemke

Date: 2013-06-25

Definition in file m_tmatrix.cc.

Function Documentation

◆ diameter_maxFromDiameter_volume_equ()

void diameter_maxFromDiameter_volume_equ	(	Numeric &	diameter_max,
		Numeric &	diameter_area_equ,
		const String &	shape,
		const Numeric &	diameter_volume_equ,
		const Numeric &	aspect_ratio,
		const Verbosity &	verbosity
	)

WORKSPACE METHOD: diameter_maxFromDiameter_volume_equ.

Calculates maximum and area equivalent diameters from volume equivalent diameter.

This is primarily a help function for using the T-matrix method and only a few particle shapes are handled. For shapes handled and further comments on the input arguments, see scat_data_singleTmatrix*.

Area equivalent diameter is the equivalent sphere diameter corresponding to the "maximum axial area". This is the largest cross-sectional area of the particle, observed either along the particle's main axis or in the perpendicular direction. That is, for a cylinder having diameter d and thickness h, this area is either (pi*d^2)/4 or (h*d).

Author: Johan Strandgren; Patrick Eriksson

Parameters

[out]	diameter_max	Generic output
[out]	diameter_area_equ	Generic output
[in]	shape	Generic Input
[in]	diameter_volume_equ	Generic Input
[in]	aspect_ratio	Generic Input

Definition at line 40 of file m_tmatrix.cc.

References max, PI, and pow().

Referenced by diameter_maxFromDiameter_volume_equ_g().

◆ diameter_volume_equFromDiameter_max()

void diameter_volume_equFromDiameter_max	(	Numeric &	diameter_volume_equ,
		Numeric &	volume,
		const String &	shape,
		const Numeric &	diameter_max,
		const Numeric &	aspect_ratio,
		const Verbosity &	verbosity
	)

WORKSPACE METHOD: diameter_volume_equFromDiameter_max.

Converts from maximum to volume equivalent diameter.

This is primarily a help function for using the T-matrix part and only a few particle shapes are handled. For shapes handled and further comments on the input arguments, see scat_data_singleTmatrix.

Also the volume is provided. It is simply sqrt(pi*dveq^3/6).

Author: Johan Strandgren; Patrick Eriksson

Parameters

[out]	diameter_volume_equ	Generic output
[out]	volume	Generic output
[in]	shape	Generic Input
[in]	diameter_max	Generic Input
[in]	aspect_ratio	Generic Input

Definition at line 82 of file m_tmatrix.cc.

References PI, and pow().

Referenced by diameter_volume_equFromDiameter_max_g().

◆ scat_data_singleTmatrix()

void scat_data_singleTmatrix	(	SingleScatteringData &	scat_data_single,
		ScatteringMetaData &	scat_meta_single,
		const GriddedField3 &	complex_refr_index,
		const String &	shape,
		const Numeric &	diameter_volume_equ,
		const Numeric &	aspect_ratio,
		const Numeric &	mass,
		const String &	ptype,
		const Vector &	data_f_grid,
		const Vector &	data_t_grid,
		const Vector &	data_za_grid,
		const Vector &	data_aa_grid,
		const Numeric &	precision,
		const String &	cri_source,
		const Index &	ndgs,
		const Index &	robust,
		const Index &	quiet,
		const Verbosity &	verbosity
	)

WORKSPACE METHOD: scat_data_singleTmatrix.

A basic interface to Mishchenko's T-matrix code linked to ARTS.

The method performs a T-matrix calculation for a single scattering element, i.e. a combination of particle shape, size, aspect ratio and orientation.

Particle shape (shape) has two options: "spheroidal" and "cylindrical"

Particle size (diameter_volume_equ) is given as the equivalent volume sphere diameter. That is, the diameter obtained if all the particle's material is rearranged into a (solid) sphere.

Particle aspect ratio ar (aspect_ratio) is a numeric value, defined according to Mishchenko's definition as ratio of horizontal axis a to vertical (rotational) axis b: ar=a/b. That is, oblates have ar>1, prolates ar<1. Perfect spheres (spheroidals with ar=1) can trigger numerical issues. To avoid these, we internally increase their aspect ratio by 1e-6, i.e. turning perfect spheres into very light oblates.

Particle type (ptype) has two options: "totally_random" and "azimuthally_random" For totally randomly oriented particles, data_aa_grid is not taken into account (but a Vector type container needs to be passed).

For further information on how aspect ratio and the different shapes and orientations are defined, see the documentation of the T-matrix code found http://www.giss.nasa.gov/staff/mmishchenko/t_matrix.html

Regarding ndgs, we refer to the this comment from the documentation: "Parameter controlling the number of division points in computing integrals over the particle surface. For compact particles, the recommended value is 2. For highly aspherical particles larger values (3, 4,...) may be necessary to obtain convergence. The code does not check convergence over this parameter. Therefore, control comparisons of results obtained with different NDGS-values are recommended."

Author: Johan Strandgren; Patrick Eriksson

Parameters

[out]	scat_data_single	WS Output
[out]	scat_meta_single	WS Output
[in]	complex_refr_index	WS Input
[in]	shape	Generic Input
[in]	diameter_volume_equ	Generic Input
[in]	aspect_ratio	Generic Input
[in]	mass	Generic Input (Default: "NaN")
[in]	ptype	Generic Input
[in]	data_f_grid	Generic Input
[in]	data_t_grid	Generic Input
[in]	data_za_grid	Generic Input
[in]	data_aa_grid	Generic Input (Default: "[]")
[in]	precision	Generic Input (Default: "0.001")
[in]	cri_source	Generic Input (Default: "Set by user, unknown source.")
[in]	ndgs	Generic Input (Default: "2")
[in]	robust	Generic Input (Default: "0")
[in]	quiet	Generic Input (Default: "1")