[Mesa-users] How to extract values
Rob Farmer
robert.j.farmer37 at gmail.com
Wed Jul 13 11:59:57 UTC 2022
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/20220713/8318175b/attachment.htm>
More information about the Mesa-users
mailing list