[Mesa-users] star pointer for bolometric corrections

[Jeremy Sakstein]

Thanks Rob,

I think I figured out the problem. The initial model didn't relax so this
is a symptom rather than a cause.

On Mon, Apr 13, 2020 at 1:45 PM Rob Farmer <r.j.farmer at uva.nl> wrote:

> Well you've tried accessing s% x(0) when the array starts counting at 1.
>
> Without seeing any code i can't say why this is happening.
>
> Rob
>
> On Mon, 13 Apr 2020 at 19:28, Jeremy Sakstein <sakstein at hawaii.edu> wrote:
>
>> Thanks Rob,
>>
>> I managed to get it to spit out sensible values but occasionally I get an
>> error:
>>
>> "Fortran runtime error: Index '0' of dimension 1 of array 's%x' below
>> lower bound of 1"
>>
>> What does this mean?
>>
>> Cheers,
>>
>> Jeremy
>>
>>
>> On Mon, Apr 13, 2020 at 10:35 AM Rob Farmer <r.j.farmer at uva.nl> wrote:
>> >
>> > Hi,
>> >
>> > You can use the function get_bc_by_name from colors_lib which takes the
>> filter name as an argument. Where the filter name is the column header for
>> the bolometric correction table you are using.
>> >
>> > Rob
>> >
>> > On Mon, 13 Apr 2020 at 16:25, Jeremy Sakstein <sakstein at hawaii.edu>
>> wrote:
>> >>
>> >> Thanks Josiah,
>> >>
>> >> What I really want to do is modify run_star_extras to spit out the BC
>> at some arbitrary point (actually as soon as RSP has loaded up).
>> >>
>> >> I went through where you suggested and I think the routine is:
>> >>
>> >> get_bc_by_id(i, safe_log10(s% photosphere_T), s% photosphere_logg,
>> m_div_h, ierr)
>> >>
>> >> where c = s% history_column_spec(j) and i = c - log_lum_band_offset.
>> Is this correct?
>> >>
>> >> How do I decide what j is to get the band I want?
>> >>
>> >> Cheers,
>> >>
>> >> Jeremy
>> >>
>> >>
>> >>
>> >>
>> >>
>> >> On Mon, Apr 6, 2020 at 6:10 PM Josiah Schwab <jwschwab at ucsc.edu>
>> wrote:
>> >> >
>> >> > Hi Jeremy,
>> >> >
>> >> >>
>> >> >> I can see from history_columns.list that it is possible to write
>> bolometric corrections straight from MESA.
>> >> >>
>> >> >> Can someone tell me what the star pointers for these are so that I
>> can use them in intermediate steps?
>> >> >
>> >> >
>> >> > I don't understand this question, but you can find the code that
>> adds this information to the history output via calls to the colors module
>> in star/private/history.f90.  (It's around line 1230 in r12778.)
>> >> >
>> >> > Josiah
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20200413/20752a34/attachment.htm>