ARTS built-in documentation server

Workspace Method yRadar

Description

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

The conversion from backscatter coefficient to Ze is:
   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:
   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).

Authors: Patrick Eriksson

Synopsis

yRadar( y, y_f, y_pol, y_pos, y_los, y_aux, y_geo, jacobian, atmgeom_checked, atmfields_checked, iy_unit_radar, iy_aux_vars, stokes_dim, f_grid, cloudbox_on, cloudbox_checked, sensor_pos, sensor_los, sensor_checked, jacobian_do, jacobian_quantities, iy_radar_agenda, instrument_pol_array, range_bins, ze_tref, k2, dbze_min )

Variables

OUTy(Vector)The measurement vector.
OUTy_f(Vector)The frequencies associated with y.
OUTy_pol(ArrayOfIndex)The polarisation states associated with y.
OUTy_pos(Matrix)The sensor positions associated with y.
OUTy_los(Matrix)The line-of-sights associated with y.
OUTy_aux(ArrayOfVector)Data auxilary to y.
OUTy_geo(Matrix)The geo-positioning associated with y.
OUTjacobian(Matrix)The Jacobian matrix.
INatmgeom_checked(Index)OK-flag for the geometry of the model atmosphere.
INatmfields_checked(Index)OK-flag for atmospheric grids and (physical) fields.
INiy_unit_radar(String)Unit for radar simulations.
INiy_aux_vars(ArrayOfString)Selection of quantities for iy_aux and when applicable also y_aux.
INstokes_dim(Index)The dimensionality of the Stokes vector (1-4).
INf_grid(Vector)The frequency grid for monochromatic pencil beam calculations.
INcloudbox_on(Index)Flag to activate the cloud box.
INcloudbox_checked(Index)OK-flag for variables associated with the cloudbox.
INsensor_pos(Matrix)The sensor position for each measurement block.
INsensor_los(Matrix)The sensor line-of-sight (LOS) for each measurement block.
INsensor_checked(Index)OK-flag for sensor related variables.
INjacobian_do(Index)Flag to activate (clear-sky) Jacobian calculations.
INjacobian_quantities(ArrayOfRetrievalQuantity)The retrieval quantities in the Jacobian matrix.
INiy_radar_agenda(Agenda)Agenda calculating pointwise backscattering.
INinstrument_pol_array(ArrayOfArrayOfIndex)Multiple definition of instrument polarisation.
INrange_bins(Vector)The range bins of an active instrument.
GINze_tref(Numeric, Default: 273.15)Reference temperature for conversion to Ze.
GINk2(Numeric, Default: -1)Reference dielectric factor.
GINdbze_min(Numeric, Default: -99)Clip value for dBZe.