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

Francis Timmes fxt44 at mac.com
Sat Feb 10 16:37:50 EST 2018


be sure to take a look at all the mlt routines in the other directory

% ls $MESA_DIR/star/other/*mlt*

fxt




> On Feb 8, 2018, at 9:03 PM, Yixiao Zhou <peter.zhou at anu.edu.au> wrote:
> 
> 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




More information about the Mesa-users mailing list