[mesa-users] calling kap_get_op_mono

[Robert FARMER]

Hi
Have you looked in star/private/opacities.f ?
This file calls kap_get_op_mono
There we see the x array is called as xa, back tracking that,
kap_get_op_mono is called by get1_kap which is called by do_kap_for_cell

There we find that xa is set as
xa(1:s% species) => s% xa(1:s% species,k)
Which is the abundances for that cell

Rob


On Tue, Nov 25, 2014 at 10:15 AM, Ehsan Moravveji <e.moravveji at gmail.com>
wrote:

> 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.
>
>
>
> ------------------------------------------------------------------------------
> Download BIRT iHub F-Type - The Free Enterprise-Grade BIRT Server
> from Actuate! Instantly Supercharge Your Business Reports and Dashboards
> with Interactivity, Sharing, Native Excel Exports, App Integration & more
> Get technology previously reserved for billion-dollar corporations, FREE
>
> http://pubads.g.doubleclick.net/gampad/clk?id=157005751&iu=/4140/ostg.clktrk
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20141125/b89b4ff1/attachment.html>