<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body bgcolor="#FFFFFF" text="#000000">
    <div class="moz-cite-prefix">Le 05/10/2013 14:18, Michel Rieutord a
      écrit :<br>
    </div>
    <blockquote cite="mid:525003A6.7040505@irap.omp.eu" type="cite">
      <meta content="text/html; charset=ISO-8859-1"
        http-equiv="Content-Type">
      <div class="moz-cite-prefix">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" class="moz-txt-link-freetext"
          href="http://code.google.com/p/ester-project/">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:1310.0908   or   arXiv:1310.1017<br>
        <br>
        Don't hesitate if you have questions.<br>
        <br>
        Best wishes,<br>
        <br>
        Michel Rieutord<br>
        <br>
      </div>
    </blockquote>
    PS third attempt; hope this is the good one.<br>
    <blockquote cite="mid:525003A6.7040505@irap.omp.eu" type="cite">
      <div class="moz-cite-prefix">
        <pre class="moz-signature" cols="72">**************************************************************************
Michel Rieutord
Institut de Recherche en Astrophysique et Planétologie
Observatoire Midi-Pyrénées      email : <a moz-do-not-send="true" 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 04/10/2013 22:25, Jeremy Jones a écrit :<br>
      </div>
      <blockquote
cite="mid:CAD9nXrZ-ZGW9Jy0DD2npW=Bmr7d_D8juMsFEHgEdRsiBMoyrbw@mail.gmail.com"
        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 class="" style="white-space:pre"> </span>History_Track1_win_flag

              = .true.</div>
            <div><span class="" style="white-space:pre"> </span>History_Track1_yname

              = 'log_R'</div>
            <div><span class="" style="white-space:pre"> </span>History_Track1_yaxis_label

              = 'log(R/R_sun)'</div>
            <div><span class="" style="white-space:pre"> </span>History_Track1_xname

              = 'log_Teff'</div>
            <div><span class="" style="white-space:pre"> </span>History_Track1_xaxis_label

              = 'log(T/K)'</div>
            <div><span class="" style="white-space:pre"> </span>History_Track1_title

              = 'Radius vs. Effective Temperature'</div>
            <div><br>
            </div>
            <div><span class="" style="white-space:pre"> </span>History_Track2_win_flag

              = .true.</div>
            <div><span class="" style="white-space:pre"> </span>History_Track2_yname

              = 'surf_avg_v_rot'</div>
            <div><span class="" style="white-space:pre"> </span>History_Track2_yaxis_label

              = 'Equatorial Velocity (km/s)'</div>
            <div><span class="" style="white-space:pre"> </span>History_Track2_xname

              = 'star_age'</div>
            <div><span class="" style="white-space:pre"> </span>History_Track2_xaxis_label

              = 'Age (yr)'</div>
            <div><span class="" style="white-space:pre"> </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">jones@chara.gsu.edu</a></div>
          </div>
        </div>
        <br>
        <fieldset class="mimeAttachmentHeader"></fieldset>
        <br>
        <pre wrap="">------------------------------------------------------------------------------
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" class="moz-txt-link-freetext" href="http://pubads.g.doubleclick.net/gampad/clk?id=60134791&iu=/4140/ostg.clktrk">http://pubads.g.doubleclick.net/gampad/clk?id=60134791&iu=/4140/ostg.clktrk</a></pre>
        <br>
        <fieldset class="mimeAttachmentHeader"></fieldset>
        <br>
        <pre wrap="">_______________________________________________
mesa-users mailing list
<a moz-do-not-send="true" class="moz-txt-link-abbreviated" href="mailto:mesa-users@lists.sourceforge.net">mesa-users@lists.sourceforge.net</a>
<a moz-do-not-send="true" class="moz-txt-link-freetext" href="https://lists.sourceforge.net/lists/listinfo/mesa-users">https://lists.sourceforge.net/lists/listinfo/mesa-users</a>
</pre>
      </blockquote>
      <br>
    </blockquote>
    <br>
  </body>
</html>