[Mesa-users] Modeling a Low-Mass Gas-Rich Planet
Brian Jackson
bjackson at boisestate.edu
Wed Oct 9 19:48:39 EDT 2019
Dear MESA Users,
I'm trying to replicate some results from Chen & Rogers (2016 -
https://ui.adsabs.harvard.edu/abs/2016ApJ...831..180C/abstract), but
I'm having trouble.
I can successfully create the initial gas ball for the planet, but
when I try to insert a core, the code runs for a while and then seems
to hang when the core reaches about 0.2 Earth mass. I've attached a
text file showing the last several lines that keep repeating over and
over again.
I am running MESA version 11701 on Mac OS 10.13.1 (17B1003), and I
believe I'm using the MESA SDK from August 30 2019 (gfortran --version
returns "gcc.exec (GCC) 8.3.0"). I have attached the .mod file
representing the initial gas ball and the inlist I'm using to insert
the core.
(And apologies in advance for forgetting something in my e-mail or
sending the wrong files -- I give it 50-50 odds.)
Thanks for your help,
Brian
---------------------------------------------------
Brian Jackson, Associate Professor
Dept. of Physics
Boise State University
1910 University Drive
Boise ID 83725-1570
E-mail: bjackson at boisestate.edu
URL: http://www.astrojack.com
Phone: (208) 426-3723
---------------------------------------------------
-------------- next part --------------
__________________________________________________________________________________________________________________________________________________
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
__________________________________________________________________________________________________________________________________________________
1410 3.836313 283.888 -25.404952 -25.404952 0.000084 0.000084 0.730000 0.001008 0.250000 0.730000 -19.999850 660 0
-1.436501 -1.110590 -0.185956 -99.000000 -26.802816 -99.000000 0.000000 0.250000 0.009360 0.020000 0.250000 15.547877 3 0
2.0946E+01 10.423647 -5.605840 0.734375 3.996636 -6.009515 0.000000 0.003442 0.002099 2.000E-02 2.000E-02 0.000E+00 varcontrol
current & target core masses 5.0118723362729166D-07 6.0024130303639664D-05 9.9165023747796144D-01
s% dt < relax_core_dt*0.9d0 1.1550760906728776D+06 5.6804669971200004D+06
s% max_timestep 3.1558149984000001D+07
s% max_years_for_timestep*secyer 6.3116299968000008D+06
1411 3.836310 276.183 -25.404974 -25.404974 0.000084 0.000084 0.730000 0.001008 0.250000 0.730000 -19.999850 660 0
-1.530693 -1.110600 -0.186127 -99.000000 -26.802838 -99.000000 0.000000 0.250000 0.009360 0.020000 0.250000 15.547866 3 0
2.0975E+01 10.423631 -5.653983 0.734718 3.999893 -5.994309 0.000000 0.003442 0.002099 2.000E-02 2.000E-02 0.000E+00 varcontrol
current & target core masses 5.0118723362729166D-07 6.0024130303639664D-05 9.9165023747796144D-01
s% dt < relax_core_dt*0.9d0 9.2986235350481141D+05 5.6804669971200004D+06
s% max_timestep 3.1558149984000001D+07
s% max_years_for_timestep*secyer 6.3116299968000008D+06
1412 3.836312 281.976 -25.404959 -25.404959 0.000084 0.000084 0.730000 0.001008 0.250000 0.730000 -19.999850 641 0
-1.429412 -1.110594 -0.186013 -99.000000 -26.802824 -99.000000 0.000000 0.250000 0.009360 0.020000 0.250000 15.547875 3 0
2.1013E+01 10.423642 -5.617696 0.734490 3.997686 -6.005530 0.000000 0.003442 0.002099 2.000E-02 2.000E-02 0.000E+00 varcontrol
current & target core masses 5.0118723362729166D-07 6.0024130303639664D-05 9.9165023747796144D-01
s% dt < relax_core_dt*0.9d0 1.1740846271683362D+06 5.6804669971200004D+06
s% max_timestep 3.1558149984000001D+07
s% max_years_for_timestep*secyer 6.3116299968000008D+06
1413 3.836312 281.754 -25.404960 -25.404960 0.000084 0.000084 0.730000 0.001008 0.250000 0.730000 -19.999850 660 0
-1.303275 -1.110594 -0.186003 -99.000000 -26.802824 -99.000000 0.000000 0.250000 0.009360 0.020000 0.250000 15.547876 3 0
2.1062E+01 10.423641 -5.619046 0.734469 3.997782 -6.005091 0.000000 0.003442 0.002099 2.000E-02 2.000E-02 0.000E+00 varcontrol
current & target core masses 5.0118723362729166D-07 6.0024130303639664D-05 9.9165023747796144D-01
s% dt < relax_core_dt*0.9d0 1.5697726383627139D+06 5.6804669971200004D+06
s% max_timestep 3.1558149984000001D+07
s% max_years_for_timestep*secyer 6.3116299968000008D+06
1414 3.836313 283.854 -25.404952 -25.404952 0.000084 0.000084 0.730000 0.001008 0.250000 0.730000 -19.999850 641 0
-1.400185 -1.110590 -0.185931 -99.000000 -26.802816 -99.000000 0.000000 0.250000 0.009360 0.020000 0.250000 15.547881 3 0
2.1102E+01 10.423647 -5.606002 0.734326 3.996629 -6.009469 0.000000 0.003442 0.002099 2.000E-02 2.000E-02 0.000E+00 max decrease
current & target core masses 5.0118723362729166D-07 6.0024130303639664D-05 9.9165023747796144D-01
s% dt < relax_core_dt*0.9d0 1.2558181106901711D+06 5.6804669971200004D+06
s% max_timestep 3.1558149984000001D+07
s% max_years_for_timestep*secyer 6.3116299968000008D+06
1415 3.836315 286.934 -25.404939 -25.404939 0.000084 0.000084 0.730000 0.001008 0.250000 0.730000 -19.999850 660 0
-1.495729 -1.110585 -0.185819 -99.000000 -26.802803 -99.000000 0.000000 0.250000 0.009360 0.020000 0.250000 15.547889 3 0
2.1134E+01 10.423656 -5.587029 0.734101 3.994708 -6.016078 0.000000 0.003442 0.002099 2.000E-02 2.000E-02 0.000E+00 varcontrol
current & target core masses 5.0118723362729166D-07 6.0024130303639664D-05 9.9165023747796144D-01
s% dt < relax_core_dt*0.9d0 1.0078193489988719D+06 5.6804669971200004D+06
s% max_timestep 3.1558149984000001D+07
s% max_years_for_timestep*secyer 6.3116299968000008D+06
1416 3.836311 279.025 -25.404968 -25.404968 0.000084 0.000084 0.730000 0.001008 0.250000 0.730000 -19.999850 660 0
-1.590104 -1.110597 -0.186014 -99.000000 -26.802832 -99.000000 0.000000 0.250000 0.009360 0.020000 0.250000 15.547874 3 0
2.1160E+01 10.423636 -5.635971 0.734491 3.998948 -5.999700 0.000000 0.003442 0.002099 2.000E-02 2.000E-02 0.000E+00 varcontrol
current & target core masses 5.0118723362729166D-07 6.0024130303639664D-05 9.9165023747796144D-01
s% dt < relax_core_dt*0.9d0 8.1097505342360109D+05 5.6804669971200004D+06
s% max_timestep 3.1558149984000001D+07
s% max_years_for_timestep*secyer 6.3116299968000008D+06
1417 3.836314 285.300 -25.404946 -25.404946 0.000084 0.000084 0.730000 0.001008 0.250000 0.730000 -19.999850 641 0
-1.494416 -1.110588 -0.185861 -99.000000 -26.802810 -99.000000 0.000000 0.250000 0.009360 0.020000 0.250000 15.547887 3 0
2.1192E+01 10.423651 -5.597032 0.734185 3.995737 -6.012569 0.000000 0.003442 0.002099 2.000E-02 2.000E-02 0.000E+00 varcontrol
current & target core masses 5.0118723362729166D-07 6.0024130303639664D-05 9.9165023747796144D-01
s% dt < relax_core_dt*0.9d0 1.0108707591488528D+06 5.6804669971200004D+06
s% max_timestep 3.1558149984000001D+07
s% max_years_for_timestep*secyer 6.3116299968000008D+06
1418 3.836314 285.419 -25.404945 -25.404945 0.000084 0.000084 0.730000 0.001008 0.250000 0.730000 -19.999850 660 0
-1.372879 -1.110587 -0.185831 -99.000000 -26.802809 -99.000000 0.000000 0.250000 0.009360 0.020000 0.250000 15.547888 3 0
2.1234E+01 10.423651 -5.596247 0.734126 3.995635 -6.012852 0.000000 0.003442 0.002099 2.000E-02 2.000E-02 0.000E+00 varcontrol
current & target core masses 5.0118723362729166D-07 6.0024130303639664D-05 9.9165023747796144D-01
s% dt < relax_core_dt*0.9d0 1.3373119213698863D+06 5.6804669971200004D+06
s% max_timestep 3.1558149984000001D+07
s% max_years_for_timestep*secyer 6.3116299968000008D+06
1419 3.836311 277.471 -25.404971 -25.404971 0.000084 0.000084 0.730000 0.001008 0.250000 0.730000 -19.999850 660 0
-1.469789 -1.110599 -0.186033 -99.000000 -26.802836 -99.000000 0.000000 0.250000 0.009360 0.020000 0.250000 15.547875 3 0
2.1268E+01 10.423633 -5.645714 0.734529 3.999473 -5.996749 0.000000 0.003442 0.002099 2.000E-02 2.000E-02 0.000E+00 max decrease
current & target core masses 5.0118723362729166D-07 6.0024130303639664D-05 9.9165023747796144D-01
s% dt < relax_core_dt*0.9d0 1.0698495370959090D+06 5.6804669971200004D+06
s% max_timestep 3.1558149984000001D+07
s% max_years_for_timestep*secyer 6.3116299968000008D+06
-------------- next part --------------
A non-text attachment was scrubbed...
Name: planet_create_28_MJ.mod
Type: application/octet-stream
Size: 266567 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20191009/e1a12ef2/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist
Type: application/octet-stream
Size: 2004 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20191009/e1a12ef2/attachment-0003.obj>
More information about the Mesa-users
mailing list