[mesa-users] Experimenting with adjust_mlt_gradT_fraction

[Ehsan Moravveji]

Dear all,
Special greetings to those of us who already started their morning in the exciting Summer School. 

I was thinking of a gedankenexperiment to modify the temperature gradient in the semi-convective layers of a 15Msun main sequence star. There is a wonderful variable already implemented to vary gradT between gradr and grad_ad, called "adjust_mlt_gradT_fraction" (hereafter called f in the plots and the following text). I have traced this variable in hydro_vars.f and mlt_info.f in /star/private. 
I took three values of f = 0.00, 0.25 and 0.50. Then, I use the other_mlt hook in my run_star_extras to change the value of f "only" in the semi-convective zones. Indeed, it works. Please see the attached gradT.png plot. 

However, the temperature and pressure profiles and all other hydro variables should be modified as well to be consistent with ad-hoc modifications in gradT. In other words, direct differentiation of temperature and pressure profiles must consistently result in gradT for any arbitrary value of f == adjust_mlt_gradT_fraction.
Currently, this is not the case. To convince myself, I did three other plots of T, P and Rho in the same semi-convective zone; see the attachments, please.

I totally understand that gradT and hydro_vars are evaluated across different modules, and that can make life complicated. But, I desire having self-consistent outcome in my profile files.
In case one likes to reproduce the results, I have attached my inlist and run_star_extras.f (compatible with v.6794).

I would be grateful if you enlighten me with your comments/feedbacks. 

Kind regards,
Ehsan.

-------------- next part --------------
A non-text attachment was scrubbed...
Name: gradT.png
Type: image/png
Size: 52332 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20140811/f72e7949/attachment.png>
-------------- next part --------------


-------------- next part --------------
A non-text attachment was scrubbed...
Name: Density.png
Type: image/png
Size: 54525 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20140811/f72e7949/attachment-0001.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Pressure.png
Type: image/png
Size: 54096 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20140811/f72e7949/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Temperature.png
Type: image/png
Size: 48388 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20140811/f72e7949/attachment-0003.png>
-------------- next part --------------

-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist_general
Type: application/octet-stream
Size: 9236 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20140811/f72e7949/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: run_star_extras.f
Type: application/octet-stream
Size: 12900 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20140811/f72e7949/attachment-0001.obj>
-------------- next part --------------