Find the corresponding species (if not in the same order)

Interpolation of ref/test data on a common altitude grid

Add nreactants to the type one_reaction

Subprogram m_reac::old_reac_rrate (rtype, ncoeff, coeff, T, P, T_o, N_o, k_r)

Find a better generic way to treat polynomials with temperature ranges

Subprogram m_reac::reac_fillmatrix (sp, re, nr, nu, T, P)

This procedure must be made only once. Create a structure and a subroutine that associates directly i,j,m,k in order to make only a multiplication!

Subprogram m_reac::reac_readreac (wkdir, reacfile, sp, re)

What happens when ions have multiple charges? e.g. NO^++ ?

Subprogram m_spec::spec_readspec (wkdir, specfile, s)

Use a generic function or subroutine to select the catalytic species.

Subprogram module_chimie::compute_rrate (rtype, ncoeff, coeff, T, P, k_r)

Find a better generic way to treat polynomials with temperature ranges

Jupiter's magnetic north pole has to be reversed (lat < 0)