psdMcFarquaharHeymsfield97

McFarquahar and Heymsfield [1997] particle size distribution for cloud ice.

This is a 1-parameter PSD, i.e. pnd_agenda_input shall have one column and pnd_agenda_input_names shall contain a single string. The input data in pnd_agenda_input shall be ice hydrometeor mass content in unit of [kg/m3]. The naming used is pnd_agenda_input_names is free but the same name must be used in particle_bulkprop_names and dpnd_data_dx_names.

psd_size_grid shall contain size in terms of volume equivalent diameter.

Derivatives are obtained by perturbation of 0.1%, but not less than 1e-9 kg/m3.

The validity range of mass content is not limited. Negative mass contents will produce negative psd values following a distribution given by abs(IWC), ie. abs(psd)=f(abs(IWC)).

If temperature is outside [ t_min , t_max ] psd=0 and dpsd=0 if picky=0, or an error is thrown if picky=1.

For temperatures below t_min_psd`, the size distribution is calculated for T = t_min_psd`. Likewise, for temperatures above t_max_psd, the distribution is derived for T = t_max_psd.

Defaults of t_min_psd and t_max_psd were set considering that the parametrization has been derived from measurements over temperatures of -70C to -20C. The noisy option can not be used together with calculation of derivatives (ie. when dpnd_data_dx_names is not empty).

Author(s): Patrick Eriksson, Jana Mendrok

Parameters:

psd_data (Matrix, optional) – Particle size distribution values for a set of points. See psd_data, defaults to self.psd_data [OUT]
dpsd_data_dx (Tensor3, optional) – Partial derivates of psd_data. See dpsd_data_dx, defaults to self.dpsd_data_dx [OUT]
psd_size_grid (Vector, optional) – The particle sizes associated with psd_data. See psd_size_grid, defaults to self.psd_size_grid [IN]
pnd_agenda_input_t (Vector, optional) – Temperature input to one element of pnd_agenda_array. See pnd_agenda_input_t, defaults to self.pnd_agenda_input_t [IN]
pnd_agenda_input (Matrix, optional) – The variable input to one element of pnd_agenda_array. See pnd_agenda_input, defaults to self.pnd_agenda_input [IN]
pnd_agenda_input_names (ArrayOfString, optional) – Naming of (existing or expected) data in pnd_agenda_input. See pnd_agenda_input_names, defaults to self.pnd_agenda_input_names [IN]
dpnd_data_dx_names (ArrayOfString, optional) – Selection of partial derivatives of pnd_data. See dpnd_data_dx_names, defaults to self.dpnd_data_dx_names [IN]
scat_species_a (Numeric, optional) – Mass-size relationship parameter, for one scattering species. See scat_species_a, defaults to self.scat_species_a [IN]
scat_species_b (Numeric, optional) – Mass-size relationship parameter, for one scattering species. See scat_species_b, defaults to self.scat_species_b [IN]
t_min (Numeric, optional) – Low temperature limit to calculate a psd. Defaults to 0 [IN]
t_max (Numeric, optional) – High temperature limit to calculate a psd. Defaults to 280. [IN]
t_min_psd (Numeric, optional) – Low temperature limit to use as paramtrization temperature. Defaults to 180 [IN]
t_max_psd (Numeric, optional) – High temperature limit to use as paramtrization temperature. Defaults to 273.15 [IN]
picky (Index, optional) – Flag whether to be strict with parametrization value checks. Defaults to 0 [IN]
noisy (Index, optional) – Distribution parameter perturbance flag. Defaults to 0 [IN]
verbosity (Verbosity) – ARTS verbosity. See verbosity, defaults to self.verbosity [IN]