ARTS 2.5.9 (git: 825fa5f2)
m_cloudradar.cc File Reference

Workspace functions related to simulation of cloud radars. More...

#include <cmath>
#include <stdexcept>
#include "arts.h"
#include "arts_constants.h"
#include "arts_omp.h"
#include "auto_md.h"
#include "logic.h"
#include "matpack_eigen.h"
#include "messages.h"
#include "montecarlo.h"
#include "propagationmatrix.h"
#include "rte.h"
#include "sensor.h"

Go to the source code of this file.

Functions

void iyRadarSingleScat (Workspace &ws, Matrix &iy, ArrayOfMatrix &iy_aux, ArrayOfTensor3 &diy_dx, Vector &ppvar_p, Vector &ppvar_t, Matrix &ppvar_vmr, Matrix &ppvar_wind, Matrix &ppvar_mag, Matrix &ppvar_pnd, Matrix &ppvar_f, const Index &stokes_dim, const Vector &f_grid, const Index &atmosphere_dim, const Vector &p_grid, const Tensor3 &t_field, const EnergyLevelMap &nlte_field, const Tensor4 &vmr_field, const ArrayOfArrayOfSpeciesTag &abs_species, const Tensor3 &wind_u_field, const Tensor3 &wind_v_field, const Tensor3 &wind_w_field, const Tensor3 &mag_u_field, const Tensor3 &mag_v_field, const Tensor3 &mag_w_field, const Index &cloudbox_on, const ArrayOfIndex &cloudbox_limits, const Tensor4 &pnd_field, const ArrayOfTensor4 &dpnd_field_dx, const ArrayOfString &scat_species, const ArrayOfArrayOfSingleScatteringData &scat_data, const Index &scat_data_checked, const ArrayOfString &iy_aux_vars, const Index &jacobian_do, const ArrayOfRetrievalQuantity &jacobian_quantities, const Ppath &ppath, const Matrix &iy_transmitter, const Agenda &propmat_clearsky_agenda, const Agenda &water_p_eq_agenda, const Numeric &rte_alonglos_v, const Index &trans_in_jacobian, const Numeric &pext_scaling, const Index &t_interp_order, const Verbosity &verbosity)
 WORKSPACE METHOD: iyRadarSingleScat. More...
 
void yRadar (Workspace &ws, Vector &y, Vector &y_f, ArrayOfIndex &y_pol, Matrix &y_pos, Matrix &y_los, ArrayOfVector &y_aux, Matrix &y_geo, Matrix &jacobian, const Index &atmgeom_checked, const Index &atmfields_checked, const String &iy_unit_radar, const ArrayOfString &iy_aux_vars, const Index &stokes_dim, const Vector &f_grid, const Index &cloudbox_on, const Index &cloudbox_checked, const Matrix &sensor_pos, const Matrix &sensor_los, const Index &sensor_checked, const Index &jacobian_do, const ArrayOfRetrievalQuantity &jacobian_quantities, const Agenda &iy_radar_agenda, const ArrayOfArrayOfIndex &instrument_pol_array, const Vector &range_bins, const Numeric &ze_tref, const Numeric &k2, const Numeric &dbze_min, const Verbosity &)
 WORKSPACE METHOD: yRadar. More...
 
void particle_bulkpropRadarOnionPeeling (Workspace &ws, Tensor4 &particle_bulkprop_field, ArrayOfString &particle_bulkprop_names, const Index &atmosphere_dim, const Vector &p_grid, const Vector &lat_grid, const Vector &lon_grid, const Tensor3 &t_field, const Tensor3 &z_field, const Tensor4 &vmr_field, const Matrix &z_surface, const Index &atmfields_checked, const Index &atmgeom_checked, const Vector &f_grid, const Agenda &propmat_clearsky_agenda, const ArrayOfString &scat_species, const ArrayOfGriddedField3 &invtable, const Matrix &incangles, const Tensor3 &dBZe, const Numeric &dbze_noise, const Matrix &h_clutter, const Index &fill_clutter, const Numeric &t_phase, const Numeric &wc_max, const Numeric &wc_clip, const Index &do_atten_abs, const Index &do_atten_hyd, const Numeric &atten_hyd_scaling, const Numeric &atten_hyd_max, const Verbosity &)
 WORKSPACE METHOD: particle_bulkpropRadarOnionPeeling. More...
 
void RadarOnionPeelingTableCalc (Workspace &ws, ArrayOfGriddedField3 &invtable, const Vector &f_grid, const ArrayOfString &scat_species, const ArrayOfArrayOfSingleScatteringData &scat_data, const ArrayOfArrayOfScatteringMetaData &scat_meta, const ArrayOfAgenda &pnd_agenda_array, const ArrayOfArrayOfString &pnd_agenda_array_input_names, const Index &i_species, const Vector &dbze_grid, const Vector &t_grid, const Numeric &wc_min, const Numeric &wc_max, const Numeric &ze_tref, const Numeric &k2, const Verbosity &verbosity)
 WORKSPACE METHOD: RadarOnionPeelingTableCalc. More...
 

Variables

constexpr Numeric PI =Constant::pi
 
constexpr Numeric SPEED_OF_LIGHT =Constant::speed_of_light
 
constexpr Numeric LOG10_EULER_NUMBER =Constant::log10_euler
 
const Index GFIELD3_DB_GRID = 1
 
const Index GFIELD3_T_GRID = 2
 

Detailed Description

Workspace functions related to simulation of cloud radars.

Author
Patrick Eriksson patri.nosp@m.ck.e.nosp@m.rikss.nosp@m.on@c.nosp@m.halme.nosp@m.rs.s.nosp@m.e
Date
2010-10-31

Definition in file m_cloudradar.cc.

Function Documentation

◆ iyRadarSingleScat()

void iyRadarSingleScat ( Workspace ws,
Matrix iy,
ArrayOfMatrix iy_aux,
ArrayOfTensor3 diy_dx,
Vector ppvar_p,
Vector ppvar_t,
Matrix ppvar_vmr,
Matrix ppvar_wind,
Matrix ppvar_mag,
Matrix ppvar_pnd,
Matrix ppvar_f,
const Index stokes_dim,
const Vector f_grid,
const Index atmosphere_dim,
const Vector p_grid,
const Tensor3 t_field,
const EnergyLevelMap nlte_field,
const Tensor4 vmr_field,
const ArrayOfArrayOfSpeciesTag abs_species,
const Tensor3 wind_u_field,
const Tensor3 wind_v_field,
const Tensor3 wind_w_field,
const Tensor3 mag_u_field,
const Tensor3 mag_v_field,
const Tensor3 mag_w_field,
const Index cloudbox_on,
const ArrayOfIndex cloudbox_limits,
const Tensor4 pnd_field,
const ArrayOfTensor4 dpnd_field_dx,
const ArrayOfString scat_species,
const ArrayOfArrayOfSingleScatteringData scat_data,
const Index scat_data_checked,
const ArrayOfString iy_aux_vars,
const Index jacobian_do,
const ArrayOfRetrievalQuantity jacobian_quantities,
const Ppath ppath,
const Matrix iy_transmitter,
const Agenda propmat_clearsky_agenda,
const Agenda water_p_eq_agenda,
const Numeric rte_alonglos_v,
const Index trans_in_jacobian,
const Numeric pext_scaling,
const Index t_interp_order,
const Verbosity verbosity 
)

WORKSPACE METHOD: iyRadarSingleScat.

Simulation of radar, restricted to single scattering.

The WSM treats e.g. radar measurements of cloud and precipitation, on the condition that multiple scattering can be ignored. Beside the direct back-scattering, the two-way attenuation by gases and particles is considered. Surface scattering/clutter is ignored.

The method could potentially be used for lidars, but multiple scattering poses here a must stronger constrain for the range of applications.

The method shall be used with yRadar, NOT with yCalc.

The ppath provided should be calculated including cloudbox interior: <br> ppathCalc( cloudbox_on=0 )

The method returns the back-scattering for each point of ppath. Several frequencies can be treated in parallel. The size of iy is [ nf*np, stokes_dim ], where nf is the length of f_grid and np is the number of path points. The data are stored in blocks of [ np, stokes_dim ]. That is, all the results for the first frequency occupy the np first rows of iy etc.

The polarisation state of the transmitted pulse is taken from iy_transmitter*. If the radar transmits several polarisations at the same frequency, you need to handle this by using two frequencies in f_grid, but these can be almost identical.

This method does not consider iy_unit_radar. Unit changes are instead applied in *yRadar. The output of this method matches the option "1".

The extinction due to particles can be scaled (by pext_scaling), which could be of interest when e.g. characterising inversions or trying to compensate for ignored multiple scattering. The later is commented further for particle_bulkpropRadarOnionPeeling.

For Jacobian calculations the default is to assume that the transmittance is unaffected by the retrieval quantities. This is done to save computational time, and should be a valid approximation for the single-scattering conditions. Set trans_in_jacobian to 1 to activate full Jacobian calculations.

Some auxiliary radiative transfer quantities can be obtained. Auxiliary quantities are selected by iy_aux_vars and returned by iy_aux. Valid choices for auxiliary data are: <br> "Radiative background": Index value flagging the radiative <br> background. The following coding is used: 0=space, 1=surface <br> and 2=cloudbox (the last case should not occur!). Only column <br> matching first Stokes element filled. Other columns are set to 0. <br> "Backscattering": The unattenuated back-scattering. That is, as <br> iy but with no attenuated applied. Here all columns are filled. <br> By combing iy and this auxiliary variable, the total two-way <br> attenuation can be derived. <br> "Abs species extinction": Extinction due to abs_species at each <br> ppath point, taken as the diagonal of the local extinction matrix. <br> "Particle extinction": Extinction due to particles at each <br> ppath point, taken as the diagonal of the local extinction matrix. <br> The retunred values includes pext_scaling

Author
Patrick Eriksson
Parameters
[in,out]wsWorkspace
[out]iyWS Output
[out]iy_auxWS Output
[out]diy_dxWS Output
[out]ppvar_pWS Output
[out]ppvar_tWS Output
[out]ppvar_vmrWS Output
[out]ppvar_windWS Output
[out]ppvar_magWS Output
[out]ppvar_pndWS Output
[out]ppvar_fWS Output
[in]stokes_dimWS Input
[in]f_gridWS Input
[in]atmosphere_dimWS Input
[in]p_gridWS Input
[in]t_fieldWS Input
[in]nlte_fieldWS Input
[in]vmr_fieldWS Input
[in]abs_speciesWS Input
[in]wind_u_fieldWS Input
[in]wind_v_fieldWS Input
[in]wind_w_fieldWS Input
[in]mag_u_fieldWS Input
[in]mag_v_fieldWS Input
[in]mag_w_fieldWS Input
[in]cloudbox_onWS Input
[in]cloudbox_limitsWS Input
[in]pnd_fieldWS Input
[in]dpnd_field_dxWS Input
[in]scat_speciesWS Input
[in]scat_dataWS Input
[in]scat_data_checkedWS Input
[in]iy_aux_varsWS Input
[in]jacobian_doWS Input
[in]jacobian_quantitiesWS Input
[in]ppathWS Input
[in]iy_transmitterWS Input
[in]propmat_clearsky_agendaWS Input
[in]water_p_eq_agendaWS Input
[in]rte_alonglos_vWS Input
[in]trans_in_jacobianGeneric Input (Default: "0")
[in]pext_scalingGeneric Input (Default: "1")
[in]t_interp_orderGeneric Input (Default: "1")

Definition at line 58 of file m_cloudradar.cc.

References a, abs(), adapt_stepwise_partial_derivatives(), ARTS_USER_ERROR, ARTS_USER_ERROR_IF, bulk_backscatter(), bulk_backscatter_derivative(), chk_if_in_range(), CommutativeTransmission, cumulative_transmission(), FOR_ANALYTICAL_JACOBIANS_DO, Forward, get_pointers_for_analytical_species(), get_pointers_for_scat_species(), get_ppath_atmvars(), get_ppath_cloudvars(), get_ppath_f(), get_standard_diy_dpath(), get_standard_starting_diy_dx(), get_stepwise_clearsky_propmat(), get_stepwise_scattersky_propmat(), joker, PropagationMatrix::Kjj(), Ppath::los, Ppath::lstep, min(), mirror_los(), ConstTensor4View::nbooks(), ConstMatrixView::ncols(), Array< base >::nelem(), ConstVectorView::nelem(), Ppath::np, ConstMatrixView::nrows(), pha_mat_1ScatElem(), ppath_what_background(), Vector::resize(), Matrix::resize(), Reverse, matpack::eigen::row_vec(), rtmethods_jacobian_finalisation(), set_backscatter_radiation_vector(), stepwise_transmission(), swap(), TotalNumberOfElements(), and RadiationVector::Vec().

Referenced by iyRadarSingleScat_g().

◆ particle_bulkpropRadarOnionPeeling()

void particle_bulkpropRadarOnionPeeling ( Workspace ws,
Tensor4 particle_bulkprop_field,
ArrayOfString particle_bulkprop_names,
const Index atmosphere_dim,
const Vector p_grid,
const Vector lat_grid,
const Vector lon_grid,
const Tensor3 t_field,
const Tensor3 z_field,
const Tensor4 vmr_field,
const Matrix z_surface,
const Index atmfields_checked,
const Index atmgeom_checked,
const Vector f_grid,
const Agenda propmat_clearsky_agenda,
const ArrayOfString scat_species,
const ArrayOfGriddedField3 invtable,
const Matrix incangles,
const Tensor3 dBZe,
const Numeric dbze_noise,
const Matrix h_clutter,
const Index fill_clutter,
const Numeric t_phase,
const Numeric wc_max,
const Numeric wc_clip,
const Index do_atten_abs,
const Index do_atten_hyd,
const Numeric atten_hyd_scaling,
const Numeric atten_hyd_max,
const Verbosity verbosity 
)

WORKSPACE METHOD: particle_bulkpropRadarOnionPeeling.

Inverts radar reflectivities by in an onion peeling manner.

The method assumes space-based measurements and invert one altitude at the time, based on a pre-calculated inversion table (invtable) and starting at the top of the atmosphere. If attenuation is completely ignored, the table is effectively used as a look-up table to map dBZe to hydrometeor values. The method considers attenuation by default, where extinction due to hydrometeors is taken from the table and the one due to abs_species is obtained by propmat_clearsky_agenda*.

The inversion table consists of two GriddedField3. The first field shall match liquid hydrometeors and is applied for temperatures above t_phase*. The second field is applied for lower temperatures and shall thus correspond to ice hydrometeors.

The size of each field is (2,ndb,nt). The two page dimensions match the hydrometeor property to retrieve and extinction, respectively. The table shall hold the 10-logarithm of the property, such as log10(IWC). ndb is the number of dBZe values in the table and nt the number of temperatures. The table is interpolated in temperature in a nearest neighbour fashion, while in a linear interpolation is applied in the dBZe dimension.

The field of radar reflectivities (dBZe) shall cover the complete atmosphere and then match e.g. t_field in size. The observation geometry is here specified by giving the incidence angle for each profile of dBZe values (by incangles). A flat Earth approximation is applied inside the method.

All values below dbze_noise are treated as pure noise and particle_bulkprop_field* is set to zero for these positions. The comparison to dbze_noise is done with uncorrected values.

Further, all values at altitudes below z_surface + h_clutter are assumed to be surface clutter and are rejected. If fill_clutter is set to 1, the retrieval just above the clutter zone is assumed valid also below and is copied to all altitudes below (also for altitudes below the surface).

Unfiltered clutter can cause extremely high retrived water contents. The GIN wc_max defines an upper limit for reasonable water contents. Retrievals ending up above this value are set to zero. Values below wc_max* but above wc_clip, are set to wc_clip.

Significant radar echos (>dbze_noise and above clutter zone) are assumed to match liquid hydrometeors for temperatures >= t_phase and ice ones for lower temperatures.

Default is to consider attenuation of both hydrometeors and absorption species. These two sources to attenuation can be ignored by setting do_atten_hyd* and do_atten_abs to zero, respectively.

Default is to consider hydrometeor attenuation, but there could be two reasons to ignore it. It can cause a "run away" effect in the retrievals. Ignoring it can also compensate for impact of multiple scattering in space-based observations, as shown by: Matrosov and Battaglia, GRL, 2009. However, ignoring the hydrometeor attenuation totally gives a too high compensating effect and the GIN atten_hyd_scaling* allows to test intermediate compensations. This GIN matches the GIN pext_scaling of iyRadarSingleScat, but they have different default values. The default in this method follows the results for CloudSat in Matrosov and Battaglia. Please note that do_atten_hyd* must be true to apply atten_hyd_scaling.

Even with atten_hyd_scaling below 1, there could be a run-away in the estimated attenuation, and atten_hyd_max stops this by setting a maximum value to the hydrometeor attenuation.

Author
Patrick Eriksson
Parameters
[in,out]wsWorkspace
[out]particle_bulkprop_fieldWS Output
[out]particle_bulkprop_namesWS Output
[in]atmosphere_dimWS Input
[in]p_gridWS Input
[in]lat_gridWS Input
[in]lon_gridWS Input
[in]t_fieldWS Input
[in]z_fieldWS Input
[in]vmr_fieldWS Input
[in]z_surfaceWS Input
[in]atmfields_checkedWS Input
[in]atmgeom_checkedWS Input
[in]f_gridWS Input
[in]propmat_clearsky_agendaWS Input
[in]scat_speciesWS Input
[in]invtableGeneric Input
[in]incanglesGeneric Input
[in]dBZeGeneric Input
[in]dbze_noiseGeneric Input (Default: "-99")
[in]h_clutterGeneric Input
[in]fill_clutterGeneric Input (Default: "0")
[in]t_phaseGeneric Input (Default: "273.15")
[in]wc_maxGeneric Input (Default: "10e-3")
[in]wc_clipGeneric Input (Default: "5e-3")
[in]do_atten_absGeneric Input (Default: "1")
[in]do_atten_hydGeneric Input (Default: "1")
[in]atten_hyd_scalingGeneric Input (Default: "0.5")
[in]atten_hyd_maxGeneric Input (Default: "3")

Definition at line 875 of file m_cloudradar.cc.

References abs(), ARTS_USER_ERROR, ARTS_USER_ERROR_IF, chk_atm_field(), chk_atm_surface(), chk_if_in_range(), Conversion::cosd(), GridPos::fd, GFIELD3_DB_GRID, GFIELD3_T_GRID, gridpos(), GridPos::idx, interp(), interpweights(), joker, PropagationMatrix::Kjj(), last(), LOG10_EULER_NUMBER, max(), min(), ConstTensor3View::ncols(), ConstMatrixView::ncols(), Array< base >::nelem(), ConstTensor3View::npages(), ConstTensor3View::nrows(), ConstMatrixView::nrows(), pow(), propmat_clearsky_agendaExecute(), Tensor4::resize(), and Matrix::resize().

Referenced by particle_bulkpropRadarOnionPeeling_g().

◆ RadarOnionPeelingTableCalc()

void RadarOnionPeelingTableCalc ( Workspace ws,
ArrayOfGriddedField3 invtable,
const Vector f_grid,
const ArrayOfString scat_species,
const ArrayOfArrayOfSingleScatteringData scat_data,
const ArrayOfArrayOfScatteringMetaData scat_meta,
const ArrayOfAgenda pnd_agenda_array,
const ArrayOfArrayOfString pnd_agenda_array_input_names,
const Index i_species,
const Vector dbze_grid,
const Vector t_grid,
const Numeric wc_min,
const Numeric wc_max,
const Numeric ze_tref,
const Numeric k2,
const Verbosity verbosity 
)

WORKSPACE METHOD: RadarOnionPeelingTableCalc.

Creates a radar inversion table.

This method is tailored to make inversion tables that fit particle_bulkpropRadarOnionPeeling*. See that method for format of the table.

The method needs to be called twice to form a complete table, once for liquid and ice hydrometeors. The table can be empty at the first call.

The input data (scat_data etc.) must match two scattering species and a single frequency (the one of the radar).

Author
Patrick Eriksson
Parameters
[in,out]wsWorkspace
[out]invtableGeneric output
[in]f_gridWS Input
[in]scat_speciesWS Input
[in]scat_dataWS Input
[in]scat_metaWS Input
[in]pnd_agenda_arrayWS Input
[in]pnd_agenda_array_input_namesWS Input
[in]i_speciesGeneric Input
[in]dbze_gridGeneric Input
[in]t_gridGeneric Input
[in]wc_minGeneric Input (Default: "1e-8")
[in]wc_maxGeneric Input (Default: "2e-2")
[in]ze_trefGeneric Input (Default: "273.15")
[in]k2Generic Input (Default: "-1")

Definition at line 1133 of file m_cloudradar.cc.

References ARTS_USER_ERROR, ARTS_USER_ERROR_IF, b, gridpos(), interp(), interpweights(), is_increasing(), joker, Array< base >::nelem(), ConstVectorView::nelem(), pnd_agenda_arrayExecute(), PTYPE_TOTAL_RND, transform(), VectorLogSpace(), w, and ze_cfac().

Referenced by RadarOnionPeelingTableCalc_g().

◆ yRadar()

void yRadar ( Workspace ws,
Vector y,
Vector y_f,
ArrayOfIndex y_pol,
Matrix y_pos,
Matrix y_los,
ArrayOfVector y_aux,
Matrix y_geo,
Matrix jacobian,
const Index atmgeom_checked,
const Index atmfields_checked,
const String iy_unit_radar,
const ArrayOfString iy_aux_vars,
const Index stokes_dim,
const Vector f_grid,
const Index cloudbox_on,
const Index cloudbox_checked,
const Matrix sensor_pos,
const Matrix sensor_los,
const Index sensor_checked,
const Index jacobian_do,
const ArrayOfRetrievalQuantity jacobian_quantities,
const Agenda iy_radar_agenda,
const ArrayOfArrayOfIndex instrument_pol_array,
const Vector range_bins,
const Numeric ze_tref,
const Numeric k2,
const Numeric dbze_min,
const Verbosity verbosity 
)

WORKSPACE METHOD: yRadar.

Replaces yCalc for radar/lidar calculations.

The output format for iy when simulating radars and lidars differs from the standard one, and yCalc can not be used for such simulations. This method works largely as yCalc, but is tailored to handle the output from iyRadarSingleScat. Note that iy_radar_agenda replaces iy_main_agenda*.

The method requires additional information about the sensor, regarding its recieving properties. First of all, recieved polarisation states are taken from instrument_pol_array. Note that this WSV allows to define several measured polarisations for each transmitted signal. For example, it is possible to simulate transmittance of V and measuring backsacttered V and H.

Secondly, the range averaging is described by range_bins. These bins can either be specified in altitude or two-way travel time. In both case, the edges of the range bins shall be specified. All data (including auxiliary variables) are returned as the average inside the bins. If a bin is totally outside the model atmosphere, NaN is returned.

The options for iy_unit_radar are: <br> "1" : Backscatter coefficient. Unit is 1/(m*sr). At zero <br> attenuation, this equals the scattering matrix value for <br> the backward direction. See further AUG. <br> "Ze" : Equivalent reflectivity. Unit is mm^6/m^3. Conversion <br> formula is given below. <br> "dBZe": 10*log10(Ze/Z0), where Z0 is 1 mm^6/m^3.

The conversion from backscatter coefficient to Ze is: <br> Ze = 1e18 * lambda^4 / (k2 * pi^5) * sum(sigma), where sum(sigma) = 4 * pi * b, and b is the backscatter coefficient.

The reference dielectric factor can either specified directly by the argument k2. For example, to mimic the CloudSat data, k2 shall be set to 0.75 (citaion needed). If k2 is set to be negative (which is defualt), k2 is calculated as: <br> k2 = abs( (n^2-1)/(n^2+2) )^2, where n is the refractive index of liquid water at temperature ze_tref* and the frequency of the radar, calculated by the MPM93 parameterization.

A lower limit for dBZe is applied (dbze_min). The main reason is to handle the fact that dBZe is not defined for Ze=0, and dBZe is set to the clip value when Ze < 10^(dbze_min/10).

Author
Patrick Eriksson
Parameters
[in,out]wsWorkspace
[out]yWS Output
[out]y_fWS Output
[out]y_polWS Output
[out]y_posWS Output
[out]y_losWS Output
[out]y_auxWS Output
[out]y_geoWS Output
[out]jacobianWS Output
[in]atmgeom_checkedWS Input
[in]atmfields_checkedWS Input
[in]iy_unit_radarWS Input
[in]iy_aux_varsWS Input
[in]stokes_dimWS Input
[in]f_gridWS Input
[in]cloudbox_onWS Input
[in]cloudbox_checkedWS Input
[in]sensor_posWS Input
[in]sensor_losWS Input
[in]sensor_checkedWS Input
[in]jacobian_doWS Input
[in]jacobian_quantitiesWS Input
[in]iy_radar_agendaWS Input
[in]instrument_pol_arrayWS Input
[in]range_binsWS Input
[in]ze_trefGeneric Input (Default: "273.15")
[in]k2Generic Input (Default: "-1")
[in]dbze_minGeneric Input (Default: "-99")

Definition at line 556 of file m_cloudradar.cc.

References a, ARTS_USER_ERROR, ARTS_USER_ERROR_IF, b, chk_if_in_range(), chk_if_increasing(), ConstVectorView::empty(), Ppath::end_lstep, error_if_limb_ppath(), FOR_ANALYTICAL_JACOBIANS_DO, integration_bin_by_vecmult(), is_increasing(), iy_radar_agendaExecute(), jac_ranges_indices(), joker, Ppath::lstep, max(), min(), mult(), ConstMatrixView::ncols(), Array< base >::nelem(), ConstVectorView::nelem(), Ppath::ngroup, Ppath::np, ConstMatrixView::nrows(), Ppath::pos, pow(), Vector::resize(), Matrix::resize(), SPEED_OF_LIGHT, stokes2pol(), and ze_cfac().

Referenced by yRadar_g().

Variable Documentation

◆ GFIELD3_DB_GRID

const Index GFIELD3_DB_GRID = 1

Definition at line 53 of file m_cloudradar.cc.

Referenced by particle_bulkpropRadarOnionPeeling().

◆ GFIELD3_T_GRID

const Index GFIELD3_T_GRID = 2

Definition at line 54 of file m_cloudradar.cc.

Referenced by particle_bulkpropRadarOnionPeeling().

◆ LOG10_EULER_NUMBER

constexpr Numeric LOG10_EULER_NUMBER =Constant::log10_euler
inlineconstexpr

Definition at line 50 of file m_cloudradar.cc.

Referenced by particle_bulkpropRadarOnionPeeling().

◆ PI

constexpr Numeric PI =Constant::pi
inlineconstexpr

Definition at line 48 of file m_cloudradar.cc.

◆ SPEED_OF_LIGHT

constexpr Numeric SPEED_OF_LIGHT =Constant::speed_of_light
inlineconstexpr

Definition at line 49 of file m_cloudradar.cc.

Referenced by yRadar().