[Mesa-users] "Mirror star nuggets" in MESA with modified gravity and heating rate
Franco Cabral
franco.cabral at mail.utoronto.ca
Mon Jan 6 00:26:03 UTC 2025
Dear all,
I am reaching out for advice on an issue I am encountering while trying to model optically thick "mirror star nuggets" in MESA (version 24.03.1) as part of a dark matter study. These nuggets are intended to represent compact objects with around the following specifications:
1.
Mass: 10^20 g
2.
Radius: 10^9 cm
3.
Heating Rate: 10^10 erg/g/s (implemented via extra_power_source in the inlist), driven by suppressed interactions with dark matter.
4.
Modified gravity: Gravitational acceleration set by a mirror star core density of around 150 g/cm^3.
I have modified create_initial_model.f90 to change the local gravity definition there and its associated differential equations. With these changes, MESA is able to converge on an initial model with specified mass and radius. However, once this initial model is passed to the main MESA solver, it encounters tolerance and convergence errors in the first simulation step.
Here are my attempts so far in changing the gravity which have been unsuccessful:
1.
Modifying MESA's global m_grav definition in star_utils.f90.
2.
Using the other_cgrav hook.
3.
Using the other_momentum_implicit hook.
I suspect that the nugget is evolving too rapidly from its initial model. Adjusting certain tolerance switches may help, but I’m unsure which ones are most relevant or if the issue lies elsewhere. I would greatly appreciate any advice on this matter.
I am attaching a GitHub link with my code changes and a sample run output using the other_cgrav hook to show the errors.
https://github.com/MESAHub/mesa/compare/main...fmcabral123:mesa:main<Https://github.com/MESAHub/mesa/compare/main...fmcabral123:mesa:main>
Sincerely,
Franco Cabral
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20250106/9748fd1e/attachment.htm>
-------------- next part --------------
reading user weak rate file /home/fcabral/mesa-r24.03.1/data/rates_data/rate_tables/S13_r_be7_wk_li7.h5
version_number r24.03.1
read inlist_make_brown_dwarf
create initial model
create initial model iteration 1
log10(L/Lsun)= -5.58165088E+00
Mass= 4.81409497E+24
Radius= 4.68475525E+08
log10(L/Lsun)= -5.62253308E+00
Mass= 4.72287928E+24
Radius= 4.59541714E+08
log10(L/Lsun)= -5.62120232E+00
Mass= 4.77768093E+24
Radius= 4.60208489E+08
create initial model iteration 2
log10(L/Lsun)= -5.47258375E+00
Mass= 4.57159475E+24
Radius= 4.92697718E+08
log10(L/Lsun)= -5.50731921E+00
Mass= 4.48632513E+24
Radius= 4.83309203E+08
log10(L/Lsun)= -5.50594787E+00
Mass= 4.53789437E+24
Radius= 4.83997965E+08
create initial model iteration 3
log10(L/Lsun)= -5.40997661E+00
Mass= 4.25334303E+24
Radius= 5.07526152E+08
log10(L/Lsun)= -5.37662467E+00
Mass= 4.33386527E+24
Radius= 5.17675263E+08
log10(L/Lsun)= -5.40863640E+00
Mass= 4.30182074E+24
Radius= 5.08237913E+08
create initial model iteration 4
log10(L/Lsun)= -5.28898500E+00
Mass= 4.10052475E+24
Radius= 5.43266215E+08
log10(L/Lsun)= -5.32201600E+00
Mass= 4.02346066E+24
Radius= 5.32174341E+08
log10(L/Lsun)= -5.32070933E+00
Mass= 4.06893993E+24
Radius= 5.32907121E+08
create initial model iteration 5
log10(L/Lsun)= -5.20869234E+00
Mass= 3.86945726E+24
Radius= 5.68917871E+08
log10(L/Lsun)= -5.24181576E+00
Mass= 3.79586411E+24
Radius= 5.56853240E+08
log10(L/Lsun)= -5.24052008E+00
Mass= 3.83849412E+24
Radius= 5.57605917E+08
create initial model iteration 6
log10(L/Lsun)= -5.16846417E+00
Mass= 3.57167536E+24
Radius= 5.81325879E+08
log10(L/Lsun)= -5.13525852E+00
Mass= 3.64209331E+24
Radius= 5.94308232E+08
log10(L/Lsun)= -5.16717217E+00
Mass= 3.61162305E+24
Radius= 5.82099019E+08
create initial model iteration 7
log10(L/Lsun)= -5.10122634E+00
Mass= 3.35245934E+24
Radius= 6.05181597E+08
log10(L/Lsun)= -5.06806641E+00
Mass= 3.41941645E+24
Radius= 6.18863735E+08
log10(L/Lsun)= -5.09994526E+00
Mass= 3.38987112E+24
Radius= 6.05976206E+08
create initial model iteration 8
log10(L/Lsun)= -5.03912315E+00
Mass= 3.14017538E+24
Radius= 6.28262895E+08
log10(L/Lsun)= -5.00640822E+00
Mass= 3.20368906E+24
Radius= 6.42406231E+08
log10(L/Lsun)= -5.03793613E+00
Mass= 3.17519832E+24
Radius= 6.29083879E+08
create initial model iteration 9
log10(L/Lsun)= -4.98144649E+00
Mass= 2.93685402E+24
Radius= 6.50472925E+08
log10(L/Lsun)= -4.94905822E+00
Mass= 2.99710999E+24
Radius= 6.64975166E+08
log10(L/Lsun)= -4.98024992E+00
Mass= 2.96963579E+24
Radius= 6.51317956E+08
create initial model iteration 10
log10(L/Lsun)= -4.92766260E+00
Mass= 2.74270875E+24
Radius= 6.71628969E+08
log10(L/Lsun)= -4.89566476E+00
Mass= 2.79937735E+24
Radius= 6.86299615E+08
log10(L/Lsun)= -4.92646009E+00
Mass= 2.77338262E+24
Radius= 6.72499912E+08
create initial model iteration 11
log10(L/Lsun)= -4.87688311E+00
Mass= 2.55817787E+24
Radius= 6.91768668E+08
log10(L/Lsun)= -4.84518918E+00
Mass= 2.61162433E+24
Radius= 7.06524024E+08
log10(L/Lsun)= -4.87567710E+00
Mass= 2.58686169E+24
Radius= 6.92667488E+08
create initial model iteration 12
log10(L/Lsun)= -4.79698316E+00
Mass= 2.43340349E+24
Radius= 7.25719115E+08
log10(L/Lsun)= -4.82849461E+00
Mass= 2.38357531E+24
Radius= 7.10935040E+08
log10(L/Lsun)= -4.82728980E+00
Mass= 2.41037880E+24
Radius= 7.11863448E+08
create initial model iteration 13
log10(L/Lsun)= -4.78198291E+00
Mass= 2.21790523E+24
Radius= 7.29095369E+08
log10(L/Lsun)= -4.75065185E+00
Mass= 2.26335723E+24
Radius= 7.43816272E+08
log10(L/Lsun)= -4.78079502E+00
Mass= 2.24291325E+24
Radius= 7.30054078E+08
create initial model iteration 14
log10(L/Lsun)= -4.73658548E+00
Mass= 2.06097062E+24
Radius= 7.46544168E+08
log10(L/Lsun)= -4.70499053E+00
Mass= 2.10295033E+24
Radius= 7.61392537E+08
log10(L/Lsun)= -4.73532628E+00
Mass= 2.08428264E+24
Radius= 7.47529167E+08
create initial model iteration 15
log10(L/Lsun)= -4.66049796E+00
Mass= 1.95037548E+24
Radius= 7.78149958E+08
log10(L/Lsun)= -4.69249412E+00
Mass= 1.91224412E+24
Radius= 7.63185272E+08
log10(L/Lsun)= -4.69115898E+00
Mass= 1.93393669E+24
Radius= 7.64195604E+08
create initial model iteration 16
log10(L/Lsun)= -4.64935263E+00
Mass= 1.77112753E+24
Radius= 7.79071243E+08
log10(L/Lsun)= -4.61694553E+00
Mass= 1.80577171E+24
Radius= 7.94122319E+08
log10(L/Lsun)= -4.64806144E+00
Mass= 1.79125482E+24
Radius= 7.80108107E+08
create initial model iteration 17
log10(L/Lsun)= -4.60683021E+00
Mass= 1.63781704E+24
Radius= 7.94201938E+08
log10(L/Lsun)= -4.57349053E+00
Mass= 1.66898085E+24
Radius= 8.09329250E+08
log10(L/Lsun)= -4.60556932E+00
Mass= 1.65649446E+24
Radius= 7.95265499E+08
create initial model iteration 18
log10(L/Lsun)= -4.52822492E+00
Mass= 1.54000511E+24
Radius= 8.23831123E+08
log10(L/Lsun)= -4.56368392E+00
Mass= 1.51197710E+24
Radius= 8.08638064E+08
log10(L/Lsun)= -4.56246124E+00
Mass= 1.52924707E+24
Radius= 8.09724099E+08
create initial model iteration 19
log10(L/Lsun)= -4.47748749E+00
Mass= 1.41892299E+24
Radius= 8.37574529E+08
log10(L/Lsun)= -4.51812592E+00
Mass= 1.39359657E+24
Radius= 8.22343314E+08
log10(L/Lsun)= -4.51677954E+00
Mass= 1.40955481E+24
Radius= 8.23443829E+08
create initial model iteration 20
log10(L/Lsun)= -4.46491925E+00
Mass= 1.28278699E+24
Radius= 8.35243179E+08
log10(L/Lsun)= -4.41179602E+00
Mass= 1.30569221E+24
Radius= 8.50385938E+08
log10(L/Lsun)= -4.46357068E+00
Mass= 1.29748188E+24
Radius= 8.36359036E+08
create initial model iteration 21
log10(L/Lsun)= -4.32074917E+00
Mass= 1.20075450E+24
Radius= 8.62467297E+08
log10(L/Lsun)= -4.39247328E+00
Mass= 1.17971485E+24
Radius= 8.47337847E+08
log10(L/Lsun)= -4.39105717E+00
Mass= 1.19323243E+24
Radius= 8.48465526E+08
create initial model iteration 22
log10(L/Lsun)= -4.29609693E+00
Mass= 1.08410955E+24
Radius= 8.58676917E+08
log10(L/Lsun)= -4.19872069E+00
Mass= 1.10353190E+24
Radius= 8.73926495E+08
log10(L/Lsun)= -4.29460543E+00
Mass= 1.09653385E+24
Radius= 8.59816442E+08
create initial model iteration 23
log10(L/Lsun)= -4.16307226E+00
Mass= 9.95477070E+23
Radius= 8.69171467E+08
log10(L/Lsun)= -4.01674430E+00
Mass= 1.01346590E+24
Radius= 8.84330940E+08
log10(L/Lsun)= -4.16144783E+00
Mass= 1.00689073E+24
Radius= 8.70322328E+08
create initial model iteration 24
log10(L/Lsun)= -3.96752339E+00
Mass= 9.13665542E+23
Radius= 8.78866912E+08
log10(L/Lsun)= -3.78044270E+00
Mass= 9.30265059E+23
Radius= 8.93865839E+08
log10(L/Lsun)= -3.96576839E+00
Mass= 9.24148043E+23
Radius= 8.80030602E+08
create initial model iteration 25
log10(L/Lsun)= -3.72705671E+00
Mass= 8.38334508E+23
Radius= 8.87831672E+08
log10(L/Lsun)= -3.57213109E+00
Mass= 8.53651445E+23
Radius= 9.03390159E+08
log10(L/Lsun)= -3.72535917E+00
Mass= 8.47951027E+23
Radius= 8.89008667E+08
create initial model iteration 26
log10(L/Lsun)= -3.56416936E+00
Mass= 7.68306156E+23
Radius= 8.96582065E+08
log10(L/Lsun)= -3.51741216E+00
Mass= 7.82466537E+23
Radius= 9.14089104E+08
log10(L/Lsun)= -3.56299775E+00
Mass= 7.77121807E+23
Radius= 8.97778052E+08
create initial model iteration 27
log10(L/Lsun)= -3.52420046E+00
Mass= 7.02403754E+23
Radius= 9.04606067E+08
log10(L/Lsun)= -3.48599791E+00
Mass= 7.15426063E+23
Radius= 9.22349057E+08
log10(L/Lsun)= -3.52256305E+00
Mass= 7.10463801E+23
Radius= 9.05798157E+08
create initial model iteration 28
log10(L/Lsun)= -3.49145280E+00
Mass= 6.41681498E+23
Radius= 9.11873601E+08
log10(L/Lsun)= -3.45944029E+00
Mass= 6.53625461E+23
Radius= 9.29795534E+08
log10(L/Lsun)= -3.49345146E+00
Mass= 6.49045387E+23
Radius= 9.13135543E+08
create initial model iteration 29
log10(L/Lsun)= -3.46903985E+00
Mass= 5.86158414E+23
Radius= 9.19197446E+08
log10(L/Lsun)= -3.43458654E+00
Mass= 5.97123529E+23
Radius= 9.37142585E+08
log10(L/Lsun)= -3.46577044E+00
Mass= 5.92885874E+23
Radius= 9.20366162E+08
create initial model iteration 30
log10(L/Lsun)= -3.44823698E+00
Mass= 5.34792553E+23
Radius= 9.25152038E+08
log10(L/Lsun)= -3.41788283E+00
Mass= 5.44840543E+23
Radius= 9.43220077E+08
log10(L/Lsun)= -3.44828884E+00
Mass= 5.40930508E+23
Radius= 9.26384997E+08
create initial model iteration 31
log10(L/Lsun)= -3.42940423E+00
Mass= 4.87934835E+23
Radius= 9.31114622E+08
log10(L/Lsun)= -3.39684531E+00
Mass= 4.97140727E+23
Radius= 9.49182410E+08
log10(L/Lsun)= -3.42601635E+00
Mass= 4.93535198E+23
Radius= 9.32283443E+08
create initial model iteration 32
log10(L/Lsun)= -3.41082717E+00
Mass= 4.44782023E+23
Radius= 9.36037750E+08
log10(L/Lsun)= -3.38192625E+00
Mass= 4.53201277E+23
Radius= 9.54201829E+08
log10(L/Lsun)= -3.41031213E+00
Mass= 4.49887309E+23
Radius= 9.37263119E+08
create initial model iteration 33
log10(L/Lsun)= -3.36516778E+00
Mass= 4.13162677E+23
Radius= 9.59189925E+08
log10(L/Lsun)= -3.39426273E+00
Mass= 4.05460953E+23
Radius= 9.41004499E+08
log10(L/Lsun)= -3.39457962E+00
Mass= 4.10115437E+23
Radius= 9.42248716E+08
create initial model iteration 34
log10(L/Lsun)= -3.37905904E+00
Mass= 3.69540438E+23
Radius= 9.45733286E+08
log10(L/Lsun)= -3.34796566E+00
Mass= 3.76583810E+23
Radius= 9.63895445E+08
log10(L/Lsun)= -3.37797856E+00
Mass= 3.73782931E+23
Radius= 9.46950543E+08
create initial model iteration 35
log10(L/Lsun)= -3.36349541E+00
Mass= 3.36625934E+23
Radius= 9.49835094E+08
log10(L/Lsun)= -3.33210773E+00
Mass= 3.43062733E+23
Radius= 9.67996606E+08
log10(L/Lsun)= -3.36238174E+00
Mass= 3.40490712E+23
Radius= 9.51051930E+08
create initial model iteration 36
log10(L/Lsun)= -3.31644100E+00
Mass= 3.12436703E+23
Radius= 9.71734423E+08
log10(L/Lsun)= -3.34778726E+00
Mass= 3.06555330E+23
Radius= 9.53574751E+08
log10(L/Lsun)= -3.34715612E+00
Mass= 3.10074999E+23
Radius= 9.54802011E+08
create initial model iteration 37
log10(L/Lsun)= -3.33163780E+00
Mass= 2.79123224E+23
Radius= 9.57068132E+08
log10(L/Lsun)= -3.29987839E+00
Mass= 2.84497321E+23
Radius= 9.75207769E+08
log10(L/Lsun)= -3.33143439E+00
Mass= 2.82328067E+23
Radius= 9.58305050E+08
create initial model iteration 38
log10(L/Lsun)= -3.31656765E+00
Mass= 2.54109289E+23
Radius= 9.60371643E+08
log10(L/Lsun)= -3.28190095E+00
Mass= 2.59020268E+23
Radius= 9.78422719E+08
log10(L/Lsun)= -3.31488960E+00
Mass= 2.57027044E+23
Radius= 9.61575607E+08
create initial model iteration 39
log10(L/Lsun)= -3.30016745E+00
Mass= 2.31237687E+23
Radius= 9.63110441E+08
log10(L/Lsun)= -3.26306780E+00
Mass= 2.35722941E+23
Radius= 9.81056058E+08
log10(L/Lsun)= -3.29921369E+00
Mass= 2.33892876E+23
Radius= 9.64328978E+08
create initial model iteration 40
log10(L/Lsun)= -3.28210903E+00
Mass= 2.10428162E+23
Radius= 9.65748288E+08
log10(L/Lsun)= -3.24041753E+00
Mass= 2.14526320E+23
Radius= 9.83504595E+08
log10(L/Lsun)= -3.28086626E+00
Mass= 2.12844504E+23
Radius= 9.66960353E+08
create initial model iteration 41
log10(L/Lsun)= -3.26103601E+00
Mass= 1.91474306E+23
Radius= 9.68112852E+08
log10(L/Lsun)= -3.21507651E+00
Mass= 1.95219544E+23
Radius= 9.85645369E+08
log10(L/Lsun)= -3.25988519E+00
Mass= 1.93673120E+23
Radius= 9.69327470E+08
create initial model iteration 42
log10(L/Lsun)= -3.23633589E+00
Mass= 1.74233958E+23
Radius= 9.70285762E+08
log10(L/Lsun)= -3.19046871E+00
Mass= 1.77658364E+23
Radius= 9.87703617E+08
log10(L/Lsun)= -3.23531186E+00
Mass= 1.76234982E+23
Radius= 9.71503106E+08
create initial model iteration 43
log10(L/Lsun)= -3.21140822E+00
Mass= 1.58541796E+23
Radius= 9.72274093E+08
log10(L/Lsun)= -3.16708820E+00
Mass= 1.61673977E+23
Radius= 9.89621657E+08
log10(L/Lsun)= -3.21005007E+00
Mass= 1.60362907E+23
Radius= 9.73481723E+08
create initial model iteration 44
log10(L/Lsun)= -3.18831541E+00
Mass= 1.44241108E+23
Radius= 9.74003410E+08
log10(L/Lsun)= -3.14574381E+00
Mass= 1.47106056E+23
Radius= 9.91287563E+08
log10(L/Lsun)= -3.18722785E+00
Mass= 1.45898146E+23
Radius= 9.75219814E+08
create initial model iteration 45
log10(L/Lsun)= -3.16665659E+00
Mass= 1.31234123E+23
Radius= 9.75655183E+08
log10(L/Lsun)= -3.12579431E+00
Mass= 1.33855914E+23
Radius= 9.92880984E+08
log10(L/Lsun)= -3.16530900E+00
Mass= 1.32741938E+23
Radius= 9.76862259E+08
create initial model iteration 46
log10(L/Lsun)= -3.14651268E+00
Mass= 1.19384612E+23
Radius= 9.77064909E+08
log10(L/Lsun)= -3.10734308E+00
Mass= 1.21784083E+23
Radius= 9.94238073E+08
log10(L/Lsun)= -3.14521356E+00
Mass= 1.20756481E+23
Radius= 9.78273348E+08
create initial model iteration 47
log10(L/Lsun)= -3.12774525E+00
Mass= 1.08603333E+23
Radius= 9.78352015E+08
log10(L/Lsun)= -3.09012050E+00
Mass= 1.10800156E+23
Radius= 9.95476527E+08
log10(L/Lsun)= -3.12663765E+00
Mass= 1.09851522E+23
Radius= 9.79568866E+08
create initial model iteration 48
log10(L/Lsun)= -3.11014976E+00
Mass= 9.88015064E+22
Radius= 9.79589824E+08
log10(L/Lsun)= -3.07388155E+00
Mass= 1.00813967E+23
Radius= 9.96669567E+08
log10(L/Lsun)= -3.10908633E+00
Mass= 9.99372743E+22
Radius= 9.80809719E+08
create initial model iteration 49
log10(L/Lsun)= -3.09373479E+00
Mass= 8.98855906E+22
Radius= 9.80735417E+08
log10(L/Lsun)= -3.05861322E+00
Mass= 9.17300081E+22
Radius= 9.97772089E+08
log10(L/Lsun)= -3.09259983E+00
Mass= 9.09191088E+22
Radius= 9.81952665E+08
create initial model iteration 50
log10(L/Lsun)= -3.07846182E+00
Mass= 8.17718001E+22
Radius= 9.81747025E+08
log10(L/Lsun)= -3.04426787E+00
Mass= 8.34627982E+22
Radius= 9.98740054E+08
log10(L/Lsun)= -3.07716055E+00
Mass= 8.27122689E+22
Radius= 9.82955700E+08
create initial model iteration 51
log10(L/Lsun)= -3.06425590E+00
Mass= 7.43850132E+22
Radius= 9.82584139E+08
log10(L/Lsun)= -3.03112751E+00
Mass= 7.59357355E+22
Radius= 9.99554891E+08
log10(L/Lsun)= -3.06303120E+00
Mass= 7.52407654E+22
Radius= 9.83798056E+08
create initial model iteration 52
log10(L/Lsun)= -3.05071155E+00
Mass= 6.76695840E+22
Radius= 9.83366263E+08
log10(L/Lsun)= -3.01857544E+00
Mass= 6.90925221E+22
Radius= 1.00032535E+09
log10(L/Lsun)= -3.04960995E+00
Mass= 6.84483303E+22
Radius= 9.84589520E+08
create initial model iteration 53
log10(L/Lsun)= -3.03789206E+00
Mass= 6.15666032E+22
Radius= 9.84142640E+08
log10(L/Lsun)= -3.00650710E+00
Mass= 6.28732663E+22
Radius= 1.00109003E+09
log10(L/Lsun)= -3.03667135E+00
Mass= 6.22753772E+22
Radius= 9.85359339E+08
create initial model iteration 54
log10(L/Lsun)= -3.02593919E+00
Mass= 5.60112406E+22
Radius= 9.84774539E+08
log10(L/Lsun)= -2.99530169E+00
Mass= 5.72115379E+22
Radius= 1.00172413E+09
log10(L/Lsun)= -3.02477725E+00
Mass= 5.66563161E+22
Radius= 9.85997389E+08
create initial model iteration 55
log10(L/Lsun)= -2.98450452E+00
Mass= 5.20645685E+22
Radius= 1.00233445E+09
log10(L/Lsun)= -3.01446415E+00
Mass= 5.09612100E+22
Radius= 9.85371938E+08
log10(L/Lsun)= -3.01326675E+00
Mass= 5.15483839E+22
Radius= 9.86594136E+08
create initial model iteration 56
log10(L/Lsun)= -2.97424645E+00
Mass= 4.73817384E+22
Radius= 1.00286349E+09
log10(L/Lsun)= -3.00367319E+00
Mass= 4.63669978E+22
Radius= 9.85885412E+08
log10(L/Lsun)= -3.00253835E+00
Mass= 4.69014948E+22
Radius= 9.87115304E+08
create initial model iteration 57
log10(L/Lsun)= -2.96447002E+00
Mass= 4.31253290E+22
Radius= 1.00340859E+09
log10(L/Lsun)= -2.99332399E+00
Mass= 4.21913508E+22
Radius= 9.86394213E+08
log10(L/Lsun)= -2.99208738E+00
Mass= 4.26779748E+22
Radius= 9.87618134E+08
create initial model iteration 58
log10(L/Lsun)= -2.95518222E+00
Mass= 3.92498633E+22
Radius= 1.00381854E+09
log10(L/Lsun)= -2.98350355E+00
Mass= 3.83899751E+22
Radius= 9.86758711E+08
log10(L/Lsun)= -2.98237005E+00
Mass= 3.88330389E+22
Radius= 9.87995603E+08
create initial model iteration 59
log10(L/Lsun)= -2.97408417E+00
Mass= 3.49366476E+22
Radius= 9.87175710E+08
log10(L/Lsun)= -2.94609113E+00
Mass= 3.57290562E+22
Radius= 1.00427970E+09
log10(L/Lsun)= -2.97295030E+00
Mass= 3.53401161E+22
Radius= 9.88416467E+08
create initial model iteration 60
log10(L/Lsun)= -2.96506654E+00
Mass= 3.17960618E+22
Radius= 9.87556154E+08
log10(L/Lsun)= -2.93743676E+00
Mass= 3.25267513E+22
Radius= 1.00472335E+09
log10(L/Lsun)= -2.96387083E+00
Mass= 3.21635111E+22
Radius= 9.88794805E+08
create initial model iteration 61
log10(L/Lsun)= -2.95647451E+00
Mass= 2.89376704E+22
Radius= 9.87836977E+08
log10(L/Lsun)= -2.92921295E+00
Mass= 2.96117029E+22
Radius= 1.00508150E+09
log10(L/Lsun)= -2.95523366E+00
Mass= 2.92723264E+22
Radius= 9.89075031E+08
create initial model iteration 62
log10(L/Lsun)= -2.92133370E+00
Mass= 2.69586016E+22
Radius= 1.00536701E+09
log10(L/Lsun)= -2.94827730E+00
Mass= 2.63365825E+22
Radius= 9.88036715E+08
log10(L/Lsun)= -2.94704759E+00
Mass= 2.66413879E+22
Radius= 9.89280232E+08
create initial model iteration 63
log10(L/Lsun)= -2.91369418E+00
Mass= 2.45456690E+22
Radius= 1.00564253E+09
log10(L/Lsun)= -2.94036818E+00
Mass= 2.39712685E+22
Radius= 9.88215034E+08
log10(L/Lsun)= -2.93919790E+00
Mass= 2.42489235E+22
Radius= 9.89470007E+08
create initial model iteration 64
log10(L/Lsun)= -2.93265176E+00
Mass= 2.18215454E+22
Radius= 9.88423225E+08
log10(L/Lsun)= -2.90620138E+00
Mass= 2.23524064E+22
Radius= 1.00596120E+09
log10(L/Lsun)= -2.93152434E+00
Mass= 2.20745215E+22
Radius= 9.89688755E+08
create initial model iteration 65
log10(L/Lsun)= -2.92517029E+00
Mass= 1.98671503E+22
Radius= 9.88650385E+08
log10(L/Lsun)= -2.89886771E+00
Mass= 2.03581140E+22
Radius= 1.00630744E+09
log10(L/Lsun)= -2.92403951E+00
Mass= 2.00976845E+22
Radius= 9.89921778E+08
create initial model iteration 66
log10(L/Lsun)= -2.89174979E+00
Mass= 1.85434008E+22
Radius= 1.00663521E+09
log10(L/Lsun)= -2.91794382E+00
Mass= 1.80891355E+22
Radius= 9.88848630E+08
log10(L/Lsun)= -2.91679836E+00
Mass= 1.82992453E+22
Radius= 9.90124950E+08
create initial model iteration 67
log10(L/Lsun)= -2.91095440E+00
Mass= 1.64712922E+22
Radius= 9.89007948E+08
log10(L/Lsun)= -2.88484729E+00
Mass= 1.68917344E+22
Radius= 1.00693547E+09
log10(L/Lsun)= -2.90979937E+00
Mass= 1.66628222E+22
Radius= 9.90290823E+08
create initial model iteration 68
log10(L/Lsun)= -2.87812490E+00
Mass= 1.53886281E+22
Radius= 1.00722289E+09
log10(L/Lsun)= -2.90416815E+00
Mass= 1.49993797E+22
Radius= 9.89143188E+08
log10(L/Lsun)= -2.90300694E+00
Mass= 1.51739782E+22
Radius= 9.90432665E+08
create initial model iteration 69
log10(L/Lsun)= -2.89757161E+00
Mass= 1.36600771E+22
Radius= 9.89252327E+08
log10(L/Lsun)= -2.87156940E+00
Mass= 1.40205512E+22
Radius= 1.00749686E+09
log10(L/Lsun)= -2.89640866E+00
Mass= 1.38192582E+22
Radius= 9.90549414E+08
create initial model iteration 70
log10(L/Lsun)= -2.89114188E+00
Mass= 1.24414401E+22
Radius= 9.89341990E+08
log10(L/Lsun)= -2.86515930E+00
Mass= 1.27753362E+22
Radius= 1.00776374E+09
log10(L/Lsun)= -2.88998099E+00
Mass= 1.25865821E+22
Radius= 9.90647644E+08
create initial model iteration 71
log10(L/Lsun)= -2.88489348E+00
Mass= 1.13325825E+22
Radius= 9.89423031E+08
log10(L/Lsun)= -2.85887515E+00
Mass= 1.16418971E+22
Radius= 1.00802892E+09
log10(L/Lsun)= -2.88370197E+00
Mass= 1.14649383E+22
Radius= 9.90733154E+08
create initial model iteration 72
log10(L/Lsun)= -2.87885030E+00
Mass= 1.03230294E+22
Radius= 9.89457780E+08
log10(L/Lsun)= -2.85275614E+00
Mass= 1.06095362E+22
Radius= 1.00825397E+09
log10(L/Lsun)= -2.87760855E+00
Mass= 1.04437314E+22
Radius= 9.90769710E+08
create initial model iteration 73
log10(L/Lsun)= -2.87301644E+00
Mass= 9.40358630E+21
Radius= 9.89424977E+08
log10(L/Lsun)= -2.84681637E+00
Mass= 9.66888396E+21
Radius= 1.00841679E+09
log10(L/Lsun)= -2.87171712E+00
Mass= 9.51367212E+21
Radius= 9.90738400E+08
create initial model iteration 74
log10(L/Lsun)= -2.86736928E+00
Mass= 8.56621950E+21
Radius= 9.89327628E+08
log10(L/Lsun)= -2.84103926E+00
Mass= 8.81179273E+21
Radius= 1.00852036E+09
log10(L/Lsun)= -2.86607216E+00
Mass= 8.66661325E+21
Radius= 9.90650508E+08
create initial model iteration 75
log10(L/Lsun)= -2.86179296E+00
Mass= 7.80432644E+21
Radius= 9.89239290E+08
log10(L/Lsun)= -2.83535004E+00
Mass= 8.03162190E+21
Radius= 1.00864959E+09
log10(L/Lsun)= -2.86051430E+00
Mass= 7.89589183E+21
Radius= 9.90574560E+08
create initial model iteration 76
log10(L/Lsun)= -2.85625234E+00
Mass= 7.11118588E+21
Radius= 9.89177518E+08
log10(L/Lsun)= -2.82972634E+00
Mass= 7.32153714E+21
Radius= 1.00882437E+09
log10(L/Lsun)= -2.85499346E+00
Mass= 7.19471017E+21
Radius= 9.90525978E+08
create initial model iteration 77
log10(L/Lsun)= -2.82411729E+00
Mass= 6.67509113E+21
Radius= 1.00903273E+09
log10(L/Lsun)= -2.85074287E+00
Mass= 6.48046644E+21
Radius= 9.89138161E+08
log10(L/Lsun)= -2.84950519E+00
Mass= 6.55666426E+21
Radius= 9.90500651E+08
create initial model iteration 78
log10(L/Lsun)= -2.81850996E+00
Mass= 6.08649484E+21
Radius= 1.00927587E+09
log10(L/Lsun)= -2.84530070E+00
Mass= 5.90647891E+21
Radius= 9.89130957E+08
log10(L/Lsun)= -2.84403753E+00
Mass= 5.97600702E+21
Radius= 9.90501417E+08
create initial model iteration 79
log10(L/Lsun)= -2.81296993E+00
Mass= 5.55006765E+21
Radius= 1.00949689E+09
log10(L/Lsun)= -2.84000443E+00
Mass= 5.38366203E+21
Radius= 9.89105603E+08
log10(L/Lsun)= -2.83865526E+00
Mass= 5.44709827E+21
Radius= 9.90473348E+08
create initial model iteration 80
log10(L/Lsun)= -2.80763034E+00
Mass= 5.06049544E+21
Radius= 1.00960355E+09
log10(L/Lsun)= -2.83490491E+00
Mass= 4.90681693E+21
Radius= 9.88964940E+08
log10(L/Lsun)= -2.83357102E+00
Mass= 4.96468731E+21
Radius= 9.90346213E+08
create initial model iteration 81
log10(L/Lsun)= -2.80235119E+00
Mass= 4.61470187E+21
Radius= 1.00977220E+09
log10(L/Lsun)= -2.82977957E+00
Mass= 4.47284445E+21
Radius= 9.88860692E+08
log10(L/Lsun)= -2.82847440E+00
Mass= 4.52564236E+21
Radius= 9.90258247E+08
create initial model iteration 82
log10(L/Lsun)= -2.79700960E+00
Mass= 4.20874284E+21
Radius= 1.00999314E+09
log10(L/Lsun)= -2.82460810E+00
Mass= 4.07786972E+21
Radius= 9.88800328E+08
log10(L/Lsun)= -2.82333272E+00
Mass= 4.12604400E+21
Radius= 9.90214912E+08
create initial model iteration 83
log10(L/Lsun)= -2.81950286E+00
Mass= 3.71833144E+21
Radius= 9.88805999E+08
log10(L/Lsun)= -2.79159578E+00
Mass= 3.83899524E+21
Radius= 1.01026695E+09
log10(L/Lsun)= -2.81813833E+00
Mass= 3.76229454E+21
Radius= 9.90218352E+08
create initial model iteration 84
log10(L/Lsun)= -2.81457006E+00
Mass= 3.39024466E+21
Radius= 9.88697665E+08
log10(L/Lsun)= -2.78636915E+00
Mass= 3.50136541E+21
Radius= 1.01042524E+09
log10(L/Lsun)= -2.81319447E+00
Mass= 3.43035390E+21
Radius= 9.90119957E+08
create initial model iteration 85
log10(L/Lsun)= -2.80962450E+00
Mass= 3.09136053E+21
Radius= 9.88588469E+08
log10(L/Lsun)= -2.78122102E+00
Mass= 3.19360249E+21
Radius= 1.01060289E+09
log10(L/Lsun)= -2.80828361E+00
Mass= 3.12795355E+21
Radius= 9.90028960E+08
create initial model iteration 86
log10(L/Lsun)= -2.80460352E+00
Mass= 2.81928624E+21
Radius= 9.88544965E+08
log10(L/Lsun)= -2.77594925E+00
Mass= 2.91328018E+21
Radius= 1.01084563E+09
log10(L/Lsun)= -2.80326789E+00
Mass= 2.85267363E+21
Radius= 9.89998814E+08
create initial model iteration 87
log10(L/Lsun)= -2.79965814E+00
Mass= 2.57141225E+21
Radius= 9.88539751E+08
log10(L/Lsun)= -2.77062536E+00
Mass= 2.65773558E+21
Radius= 1.01110706E+09
log10(L/Lsun)= -2.79823289E+00
Mass= 2.60187695E+21
Radius= 9.89990530E+08
create initial model iteration 88
log10(L/Lsun)= -2.79482362E+00
Mass= 2.34514405E+21
Radius= 9.88418060E+08
log10(L/Lsun)= -2.76555775E+00
Mass= 2.42431753E+21
Radius= 1.01126921E+09
log10(L/Lsun)= -2.79343742E+00
Mass= 2.37293413E+21
Radius= 9.89888024E+08
create initial model iteration 89
log10(L/Lsun)= -2.78984010E+00
Mass= 2.13916219E+21
Radius= 9.88371871E+08
log10(L/Lsun)= -2.76031839E+00
Mass= 2.21170214E+21
Radius= 1.01150776E+09
log10(L/Lsun)= -2.78849088E+00
Mass= 2.16451350E+21
Radius= 9.89860967E+08
create initial model iteration 90
log10(L/Lsun)= -2.78487298E+00
Mass= 1.95159286E+21
Radius= 9.88426881E+08
log10(L/Lsun)= -2.75496549E+00
Mass= 2.01797770E+21
Radius= 1.01182556E+09
log10(L/Lsun)= -2.78341719E+00
Mass= 1.97471992E+21
Radius= 9.89908075E+08
create initial model iteration 91
log10(L/Lsun)= -2.78004427E+00
Mass= 1.78021512E+21
Radius= 9.88320700E+08
log10(L/Lsun)= -2.74987009E+00
Mass= 1.84088121E+21
Radius= 1.01199005E+09
log10(L/Lsun)= -2.77863234E+00
Mass= 1.80130740E+21
Radius= 9.89822011E+08
create initial model iteration 92
log10(L/Lsun)= -2.77504958E+00
Mass= 1.62419591E+21
Radius= 9.88312749E+08
log10(L/Lsun)= -2.74454455E+00
Mass= 1.67956903E+21
Radius= 1.01224208E+09
log10(L/Lsun)= -2.77364194E+00
Mass= 1.64343461E+21
Radius= 9.89827376E+08
create initial model iteration 93
log10(L/Lsun)= -2.77006740E+00
Mass= 1.48201126E+21
Radius= 9.88364315E+08
log10(L/Lsun)= -2.73924120E+00
Mass= 1.53248853E+21
Radius= 1.01255102E+09
log10(L/Lsun)= -2.76859885E+00
Mass= 1.49955677E+21
Radius= 9.89878377E+08
create initial model iteration 94
log10(L/Lsun)= -2.76510302E+00
Mass= 1.35219787E+21
Radius= 9.88319561E+08
log10(L/Lsun)= -2.73396627E+00
Mass= 1.39815363E+21
Radius= 1.01275363E+09
log10(L/Lsun)= -2.76368291E+00
Mass= 1.36819628E+21
Radius= 9.89854168E+08
create initial model iteration 95
log10(L/Lsun)= -2.76001519E+00
Mass= 1.23399850E+21
Radius= 9.88401112E+08
log10(L/Lsun)= -2.72854021E+00
Mass= 1.27578618E+21
Radius= 1.01307545E+09
log10(L/Lsun)= -2.75851088E+00
Mass= 1.24858629E+21
Radius= 9.89929593E+08
create initial model iteration 96
log10(L/Lsun)= -2.75504559E+00
Mass= 1.12601085E+21
Radius= 9.88356575E+08
log10(L/Lsun)= -2.72323482E+00
Mass= 1.16396664E+21
Radius= 1.01325601E+09
log10(L/Lsun)= -2.75358612E+00
Mass= 1.13930981E+21
Radius= 9.89904363E+08
create initial model iteration 97
log10(L/Lsun)= -2.74986518E+00
Mass= 1.02767652E+21
Radius= 9.88425267E+08
log10(L/Lsun)= -2.71773463E+00
Mass= 1.06211246E+21
Radius= 1.01355010E+09
log10(L/Lsun)= -2.74838617E+00
Mass= 1.03980153E+21
Radius= 9.89979775E+08
create initial model iteration 98
log10(L/Lsun)= -2.74469407E+00
Mass= 9.37976126E+20
Radius= 9.88501667E+08
log10(L/Lsun)= -2.71219147E+00
Mass= 9.69188804E+20
Radius= 1.01383203E+09
log10(L/Lsun)= -2.74320270E+00
Mass= 9.49029157E+20
Radius= 9.90063046E+08
create initial model iteration 99
log10(L/Lsun)= -2.70658267E+00
Mass= 8.84428103E+20
Radius= 1.01411224E+09
log10(L/Lsun)= -2.73939185E+00
Mass= 8.56160926E+20
Radius= 9.88571074E+08
log10(L/Lsun)= -2.73792224E+00
Mass= 8.66236389E+20
Radius= 9.90146109E+08
create initial model iteration 100
log10(L/Lsun)= -2.70081549E+00
Mass= 8.07169064E+20
Radius= 1.01444679E+09
log10(L/Lsun)= -2.73403410E+00
Mass= 7.81588877E+20
Radius= 9.88721867E+08
log10(L/Lsun)= -2.73252199E+00
Mass= 7.90773334E+20
Radius= 9.90296751E+08
create initial model iteration 101
log10(L/Lsun)= -2.72858548E+00
Mass= 7.13519269E+20
Radius= 9.88803434E+08
log10(L/Lsun)= -2.69506343E+00
Mass= 7.36656038E+20
Radius= 1.01472396E+09
log10(L/Lsun)= -2.72710363E+00
Mass= 7.21891137E+20
Radius= 9.90392992E+08
create initial model iteration 102
log10(L/Lsun)= -2.72305071E+00
Mass= 6.51487999E+20
Radius= 9.88981603E+08
log10(L/Lsun)= -2.68907761E+00
Mass= 6.72403464E+20
Radius= 1.01506178E+09
log10(L/Lsun)= -2.72152672E+00
Mass= 6.59120202E+20
Radius= 9.90570764E+08
create initial model iteration 103
log10(L/Lsun)= -2.68312788E+00
Mass= 6.13779054E+20
Radius= 1.01535965E+09
log10(L/Lsun)= -2.71740844E+00
Mass= 5.94876349E+20
Radius= 9.89099217E+08
log10(L/Lsun)= -2.71589724E+00
Mass= 6.01834209E+20
Radius= 9.90699080E+08
create initial model iteration 104
log10(L/Lsun)= -2.67692767E+00
Mass= 5.60351996E+20
Radius= 1.01567719E+09
log10(L/Lsun)= -2.71168211E+00
Mass= 5.43271575E+20
Radius= 9.89277585E+08
log10(L/Lsun)= -2.71015255E+00
Mass= 5.49615307E+20
Radius= 9.90880585E+08
create initial model iteration 105
log10(L/Lsun)= -2.70580603E+00
Mass= 4.96205739E+20
Radius= 9.89439912E+08
log10(L/Lsun)= -2.67070646E+00
Mass= 5.11640924E+20
Radius= 1.01600173E+09
log10(L/Lsun)= -2.70425200E+00
Mass= 5.01990293E+20
Radius= 9.91044656E+08
create initial model iteration 106
log10(L/Lsun)= -2.66434054E+00
Mass= 4.67225996E+20
Radius= 1.01627606E+09
log10(L/Lsun)= -2.69988594E+00
Mass= 4.53274444E+20
Radius= 9.89579482E+08
log10(L/Lsun)= -2.69835830E+00
Mass= 4.58549965E+20
Radius= 9.91196613E+08
create initial model iteration 107
log10(L/Lsun)= -2.65772594E+00
Mass= 4.26783146E+20
Radius= 1.01662575E+09
log10(L/Lsun)= -2.69373126E+00
Mass= 4.14169349E+20
Radius= 9.89798381E+08
log10(L/Lsun)= -2.69214438E+00
Mass= 4.18981899E+20
Radius= 9.91408956E+08
create initial model iteration 108
log10(L/Lsun)= -2.65111421E+00
Mass= 3.89898419E+20
Radius= 1.01688156E+09
log10(L/Lsun)= -2.68752726E+00
Mass= 3.78487078E+20
Radius= 9.89918338E+08
log10(L/Lsun)= -2.68600701E+00
Mass= 3.82878392E+20
Radius= 9.91550180E+08
create initial model iteration 109
log10(L/Lsun)= -2.64403629E+00
Mass= 3.56354505E+20
Radius= 1.01727531E+09
log10(L/Lsun)= -2.68103086E+00
Mass= 3.46025826E+20
Radius= 9.90210577E+08
log10(L/Lsun)= -2.67945238E+00
Mass= 3.50034568E+20
Radius= 9.91836270E+08
create initial model iteration 110
log10(L/Lsun)= -2.67440021E+00
Mass= 3.16424709E+20
Radius= 9.90406333E+08
log10(L/Lsun)= -2.63696923E+00
Mass= 3.25781892E+20
Radius= 1.01759998E+09
log10(L/Lsun)= -2.67285094E+00
Mass= 3.20085551E+20
Radius= 9.92045303E+08
create initial model iteration 111
log10(L/Lsun)= -2.66755410E+00
Mass= 2.89490062E+20
Radius= 9.90703100E+08
log10(L/Lsun)= -2.62952345E+00
Mass= 2.97974871E+20
Radius= 1.01799235E+09
log10(L/Lsun)= -2.66598526E+00
Mass= 2.92834843E+20
Radius= 9.92344002E+08
create initial model iteration 112
log10(L/Lsun)= -2.66048304E+00
Mass= 2.64967995E+20
Radius= 9.90983934E+08
log10(L/Lsun)= -2.62183319E+00
Mass= 2.72670852E+20
Radius= 1.01836567E+09
log10(L/Lsun)= -2.65885918E+00
Mass= 2.68025659E+20
Radius= 9.92618265E+08
create initial model iteration 113
log10(L/Lsun)= -2.61406324E+00
Mass= 2.49644509E+20
Radius= 1.01867123E+09
log10(L/Lsun)= -2.65329890E+00
Mass= 2.42641665E+20
Radius= 9.91189323E+08
log10(L/Lsun)= -2.65174118E+00
Mass= 2.45438783E+20
Radius= 9.92845770E+08
create initial model iteration 114
log10(L/Lsun)= -2.60567958E+00
Mass= 2.28758576E+20
Radius= 1.01913381E+09
log10(L/Lsun)= -2.64568319E+00
Mass= 2.22382569E+20
Radius= 9.91582730E+08
log10(L/Lsun)= -2.64406737E+00
Mass= 2.24943578E+20
Radius= 9.93232503E+08
create initial model iteration 115
log10(L/Lsun)= -2.63784666E+00
Mass= 2.03970565E+20
Radius= 9.91888569E+08
log10(L/Lsun)= -2.59716441E+00
Mass= 2.09786130E+20
Radius= 1.01953046E+09
log10(L/Lsun)= -2.63627844E+00
Mass= 2.06317682E+20
Radius= 9.93556337E+08
create initial model iteration 116
log10(L/Lsun)= -2.62964554E+00
Mass= 1.87297114E+20
Radius= 9.92358603E+08
log10(L/Lsun)= -2.58807393E+00
Mass= 1.92611907E+20
Radius= 1.02005461E+09
log10(L/Lsun)= -2.62803510E+00
Mass= 1.89450837E+20
Radius= 9.94023890E+08
create initial model iteration 117
log10(L/Lsun)= -2.62114305E+00
Mass= 1.72202153E+20
Radius= 9.92786827E+08
log10(L/Lsun)= -2.57864714E+00
Mass= 1.77070060E+20
Radius= 1.02053530E+09
log10(L/Lsun)= -2.61951013E+00
Mass= 1.74181363E+20
Radius= 9.94455217E+08
create initial model iteration 118
log10(L/Lsun)= -2.61234441E+00
Mass= 1.58585861E+20
Radius= 9.93262555E+08
log10(L/Lsun)= -2.56882932E+00
Mass= 1.63055593E+20
Radius= 1.02105093E+09
log10(L/Lsun)= -2.61072855E+00
Mass= 1.60407961E+20
Radius= 9.94942974E+08
create initial model iteration 119
log10(L/Lsun)= -2.55834402E+00
Mass= 1.50482887E+20
Radius= 1.02165553E+09
log10(L/Lsun)= -2.60306119E+00
Mass= 1.46366418E+20
Radius= 9.93852244E+08
log10(L/Lsun)= -2.60139782E+00
Mass= 1.48047715E+20
Radius= 9.95529961E+08
create initial model iteration 120
log10(L/Lsun)= -2.54753264E+00
Mass= 1.39268884E+20
Radius= 1.02227759E+09
log10(L/Lsun)= -2.59353181E+00
Mass= 1.35464226E+20
Radius= 9.94462486E+08
log10(L/Lsun)= -2.59189211E+00
Mass= 1.37020266E+20
Radius= 9.96154528E+08
create initial model iteration 121
log10(L/Lsun)= -2.53603287E+00
Mass= 1.29393100E+20
Radius= 1.02302969E+09
log10(L/Lsun)= -2.58351369E+00
Mass= 1.25860669E+20
Radius= 9.95225377E+08
log10(L/Lsun)= -2.58187734E+00
Mass= 1.27306568E+20
Radius= 9.96927454E+08
create initial model iteration 122
log10(L/Lsun)= -2.57333590E+00
Mass= 1.17563172E+20
Radius= 9.96136404E+08
log10(L/Lsun)= -2.52416572E+00
Mass= 1.20862012E+20
Radius= 1.02390047E+09
log10(L/Lsun)= -2.57164028E+00
Mass= 1.18914030E+20
Radius= 9.97833935E+08
create initial model iteration 123
log10(L/Lsun)= -2.51247334E+00
Mass= 1.13752771E+20
Radius= 1.02478387E+09
log10(L/Lsun)= -2.56339331E+00
Mass= 1.10647800E+20
Radius= 9.97070410E+08
log10(L/Lsun)= -2.56170748E+00
Mass= 1.11919674E+20
Radius= 9.98779754E+08
create initial model iteration 124
log10(L/Lsun)= -2.55417002E+00
Mass= 1.05309781E+20
Radius= 9.98133239E+08
log10(L/Lsun)= -2.50149272E+00
Mass= 1.08265881E+20
Radius= 1.02578316E+09
log10(L/Lsun)= -2.55246866E+00
Mass= 1.06520896E+20
Radius= 9.99848868E+08
create initial model iteration 125
log10(L/Lsun)= -2.54696498E+00
Mass= 1.01821645E+20
Radius= 9.99133495E+08
log10(L/Lsun)= -2.49284262E+00
Mass= 1.04680929E+20
Radius= 1.02672648E+09
log10(L/Lsun)= -2.54524846E+00
Mass= 1.02993067E+20
Radius= 1.00085239E+09
create initial model iteration 126
log10(L/Lsun)= -2.48831306E+00
Mass= 1.03100475E+20
Radius= 1.02736462E+09
log10(L/Lsun)= -2.54320468E+00
Mass= 1.00283358E+20
Radius= 9.99798652E+08
log10(L/Lsun)= -2.54148607E+00
Mass= 1.01437333E+20
Radius= 1.00152082E+09
create initial model iteration 127
log10(L/Lsun)= -2.48731461E+00
Mass= 1.02816772E+20
Radius= 1.02754767E+09
log10(L/Lsun)= -2.54237069E+00
Mass= 1.00007014E+20
Radius= 9.99985840E+08
log10(L/Lsun)= -2.54065705E+00
Mass= 1.01157872E+20
Radius= 1.00171007E+09
Model data written to model_output.txt
done build_initial_model
relax_num_steps 10
hydro_solver_step returned ierr -1
s% model_number 2
nz 403
s% num_retries 0
retry: logT < hydro_mtx_min_allowed_logT 374 2
solver rejected trial model
s% model_number 2
s% solver_call_number 3
nz 403
s% num_retries 1
dt 1.2623040000000002D+02
log dt/secyer -5.3979400086720375D+00
retry: avg corr too large -- give up in solver 2
solver rejected trial model
s% model_number 2
s% solver_call_number 4
nz 403
s% num_retries 2
dt 6.3115200000000009D+01
log dt/secyer -5.6989700043360179D+00
retry: avg corr too large -- give up in solver 2
solver rejected trial model
s% model_number 2
s% solver_call_number 5
nz 403
s% num_retries 3
dt 3.1557600000000004D+01
log dt/secyer -5.9999999999999991D+00
retry: max residual jumped -- give up in solver 2
solver rejected trial model
s% model_number 2
s% solver_call_number 6
nz 403
s% num_retries 4
dt 1.5778800000000002D+01
log dt/secyer -6.3010299956639804D+00
retry: avg corr too large -- give up in solver 2
hydro_solver_step returned ierr -1
s% model_number 2
nz 403
s% num_retries 5
retry: logRho > hydro_mtx_max_allowed_logRho 390 2
solver rejected trial model
s% model_number 2
s% solver_call_number 8
nz 403
s% num_retries 6
dt 3.9447000000000005D+00
log dt/secyer -6.9030899869919429D+00
retry: avg corr too large -- give up in solver 2
solver rejected trial model
s% model_number 2
s% solver_call_number 9
nz 403
s% num_retries 7
dt 1.9723500000000003D+00
log dt/secyer -7.2041199826559241D+00
retry: avg corr too large -- give up in solver 2
solver rejected trial model
s% model_number 2
s% solver_call_number 10
nz 403
s% num_retries 8
dt 9.8617500000000013D-01
log dt/secyer -7.5051499783199054D+00
retry: avg corr too large -- give up in solver 2
solver rejected trial model
s% model_number 2
s% solver_call_number 11
nz 403
s% num_retries 9
dt 4.9308750000000007D-01
log dt/secyer -7.8061799739838866D+00
retry: avg corr too large -- give up in solver 2
solver rejected trial model
s% model_number 2
s% solver_call_number 12
nz 403
s% num_retries 10
dt 2.4654375000000003D-01
log dt/secyer -8.1072099696478670D+00
retry: avg corr too large -- give up in solver 2
solver rejected trial model
s% model_number 2
s% solver_call_number 13
nz 403
s% num_retries 11
dt 1.2327187500000002D-01
log dt/secyer -8.4082399653118483D+00
retry: avg corr too large -- give up in solver 2
hydro_solver_step returned ierr -1
s% model_number 2
nz 403
s% num_retries 12
retry: get_T_tau -- L <= 0 2
hydro_solver_step returned ierr -1
s% model_number 2
nz 403
s% num_retries 13
retry: get_T_tau -- L <= 0 2
hydro_solver_step returned ierr -1
s% model_number 2
nz 403
s% num_retries 14
retry: get_T_tau -- L <= 0 2
hydro_solver_step returned ierr -1
s% model_number 2
nz 403
s% num_retries 15
retry: get_T_tau -- L <= 0 2
hydro_solver_step returned ierr -1
s% model_number 2
nz 403
s% num_retries 16
retry: get_T_tau -- L <= 0 2
hydro_solver_step returned ierr -1
s% model_number 2
nz 403
s% num_retries 17
retry: get_T_tau -- L <= 0 2
hydro_solver_step returned ierr -1
s% model_number 2
nz 403
s% num_retries 18
retry: get_T_tau -- L <= 0 2
solver rejected trial model
s% model_number 2
s% solver_call_number 21
nz 403
s% num_retries 19
dt 4.8153076171875007D-04
log dt/secyer -1.0816479930623698D+01
retry: avg corr too large -- give up in solver 2
solver rejected trial model
s% model_number 2
s% solver_call_number 22
nz 403
s% num_retries 20
dt 2.4076538085937503D-04
log dt/secyer -1.1117509926287680D+01
retry: avg corr too large -- give up in solver 2
solver rejected trial model
s% model_number 2
s% solver_call_number 23
nz 403
s% num_retries 21
dt 1.2038269042968752D-04
log dt/secyer -1.1418539921951661D+01
retry: avg corr too large -- give up in solver 2
solver rejected trial model
s% model_number 2
s% solver_call_number 24
nz 403
s% num_retries 22
dt 6.0191345214843758D-05
log dt/secyer -1.1719569917615642D+01
retry: avg corr too large -- give up in solver 2
solver rejected trial model
s% model_number 2
s% solver_call_number 25
nz 403
s% num_retries 23
dt 3.0095672607421879D-05
log dt/secyer -1.2020599913279623D+01
retry: avg corr too large -- give up in solver 2
solver rejected trial model
s% model_number 2
s% solver_call_number 26
nz 403
s% num_retries 24
dt 1.5047836303710940D-05
log dt/secyer -1.2321629908943605D+01
retry: avg corr too large -- give up in solver 2
solver rejected trial model
s% model_number 2
s% solver_call_number 27
nz 403
s% num_retries 25
dt 7.5239181518554698D-06
log dt/secyer -1.2622659904607586D+01
retry: avg corr too large -- give up in solver 2
solver rejected trial model
s% model_number 2
s% solver_call_number 28
nz 403
s% num_retries 26
dt 3.7619590759277349D-06
log dt/secyer -1.2923689900271567D+01
retry: avg corr too large -- give up in solver 2
solver rejected trial model
s% model_number 2
s% solver_call_number 29
nz 403
s% num_retries 27
dt 1.8809795379638674D-06
log dt/secyer -1.3224719895935548D+01
retry: avg corr too large -- give up in solver 2
dt 9.4048976898193372D-07
min_timestep_limit 9.9999999999999995D-07
stopping because of problems dt < min_timestep_limit
failed in relax num steps
failed in do_relax_num_steps
star_create_initial_model ierr -1
do_load1_star ierr -1
do_before_evolve_loop ierr -1
do_before_evolve_loop ierr -1
DATE: 2024-12-05
TIME: 14:48:53
******************************************************************
failed to create final.mod when running inlist_make_brown_dwarf_header
******************************************************************
More information about the Mesa-users
mailing list