[Mesa-users] Unable to add extra heating to stars

[Thomas Steindl]

Hi Shafaat, Hi Warrick, 

With the change to autodiff, is there a need to set the partials at all?
And if you are setting the partial derivative to zero, is there at all a need to use the other_energy_implicit instead of the other_energy routine?

Best,
Thomas


__________________________________________________________________

Thomas Steindl, MSc
University of Innsbruck | Institute for Astro- and Particle Physics
Tel: +43 512507 52031
Mail: thomas.steindl at uibk.ac.at
Homepage: www.thomassteindl.com

> Am 06.12.2022 um 18:18 schrieb Warrick Ball <W.H.Ball at bham.ac.uk>:
> 
> Hi Shafaat,
> 
> I think this is closely related to the question you asked about a month ago [1] to which I replied [2] explaining that the `extra_heat` variable is now of a type from the automatic differentiation module.  I don't have a complete solution but in principle think what you need to do is assign the derivatives to things like `s% extra_heat(k)% d1Array(j)` where `j` is the index of the relevant variable.
> 
> The indices are defined from line 259 onward in `star_data/public/star_data_def.inc` and you can use those index names (you might need `use star_data_def` in your `run_star_extras.f90`). e.g., instead of
> 
>    s% d_extra_heat_dlnR00(k) = 0d0
> 
> I think you'd set
> 
>    s% extra_heat% d1Array(i_lnR_00) = 0d0
> 
> I stand to be corrected though, since I haven't used these `auto_diff` types myself.
> 
> That said, all the derivatives you're getting errors on so far are being set to zero, in which case I think you can just use something like
> 
>    s% extra_heat(k)% d1Array(:) = 0
> 
> This is a good question though, since I don't think this is adequately documented.  If nothing else, I'll update the module documentation for `auto_diff` to point to where the indices of the `auto_diff_real_star_order1` are found.
> 
> Cheers,
> Warrick
> 
> [1] https://lists.mesastar.org/pipermail/mesa-users/2022-November/014153.html
> [2] https://lists.mesastar.org/pipermail/mesa-users/2022-November/014154.html
> 
> ___________
> 
> Warrick Ball
> Postdoc, School of Physics and Astronomy
> University of Birmingham, Edgbaston, Birmingham B15 2TT
> W.H.Ball at bham.ac.uk
> +44 (0)121 414 4552
> 
> On Tue, 6 Dec 2022, mesa-users at lists.mesastar.org wrote:
> 
>> Hello,
>> I am trying to adapt codes from a paper and add extra heating to my stellar evolution model. However, whenever I try to make the file (./mk), I run into the following problem:
>> ➜  work ./mk
>> gfortran -Wno-uninitialized -fno-range-check -fmax-errors=7  -fprotect-parens -fno-sign-zero -fbacktrace -ggdb -finit-real=snan -fopenmp -fbounds-check -Wuninitialized -Warray-bounds -ggdb -ffree-form -ffree-line-length-none -x f95-cpp-input
>> -I/Users/jareddiks/Downloads/mesa-r21.12.1/include -I../src -c ../src/run_star_extras.f90
>> ../src/run_star_extras.f90:103:36:
>>   103 |              s% d_extra_heat_dlnR00(k) = 0d0
>>       |                                    1
>> Error: 'd_extra_heat_dlnr00' at (1) is not a member of the 'star_info' structure; did you mean 'extra_heat'?
>> ../src/run_star_extras.f90:104:36:
>>   104 |              s% d_extra_heat_dlnRp1(k) = 0.
>>       |                                    1
>> Error: 'd_extra_heat_dlnrp1' at (1) is not a member of the 'star_info' structure; did you mean 'extra_heat'?
>> ../src/run_star_extras.f90:105:36:
>>   105 |              s% d_extra_heat_dlndm1(k) = 0.
>>       |                                    1
>> Error: 'd_extra_heat_dlndm1' at (1) is not a member of the 'star_info' structure; did you mean 'extra_heat'?
>> ../src/run_star_extras.f90:106:36:
>>   106 |              s% d_extra_heat_dlnd00(k) = 0.
>>       |                                    1
>> Error: 'd_extra_heat_dlnd00' at (1) is not a member of the 'star_info' structure; did you mean 'extra_heat'?
>> ../src/run_star_extras.f90:107:36:
>>   107 |              s% d_extra_heat_dlndp1(k) = 0.
>>       |                                    1
>> Error: 'd_extra_heat_dlndp1' at (1) is not a member of the 'star_info' structure; did you mean 'extra_heat'?
>> ../src/run_star_extras.f90:109:36:
>>   109 |              s% d_extra_heat_dlnTm1(k) = 0.
>>       |                                    1
>> Error: 'd_extra_heat_dlntm1' at (1) is not a member of the 'star_info' structure; did you mean 'extra_heat'?
>> ../src/run_star_extras.f90:110:36:
>>   110 |              s% d_extra_heat_dlnT00(k) = 0.
>>       |                                    1
>> Error: 'd_extra_heat_dlnt00' at (1) is not a member of the 'star_info' structure; did you mean 'extra_heat'?
>> compilation terminated due to -fmax-errors=7.
>> make: *** [run_star_extras.o] Error 1
>> FAILED
>> I'm guessing the following problem may be due to a version error. I've looked into the MESA change log but I'm not sure how to get rid of this error. 
>> Is there any way I might just find a workaround/avoid it?  That would be really helpful. 
>> Thank you!
>> Regards,
>> Shafaat
>> 
> _______________________________________________
> 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/20221206/1f207f91/attachment.htm>