[mesa-users] colors module

[Ehsan Moravveji]

Hello to Astrophysics!
Here is the way I do it in my run_star_extras:
1. before *implicit none*, add the following
      use colors_def, only : n_colors, init_colors_names
      use colors_lib

2. after *implicit none*, add the following
      real*8 :: colors_out(n_colors), log_g              ! for color output

3. in your *extras_finish_step*:
      log_Teff = dlog10(s% Teff)
      log_L = dlog10(s% photosphere_L)
      mass = s% star_mass
      Z_surf = 1d0 - s% surface_h1 - s% surface_he4
      denumerator = 1d0 - s% surface_he4 - (1d0+1.2857)*Z_surf
      Fe_H = dlog10(Z_surf / denumerator) - dlog10(0.0174d0)  ! See Eq. A.4
in Mowlavi et al. (2012)
      call colors_get(log_Teff, log_L, mass, Fe_H, colors_out, log_g, ierr)

4. in your *how_many_extra_history_columns*:
         how_many_extra_history_columns = n_colors

5. and finally in your *data_for_extra_history_columns*
         do cnt = 1, n_colors
            names(cnt) = colors_name(cnt)
            vals(cnt) = colors_out(cnt)
         enddo

I hope this helps.
Cheers.

Ehsan.


On Thu, Jun 20, 2013 at 8:06 AM, astrophysics <astrophysics at go2.pl> wrote:

> Just missed the title.
> ;-)
> BR
> Marcin
>
>
>
>
> -------- Oryginalna wiadomość --------
> Temat:
> Od:astrophysics <astrophysics at go2.pl>
> Do:mesa-users at lists.sourceforge.net
> DW:
>
>
> Hi,
>
> MESA has the 'colors' module. Among other values it calculates the color
> indices during the run. Unfortunately I lost track how to enable logging to
> an output file for those.
> Does anyone remember how to enable it and could remind me?
>
> Best regards
> Marcin
>
>
>
>
> ------------------------------------------------------------------------------
> This SF.net email is sponsored by Windows:
>
> Build for Windows Store.
>
> http://p.sf.net/sfu/windows-dev2dev
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users
>
>


-- 
Kind regards

Ehsan Moravveji
http://moravveji.webs.com/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20130620/e1edd1ab/attachment.html>