[Mesa-users] CO enhancement opacity table

[孙天成]

Dear Aaron,
       Thank you for your help!
       I find a way to add a new opacity table in the folder mesa/kap/preprocessor/, which was introduced in the README file. I imported my opacity table into mesa/kap/preprocessor/kap_input_data, and through the command: rebulid_all, I can change the newly added opacity table to the format that mesa expects and import it into mesa/data/kap_data.
       After testing, my opacity table can be successfully read and run by mesa without errors.
       The method you mentioned gave me a better understanding of the internal format of Mesa's opacity table. Thanks again!
Tiancheng
| |
孙天成
|
|
isaac0619 at 163.com
|
签名由网易邮箱大师定制


On 09/9/2020 22:04，Aaron Dotter<aaron.dotter at gmail.com> wrote：
Hi Tiancheng,
 
       The following attachments are my inlist, an example of the high-temperature opacity table, and the evolution trajectory from the calculation.

       I want to know how to use my opacity table correctly, thank you very much for your help.


I've looked at your inlists and the opacity table you provided.  You sent a pure helium (Y=1.0) table so I can't tell much about your new abundance pattern from that.


However, the format of your table appears to be different from what MESA is expecting.  Please compare the header of your table to one of those in mesa/data/kap_data/.  There are blank lines in the MESA file headers that are not in your table.  There are also blank fields in the high-T, high-R corner of your table that MESA might misinterpret.


It appears in your inlist that you are also using a different set of low-T opacity tables that correspond to your new high-T tables.  You should also check those to make sure they are in the correct format.



Unfortunately it IS possible for the opacity table file to be in the wrong format and yet not cause an error when the file is read in by MESA.


The file mesa/kap/preprocessor/src/fixed_metal.f includes a subroutine called Do_Table that shows how to correctly format the table (this is where the MESA fixed-metal opacity tables are created).


I hope this is helpful!


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