<html> <head> <meta content="text/html; charset=ISO-8859-1" http-equiv="Content-Type"> </head> <body text="#000000" bgcolor="#FFFFFF"> <div class="moz-cite-prefix">Hi,<br> Just to have in mind orders of magnitude for a 2.5Msun star (early A-type, solar Z), the 1D model gives a good estimate of the polar radius (2% err), it underestimates polar temperature (500K over 11000K) and overestimates luminosity (10%), when the rotation of the star is 95% critical; if rotation is less errors are less.<br> <br> best,<br> <br> Michel<br> <pre class="moz-signature" cols="72">***************************************************************************** Michel Rieutord Institut de Recherche en Astrophysique et Planétologie Observatoire Midi-Pyrénées email : <a class="moz-txt-link-abbreviated" href="mailto:Michel.Rieutord@irap.omp.eu">Michel.Rieutord@irap.omp.eu</a> 14 avenue Edouard Belin Tel +33 -(0)5 61 33 29 49 31400 Toulouse, France Fax +33 -(0)5 61 33 28 40 *****************************************************************************</pre> Le 06/10/2013 22:06, Jeremy Jones a écrit :<br> </div> <blockquote cite="mid:CAD9nXrZ2cHc+D0J3F=9aCYTg_VQ+rRKxpabZcUv-nZKuwC6ekA@mail.gmail.com" type="cite"> <div dir="ltr">Thanks for all of your responses. I'll definitely take a look at the ESTER code. <div><br> </div> <div>In the meanwhile, I was thinking I would compare the polar radius, equatorial velocity, and luminosity calculated by my model fitting to the interferometric data to those same outputs from mesa (now that I know the polar radius parameter exists... I must have missed it when I was looking through all the parameters) in order to get the age/mass of the stars I observed. Since I've got a handful of stars that should be coeval (UMa moving group, some slow rotators (v_eq <100 km/s) and some rapid rotators), it'll be interesting to see the results I get out of mesa even if it is a 1-D code. Ideally, I'll be using as many evolution models as I can get my hands on for my thesis, but there are only so many hours in the day...</div> <div><br> </div> <div>~Jeremy</div> <div><br> </div> <div> <div style="color:rgb(136,136,136);font-family:arial,sans-serif;font-size:13.333333015441895px">-- </div> <div style="color:rgb(136,136,136);font-family:arial,sans-serif;font-size:13.333333015441895px">Jeremy Jones</div> <div style="color:rgb(136,136,136);font-family:arial,sans-serif;font-size:13.333333015441895px">Graduate Student Assistant</div> <div style="color:rgb(136,136,136);font-family:arial,sans-serif;font-size:13.333333015441895px">Georgia State University</div> <div style="color:rgb(136,136,136);font-family:arial,sans-serif;font-size:13.333333015441895px">Department of Physics & Astronomy</div> <div style="color:rgb(136,136,136);font-family:arial,sans-serif;font-size:13.333333015441895px"><a moz-do-not-send="true" href="mailto:jones@chara.gsu.edu" target="_blank">jones@chara.gsu.edu</a></div> </div> </div> <div class="gmail_extra"><br> <br> <div class="gmail_quote">On Sun, Oct 6, 2013 at 7:16 AM, Michel Rieutord <span dir="ltr"><<a moz-do-not-send="true" href="mailto:mrieutord@irap.omp.eu" target="_blank">mrieutord@irap.omp.eu</a>></span> wrote:<br> <blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"> <div bgcolor="#FFFFFF" text="#000000"> <div>Le 05/10/2013 14:18, Michel Rieutord a écrit :<br> </div> <div class="im"> <blockquote type="cite"> <div>Hi Jeremy,<br> <br> I think the best answer to your question is to be found with the ESTER <br> code that you can freely download at<br> <a moz-do-not-send="true" href="http://code.google.com/p/ester-project/" target="_blank">http://code.google.com/p/ester-project/</a><br> <br> This is a code that computes the 2D structure, including differential<br> rotation and meridional circulation, of rotating early-type stars at <br> any rotation rate (below critical of course). With this code you'll be <br> able to compare your data with those of the model including the <br> temperature dependence with latitude; by adjusting the hyrdogen content<br> of the core you may get an idea of the evolution status of your star on<br> the main sequence. This won't give you an age in years but this can <br> obtained with a 1D code afterwards. We are working out the time <br> evolution but this is not yet ready; this a tough problem all the more<br> that we need to also evolve the dynamics...<br> <br> You may find examples of the results in<br> Espinosa Lara & Rieutord 2013, A&A, 552, A35<br> Rieutord & Espinosa Lara 2013, LNP 865, 49, Arxive 1208.4929<br> or 2013arXiv:<a moz-do-not-send="true" href="tel:1310.0908" value="+13100908" target="_blank">1310.0908</a> or arXiv:<a moz-do-not-send="true" href="tel:1310.1017" value="+13101017" target="_blank">1310.1017</a><br> <br> Don't hesitate if you have questions.<br> <br> Best wishes,<br> <br> Michel Rieutord<br> <br> </div> </blockquote> </div> PS third attempt; hope this is the good one. <div> <div class="h5"><br> <blockquote type="cite"> <div> <pre cols="72">************************************************************************** Michel Rieutord Institut de Recherche en Astrophysique et Planétologie Observatoire Midi-Pyrénées email : <a moz-do-not-send="true" href="mailto:Michel.Rieutord@irap.omp.eu" target="_blank">Michel.Rieutord@irap.omp.eu</a> 14 avenue Edouard Belin Tel <a moz-do-not-send="true" href="tel:%2B33%20-%280%295%2061%2033%2029%2049" value="+33561332949" target="_blank">+33 -(0)5 61 33 29 49</a> 31400 Toulouse, France Fax <a moz-do-not-send="true" href="tel:%2B33%20-%280%295%2061%2033%2028%2040" value="+33561332840" target="_blank">+33 -(0)5 61 33 28 40</a> **************************************************************************</pre> Le 04/10/2013 22:25, Jeremy Jones a écrit :<br> </div> <blockquote type="cite"> <div dir="ltr">Hey mesa users, <div><br> </div> <div>I'm working on a project concerning interferometry of A-stars. Many of my targets are rapid rotators, and I'm looking into using mesa to determine their ages. I've played around with mesa a bit (I've included one of my inlists below). I have a question about how to compare my observables with the mesa outputs, in particular radius and temperature.</div> <div><br> </div> <div>From the modeling I do of the inerferometric data, I get the full radius profile of the star (i.e. radius as a function of latitude). Likewise, I get the full temperature profile of the star.</div> <div> <br> </div> <div>My question is for a rapid rotator, which radius/temperature should I compare with the mesa outputs to get the age of the star? Polar radius/temperature, equatorial radius/temperature, something more involved?</div> <div><br> </div> <div>Thanks for any insight!</div> <div><br> </div> <div>~Jeremy</div> <div><br> </div> <div>The aforementioned inlist:</div> <div>================</div> <div> <div>&star_job</div> <div><br> </div> <div> mesa_dir = '../..'</div> <div><br> </div> <div> set_initial_age = .true.</div> <div> initial_age = 0 ! in years</div> <div><br> </div> <div> set_initial_model_number = .true.</div> <div> initial_model_number = 0</div> <div><br> </div> <div> create_pre_main_sequence_model = .true.</div> <div> </div> <div> new_omega_div_omega_crit = 0.5 </div> <div> near_zams_relax_omega_div_omega_crit = .true.</div> <div><br> </div> <div> pgstar_flag = .true.</div> <div><br> </div> <div><br> </div> <div>/ ! end of star_job namelist</div> <div><br> </div> <div><br> </div> <div>&controls</div> <div><br> </div> <div> min_timestep_limit = 1d-12 ! (seconds)</div> <div><br> </div> <div> ! check for retries and backups as part of test_suite</div> <div> ! you can delete this for use outside of test_suite</div> <div> max_number_backups = 400</div> <div> max_number_retries = 900</div> <div> max_model_number = 1500</div> <div><br> </div> <div> !max_allowed_nz = 20000 ! maximum number of grid points allowed</div> <div><br> </div> <div> initial_mass = 2</div> <div> initial_z = 0.016</div> <div><br> </div> <div> photostep = 10</div> <div> profile_interval = 10</div> <div> history_interval = 1</div> <div> terminal_interval = 2</div> <div> write_header_frequency = 10</div> <div><br> </div> <div>/ ! end of controls namelist</div> <div><br> </div> <div><br> </div> <div>&pgstar</div> <div><br> </div> <div><span style="white-space:pre-wrap"> </span>History_Track1_win_flag = .true.</div> <div><span style="white-space:pre-wrap"> </span>History_Track1_yname = 'log_R'</div> <div><span style="white-space:pre-wrap"> </span>History_Track1_yaxis_label = 'log(R/R_sun)'</div> <div><span style="white-space:pre-wrap"> </span>History_Track1_xname = 'log_Teff'</div> <div><span style="white-space:pre-wrap"> </span>History_Track1_xaxis_label = 'log(T/K)'</div> <div><span style="white-space:pre-wrap"> </span>History_Track1_title = 'Radius vs. Effective Temperature'</div> <div><br> </div> <div><span style="white-space:pre-wrap"> </span>History_Track2_win_flag = .true.</div> <div><span style="white-space:pre-wrap"> </span>History_Track2_yname = 'surf_avg_v_rot'</div> <div><span style="white-space:pre-wrap"> </span>History_Track2_yaxis_label = 'Equatorial Velocity (km/s)'</div> <div><span style="white-space:pre-wrap"> </span>History_Track2_xname = 'star_age'</div> <div><span style="white-space:pre-wrap"> </span>History_Track2_xaxis_label = 'Age (yr)'</div> <div><span style="white-space:pre-wrap"> </span>History_Track2_title = 'Surface Velocity vs. Age'</div> <div><br> </div> <div>/ ! end of pgstar namelist</div> </div> <div>================<br> </div> <div><br> </div> <div> <div>-- </div> <div>Jeremy Jones</div> <div>Graduate Student Assistant</div> <div>Georgia State University</div> <div> Department of Physics & Astronomy</div> <div><a moz-do-not-send="true" href="mailto:jones@chara.gsu.edu" target="_blank">jones@chara.gsu.edu</a></div> </div> </div> <br> <fieldset></fieldset> <br> <pre>------------------------------------------------------------------------------ October Webinars: Code for Performance Free Intel webinars can help you accelerate application performance. Explore tips for MPI, OpenMP, advanced profiling, and more. Get the most from the latest Intel processors and coprocessors. See abstracts and register > <a moz-do-not-send="true" href="http://pubads.g.doubleclick.net/gampad/clk?id=60134791&iu=/4140/ostg.clktrk" target="_blank">http://pubads.g.doubleclick.net/gampad/clk?id=60134791&iu=/4140/ostg.clktrk</a></pre> <br> <fieldset></fieldset> <br> <pre>_______________________________________________ mesa-users mailing list <a moz-do-not-send="true" href="mailto:mesa-users@lists.sourceforge.net" target="_blank">mesa-users@lists.sourceforge.net</a> <a moz-do-not-send="true" href="https://lists.sourceforge.net/lists/listinfo/mesa-users" target="_blank">https://lists.sourceforge.net/lists/listinfo/mesa-users</a> </pre> </blockquote> <br> </blockquote> <br> </div> </div> </div> </blockquote> </div> <br> </div> </blockquote> <br> </body> </html>