[Mesa-users] Modify radiative temperature gradient "gradr"

[Yixiao Zhou]

Hi Francis,


Thank you very much for the fast reply, it is very helpful! I didn't realize that "gradr" is actually included in the "mlt_basics" variable in other_mlt.f90


So what I should do is write my own MLT routine to change gradr and let it overwrite the "mlt_eval" function in MESA, then the new gradr will be used to determine the region of convection zone through the Schwarzschild criterion. Am I correct?


Thanks again,

Yixiao Zhou

________________________________
From: Francis Timmes <fxt44 at mac.com>
Sent: Friday, 9 February 2018 5:56:09 AM
To: Yixiao Zhou
Cc: Frank Timmes; mesa-users at lists.mesastar.org
Subject: Re: [Mesa-users] Modify radiative temperature gradient "gradr"

> can this be done in run_star_extras.f?

yes.

the star_info structure contains the stellar model and the
components are listed in the file $MESA_DIR/star/public/star_data.inc ,
which includes gradr. do read the header of this file where its pointed
out that gradr is defined at cell boundaries and depend on its neighbors.

$MESA_DIR/star/other/other_mlt.f90 contains the hook for the physics or numerics
that you might want to alter.

if any of this is mysterious, see josiah schwab's excellent 2017 mesa summer school
material on "beyond inlists", http://cococubed.asu.edu/mesa_market/education.html

fxt







> On Feb 8, 2018, at 4:23 AM, Yixiao Zhou via Mesa-users <mesa-users at lists.mesastar.org> wrote:
>
> Dear MESA users,
>
> I'm looking for a method to adjust the "gradr" value in every step:
> gradr_new = gradr_MESA * f_corr,
> where the "correction function" f_corr depends on radius.
>
> This modification should be done before MESA enters into MLT section so that I can determine convective boundaries and other convective-related quantities using the new gradr data:
> Step (i) begin ... --> MESA calculate gradr --> modify gradr --> use the new value for MLT part --> ... step (i) end goto step (i+1)
>
> I am wondering if there is a way to achieve this without changing the source code? For example, can this be done in run_star_extras.f?
>
> Many thanks in advance,
> Yixiao Zhou
> PhD student
> Australian National University
> _______________________________________________
> 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/20180209/ce2b2e00/attachment.html>