psdFieldEtAl07

Workspace.psdFieldEtAl07(self: pyarts.arts._Workspace, psd_data: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Matrix]] = self.psd_data, dpsd_data_dx: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Tensor3]] = self.dpsd_data_dx, psd_size_grid: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Vector]] = self.psd_size_grid, pnd_agenda_input_t: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Vector]] = self.pnd_agenda_input_t, pnd_agenda_input: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Matrix]] = self.pnd_agenda_input, pnd_agenda_input_names: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.ArrayOfString]] = self.pnd_agenda_input_names, dpnd_data_dx_names: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.ArrayOfString]] = self.dpnd_data_dx_names, scat_species_a: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Numeric]] = self.scat_species_a, scat_species_b: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Numeric]] = self.scat_species_b, regime: Union[pyarts.arts.WorkspaceVariable, pyarts.arts.String], t_min: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Numeric]] = 0, t_max: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Numeric]] = 290., t_min_psd: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Numeric]] = 200., t_max_psd: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Numeric]] = 273.15, beta_min: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Numeric]] = 1.01, beta_max: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Numeric]] = 4, picky: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Index]] = 0, verbosity: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Verbosity]] = self.verbosity) → None

The Field et al. [2007] particle size distribution for snow and 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 maximum diameter.

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

Both parametrization for tropics and midlatitudes are handled, governed by setting of regime, where “TR” selectes the tropical case, and “ML” the midlatitude one.

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 -60C to 0C. Checks of the sanity of the mass-dimension relationship are performed Errors are thrown if:

  • Mass-dimension relation exponent scat_species_b is outside [beta_min, beta_max].

Author(s): 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]

  • regime (String) – Parametrization regime (“TR”=tropical or “ML”=midlatitude). [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 290. [IN]

  • t_min_psd (Numeric, optional) – Low temperature limit to use as paramtrization temperature. Defaults to 200. [IN]

  • t_max_psd (Numeric, optional) – High temperature limit to use as paramtrization temperature. Defaults to 273.15 [IN]

  • beta_min (Numeric, optional) – Low b limit (only if picky). Defaults to 1.01 [IN]

  • beta_max (Numeric, optional) – High b limit (only if picky). Defaults to 4 [IN]

  • picky (Index, optional) – Flag whether to be strict with parametrization value checks. Defaults to 0 [IN]

  • verbosity (Verbosity) – ARTS verbosity. See verbosity, defaults to self.verbosity [IN]