[mesa-users] solar_calibration
Bill Paxton
paxton at kitp.ucsb.edu
Wed Aug 29 12:28:33 EDT 2012
Hi,
I just did a quick check and solar calibration seems to work on my machine.
Let me give the details of what I did in case that might help.
1) copy star/test_suite/solar_calibration to a new top level directory (not part of the mesa directory)
2) edit make/makefile to delete the line MESA_DIR =
this assumes that the environment variable MESA_DIR has been set
3) edit inlist to set mesa_dir = ''
again, this will direct the code to use the MESA_DIR environment variable
4) ./clean; ./mk; ./rn
I stopped the run by hand after it had generated 10 test cases using BOBYQA.
Here are the results for those 10:
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: bobyqa_results.txt
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120829/28ced2fc/attachment.txt>
-------------- next part --------------
Here are the inlists etc.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist
Type: application/octet-stream
Size: 294 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120829/28ced2fc/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist_calibration_controls
Type: application/octet-stream
Size: 2400 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120829/28ced2fc/attachment-0001.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist_solar
Type: application/octet-stream
Size: 3206 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120829/28ced2fc/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: solar_calibration_input.mod
Type: application/octet-stream
Size: 852256 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120829/28ced2fc/attachment-0003.obj>
-------------- next part --------------
I'm running a pre-release version (4437) -- but I doubt if that is your problem.
Good luck,
Bill
On Aug 29, 2012, at 3:52 AM, kuldeep verma wrote:
> 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
> ------------------------------------------------------------------------------
> Live Security Virtual Conference
> Exclusive live event will cover all the ways today's security and
> threat landscape has changed and how IT managers can respond. Discussions
> will include endpoint security, mobile security and the latest in malware
> threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/_______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users
More information about the Mesa-users
mailing list