[mesa-users] Additional chi2 variables in astero

[Bill Paxton]

On Feb 21, 2014, at 7:29 AM, Warrick Ball wrote:

> Bill, João,
> 
> I happen to be looking at a similar idea (i.e. adding variables for the fit), and I'm wondering if this fix will work for me too.  In particular, I wanted to add the large and small frequency separations, computed from the power spectrum via ADIPLS or GYRE.  Can I use the same framework?
> 
> I presume I can look at the "pulse" test case to see how to incorporate ADIPLS into run_star_extras, then use the frequencies to compute the large and small separations and assign them to my_var?.  Does that make sense?
> 
> Cheers,
> Warrick



Hi,

You should be able to do that.  The delta_nu value is available without needing to call adipls.  But calling adipls or gyre from run_star_extras is of course an option.

Cheers,
Bill







> 
> On Wed, 19 Feb 2014, Bill Paxton wrote:
> 
>> Hi João,
>> I've added user defined chi2 variables for you.
>> They'll show up in the next release, but you can have a pre-release if you'd like.
>> You could try 5995 -- let me know if you run into problems with it.
>> These are the new controls in the astero_search inlist:
>>       ! to include user-defined variables in chi^2 spectro
>>       ! set the "my_var" variable in your extras_check_model routine.
>>       !     use astero_data, only: my_var1, my_var2, my_var3
>>       !     my_var1 = ......
>>       include_my_var1_in_chi2_spectro = .false.
>>       my_var1_target = 0
>>       my_var1_sigma = 0
>>       my_var1_name = 'my_var1'
>>       
>>       include_my_var2_in_chi2_spectro = .false.
>>       my_var2_target = 0
>>       my_var2_sigma = 0
>>       my_var2_name = 'my_var2'
>>       
>>       include_my_var3_in_chi2_spectro = .false.
>>       my_var3_target = 0
>>       my_var3_sigma = 0
>>       my_var3_name = 'my_var3'
>> Cheers,
>> Bill
>> On Feb 19, 2014, at 8:35 AM, João Faria wrote:
>> 
>>      Hi Dennis,
>> 
>>      Sorry, I had to check with my colleague. We know the mean density independently from transit analysis
>>      and have spectroscopic parameters whose determination we trust like the effective temperature and
>>      metallicity. For the surface gravity, not so much.
>>      The goal would be to find the best model and determine its log(g). There are no seismic observations
>>      unfortunately...
>> 
>>      Cheers,
>>      João
>> 
>>      On 02/13/2014 02:41 AM, Dennis Stello wrote:
>>      Hi Jo�o
>> 
>>      In regards to the mean density you were requesting, is the delta_nu not good enough, or are you
>>      thinking of independent density input from a
>>      binary system?
>> 
>>      Cheers
>>      Dennis/
>> 
>>      On Wed, 12 Feb 2014, João Faria wrote:
>> 
>>            Hi everyone,
>> 
>>            Here I go again bugging people with requests on the astero "package"...
>> 
>>            Recently I've been confronted with the question of where to start if we want to do a
>>            spectroscopic-only chi2 model fit in MESA. Note that spectroscopic-only might mean
>>            including surface He, convection zone radius, etc, besides the usual spectroscopic
>>            observables.
>>            In the old days (like last year), I would say start from the solar_calibration test.
>>            But now (v5819) every spectroscopic variable is also in the example_astero test, so
>>            the answer might be this test with /chi2_seismo_fraction=0./
>> 
>>            On to the request: I've been asked if we can include the stellar mean density as a
>>            constraint in the chi2. So something like /include_rhoMean_in_chi2_spectro = .true.
>>            /Also, and this might be more difficult in practice, how would I go about including
>>            a new variable in the chi2 from run_star_extras.f instead of sending an email asking
>>            for it to be changed in the guts of MESA?
>>            Is it feasible something like
>>            /    include_newVar1_in_chi2_spectro = .true./
>>                newVar1_target = -1
>>                newVar1_sigma = -1
>> 
>>               ... newVar2 ...
>> 
>>            in the inlist and then somehow calculating the model's newVar1, newVar2 in
>>            run_star_extras.f ?
>> 
>>            Thanks in advance for your help.
>> 
>>            Cheers,
>>            João Faria
>> ------------------------------------------------------------------------------
>> Android apps run on BlackBerry 10
>> Introducing the new BlackBerry 10.2.1 Runtime for Android apps.
>> Now with support for Jelly Bean, Bluetooth, Mapview and more.
>> Get your Android app in front of a whole new audience.  Start now.
>> http://pubads.g.doubleclick.net/gampad/clk?id=124407151&iu=/4140/ostg.clktrk
>> _______________________________________________
>> mesa-users mailing list
>> mesa-users at lists.sourceforge.net
>> https://lists.sourceforge.net/lists/listinfo/mesa-users
>> ------------------------------------------------------------------------------
>> Managing the Performance of Cloud-Based Applications
>> Take advantage of what the Cloud has to offer - Avoid Common Pitfalls.
>> Read the Whitepaper.
>> http://pubads.g.doubleclick.net/gampad/clk?id=121054471&iu=/4140/ostg.clktrk_____________________________________
>> __________
>> mesa-users mailing list
>> mesa-users at lists.sourceforge.net
>> https://lists.sourceforge.net/lists/listinfo/mesa-users
>> 
> 
> 
> ------------
> Warrick Ball
> Postdoc, Institut für Astrophysik Göttingen
> wball at astro.physik.uni-goettingen.de
> +49 (0) 551 39 5069