[mesa-users] solar_calibration

[kuldeep verma]

hi all,
   I'm a new user of mesa star. I have been trying to get the solar
parameter (e.g. mixing length, Y and  Z/X), by running the test program
solar_calibration (mesa/star/test_suite/solar_calibration) to construct a
solar model. The problem is, it stops the run without any error even if the
value of surf_Z_div_X has not converged to the specified accuracy (the
value I get surf_Z_div_X=0.02505594313202 ). Luminosity and radius have
converged to the specified accuracy. The contribution to the chi2
(chi2=1.52168809042131) is mostly from surf_Z_div_X. Edited
inlist_calibration_controls (inlist_solar is untouched) is given below. Any
suggestion or comment will be appreciated.

&calibration_controls

   !search_type = 'just_do_first_values'
      ! run 1 case only using "first" values defined below
   search_type = 'bobyqa'
      ! search for combination minimizing chi^2 using BOBYQA algorithm
   !search_type = 'hooke'
      ! minimize chi^2 using Hooke-Jeeves algorithm

   ! targets

   target_log_luminosity = 0
   sigma_log_luminosity = 1d-5

   target_log_radius = 0
   sigma_log_radius = 1d-5

   target_surface_Z_div_X = 2.292d-2 ! GS98 value
   !target_surface_Z_div_X = 1.81d-2 ! Asplund 09 value
   sigma_surface_Z_div_X = 1d-3

   ! OPTIONAL EXTRA TERMS FOR CHI^2
   target_Rcz = 0.713d0 ! Bahcall, Serenelli, Basu, 2005
   !sigma_Rcz = 1d-3 ! use this if no overshooting
   sigma_Rcz = 0.003d0 ! use this if are using overshooting
   use_Rcz_in_search = .false.

   target_surf_He = 0.2485d0 ! Bahcall, Serenelli, Basu, 2005
   sigma_surf_He = 0.0034
   use_sigma_surf_He_in_search = .false.

   target_csound_rms = 0
   sigma_csound_rms = 1d-5
   use_csound_rms_in_search = .false.


   ! parameters
   first_Y = 0.27992069626717d0
   min_Y = 0.2784 ! 0.277d0
   max_Y = 0.2804 ! 0.283d0
   vary_Y = .true. ! else use first_Y

   first_Z_div_X = 0.02895042291518d0
   min_Z_div_X = 0.0275 ! 0.025d0
   max_Z_div_X = 0.0296 ! 0.033d0
   vary_Z_div_X = .true. ! else use first_Z_div_X

   first_alpha = 1.91053122947443d0
   min_alpha = 1.90 ! 1.84d0
   max_alpha = 1.94 ! 1.96d0
   vary_alpha = .true. ! else use first_alpha

   first_f_ov = 0.01157256164430d0
   min_f_ov = 0.0107 ! 0.000
   max_f_ov = 0.0127 ! 0.015
   ! set sigma_Rcz according to whether or not using overshooting
   !vary_f_ov = .true. ! else use first_f_ov
   vary_f_ov = .false.

   first_f0_ov_fraction = -1 ! disable this
   min_f0_ov_fraction = 0.0
   max_f0_ov_fraction = 0.9
   !vary_f0_ov_fraction = .true. ! else use first_f0_ov_fraction
   vary_f0_ov_fraction = .false.
      ! if true, then change overshoot_f0_below_nonburn along with other
parameters
      ! if false, then set overshoot_f0_below_nonburn to
first_f0_ov_fraction*first_f_ov

   Y_frac_he3 = 1d-4


   max_tries = 500 ! limit for BOBYQA
   rhoend = 1d-7 ! convergence parameter for BOBYQA

   hooke_itermax = 100
   hooke_rho = 0.50
   hooke_eps = 0.05
   hooke_xscale = 0.75

/ ! end of calibration_controls



cheers,
kuldeep verma
department of astronomy and astrophysics
tata institute of fundamental research
mumbai, india
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120829/cab9a064/attachment.html>