ARTS built-in documentation server

Workspace Method DisortCalcWithARTSSurface


Interface to the DISORT scattering solver (by Stamnes et al.).

As DisortCalc but uses surface_rtprop_agenda.

The Lambertian surface reflection is set by surface_rtprop_agenda.
If the GIN inc_angle is inside of the range [0,90], the reflection is
set according to the result of surface_rtprop_agenda for this incidence
angle. Otherwise (default) is to call surface_rtprop_agenda for
multiple angles, to estimate the hemispheric mean value.

Some auxiliary quantities can be obtained. Auxiliary
quantities are selected by disort_aux_vars and returned by disort_aux.
Valid choices for auxiliary data are:
 "Layer optical thickness": Matrix [f_grid, size of p_grid - 1] layer optical                              thickness.
 "Single scattering albedo": Matrix [f_grid, size of p_grid - 1] layer single"
                               scattering albedo.
 "Direct beam": Matrix [f_grid, p_grid]. Attenuated direct at level.
                               Zero, if no sun is present 

Authors: Claudia Emde, Jana Mendrok, Manfred Brath


DisortCalcWithARTSSurface( cloudbox_field, disort_aux, atmfields_checked, atmgeom_checked, scat_data_checked, cloudbox_checked, cloudbox_on, cloudbox_limits, propmat_clearsky_agenda, surface_rtprop_agenda, gas_scattering_agenda, atmosphere_dim, pnd_field, t_field, z_field, vmr_field, p_grid, lat_true, lon_true, refellipsoid, scat_data, suns, f_grid, za_grid, aa_grid, stokes_dim, z_surface, gas_scattering_do, suns_do, disort_aux_vars, nstreams, Npfct, only_tro, quiet, emission, intensity_correction, inc_angle )


OUTcloudbox_field(Tensor7)The spectral radiance field inside the cloudbx.
OUTdisort_aux(ArrayOfMatrix)Auxilary data to the output of the DisortCalc-Methods.
INatmfields_checked(Index)OK-flag for atmospheric grids and (physical) fields.
INatmgeom_checked(Index)OK-flag for the geometry of the model atmosphere.
INscat_data_checked(Index)OK-flag for scat_data.
INcloudbox_checked(Index)OK-flag for variables associated with the cloudbox.
INcloudbox_on(Index)Flag to activate the cloud box.
INcloudbox_limits(ArrayOfIndex)The limits of the cloud box.
INpropmat_clearsky_agenda(Agenda)Agenda calculating the absorption coefficient matrices.
INsurface_rtprop_agenda(Agenda)Agenda providing radiative properties of the surface.
INgas_scattering_agenda(Agenda)Agenda calculating gas scattering extinction and phase matrix.
INatmosphere_dim(Index)The atmospheric dimensionality (1-3).
INpnd_field(Tensor4)Particle number density field.
INt_field(Tensor3)The field of atmospheric temperatures.
INz_field(Tensor3)The field of geometrical altitudes.
INvmr_field(Tensor4)VMR field.
INp_grid(Vector)The pressure grid.
INlat_true(Vector)Latitudinal geolocation for 1D and 2D data.
INlon_true(Vector)Longitudinal geolocation for 1D and 2D data.
INrefellipsoid(Vector)Reference ellipsoid.
INscat_data(ArrayOfArrayOfSingleScatteringData)Array of single scattering data.
INsuns(ArrayOfSun)Array of Sun.
INf_grid(Vector)The frequency grid for monochromatic pencil beam calculations.
INza_grid(Vector)Zenith angle grid.
INaa_grid(Vector)Azimuthal angle grid.
INstokes_dim(Index)The dimensionality of the Stokes vector (1-4).
INz_surface(Matrix)The surface altitude.
INgas_scattering_do(Index)Flag to activate gas scattering.
INsuns_do(Index)Flag to activate the sun(s).
INdisort_aux_vars(ArrayOfString)Selection of quantities for disort_aux .
GINnstreams(Index, Default: 8)Number of polar angle directions (streams) in DISORT solution (must be an even number).
GINNpfct(Index, Default: 181)Number of angular grid points to calculate bulk phase function on (and derive Legendre polynomials from). If <0, the finest za_grid from scat_data will be used.
GINonly_tro(Index, Default: 0)Set to 1 if the scattering data is just of TRO type. Has effect only if Npfct > 3 or Npfct<0, but then leads to much faster calculations.
GINquiet(Index, Default: 0)Silence C Disort warnings.
GINemission(Index, Default: 1)Enables blackbody emission. Set to zero, if no Emission e. g. like in visible regime for earth is needed
GINintensity_correction(Index, Default: 1)Enables intensity correction. Importantant for low number of streams. Set to zero, if problems encounter or using a high number of streams (>30)
GINinc_angle(Numeric, Default: -1)Incidence angle, see above.