[mesa-users] Mesa sound speeds

[Robert FARMER]

Hi
What version of MESA are you using?

What values are you getting and what values do you expect to get for the
csound (plots are helpful here)? Have you checked the units are the same?
Can you show a plot of what you think is diverging and what you expect the
data to show?

Rob

On Thu, Jul 30, 2015 at 2:42 PM, Farag, Ebraheem Khaled <farag.11 at lanl.gov>
wrote:

>  Hello, I'm running the mesa example_solar_model in the test_suite file
> and my sound speed plots are more than an order of magnitude greater than
> the data I have from AGS09. Is there something specific I'm doing wrong? o
> haven't changed anything, but the data looks completely divergent in the
> csound department. This doesn't seem correct because the mesa power point
> files indicate they closely match the sun helioseismology inferences, yet
> I'm not seeing that data when i run the simulation. could someone help. the
> model runs diverges regardless of opacities. I've change them to check and
> they just diverge in the same manner. especially at r>0.5.
>
> thanks,
> Ebraheem
>  ------------------------------
> *From:* Ken Shen [kenshen at astro.berkeley.edu]
> *Sent:* Thursday, July 30, 2015 12:21 PM
> *To:* Josiah Schwab
> *Cc:* mesa-users at lists.sourceforge.net
> *Subject:* Re: [mesa-users] Trouble relaxing composition
>
>   That worked perfectly, thanks!
>
> On Thu, Jul 30, 2015 at 12:37 AM, Josiah Schwab <jwschwab at berkeley.edu>
> wrote:
>
>> Hi Ken,
>>
>> > The issue is that the run relaxes to a composition that is not what I've
>> > specified.  In particular, the inner 98% or so should be nearly pure
>> 4He,
>> > but instead relaxes to exactly 50/50 4He/28Si.  Similar issues occur
>> with
>> > other compositions I've tried; in one of them, the inner region relaxed
>> to
>> > exactly 1/3 each of C/O/Si instead of 50/50 C/O.  Any insight out there?
>>
>> Pretty sure that your composition file is backwards.  That is, the first
>> column (xq) should be monotonically increasing.
>>
>> If I reverse the lines (i.e., select region, M-x reverse-region), I get
>> the the attached abundance profile.
>>
>> Josiah
>>
>> P.S.  Also, you should probably have a point at xq = 1.0.
>>
>>
>
>
> ------------------------------------------------------------------------------
>
> _______________________________________________
> 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/20150730/96bed50d/attachment.html>