[Mesa-users] How to extract values

[Rob Farmer]

You want to pass the string "burn_c" not a variable called burn_c.

Rob


On Thu, 14 Jul 2022 at 18:09, harish sorrout <harishsorrout at gmail.com>
wrote:

> 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/ed82cfb5/attachment.htm>