mesa-users] problem with running MESA

Miroslav Mocak mmocak at
Wed Sep 9 09:37:13 EDT 2009

Dear MESA users,

  I have a problem with running the latest version of
MESA compiled by ifort (version 11) compiler on my machine which has
the Intel I7 processor and runs Fedora 8. After
instalation (which went ok), I tried to run the default
problem in the /mesa/star/test directory by ./rn, but the
code stopped by giving me the following message 
nvar /= s% nvar eval_equ . Could you help me to resolve
the issue? Thanks. Below you'll find the whole text output
which I get until the code stops.


Wed Sep  9 15:29:15 CEST 2009
 loading eos data
 loading kap data
 finished loading

 The terminal output contains the following information:

      'step' is the number of steps since the start of the run,
      'lg_dt_yr' 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_Pcntr' is log10 center pressure (ergs/cm^3),
      'lg_Teff' is log10 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_LZ' is log10 total burning power excluding LH and L3a and 
photodisintegrations (Lsun),
      'lg_LNuc' is log10 nuclear power excluding photodisintegration and 
neutrinos (Lsun),
      'lg_LNeu' is log10 total neutrino power (Lsun),
      'lg_Psurf' is log10 surface pressure (gas + radiation),
      'Mass' is the total stellar 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_rich' is the remaining mass outside of the hydrogen poor core,
      'H_poor' is the core mass where hydrogen abundance is <= 0.37E+00
      'He_poor' is the core mass where helium abundance is <= 0.47E+00
      '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,
      'X_avg' is the star average hydrogen mass fraction,
      'Y_avg' is the star average helium mass fraction,
      'Z_avg' is the star average metallicity,
      'max_avg' is the star total mass fraction of the heaviest isotope,
      'gam_cntr' is the center plasma interaction parameter,
      'eta_cntr' is the center electron degeneracy parameter,
      'pts' is the number of grid points in the current model,
      'jacs' is the number of jacobians created in the current step,
      'retry' is the number of step retries required during the run,
      'bckup' is the number of step backups required during the run,
      'dt_limit' is an indication of what limited the timestep.

 All this and more are saved in 'LOGS/star.log' during the run.

 load starting model from ../starting_models


       step    lg_Tcntr    lg_Teff    lg_LH     lg_Lnuc     Mass       
H_rich     H_cntr     N_cntr     X_surf     X_avg     eta_cntr   pts  retry
   lg_dt_yr    lg_Dcntr   lg_Rsurf    lg_L3a    lg_Lneu     lg_Mdot    
H_poor     He_cntr    O_cntr     Z_surf     Y_avg     gam_cntr  jacs  bckup
     age_yr    lg_Pcntr    lg_L       lg_LZ     lg_Psurf    lg_Dsurf   
He_poor    C_cntr     Ne_cntr    Z_cntr     Z_avg     max_avg      dt_limit

 nvar /= s% nvar



      Miroslav Mocak

      Institut d'Astronomie et d'Astrophysique
      Universit\E9 Libre de Bruxelles
      Campus De La Plaine
      CP 226 - Boulevard Du Triomphe
      B-1050 Bruxelles - BELGIUM

      -- -- -- --  -- -- -- -- -- -- -- -- --  -- -- -- -- -- -- -- --

      -- -- -- --  -- -- -- -- -- -- -- -- --  -- -- -- -- -- -- -- --

       Understanding is, after all, what science is all about - and science
      is a great deal more than mere mindless computation (Roger Penrose)


More information about the Mesa-users mailing list