Data and code for "Vertical and temporal H3+ structure at the auroral footprint of Io"

Data and code for "Vertical and temporal H3+ structure at the auroral footprint of Io".

Data  used to obtain the Figure 2 is in ascii format and has the following columns:

Date
Time
File name in PDS
Average radiance in the ROI (mW m^-2 sr^-1)

Code can be executed in the MATLAB environment. The code can be executed as:

v=myfunc(D0,F0,t)

given the time t, the peak density of H3+ D0 at time t=0,  and the conversion factor from column density to radiance F0, gives the vertically-emitted radiance v (mW m^-2 sr^-1). Note that, since the step is 10 km and the density is in cm^-3, there is a factor 10^6 to take into account. To obtain radiances in mW m^-2 sr^-1 instead of W m^-2 sr^-1, another factor 10^3 has to be included.

For example, if one calls:

>> myfunc(366362,  38.26e-18,0)*1d6*1d3

     ans =      0.2993

The output is the radiance at time=0 for a peak density of 366362 cm^-3 and a conversion factor of 38.26e-18 (these are the exact output from the best-fit routine, and in the paper they have been rounded to 3.7E5 and 40E-18). 

Version with myfunc_1 is with variable destruction rate (k) but uniform emission rate (f); version myfunc_2 has also variable emission rate.

The result may need to be scaled depending on the units used for F0.