[mesa-users] astero...

Warrick Ball wball at bison.ph.bham.ac.uk
Tue May 16 11:43:01 EDT 2017


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
wball at bison.ph.bham.ac.uk
+44 (0)121 414 4552

On Tue, 16 May 2017, bruno.moura.pesquisa at gmail.com 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 <wball at bison.ph.bham.ac.uk>:
>       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
>       wball at bison.ph.bham.ac.uk
>       +44 (0)121 414 4552
>
>       On Tue, 16 May 2017, bruno.moura.pesquisa at gmail.com 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 <wball at bison.ph.bham.ac.uk>:
>                   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
>                   wball at bison.ph.bham.ac.uk
>                   +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
> 
> 
> 
> 
> 
>


More information about the Mesa-users mailing list