290 likes | 455 Views
|> cout << ”Last updated :” << last_update << “<br>”; |> Last updated : January 24, 2007 |>. January 2007 Bogdan Popescu. MASSCLEAN 1.00.
E N D
|> cout << ”Last updated :” << last_update << “\n”; |> Last updated : January 24, 2007 |> January 2007 Bogdan Popescu
MASSCLEAN 1.00 MASSive CLuster Evolution and ANalysis Package January 2007 Bogdan Popescu
Directory Structure install Makefile writeindex generator.cc index_trek.cc n_distribution.cc trek.cc generator.in start.trek COPYING generator n_distribution index_trek trek massclean (main directory) bin (copies of the compiled binaries) data indexiso mass_bins.in n_distribution.out trek.in *.eps gnuplot_n.gpl gnuplot_n.out iso.trek *.trek files docs (soon) examples (soon) tracksisoUBVRIJHKLM *. UBVRIJHKLM files (Geneva Database) iso.UBVRIJHKLM tracks.UBVRIJHKLM January 2007 Bogdan Popescu
Data Flow input output mass_bins.in standard input (or generator.in script file) generator n_distribution n_distribution.out gnuplot files mass_bins.in standard input (mass of the cluster) writeindex (script) files in the directory : /tracksisoUBVRIJHKLM/iso.UBVRIJHKLM/ indexiso index_trek indexiso trek.in trek trek.in *.trek files n_distribution.out *. UBVRIJHKLM files(Geneva database) ??? *.trek files *.FITS images January 2007 Bogdan Popescu
Data Flow input output mass_bins.in standard input (or generator.in script file) generator n_distribution n_distribution.out gnuplot files mass_bins.in standard input (mass of the cluster) writeindex (script) files in the directory : /tracksisoUBVRIJHKLM/iso.UBVRIJHKLM/ indexiso index_trek indexiso trek.in trek trek.in *.trek files n_distribution.out *. UBVRIJHKLM files(Geneva database) ??? *.trek files *.FITS images January 2007 Bogdan Popescu
generator.cc input output generator.in script : ./generator << ! 0.8 [Mmin] 130 [Mmax] 0.1 [bin1] 0.05 [bin2] 200 [max no of bins] ! mass_bins.in file : ## Mass Bin (DeltaM/2) ## (in solar units) ## 0.8 0.1 [mass interval 0.7-0.9] 1.05 0.15 [mass interval 0.9-1.2] 1.4 0.2 [mass interval 1.2-1.4] 1.85 0.25 2.4 0.3 3.05 0.35 3.8 0.4 4.65 0.45 5.6 0.5 ........... ........... ........... 115.8 2.4 [mass interval 113.4-118.2] 120.65 2.45 [mass interval 118.2-123.1] or standard input January 2007 Bogdan Popescu
generator.cc output mass_bins.in file : ## Mass Bin (DeltaM/2) ## (in solar units) ## 0.8 0.1 1.05 0.15 1.4 0.2 1.85 0.25 2.4 0.3 3.05 0.35 3.8 0.4 4.65 0.45 5.6 0.5 ........... ........... ........... 115.8 2.4 120.65 2.45 Mass interval [M - bin1 - kx bin2, M + bin1 + kx bin2] January 2007 Bogdan Popescu
generator.cc - screenshot January 2007 Bogdan Popescu
Data Flow input output mass_bins.in standard input (or generator.in script file) generator n_distribution n_distribution.out gnuplot files mass_bins.in standard input (mass of the cluster) writeindex (script) files in the directory : /tracksisoUBVRIJHKLM/iso.UBVRIJHKLM/ indexiso index_trek indexiso trek.in trek trek.in *.trek files n_distribution.out *. UBVRIJHKLM files(Geneva database) ??? *.trek files *.FITS images January 2007 Bogdan Popescu
n_distribution.cc input output mass_bins.in file standard input (mass of the cluster) n_distribution.out file : # The number of bins : 48 # The total mass : 100000 # The mass in the distribution is : 99820.7 ### Mass Number # 0.8 12629 1.05 10063 1.4 6824 1.85 4418 ........... ........... ........... 111.05 3 115.8 2 120.65 2 ## Mass Bin (DeltaM/2) ## (in solar units) ## 0.8 0.1 1.05 0.15 1.4 0.2 1.85 0.25 2.4 0.3 3.05 0.35 3.8 0.4 4.65 0.45 5.6 0.5 ........... ........... ........... 115.8 2.4 120.65 2.45 January 2007 Bogdan Popescu
n_distribution.cc output January 2007 Bogdan Popescu
n_distribution.cc output January 2007 Bogdan Popescu
n_distribution.cc – screenshot January 2007 Bogdan Popescu
Data Flow input output mass_bins.in standard input (or generator.in script file) generator n_distribution n_distribution.out gnuplot files mass_bins.in standard input (mass of the cluster) writeindex (script) files in the directory : /tracksisoUBVRIJHKLM/iso.UBVRIJHKLM/ indexiso index_trek indexiso trek.in trek trek.in *.trek files n_distribution.out *. UBVRIJHKLM files(Geneva database) ??? *.trek files *.FITS images January 2007 Bogdan Popescu
writeindex (script) input output files in the directory : /tracksisoUBVRIJHKLM/iso.UBVRIJHKLM/ (Geneva Database) indexiso file : iso_c020_0300.UBVRIJHKLM iso_c020_0500.UBVRIJHKLM iso_c020_0530.UBVRIJHKLM iso_c020_0559.UBVRIJHKLM iso_c020_0580.UBVRIJHKLM iso_c020_0590.UBVRIJHKLM iso_c020_0600.UBVRIJHKLM ........................ iso_e020_0744.UBVRIJHKLM iso_e020_0750.UBVRIJHKLM writeindex script : ls tracksisoUBVRIJHKLM/iso.UBVRIJHKLM/ > ././data/indexiso January 2007 Bogdan Popescu
Data Flow input output mass_bins.in standard input (or generator.in script file) generator n_distribution n_distribution.out gnuplot files mass_bins.in standard input (mass of the cluster) writeindex (script) files in the directory : /tracksisoUBVRIJHKLM/iso.UBVRIJHKLM/ indexiso index_trek indexiso trek.in trek trek.in *.trek files n_distribution.out *. UBVRIJHKLM files(Geneva database) ??? *.trek files *.FITS images January 2007 Bogdan Popescu
index_trek.cc input output indexiso file trek.in file : # # B.Popescu and M.M.Hanson,2007 # trek.in # # Note : EXACTLY 7 extra lines required # ## Isochrone file Output file iso_c020_0300.UBVRIJHKLM iso_c020_0300.trek iso_c020_0500.UBVRIJHKLM iso_c020_0500.trek iso_c020_0530.UBVRIJHKLM iso_c020_0530.trek ............................................... ............................................... ............................................... iso_e020_0740.UBVRIJHKLM iso_e020_0740.trek iso_e020_0744.UBVRIJHKLM iso_e020_0744.trek iso_e020_0750.UBVRIJHKLM iso_e020_0750.trek January 2007 Bogdan Popescu
index_trek.cc – screenshot January 2007 Bogdan Popescu
Data Flow input output mass_bins.in standard input (or generator.in script file) generator n_distribution n_distribution.out gnuplot files mass_bins.in standard input (mass of the cluster) writeindex (script) files in the directory : /tracksisoUBVRIJHKLM/iso.UBVRIJHKLM/ indexiso index_trek indexiso trek.in trek trek.in *.trek files n_distribution.out *. UBVRIJHKLM files(Geneva database) ??? *.trek files *.FITS images January 2007 Bogdan Popescu
trek.cc input output trek.in n_distribution.out *. UBVRIJHKLM files(Geneva database) *.trek files (the variation of the mass in a few Geneva *. UBVRIJHKLM isochrone files) January 2007 Bogdan Popescu
trek.cc – screenshot January 2007 Bogdan Popescu
trek.cc – screenshot January 2007 Bogdan Popescu
trek.cc output *.trek file example : ## ## B.Popescu and M.M.Hanson,2007 ## ## # iso_c020_0885.trek # # N M M_init M_act logTeff logg logL M( V) U-B B-V V-R V-I V-K R-I I-K J-H H-K K-L J-K J-L J-L2 K-M # # 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 12629 0.8 0.8 0.8 3.687 4.651 -0.606 6.613 0.847 0.96 0.526 0.949 2.274 0.423 1.319 0.513 0.085 0.076 0.598 0.674 0.676 0.552 10063 1.05 1.05 1.05 3.766 4.506 -0.029 4.932 0.125 0.604 0.354 0.674 1.494 0.32 0.815 0.331 0.051 0.04 0.383 0.422 0.421 0.37 6824 1.4 1.4 1.4 3.835 4.291 0.588 3.28 0 0.38 0.215 0.419 0.841 0.204 0.417 0.162 0.034 0.007 0.196 0.203 0.204 0.197 4418 1.85 1.85 1.85 3.902 4.102 1.165 1.821 0.119 0.173 0.09 0.202 0.363 0.112 0.156 0.074 0.029 -0.02 0.103 0.083 0.083 0.081 2866 2.4 2.4034 2.3982 3.688 2.749 1.776 0.641 0.753 0.94 0.525 0.962 2.224 0.437 1.256 0.516 0.085 0.068 0.601 0.669 0.671 0.506 # iso_c020_0885.trek c = basic grid log(t) range : [3.00,10.20] 0.020 = metalicity Z 0885 - log(t)=8.85 # iso_e020_0750.trek e = extended grid (high mass loss for massive stars) log(t) range : [3.00,7.50] 0.020 = metalicity Z 0750 - log(t)=7.50 January 2007 Bogdan Popescu
trek.cc output *.trek file example : ## ## B.Popescu and M.M.Hanson,2007 ## ## # iso_c020_0600.trek # # N M M_init M_act logTeff logg logL M( V) U-B B-V # # 1 2 3 4 5 6 7 8 9 10 12629 0.8 0.8 0.8 3.687 4.656 -0.612 6.628 0.847 0.96 10063 1.05 1.05 1.05 3.763 4.523 -0.058 5.009 0.139 0.615 6824 1.4 1.4 1.4 3.834 4.336 0.54 3.402 -0.003 0.385 4418 1.85 1.85 1.85 3.931 4.314 1.068 2.13 0.083 0.103 2866 2.4 2.4 2.4 4.018 4.315 1.528 1.31 -0.121 -0.042 1898 3.05 3.05 3.05 4.093 4.312 1.936 0.696 -0.362 -0.114 1290 3.8 3.82 3.82 4.159 4.304 2.309 0.165 -0.538 -0.157 .................................................................................................................................. 3 97.4 97.02 93.3075 4.677 3.918 6.152 -6.754 -1.203 -0.322 3 101.85 102.02 97.5413 4.678 3.91 6.185 -6.838 -1.205 -0.321 3 106.4 106.02 100.8958 4.679 3.903 6.211 -6.903 -1.207 -0.32 3 111.05 111.02 105.0492 4.681 3.895 6.241 -6.976 -1.211 -0.32 2 115.8 116.02 109.1592 4.682 3.888 6.27 -7.049 -1.213 -0.319 January 2007 Bogdan Popescu
Data Flow input output mass_bins.in standard input (or generator.in script file) generator n_distribution n_distribution.out gnuplot files mass_bins.in standard input (mass of the cluster) writeindex (script) files in the directory : /tracksisoUBVRIJHKLM/iso.UBVRIJHKLM/ indexiso index_trek indexiso trek.in trek trek.in *.trek files n_distribution.out *. UBVRIJHKLM files(Geneva database) ??? *.trek files *.FITS images January 2007 Bogdan Popescu
Near Future ???.cc Combines c and e isochrones, to include high mass loss for massive stars. Not necessary. ???.cc Interface to FITS images. Maybe CFITSIO, maybe input for a photometry program. ???.cc Random coordinate (position) generator – possibly, depending on previous step. ???.cc Red giants, extinction, galactic plane... January 2007 Bogdan Popescu
Miscellaneous All programs have been intense tested, using : * 6-506 bins, * 1-128 isochrone files, * 103-107 for the mass of the cluster (in solar units). gcc 3.3.2 (Fedora Core 1) and gcc 3.4.4 (Fedora Core 4) are recommended. gcc 4.1.1(Fedora Core 6)IS NOT RECOMMENDED. (precompiled binaries are available in the “/bin” directory.) January 2007 Bogdan Popescu
References 1. Geneva Database http://webast.ast.obs-mip.fr/equipe/stellar/ http://obswww.unige.ch/~mowlavi/evol/stev_database.html 2. Schaller, G., Schaerer, D., Meynet, G., Maeder, A., 1992, A&AS, 96, 269 3. Schaerer, D., Meynet, G., Maeder, A., Schaller, G., 1993, A&AS, 98, 523 4. Charbonnel, C., Meynet, G., Maeder, A., Schaller, G., Schaerer, D., 1993, A&AS, 101, 415 5. Schaerer, D., Charbonnel, C., Meynet, G., Maeder, A., Schaller, G., 1993, A&AS, 102, 339 6. Meynet, G., Maeder, A., Schaller, G., Schaerer, D., Charbonnel, C.,1994, A&AS, 103, 97 7. Charbonnel, C., Meynet, G., Maeder, A., Schaerer, D., 1996, A&AS, 115, 339 8. Mowlavi, N., Schaerer, D., Meynet, G., Bernasconi, P. A., Charbonnel, C., Maeder, A., 1998a, A&AS, 128, 471 9. Charbonnel, C., Dappen, W., Schaerer, D., Bernasconi, P.A., Maeder, A., Meynet, G., Mowlavi, N., 1999, A&AS, 135, 405 10. Lejeune, T. and Schaerer, D., A&A 366, 538-546 (2001) 11. Salpeter, E.E. 1955, ApJ 121,161S January 2007 Bogdan Popescu
January 2007 Bogdan Popescu