Ozone

class retrieval.species_o3.OzoneLimbProfile(atmosphere)

Handles the retrieval specifics for the ozone species including configuration of retrieval method, measurement vector and jacobian calculations and updating of atmoshperic state with new ozone profile.

The ozone measurement vector are triplets normalized by high altitude to help decrease senstivity to albedo and instrument calibration.

Methods

check_convergence
configure_instrument_and_method
get_required_wavelengths
get_retrieval_alts
get_x
lm_weighting_matrix
mart_weighting_matrix
measurement_vector
measurement_vector_jacobian
measurement_vector_jacobian_num
measurement_vector_jacobian_wf
output_data
output_data_bounded
set_x
update_profile
configure_instrument_and_method(ret_inst, meas_inst, method)

Sets up the retrieval method, and retrieval instrument. Uses wavelengths from the measurement instrument.

get_required_wavelengths()

Calculates the suggested wavelengths for the retrieval

Returns:

wavel : numpy array

Wavelengths suggested for the retrieval

lm_weighting_matrix(inst, y_meas)

Compute the relative weights of the measurements for least squares fitting.

mart_weighting_matrix(inst, y_meas)

Compute the relative weights of the measurements for MART fitting.

measurement_vector(inst, do_error=False)

Compute the ozone measurement vector

Parameters:

inst : instrument class

Instance of an instrument, either real or simulated that contains radiances.

measurement_vector_jacobian_num(inst)

Computes the measurement vector Jacobian through numerical perturbations. This is both slower and less accurate (assuming single scatter) than the analytic method, but is included incase the radiative transfer model does not support analytic jacobians.

measurement_vector_jacobian_wf(inst)

Computes the measurement vector Jacobian using weighting functions from the HR model

update_profile(delta, mode='additive')

Updates the ozone profile used by the radiative transfer model based on recommended step size, delta. Additional constraints can be imposed on delta, such as only using certain altitudes, what to do above/below retrieval limits, and maximum step size.

Parameters:

delta : 1d array

An array of numbers on the retrieval grid (defined by the jacobian calculation) in units of particles per cm^3.