[mesa-users] calling kap_get_op_mono

[Bill Paxton]

Hi Ehsan,

Did you check the mesa/kap/test/src files?  Look for a call on 'op_mono_get_kap'.   Perhaps you can use that to figure out how to use it.  

Since we don't distribute the large OP_mono database with mesa, it doesn't get tested as part of the installation.  So the test case in kap may have gotten "software rot" from lack of use.  If so, let me know and we'll fix it.   The OP mono opacities do get tested regularly as part of the mesa/star test suite (in the radiative_levitation test case), but that only is run "on demand", not as part of the usual mesa installation.    If you've run that test, you'll know that the mono opacities are an expensive luxury -- it takes lots of processing to get all of those separate opacities.

Cheers,
Bill

p.s.  I just saw Rob's response.  (he beat me -- good work Rob!)   



On Nov 25, 2014, at 9:15 AM, Ehsan Moravveji 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