ARTS built-in documentation server

Workspace Method OptimizeDoitPressureGrid


Optimization of the pressure grid for RT calculation.
The methods consists of three parts:
1) Calculate the single scattering albedo and the scattering opticalthickness 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 .

Authors: Jakob Doerr


OptimizeDoitPressureGrid( p_grid, pnd_field, t_field, scat_data_mono, z_field, cloudbox_limits, cloudbox_field_mono, pha_mat_doit, vmr_field, p_grid_orig, f_grid, f_index, propmat_clearsky_agenda, tau_scat_max, sgl_alb_max, cloudbox_size_max )


OUT+INp_grid(Vector)The pressure grid.
OUT+INpnd_field(Tensor4)Particle number density field.
OUT+INt_field(Tensor3)The field of atmospheric temperatures.
OUT+INscat_data_mono(ArrayOfArrayOfSingleScatteringData)Monochromatic single scattering data.
OUT+INz_field(Tensor3)The field of geometrical altitudes.
OUT+INcloudbox_limits(ArrayOfIndex)The limits of the cloud box.
OUT+INcloudbox_field_mono(Tensor6)Monochromatic radiation field inside the cloudbox.
OUT+INpha_mat_doit(Tensor7)Ensemble averaged phase matrix for DOIT calculation.
OUT+INvmr_field(Tensor4)VMR field.
OUTp_grid_orig(Vector)The original pressure grid before optimization.
INf_grid(Vector)The frequency grid for monochromatic pencil beam calculations.
INf_index(Index)Frequency index.
INpropmat_clearsky_agenda(Agenda)Agenda calculating the absorption coefficient matrices.
GINtau_scat_max(Numeric, Default: 0.1)Maximum scattering optical thickness
GINsgl_alb_max(Numeric, Default: 0.9)Maximum single scattering albedo
GINcloudbox_size_max(Index, Default: 200)Maximum cloudbox size