[Mesa-users] modeling a planet

[Nuno Moedas]

Hi,

I am trying to model a planet of about 2 Jupiter masses using 
mesa-r12778. I used as starting point the test suit that is available in 
MESA. I was able to create the model and add the core and evolve it. Now 
I am trying to create a grid where a change the mass of the core from 0 
until reaches the total mass of the planet. However when I reach more or 
less half the mass of the planet in the core the model does not converge 
in the step of adding the core. I was trying to change the parameters of 
dlg_core_mass_per_step and relax_core_years_for_dt, but it seems that 
still gives me the same error.
  stopping because of problems num_backups > max_number_backups


This is my first time modelling a planet. I want to reach a point where 
the planet only is a core and does not present any mass in the envelope. 
Now I don't know what to do to pass half of the mass of the planet on 
the core, or if I missing any command in it.
I send also the inlists that I use to create the planet, add the core of 
the planet, and evolve the planet.


Thank you.

Sincerely,
Nuno Moedas
-------------- next part --------------
&star_job
	save_model_when_terminate = .true.
	save_model_filename = "creat_planet_Mp_2.024_mc_340_rc_10_Z_0.0144072_Y_0.2637283.mod"

	create_initial_model = .true.
	radius_in_cm_for_create_initial_model = 20414012000.0
	mass_in_gm_for_create_initial_model = 1.8107904d+30
	initial_model_relax_num_steps = 0

	set_initial_age = .true.
	initial_age = 0.d0 ! in years
	eos_file_prefix = "mesa"
	initial_zfracs = 6
	kappa_file_prefix = "a09"
	kappa_lowT_prefix = "lowT_Freedman11"

	change_lnPgas_flag = .true.
	new_lnPgas_flag = .true.

	profile_columns_file = "profile.list"
	history_columns_file = "history.list"

/ !end of star_job

&controls
	steps_before_use_gold_tolerances = 50
	warn_when_stop_checking_residuals = .false. ! off until evolve
	warn_when_large_virial_thm_rel_err = 1d99 ! off until evolve

! check for retries and backups as part of test_suite
	max_number_backups = 0
	max_number_retries = 0
	max_model_number = 100

! for create_initial_model
	initial_Z = 0.0144072
	initial_Y = 0.2637283

! when to stop
	max_age = 1.d3
	max_model_number = 1000 

! knobs
	T_mix_limit = 0
	max_resid_jump_limit = 1d18 
	max_corr_jump_limit = 1d18 
	photo_interval = 50
	max_num_profile_models=2000
	history_interval = 2
	terminal_interval = 10
	write_header_frequency = 10
/
-------------- next part --------------
&star_job
	load_saved_model = .true.
	saved_model_name = "creat_planet_Mp_2.024_mc_340_rc_10_Z_0.0144072_Y_0.2637283.mod"

	save_model_when_terminate = .true.
	save_model_filename = "core_planet_Mp_2.024_mc_340_rc_10_Z_0.0144072_Y_0.2637283.mod"

	relax_core = .true.
	new_core_mass = 0.0010215041222602051
	core_avg_rho = 10
	core_avg_eps = 0
	dlg_core_mass_per_step = 0.001
	relax_core_years_for_dt = 0.1

	eos_file_prefix = "mesa"
	initial_zfracs = 6
	kappa_file_prefix = "a09"
	kappa_lowT_prefix = "lowT_Freedman11"

	change_lnPgas_flag = .true.
	new_lnPgas_flag = .true.

	profile_columns_file = "profile.list"
	history_columns_file = "history.list"

/ !end of star_job

&controls
	initial_Z = 0.0144072
	initial_Y = 0.2637283

! when to stop
	max_age = 2.d3
	max_model_number = 1000 

! knobs
	T_mix_limit = 0
	photo_interval = 50
	max_num_profile_models=2000
	history_interval = 1
	terminal_interval = 10
	write_header_frequency = 10
/
-------------- next part --------------
&star_job
	load_saved_model = .true.
	saved_model_name = "core_planet_Mp_2.024_mc_340_rc_10_Z_0.0144072_Y_0.2637283.mod"

	save_model_when_terminate = .true.
	save_model_filename = "planet_Mp_2.024_mc_340_rc_10_Z_0.0144072_Y_0.2637283.mod"

	eos_file_prefix = "mesa"
	initial_zfracs = 6
	kappa_file_prefix = "a09"
	kappa_lowT_prefix = "lowT_Freedman11"

	change_lnPgas_flag = .true.
	new_lnPgas_flag = .true.

	profile_columns_file = "profile.list"
	history_columns_file = "history.list"

/ !end of star_job

&controls
	initial_Z = 0.0144072
	initial_Y = 0.2637283

! surface heating
	column_depth_for_irradiation = 300
	irradiation_flux = 8236372915.813784
	! when to stop
	max_age = 4000000000.0
	max_years_for_timestep = 4d7
	max_model_number = 2000 
	Teff_lower_limit = 100 

! knobs
	T_mix_limit = 0
	photo_interval = 50
	profile_interval = 5
	max_num_profile_models=2000
	history_interval = 2
	terminal_interval = 10
	write_header_frequency = 10
/