Utility

radtran.climatologies.utility.extinction_to_number_density(heights_m, kext, wavel, optprop, neutral=None, lat=0, lon=0, mjd=54372)

Converts a (height,extinction) profile to (height,numberdensity)

Parameters:

heights_m : numpy array

Heights in m for the extinction profile

kext : numpy array

Same shape as heights_m, extinction in per cm

wavel : scalar

Wavelength in which the extinction is specified

optprop : ISKOpticalProperty

Optical property to use to calculate cross sections

neutral : ISKClimatology (optional)

Climatology to use for the neutral background, if not specified ECMWF is used

lat : scalar, optional

Latitude, defaults to 0

lon : scalar, optional

Longitude, defaults to 0

mjd : scalar, optional

modified julian date, defaults to 54372

Returns:

n : numpy array

Number density specified at heights_m in /cm3

radtran.climatologies.utility.height_profile_from_clim(climatology, guid, heights, lat, lon, mjd)

Get the height profile from a climatology at the specified location

Parameters:

climatology: ISKClimatology object

guid: ISKGuid object

heights: array

an array of heights in meters

lat: double

latitude

lon: double

longitude

mjd: double

modified julian date

Returns:

profile: numpy.1darray

array the same size as heights containing the profile values.

radtran.climatologies.utility.number_density_to_extinction(heights_m, n_cm3, wavel, optprop, neutral=None, lat=0, lon=0, mjd=54372)

Converts a (height,number density) profile to (height,extinction)

Parameters:

heights_m : numpy array

Heights in m for the extinction profile

n_cm3 : numpy array

Same shape as heights_m, aerosol number density in particles per cm^3

wavel : scalar

Wavelength in which the extinction is specified

optprop : ISKOpticalProperty

Optical property to use to calculate cross sections

neutral : ISKClimatology (optional)

Climatology to use for the neutral background, if not specified ECMWF is used

lat : scalar, optional

Latitude, defaults to 0

lon : scalar, optional

Longitude, defaults to 0

mjd : scalar, optional

modified julian date, defaults to 54372

Returns:

k : numpy array

extinction specified in per km