[Mesa-users] Fwd: How to extract values
harish sorrout
harishsorrout at gmail.com
Fri Jul 22 16:27:40 UTC 2022
Thank you Rob,
Now I have added the extra profile columns in run_star_extras.f90. By
using mesa 12778, I have accreted the carbon and oxygen on to CO white
dwarf with accretion rate 2.3d-6. Although white dwarf accreted 1.369 solar
mass and central temperature log 8.69 and then simulation terminate with
following errors:
num_newton_iterations dt_limit 1.1642857143 17 14
log eps nuc ratio 368
764 164 5.0531126467963694D-09 -5.0531126467963694D-09 1.
5160424637814009D+01 1.5160424814208914D+01
timestep_hold > model_number, so no timestep increase 165
164
WARNING: rel_run_E_err 164
4.1014717634964032D-02
retry 165
retry 165
retry 165
retry 165
retry 165
retry 165
retry 165
num_newton_iterations dt_limit 1.2357142857 18 14
log eps nuc ratio 362
775 165 -1.6962520136639194D-12 1.6962520136639194D-12 8.
3615846268716751D+00 8.3615846268390168D+00
timestep_hold > model_number, so no timestep increase 166
165
WARNING: rel_run_E_err 165
4.1014909100627141D-02
png/trho_000165.png/png
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
retry 166
dt
8.4894335288525781D-16
min_timestep_limit
1.0000000000000001D-15
stopping because of problems dt < min_timestep_limit
WARNING: rel_run_E_err 166
4.1014909100627141D-02
terminated evolution: cannot find acceptable model
termination code: min_timestep_limit
DATE: 2022-07-22
TIME: 21:46:20
******************************************************************
failed to create lgT_8.45.mod when running inlist_to_c_ignition_header
******************************************************************
SIr, I know it is difficult to further simulate upto log Tc = 9.5. but I
want to recreate the results of paper "The outcomes of carbon–oxygen white
dwarfs accreting CO-rich material " https://doi.org/10.1093/mnras/sty3176 in
which they achieved that temperature ranges. Also I have checked their
published inlist.
So, I am requesting that you and all mesa users please suggest what
parameters I should use in my inlist and what physics I am missing to
achieve that temperature range. Please do me a favour.
Thankyou
with regards
Harish Kumar
On Tue, Jul 19, 2022 at 8:11 PM Rob Farmer <robert.j.farmer37 at gmail.com>
wrote:
> > By doing so, I did not get any errors while compiling. Sir Is there
> anything wrong with the above command.
> Sir I could not understand your words "You want to pass the string
> "burn_c" not a variable called burn_c".
> Please enlighten me.
>
> Yes your still passing a character called burn_c and not the string
> 'burn_c' (with the quotes)
>
> E_cburn(nz) = star_get_profile_output(s,'burn_c',k)
>
> Rob
>
>
> On Fri, 15 Jul 2022 at 17:38, harish sorrout via Mesa-users <
> mesa-users at lists.mesastar.org> wrote:
>
>>
>> ---------- Forwarded message ---------
>> From: harish sorrout <harishsorrout at gmail.com>
>> Date: Fri, 15 Jul 2022, 21:02
>> Subject: Re: [Mesa-users] How to extract values
>> To: Rob Farmer <robert.j.farmer37 at gmail.com>
>>
>>
>> Sir
>> As you told to use
>> real(dp) function star_get_profile_output(s,name,k) to get
>> the value of 'burn_c' . So for this I need to put 'burn_c' at the place of
>> 'name' argument.
>> So i select character (len=maxlen_profile_column_name) :: burn_c. And use
>> E_cburn(nz) = star_get_profile_output(s,burn_c,k) where real(dp) ::
>> E_cburn(nz). and
>> then used burn_timescales = (s% cp(k) * s% T) / E_cburn(k)
>>
>> By doing so, I did not get any errors while compiling. Sir Is there
>> anything wrong with the above command.
>> Sir I could not understand your words "You want to pass the string
>> "burn_c" not a variable called burn_c".
>> Please enlighten me.
>>
>> 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/20220722/52dfd6ce/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist_to_lgT_8.45
Type: application/octet-stream
Size: 8501 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20220722/52dfd6ce/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 0.85Mwd_0.5crelax.mod
Type: audio/x-mod
Size: 444151 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20220722/52dfd6ce/attachment.bin>
More information about the Mesa-users
mailing list