[Mesa-users] error when accreting material similar to earth(high z)
李宗恒
lzh17513652589 at 163.com
Thu Aug 7 14:35:16 UTC 2025
hello mesa_users:
following is my issue
mesa-r23051
mesasdk-x86_64-linux-23.7.3.tar
when I try to make the sun-like star to accrete earth-like materials(several earth-masses),some error come out
retry: adjust_correction failed in eval_equations -- give up in solver 159
retry: get_T_tau -- L <= 0 226
the following is the composition I used
if (s% star_age >1d7 .and. s% mstar < m_limit .and. s% x_logical_ctrl(3))then
write(,)"accreting**"
s% accrete_same_as_surface = .false.
s% accrete_given_mass_fractions = .true.
s% num_accretion_species = 4
! s% accretion_species_id(1) = 'li7'
! s% accretion_species_xa(1) = 0.5
s% accretion_species_id(1) = 'o16'
s% accretion_species_xa(1) = 0.3
s% accretion_species_id(2) = 'mg24'
s% accretion_species_xa(2) = 0.2
s% accretion_species_id(3) = 'si28'
s% accretion_species_xa(3) = 0.2
s% accretion_species_id(4) = 'fe56'
s% accretion_species_xa(4) = 0.3
s% mass_change = 1d-11
write(,)"mstar=,m_limit=",s% mstar,m_limit
write(,)"accreting*"
else
write(,)"mstar=",s% mstar
s% mass_change = 0
endif
But if I use the abundance below ,everything is fine
s% num_accretion_species = 4
s% accretion_species_id(1) = 'h1'
s% accretion_species_xa(1) = 0.75
s% accretion_species_id(2) = 'he4'
s% accretion_species_xa(2) = 0.2
s% accretion_species_id(3) = 'li7'
s% accretion_species_xa(3) = 1d-9
s% accretion_species_id(4) = 'fe56'
s% accretion_species_xa(4) = 0.04999999d0
s% mass_change = 1d-11
Linux LL 6.6.87.2-microsoft-standard-WSL2 #1 SMP PREEMPT_DYNAMIC Thu Jun 5 18:30:46 UTC 2025 x86_64 x86_64 x86_64 GNU/Linux
li at LL:~$ gfortran -v
Using built-in specs.
COLLECT_GCC=/home/li/mesasdk/bin/gfortran.exec
COLLECT_LTO_WRAPPER=/home/li/mesasdk/bin/../libexec/gcc/x86_64-pc-linux-gnu/13.1.0/lto-wrapper
Target: x86_64-pc-linux-gnu
Configured with: /home/user/sdk2-tmp/build/gcc/configure CC= CXX= --build=x86_64-pc-linux-gnu --host=x86_64-pc-linux-gnu --target=x86_64-pc-linux-gnu --prefix=/home/user/sdk2-tmp/mesasdk --with-gmp=/home/user/sdk2-tmp/mesasdk --with-mpfr=/home/user/sdk2-tmp/mesasdk --with-mpc=/home/user/sdk2-tmp/mesasdk --enable-languages=c,c++,fortran --disable-multilib --disable-nls --disable-libsanitizer --enable-clocale=generic
Thread model: posix
Supported LTO compression algorithms: zlib
gcc version 13.1.0 (GCC)
Thanks a lot
Jason Lee
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20250807/170ef405/attachment.htm>
More information about the Mesa-users
mailing list