[Mesa-users] Test not terminating at the required phase

Kimberly Appel kimberlyappel00 at gmail.com
Sat Jun 13 21:33:46 EDT 2020 

 1M (working directory).zip
<https://drive.google.com/file/d/14SXtM1Fwb-tziLNTDcv94hGPTSG-UzeB/view?usp=drive_web>

Hi Evan,

I have attached my working directory. Do you know why the files in the LOGS
disappear after the run terminates? I need the data in the history.data
file for later analysis.

KInd regards,
Kim.

On Sun, Jun 14, 2020 at 1:26 AM Evan Bauer <ebauer at physics.ucsb.edu> wrote:

> Hi Kim,
>
> (please keep replies on the mesa-users list so others can benefit from
> seeing a solution if they have similar problems)
>
> This png file was saved at model 5631, but your terminal output shows that
> the run went to model 5864 before terminating. My guess is that something
> about your pgstar settings caused you to stop saving pngs when
> switching from inlist_to_end_agb to inlist_to_wd. Or maybe you saved them
> in some other location? It’s hard to say any more without seeing a full
> copy of your work directory.
>
> Cheers,
> Evan
>
>
> [image: unknown.png]
>
> On Jun 13, 2020, at 4:15 PM, Kimberly Appel <kimberlyappel00 at gmail.com>
> wrote:
>
> Sorry, that's a wrong model photo. I attached the correct photo
>
>
> On Sun, Jun 14, 2020 at 1:12 AM Kimberly Appel <kimberlyappel00 at gmail.com>
> wrote:
>
>> Hi Evan,
>>
>> The HR diagram in the PGPlot seems to indicate that it has not reached
>> the WD stage. I am not sure. Maybe it's because I set static plot bounds
>> for log Teff (HR_logT_min = 3.0, HR_logT_max = 6.0)?
>>
>> I have attached the last model photo (when the run terminates).
>>
>> Kind regards,
>> KIm.
>>
>> On Sun, Jun 14, 2020 at 12:22 AM Evan Bauer <ebauer at physics.ucsb.edu>
>> wrote:
>>
>>> Hi Kim,
>>>
>>> This output indicates that the model has reached the desired stopping
>>> condition:
>>> stop because log_surface_luminosity <= log_L_lower_limit
>>>
>>> The column for luminosity in the output also shows lg_L =  -1.004672
>>>
>>> Can you clarify what seems to indicate that it has not reached the WD
>>> stage? Could this be an issue with your pgstar plots perhaps?
>>>
>>> Cheers,
>>> Evan
>>>
>>>
>>> On Jun 13, 2020, at 5:31 AM, Kimberly Appel <kimberlyappel00 at gmail.com>
>>> wrote:
>>>
>>> Hi Evan,
>>>
>>> This is the terminal output I get when the run terminates:
>>>
>>> stop because log_surface_luminosity <= log_L_lower_limit
>>>    -0.1004671865E+01   -0.1000000000E+01
>>>
>>>        5864   7.620999  2.535E+04  -1.712706  -1.712706   0.510344
>>> 0.000279   0.000000   0.000000   0.301345   0.000289  61.736875   1139
>>> 16
>>>    6.074703   6.327943  -1.786974 -18.571076  -1.239491 -25.392495
>>> 0.510065   0.000000   0.700169   0.020075   0.074805  15.622100      6
>>>  0
>>>  1.2547E+10  22.916722  -1.004672   7.719701   4.506902  -5.426933
>>> 0.446952   0.276606   0.002056  1.000E+00  9.249E-01  0.000E+00  max
>>> increase
>>>                                rel_E_err       -0.0005495591127248
>>>                        log_rel_run_E_err       -1.9417731348973926
>>>
>>> save LOGS/profile7.data for model 5864
>>> save photos/x864 for model 5864
>>>  saved to final.mod
>>> termination code: log_L_lower_limit
>>>
>>>         runtime (minutes), retries, backups, steps        2.99        16
>>>         0      5864
>>>
>>> __________________________________________________________________________________________________________________________________________________
>>>
>>>        step    lg_Tcntr    Teff     lg_LH       lg_Lnuc     Mass
>>> H_rich     H_cntr     N_cntr     Y_surf     X_avg     eta_cntr  zones retry
>>>    lg_dt_yr    lg_Dcntr    lg_R     lg_L3a      lg_Lneu     lg_Mdot
>>>  He_core    He_cntr    O_cntr     Z_surf     Y_avg     gam_cntr  iters bckup
>>>      age_yr    lg_Pcntr    lg_L     lg_gsurf    lg_Psurf    lg_Dsurf
>>> C_core     C_cntr     Ne_cntr    Z_cntr     Z_avg     v_div_cs     dt_limit
>>>
>>> __________________________________________________________________________________________________________________________________________________
>>>
>>>        5864   7.620999  2.535E+04  -1.712706  -1.712706   0.510344
>>> 0.000279   0.000000   0.000000   0.301345   0.000289  61.736875   1139
>>> 16
>>>    6.074703   6.327943  -1.786974 -18.571076  -1.239491 -25.392495
>>> 0.510065   0.000000   0.700169   0.020075   0.074805  15.622100      6
>>>  0
>>>  1.2547E+10  22.916722  -1.004672   7.719701   4.506902  -5.426933
>>> 0.446952   0.276606   0.002056  1.000E+00  9.249E-01  0.000E+00  max
>>> increase
>>>                                rel_E_err       -0.0005495591127248
>>>                        log_rel_run_E_err       -1.9417731348973926
>>>
>>> DATE: 2020-06-13
>>> TIME: 14:06:24
>>>
>>> finished inlist_to_wd_header
>>>
>>> ***
>>> ***
>>> ***
>>> ***
>>> ***
>>> ***
>>>
>>> DATE: 2020-06-13
>>> TIME: 14:06:24
>>> finished all inlists for 1M_pre_ms_to_wd
>>> (base) kim at kim-Lenovo-V130-15IKB:~/1M/1M_pre_ms_to_wd_COPY$
>>>
>>> I increased the maximum number of models to run but that did not help -
>>> the test does not get past the AGB stage, i.e., it never reaches the white
>>> dwarf stage. I did set  log_L_lower_limit = -1 in my the model
>>> (inlist_to_wd) but still, that also did not help. Also, my LOGS folder is
>>> empty. I am really worried about that.
>>>
>>> Kind regards,
>>> Kim
>>>
>>> On Thu, Jun 11, 2020 at 7:35 PM Evan Bauer <ebauer at physics.ucsb.edu>
>>> wrote:
>>>
>>>> Hi Kim,
>>>>
>>>> Can you send the terminal output that you get when the run terminates?
>>>> That should tell us what stopping condition the run has reached.
>>>>
>>>> Cheers,
>>>> Evan
>>>>
>>>>
>>>> > On Jun 10, 2020, at 1:46 AM, Kimberly Appel via Mesa-users <
>>>> mesa-users at lists.mesastar.org> wrote:
>>>> >
>>>> > Hi, I am using a Lenovo-V130-15IKB laptop. My operating system is
>>>> Ubuntu. The version of MESA I am using is r12778 and the version of MESA
>>>> SDK I am using is x86_64-Linux-20.3.2.
>>>> >
>>>> > So, I am trying to run a 1Msun star model from the pre-main sequence
>>>> to the white dwarf phase but the test is not terminating when the log of
>>>> the surface luminosity drops below -1 as it should. Instead, it terminates
>>>> at around the end of the thermally pulsing AGB phase (I think) when the log
>>>> of the surface luminosity is around 3.5. I am not sure what I am missing? I
>>>> set  log_L_lower_limit = -1 in the inlist_project file.
>>>> >
>>>> > I have attached the inlist_to_wd, inlist_pgstar, history_columns.list
>>>> files and a photo. Please let me know if you need additional files.
>>>> >
>>>> > KInd regards,
>>>> > Kim.
>>>> >
>>>> > <inlist_to_wd><inlist_pgstar><history_columns.list><HR
>>>> diagram.png>_______________________________________________
>>>> > mesa-users at lists.mesastar.org
>>>> > https://lists.mesastar.org/mailman/listinfo/mesa-users
>>>> >
>>>>
>>>>
>>> <grid1005631.png>
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20200614/028168a4/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: unknown.png
Type: image/png
Size: 45936 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20200614/028168a4/attachment.png>

More information about the Mesa-users mailing list