[Mesa-users] Test not terminating at the required phase
Evan Bauer
ebauer at physics.ucsb.edu
Sat Jun 13 22:06:28 EDT 2020
Hi Kim,
Thanks! That makes it easy to look for problems. Here’s one in inlist_to_wd_header that explains your pgstar issue:
&pgstar
read_extra_pgstar_inlist1 = .true.
extra_pgstar_inlist1_name = 'inlist_to_wd'
/ ! end of pgstar namelist
That’s pointing to inlist_to_wd, which has an empty pgstar section, so I think you want to change that to point to inlist_pgstar like inlist_to_end_agb_header does. With your current setup, there will be no pgstar output during the last phase of evolution onto the white dwarf track even though the model is successfully getting there.
As for the LOGS output, the “rn” script for this case is set up to move LOGS to a new directory at the end of each stage of the run. E.g. for the final stage the LOGS information should end up in a directory called LOGS_to_wd. Do you see that?
Cheers,
Evan
> On Jun 13, 2020, at 6:33 PM, Kimberly Appel <kimberlyappel00 at gmail.com> wrote:
>
> 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 <mailto: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
>
>
> <unknown.png>
>
>> On Jun 13, 2020, at 4:15 PM, Kimberly Appel <kimberlyappel00 at gmail.com <mailto: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 <mailto: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 <mailto: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 <mailto: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 <mailto: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 <mailto: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 <mailto:mesa-users at lists.mesastar.org>
>>> > https://lists.mesastar.org/mailman/listinfo/mesa-users <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/20200613/4505fed7/attachment.htm>
More information about the Mesa-users
mailing list