From: Warrick B. <wb...@bi...> - 2017-05-16 15:43:10
|
Hi Bruno, The astero/work folder has all the inlists and code that you need. Try copying it to a new location, then running ./mk and ./rn, and see what happens! W ------------ Warrick Ball Postdoc, School of Physics and Astronomy University of Birmingham, Edgbaston, Birmingham B15 2TT wb...@bi... +44 (0)121 414 4552 On Tue, 16 May 2017, bru...@gm... wrote: > Thanks again. It means then that all the templates I want to do I must transcribe to the astero module and not include the $ MESA_DIR / star / work inlists. > > > Bruno > > > 2017-05-16 12:34 GMT-03:00 Warrick Ball <wb...@bi...>: > Hi Bruno, > > If you want to work with the astero module, you should use the folder $MESA_DIR/star/astero/work instead of the usual work folder ($MESA_DIR/star/work). That is, whereas one would usually copy > $MESA_DIR/star/work to a new location and start editing the inlists, you should instead copy $MESA_DIR/star/astero/work. Does that make more sense? > > Cheers, > Warrick > > PS: Keep replies on the list so that they're available for future readers (and so that someone else might help you if I'm unavailable/slow to respond). > > > ------------ > Warrick Ball > Postdoc, School of Physics and Astronomy > University of Birmingham, Edgbaston, Birmingham B15 2TT > wb...@bi... > +44 (0)121 414 4552 > > On Tue, 16 May 2017, bru...@gm... wrote: > > Thanks also. One question: should I just include in the folder astro_serarch? > > Cheers > > Bruno > 2017-05-16 11:21 GMT-03:00 Warrick Ball <wb...@bi...>: > Hi Bruno, > > Not sure if someone has already replied to this but my guess is that you've added the astero stuff to &controls in the usual MESA inlist, rather than to > &astero_search_controls. The astero > module requires extra code that is compiled in addition to the usual MESA code, and that this code has its own set of extra parameters in two namelists called > &astero_search_controls and > &astero_pgstar_controls in the file inlist_astero_search_controls. > > Cheers, > Warrick > > ------------ > Warrick Ball > Postdoc, School of Physics and Astronomy > University of Birmingham, Edgbaston, Birmingham B15 2TT > wb...@bi... > +44 (0)121 414 4552 > > On Tue, 16 May 2017, Bruno Lustosa wrote: > > Hello everyone > > I have a tense question about my model with regard to astero controls. By including chi2_seismo_fraction directly in my inlist in & controls the following: > > ! Seismic data > Chi2_seismo_fraction = 0.667d0 > Include_Teff_in_chi2_spectro = .true. > Teff_target = 4400d0 > Teff_sigma = 100d0 > Include_logg_in_chi2_spectro = .true. > Logg_target = 2.5d0 > Logg_sigma = 0.02d0 > Include_logL_in_chi2_spectro = .true. > LogL_target = 1.62d0 > LogL_sigma = 0.04d0 > Include_FeH_in_chi2_spectro = .true. ! [Fe / H] > FeH_target = 0.1 > FeH_sigma = 0.05 > Include_age_in_chi2_spectro = .true. > Age_target = 8.87d9! (See Bahcall, Serenelli, and Basu, 2006) > Age_sigma = 0.0065d9 > Num_smaller_steps_before_age_target = 50 > Dt_for_smaller_steps_before_age_target = 0.0065d8 > Chi2_seismo_delta_nu_fraction = 0d0 > ! If> 0 then delta_nu and delta_nu_sigma must be set (see below) > Chi2_seismo_nu_max_fraction = 0d0 > ! If> 0 then nu_max and nu_max_sigma must be set (see below) > Chi2_seismo_r_010_fraction = 0d0 > ! If> 0, then include r_010 frequency ratios > Chi2_seismo_r_02_fraction = 0d0 > ! If> 0, then include r_02 frequency ratios > ! Fraction for frequencies = 1 - (frac_r_010_ratios + frac_r_02_ratios + frac_delta_nu + frac_nu_max) > ! Nu_max = 70 > ! Nu_max_sigma = -1 > Delta_nu = -1 > Delta_nu_sigma = -1.0 > Eval_chi2_at_target_age_only = .true. > Search_type = 'use_first_values' > Simplex_itermax = 1000! Each iteration revises the simplex > Simplex_fcn_calls_max = 10000! May use several function calls per iteration > Simplex_x_atol = 1d-10! Tolerance for absolute differences > Simplex_x_rtol = 1d-10! Tolerance for relative differences > Simplex_centroid_weight_power = 0d0 > Simplex_output_filename = '16_maio.data' > Simplex_seed = 1074698122 > Oscillation_code = 'adipls > > The following error is reported: > > At line 1700 of file ../private/ctrls_io.f90 (unit = 30, file = 'inlist_2.0') > Fortran runtime error: Can not match namelist object name chi2_seismo_fraction > > > Have I or should I change some subroutine I did not do? > > > Cheers > > > Bruno > > > > > > |