[Mesa-users] Varying radiative opacity

[Evan Bauer]

Hi Jeremy,

I believe your suspicion is right, and the opacity_factor control will scale the net opacity (including conduction), so you can’t do what you want with an inlist option.

In the most recent MESA release, it would be quite straightforward to implement what you want in your run_star_extras using the kap/other/other_radiative_opacity.f90 hook. 

The hooks back in r12778 don’t give you quite as direct of access, but I imagine you could find a roundabout way of modifying just the radiative opacity. E.g. you could use the star/other/other_kap.f90 hook along with the kap/public function kap_get_elect_cond_opacity. That would allow you to subtract out the conductive opacity, modify what’s left by the factor you want, and then add the conductive opacity back in when you’re done. Does that make sense?

Cheers,
Evan

> On Mar 29, 2022, at 3:58 PM, Jeremy Sakstein <sakstein at hawaii.edu> wrote:
> 
> Hi everyone,
> 
> I'm trying to vary the radiative opacity leaving the conduction opacity fixed in r12778.
> 
> Does anyone have any suggestions for the best way to do this?
> 
> I noticed the option opacity_factor in the controls but I am unsure what this scales. Currently, I have include_electron_conduction = .true. so my suspicion is that it will multiply the sum of the conduction and radiative opacities.
> 
> Is there any way to just scale the radiative opacities?
> 
> Cheers,
> 
> Jeremy
> 
> -- 
> Jeremy Sakstein (He/Him/His)
> Assistant Professor of Theoretical Physics
> Department of Physics & Astronomy
> University of Hawaiʻi at Mānoa
> Watanabe Hall
> 2505 Correa Road
> Honolulu, HI, 96822, USA
> http://www.jeremysakstein.com <http://www.jeremysakstein.com/>
> _______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20220330/f43b4836/attachment.htm>