[Mesa-users] Problems when running MESA program

Thu Feb 24 14:04:50 UTC 2022

hi nie,

the issue appears to be the model, not the installation of mesa,
but not enough information is provided for anyone to potentially help.
share your inlist and any custom run_star_extras.

fxt

> On Feb 24, 2022, at 1:37 AM, 聂昊明 via Mesa-users <mesa-users at lists.mesastar.org> wrote:
> 
> Greetings,
> 
> Firstly, I appreciate the wonderful program you provided to all of us, but I've met some problems running MESA.
> 
> MESA version: r15140
> MESASDK version: x86_64-macos-21.11.1
> Computer info: MacBook Pro (13-inch, 2019, Four Thunderbolt 3 ports)
> OS: macOS 10.15.7 (Catalina)
> 
> uname -a
> 
> Darwin niehaomingdeMacBook-Pro.local 19.6.0 Darwin Kernel Version 19.6.0: Mon Apr 12 20:57:45 PDT 2021; root:xnu-6153.141.28.1~1/RELEASE_X86_64 x86_64
> 
> gfortran -v
> 
> Using built-in specs.
> COLLECT_GCC=/Applications/mesasdk/bin/gfortran.exec
> COLLECT_LTO_WRAPPER=/Applications/mesasdk/bin/../libexec/gcc/x86_64-apple-darwin21.1.0/10.2.0/lto-wrapper
> Target: x86_64-apple-darwin21.1.0
> Configured with: /opt/sdk2-tmp/build/gcc/configure CC=clang CXX=clang++ --build=x86_64-apple-darwin21.1.0 --host=x86_64-apple-darwin21.1.0 --target=x86_64-apple-darwin21.1.0 --prefix=/opt/sdk2-tmp/mesasdk --with-gmp=/opt/sdk2-tmp/mesasdk --with-mpfr=/opt/sdk2-tmp/mesasdk --with-mpc=/opt/sdk2-tmp/mesasdk --enable-languages=c,c++,fortran --disable-multilib --disable-nls --disable-libsanitizer --with-sysroot=/opt/sdk2-tmp/mesasdk/sysroot
> Thread model: posix
> Supported LTO compression algorithms: zlib
> gcc version 10.2.0 (GCC) 
> 
> echo $MESASDK_ROOT
> 
> /Applications/mesasdk
> 
> echo $PATH
> 
> /Users/niehaoming/opt/anaconda3/bin:/Users/niehaoming/opt/anaconda3/condabin:/Applications/mesasdk/bin:/usr/local/opt/qt/bin:/usr/local/Cellar/geant4/geant4_install/bin:/usr/local/Cellar/root/6.22.02_2/bin:/bin:/usr/bin:/usr/local/bin
> 
> echo $MESA_DIR
> 
> /Users/niehaoming/mesa-r15140
> 
> I believe that installation is successful, with all enviroment variables set properly. Then I copied the /star/work folder to a new place and made with ./mk then ran with ./rn , the output was a little different from what was expected:
> 
> DATE: 2022-02-24
> TIME: 16:32:17
>  version_number 15140
>  read inlist_project
> 
> 
>  The terminal output contains the following information
> 
>       'step' is the number of steps since the start of the run,
>       'lg_dt' is log10 timestep in years,
>       'age_yr' is the simulated years since the start run,
>       'lg_Tcntr' is log10 center temperature (K),
>       'lg_Dcntr' is log10 center density (g/cm^3),
>       'lg_Tmax' is log10 max temperature (K),
>       'Teff' is the surface temperature (K),
>       'lg_R' is log10 surface radius (Rsun),
>       'lg_L' is log10 surface luminosity (Lsun),
>       'lg_LH' is log10 total PP and CNO hydrogen burning power (Lsun),
>       'lg_L3a' is log10 total triple-alpha helium burning power (Lsun),
>       'lg_gsurf' is log10 surface gravity,
>       'lg_LNuc' is log10 nuclear power (Lsun),
>       'lg_LNeu' is log10 total neutrino power (Lsun),
>       'lg_Lphoto' is log10 total photodisintegration (Lsun),
>       'Mass' is the total stellar baryonic mass (Msun),
>       'lg_Mdot' is log10 magnitude of rate of change of mass (Msun/year),
>       'lg_Dsurf' is log10 surface density (g/cm^3),
>       'H_env' is the amount of mass where H is the most abundant iso,
>       'He_core' is the largest mass where He is most abundant iso.
>       'C_core' is the largest mass where C is most abundant iso.
>       'H_cntr' is the center H1 mass fraction,
>       'He_cntr' is the center He4 mass fraction,
>       'C_cntr' is the center C12 mass fraction,
>       'N_cntr' is the center N14 mass fraction,
>       'O_cntr' is the center O16 mass fraction,
>       'Ne_cntr' is the center Ne20 mass fraction,
>       'gam_cntr' is the center plasma interaction parameter,
>       'eta_cntr' is the center electron degeneracy parameter,
>       'zones' is the number of zones in the current model,
>       'iters' is the number of solver iterations for the current step,
>       'retry' is the number of step retries required during the run,
>       'dt_limit' is an indication of what limited the timestep.
> 
>  All this and more are saved in the LOGS directory during the run.
>  create pre-main-sequence model
>                                           species  mass           8    1.5000000000000002D+01
> 
>                                         relax_num_steps         300
> 
>  retry: get_T_tau -- L <= 0       1
>           1   5.479738   3311.350 -21.748047 -21.748047  15.000000  15.000000   0.700000   0.001010   0.280000 -12.926136   1078      1
>   -5.301030   5.479738   2.521737 -99.000000 -23.145915 -99.000000   0.000000   0.280000   0.009381   0.020000   0.013983      8
>  5.0000E-06  -5.414988   4.078188 -99.000000 -99.000000  -8.223243   0.000000   0.003447   0.002085   0.020000  0.000E+00         retry
> 
>           2   5.478886   3347.984 -21.748142 -21.748142  15.000000  15.000000   0.700000   0.001010   0.280000 -12.926202    610      1
>   -5.397940   5.478886   2.521754 -99.000000 -23.146010 -99.000000   0.000000   0.280000   0.009381   0.020000   0.013996      3
>  9.0000E-06  -5.416294   4.097334 -99.000000 -99.000000  -8.242255   0.000000   0.003447   0.002085   0.020000  0.000E+00  max decrease
> 
> Then the program stuck with no further output, and it never stopped running after that. I'm not sure what does "max decrease" mean. I'll be thankful if you could provide any help.
> 
> Thank you.
> NIE Haoming, from Tsinghua University
> _______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
> 