[mesa-users] colors module

[astrophysics]

Thanks Ehsan.

That was quick. 
I'll try it and let you know how it worked.

Cheers,
Marcin

Ehsan Moravveji <e.moravveji at gmail.com> pisze: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/3fba3c79/attachment.html>