[Mesa-users] How to extract values

[harish sorrout]

Thank you sir,
As you suggested, I extracted the values in run_star_extras.f. and then
compiled using ./mk and it compiled. but when ran the simulation using ./rn.
I got following error :
reduce dt because of max decrease       1
 FATAL ERROR Bad value for profile name 0�{9�˖]�@D�72�P�7�!���
STOP star_get_profile_id

sir kindly suggest what I should do to resolve this error?
thankyou

 with regards
Harish Kumar


On Wed, Jul 13, 2022 at 5:30 PM Rob Farmer <robert.j.farmer37 at gmail.com>
wrote:

> So its a function that is in MESA_DIR/star/public/star_lib.f90.
>
> So first make sure you have
> use star_lib
> so you can call the function.
>
> It has a function arguments:
>
>       real(dp) function star_get_profile_output(s,name,k)
>          type (star_info), pointer :: s
>          character(len=*),intent(in) :: name
>          integer,intent(in) :: k
>
> So you just need to pass (s, the name of the column you want, the zone you
> want that value at) and you'll get back that value for the column you asked
> for.
>
> Rob
>
>
>
>
> On Wed, 13 Jul 2022 at 11:31, harish sorrout <harishsorrout at gmail.com>
> wrote:
>
>> thank you sir
>>    I have attached the file of run_star_extras.f. As expected, after ./mk
>> there is a error message
>>
>>   221 |       vals(k,1) = (s% cp(k)*s% T(k))/ burn_c
>>       |                                            1
>> Error: Symbol 'burn_c' at (1) has no IMPLICIT type. Although burn_c is in
>> profile column list.
>>
>> sir, how and where used 'function star_get_profile_output' in
>> run_star_extra.f file to get burn_c.
>> Please enlighten me sir.
>>
>> with regards
>> Harish Kumar
>>
>>
>> On Tue, Jul 12, 2022 at 11:56 PM Rob Farmer <robert.j.farmer37 at gmail.com>
>> wrote:
>>
>>> >Although I tried.
>>>
>>> Okay let's see what you have done so we can help you fix your code.
>>>
>>> Rob
>>>
>>>
>>> On Tue, 12 Jul 2022 at 18:46, harish sorrout <harishsorrout at gmail.com>
>>> wrote:
>>>
>>>> Thanks sir, I try to find out *burn_c* parameter which is given in
>>>> profile column list. I want
>>>>
>>>> T = Cp * T / burn_c ;  where burn_c is the energy generated due to
>>>> carbon burning ( erg/g/s )
>>>>
>>>>
>>>> To get  these variables, Cp = s%cp , T= s% T. But burn_c is not member
>>>> of star family. So I cannot use s%burn_c. As you suggested to use  function
>>>> star_get_profile_output() to extract profile variables.
>>>> Sir, sorry for sillyness but I don't know how to use this function.
>>>> Although I tried.
>>>> Sir can you please tell me this by giving one example.
>>>> I will be thankful to you.
>>>>
>>>>
>>>> With regards
>>>> Harish Kumar
>>>>
>>>>
>>>>
>>>> On Sun, 10 Jul 2022, 18:53 Rob Farmer, <robert.j.farmer37 at gmail.com>
>>>> wrote:
>>>>
>>>>> Hi,
>>>>>
>>>>> There is the function star_get_profile_output() in star_lib that can
>>>>> be used to get any profile output (and star_get_history_output() for
>>>>> history output)
>>>>>
>>>>> Rob
>>>>>
>>>>> On Sun, 10 Jul 2022 at 12:28, harish sorrout via Mesa-users <
>>>>> mesa-users at lists.mesastar.org> wrote:
>>>>>
>>>>>> Respected mesa users,
>>>>>>
>>>>>> I am trying to see the accretion effect on CO white dwarfs by using
>>>>>> MESA version: 12778. And
>>>>>> I want to calculate nuclear burning time scale only from carbon
>>>>>> burning where t_b = Cp T / q , where Cp: specific heat, T: temperature and
>>>>>> q: burning rates in run_star_extras.f.  To get burning time scale, I need
>>>>>> three variables from which two I have extract using star pointer s% Cp, s%
>>>>>> T. But I don't know how to extract third variable that is nuclear burning
>>>>>> rates ( burn_c and c12_c12 which are available in profile column list) due
>>>>>> to carbon burning.
>>>>>> Please advise me how can extract.
>>>>>>
>>>>>> I will be thankful to you
>>>>>>
>>>>>> With regards
>>>>>>
>>>>>> Harish Kumar
>>>>>>
>>>>>>
>>>>>>
>>>>>> _______________________________________________
>>>>>> 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/20220714/68ce6a28/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: run_star_extras.f
Type: text/x-fortran
Size: 11816 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20220714/68ce6a28/attachment.bin>