ARTS built-in documentation server
Workspace Method pndFromPsd
Calculates pnd_data from given psd_data for one scattering species.
Performs integration of the size distribution over the size grid
bin deriving pnd (units #/m3) from psd (units #/m3/m). Some checks
on the sufficiency of the size grid range and coverage are applied.
*quad_order* can be 0 for rectangular or 1 for trapezoidal
integration. The only difference is the treatment of the start and
end nodes. For trapezoidal their corresponding bins end exactly at
the nodes, while for rectangular they extend further out by the half
distance to the neighbor node (but not beyond 0).
Attempts to check that the size grids and scat_data represent the
bulk extinction sufficiently. Specifically, it is tested that
(a) psd*ext is decreasing at the small and large particle size
ends of the size grid - but only if scattering species bulk
extinction exceeds 1% of *threshold_ss_ext*.
(b) removing the smallest and largest particles changes the
resulting bulk extinction by less then a fraction of
*threshold_se_ext* - but only if scattering species bulk
extinction exceeds *threshold_ss_ext* and number density (pnd)
of the edge size point at this atmospheric level is larger
than *threshold_se_pnd* times the maximum pnd of this
scattering element over all atmospheric levels.
Skipping tests in case of low extinction is done in order to
minimize issues arising from very low mass densities,
particularly at single atmospheric levels, and very low bulk
extinctions, i.e. in cases where the effects on the radiance fields
are estimated to be low.
NOTE: The tests are only approximate and do not guarantee the
validity of the resulting bulk properties (and increasing the
thresholds will decrease the reliability of the bulk properties).
Authors: Jana Mendrok, Patrick Eriksson
|pndFromPsd( ||pnd_data, dpnd_data_dx, pnd_size_grid, psd_data, psd_size_grid, dpsd_data_dx, scat_data, f_grid, scat_data_checked, quad_order, scat_index, threshold_se_ext, threshold_ss_ext, threshold_se_pnd )|