[mesa-users] Error while running the code for 1.4M☉, 1.6M☉, 1.8M☉, 2.0M☉

Maya Mahira mayamahira89 at yahoo.com
Tue Mar 6 03:58:13 EST 2012


Dear all,



Based on the subject above, I'm having difficulty to find the right way to make correction to the code. I have set the star age as 4.65 x 10^9 but when running the code, it cannot achieve until that years then it will give the notice that I'm having a convergence problem.

Does anyone knows how to solve it?

Here, I attached the inlist that has been used for running the code.


Regard,
Mahirah Mohd Idris




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

!inlist_Mahirah

&star_job

      mesa_data_dir = '/home/teori/Mahirah/mesa3/data'
      
      create_pre_main_sequence_model = .true.

      initial_zfracs = 2
         
      set_rates_preference = .false. ! for use by net + rates modules
      new_rates_preference = 1 
         ! 1 = NACRE rates -- this is the default
         ! 2 = jina reaclib rates -- to match jina where possible
         
      set_rate_c12ag = 'NACRE' ! empty string means ignore this control
         ! one of 'NACRE', 'Buchmann', 'Kunz', or 'CF88'
      set_rate_n14pg = 'NACRE' ! empty string means ignore this control
         ! one of 'NACRE', 'Imbriani', or 'CF88'
      set_rate_3a = 'NACRE' ! empty string means ignore this control
         ! one of 'NACRE', 'Fynbo', 'CF88', or 'FL87'
      set_rate_1212_to_ne20 = '' ! empty string means ignore this control
         ! one of 'CF88' or 'G05'
    
! solar model calibration
      
         !run_calibration = .true. ! calibrate a solar model
         !calibration_inlist = 'calibration_zahra'
        
         kappa_file_prefix = 'gn93'

/ ! end of star_job namelist


&controls

! starting specifications

      initial_mass =  1.40
      initial_z = 0.01d0 
      initial_y = 0.2755

      photostep = 1
      profile_interval = 1
      log_cnt = 10
      terminal_cnt = 10
      write_header_frequency = 10

! when  to stop
     
      max_model_number = 100000    
      max_age = 4.65d9
    
! mixing parameters     

      mixing_length_alpha =2.1895
      use_Henyey_MLT = .true.
! element diffusion

      do_element_diffusion = .true. 

! timesteps
      
      max_timestep = 3.15d14  
      varcontrol_target = 1d-4 

!trace_mass_location = 0.3

      !report_hydro_solver_progress = .true. ! set true to see info about newton iterations
      !report_ierr = .false. ! if true, produce terminal output when have some internal error

      !RGB_wind_scheme = 'Reimers'    
      !AGB_wind_scheme = 'Blocker'        
      !RGB_to_AGB_wind_switch = 1d-4    
      !Reimers_wind_eta = 0.5d0      
      !Blocker_wind_eta = 0.08d0     

      !CO_full_on = 0    
      !CO_full_off = 0    
      !use_CO_enhanced_opacities = .false.    
     
/ ! end of controls namelist
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120306/8f67415f/attachment.html>


More information about the Mesa-users mailing list