[mesa-users] calling kap_get_op_mono

[Ehsan Moravveji]

Hello everybody,
I would like to play around with the public routine kap_get_op_mono() in my run_star_extras to compare the resulting opacity profile with that of e.g. OPAL type 1 tables. However, I have difficulty passing the appropriate arguments. So, let’s see first how kap_get_op_mono() is defined in mesa/kap/public/kap_lib.f:

      subroutine kap_get_op_mono( &
            handle, zbar, log10_rho, log10_T, &
            lnfree_e, d_lnfree_e_dlnRho, d_lnfree_e_dlnT, &
            ! args for op_mono
            use_op_mono_alt_get_kap, &
            nel, izzp, fap, screening, umesh, ff, rs, ss, &                  <<---
            ! output
            kap, dlnkap_dlnRho, dlnkap_dlnT, ierr)

For me understanding the following line is uneasy, and I would appreciate any help on how to initialize them properly:

            nel, izzp, fap, screening, umesh, ff, rs, ss

Additionally, I realized another public routine called get_op_mono_args(), which returns “nel”, but there I cannot figure out what the input array “X” is.
Is it hydrogen mass fraction? Or perhaps I have to loop over the present species in the net, and call get_op_mono_args() for each of them?

I would appreciate if one of the “stewards” close to the “cockpit” give me some directions ;-)

Kind regards
Ehsan.