[mesa-users] colors module

astrophysics astrophysics at go2.pl
Fri Jun 21 02:20:34 EDT 2013


Ehsan,
It worked just great. Thanks again.
The only thing one needs to add are data type declarations for the code to be compiled correctly.
These are as follows

      REAL(8) :: denumerator
      REAL(8) :: Fe_H
      REAL(8) :: log_L
      REAL(8) :: log_Teff
      REAL(8) :: mass
      DOUBLE PRECISION :: Z_surf
and
      integer :: cnt

Cheers 
Marcin


astrophysics <astrophysics at go2.pl> pisze: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/20130621/d528a43d/attachment.html>


More information about the Mesa-users mailing list