OptimizeDoitPressureGrid

Workspace.OptimizeDoitPressureGrid(self: pyarts.arts._Workspace, p_grid: pyarts.arts.WorkspaceVariable | pyarts.arts.Vector | None = self.p_grid, pnd_field: pyarts.arts.WorkspaceVariable | pyarts.arts.Tensor4 | None = self.pnd_field, t_field: pyarts.arts.WorkspaceVariable | pyarts.arts.Tensor3 | None = self.t_field, scat_data_mono: pyarts.arts.WorkspaceVariable | pyarts.arts.ArrayOfArrayOfSingleScatteringData | None = self.scat_data_mono, z_field: pyarts.arts.WorkspaceVariable | pyarts.arts.Tensor3 | None = self.z_field, cloudbox_limits: pyarts.arts.WorkspaceVariable | pyarts.arts.ArrayOfIndex | None = self.cloudbox_limits, cloudbox_field_mono: pyarts.arts.WorkspaceVariable | pyarts.arts.Tensor6 | None = self.cloudbox_field_mono, pha_mat_doit: pyarts.arts.WorkspaceVariable | pyarts.arts.Tensor7 | None = self.pha_mat_doit, vmr_field: pyarts.arts.WorkspaceVariable | pyarts.arts.Tensor4 | None = self.vmr_field, p_grid_orig: pyarts.arts.WorkspaceVariable | pyarts.arts.Vector | None = self.p_grid_orig, f_grid: pyarts.arts.WorkspaceVariable | pyarts.arts.Vector | None = self.f_grid, f_index: pyarts.arts.WorkspaceVariable | pyarts.arts.Index | None = self.f_index, propmat_clearsky_agenda: pyarts.arts.WorkspaceVariable | pyarts.arts.Agenda | None = self.propmat_clearsky_agenda, tau_scat_max: pyarts.arts.WorkspaceVariable | pyarts.arts.Numeric | None = 0.1, sgl_alb_max: pyarts.arts.WorkspaceVariable | pyarts.arts.Numeric | None = 0.9, cloudbox_size_max: pyarts.arts.WorkspaceVariable | pyarts.arts.Index | None = 200, verbosity: pyarts.arts.WorkspaceVariable | pyarts.arts.Verbosity | None = self.verbosity) → None

Optimization of the pressure grid for RT calculation.

The methods consists of three parts:

  1. Calculate the single scattering albedo and the scattering optical thickness from the scattering and absorption species.

  2. Enhance z_field according to the two thresholds sgl_alb_max and tau_scat_max. If the resulting cloudbox size is bigger than cloudbox_size_max, this step is repeated with a higher threshold of tau_scat_max.

  3. Interpolate all variables used in doit_mono_agenda to the new z_field This method should be called inside doit_mono_agenda, right before cloudbox_field_monoIterate(). It can only be used if ScatSpeciesMerge() has been called and if it is called, cloudbox_field_monoOptimizeReverse() has to be called right after cloudbox_field_monoIterate() to interpolate cloudbox_field_mono back to the original size.

Optimization currently only works with stokes_dim = 1 .

Author(s): Jakob Doerr

Parameters:

  • p_grid (Vector, optional) – The pressure grid. See p_grid, defaults to self.p_grid [INOUT]

  • pnd_field (Tensor4, optional) – Particle number density field. See pnd_field, defaults to self.pnd_field [INOUT]

  • t_field (Tensor3, optional) – The field of atmospheric temperatures. See t_field, defaults to self.t_field [INOUT]

  • scat_data_mono (ArrayOfArrayOfSingleScatteringData, optional) – Monochromatic single scattering data. See scat_data_mono, defaults to self.scat_data_mono [INOUT]

  • z_field (Tensor3, optional) – The field of geometrical altitudes. See z_field, defaults to self.z_field [INOUT]

  • cloudbox_limits (ArrayOfIndex, optional) – The limits of the cloud box. See cloudbox_limits, defaults to self.cloudbox_limits [INOUT]

  • cloudbox_field_mono (Tensor6, optional) – Monochromatic radiation field inside the cloudbox. See cloudbox_field_mono, defaults to self.cloudbox_field_mono [INOUT]

  • pha_mat_doit (Tensor7, optional) – Ensemble averaged phase matrix for DOIT calculation. See pha_mat_doit, defaults to self.pha_mat_doit [INOUT]

  • vmr_field (Tensor4, optional) – VMR field. See vmr_field, defaults to self.vmr_field [INOUT]

  • p_grid_orig (Vector, optional) – The original pressure grid before optimization. See p_grid_orig, defaults to self.p_grid_orig [OUT]

  • f_grid (Vector, optional) – The frequency grid for monochromatic pencil beam calculations. See f_grid, defaults to self.f_grid [IN]

  • f_index (Index, optional) – Frequency index. See f_index, defaults to self.f_index [IN]

  • propmat_clearsky_agenda (Agenda, optional) – Calculate the absorption coefficient matrix. See propmat_clearsky_agenda, defaults to self.propmat_clearsky_agenda [IN]

  • tau_scat_max (Numeric, optional) – Maximum scattering optical thickness. Defaults to 0.1 [IN]

  • sgl_alb_max (Numeric, optional) – Maximum single scattering albedo. Defaults to 0.9 [IN]

  • cloudbox_size_max (Index, optional) – Maximum cloudbox size. Defaults to 200 [IN]

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