Commit e3ddb2e0cc85ecb19592c59bb834ec7c0bb430c8

Authored by Jean-Philippe Bernard
1 parent 21924e89
Exists in master

improved

Showing 1 changed file with 13 additions and 20 deletions   Show diff stats
src/idl/dustem_plugin_stellar_population.pro
1   -FUNCTION dustem_plugin_stellar_population, key=key, val=val, scope=scope, paramtag=paramtag,help=help,test=test
  1 +FUNCTION dustem_plugin_stellar_population, key=key, val=val, scope=scope, paramtag=paramtag,help=help
2 2  
3 3 ;+
4 4 ; NAME:
... ... @@ -30,9 +30,19 @@ FUNCTION dustem_plugin_stellar_population, key=key, val=val, scope=scope, paramt
30 30 ; PROCEDURE:
31 31 ; This is a dustem plugin
32 32 ;-
33   -;- HISTORY: ADDED STELLAR LUMINOSITY CLASSES
34 33  
35   -IF keyword_set(test) THEN stop
  34 +;================NOTA BENE==========================================================================
  35 +;DATA IS NOW RETRIEVED FROM THIS TEXT FILE: "A Modern Mean Dwarf Stellar Color and Effective Temperature Sequence"
  36 +;NB: intermediate spectral classes '.5' are not taken into account - this will probably have to change
  37 +;NB: Also contact we need to contact the reasercher who wrote the text file because he said so in it.
  38 +;Other luminosity classes Should be included (only MS so far).
  39 +;BB approximation is a first degree 'bad' approximation because of the lack of radiative transfer especially at the
  40 +;photosphere of stars.
  41 +;REMARKS: BECAUSE WE STILL HAVEN'T SET THE DEFAULT VALUE FOR THE MAJORIY OF THE STELLAR POULATIONS,THE CORRESPONDING LINES ARE COMMENTED INSTEAD OF SETTING ARBITRARY VALUES.
  42 +;CONSIDERED STARS: SPEC_TYPE(N=5) = OBAFG, LUM_CLASS(N=10) = IA+,IA,IAB,IB,II,III,IV,V,VI,VII
  43 +;KM spectral types are not included because their UV part was not that important to excite the dust.
  44 +;This will help ease and shorten the fitting procedure
  45 +;If the user wants to use them without having to read the entirety of this plugin please contact the DustEmWrap team.
36 46  
37 47 IF keyword_set(help) THEN BEGIN
38 48 doc_library,'dustem_plugin_stellar_population'
... ... @@ -46,12 +56,9 @@ ENDIF
46 56  
47 57 IF keyword_set(scope) THEN BEGIN
48 58 message,'Scope keyword was set',/continue
49   - ;scope='STELLAR_POPULATION'
50   - ;scope='REPLACE_ISRF'
51 59 scope='ADD_ISRF'
52 60 stellar_isrf=0.
53 61 goto,the_end
54   - ;stop
55 62 ENDIF
56 63  
57 64 ;==== read stellar population parameters if not already in memory
... ... @@ -79,7 +86,6 @@ IF keyword_set(paramtag) THEN BEGIN
79 86 paramtag[ii]=cc+spectral_types[ij[0,0]]
80 87 ii=ii+1
81 88 ENDFOR
82   - ;stop
83 89 stellar_isrf=0.
84 90 goto, the_end
85 91 ENDIF
... ... @@ -92,19 +98,6 @@ rsun2cm = 6.957e10 ;Rsun in cm
92 98 c2a = 3e18 ;speed of light in ansgtroms/s (because of the Astron's PLANCK function)
93 99 pc2cm = 3.086e18 ;cm (cgs)
94 100  
95   -;================NOTA BENE==========================================================================
96   -;DATA IS NOW RETRIEVED FROM THIS TEXT FILE: "A Modern Mean Dwarf Stellar Color and Effective Temperature Sequence"
97   -;NB: intermediate spectral classes '.5' are not taken into account - this will probably have to change
98   -;NB: Also contact we need to contact the reasercher who wrote the text file because he said so in it.
99   -;Other luminosity classes Should be included (only MS so far).
100   -;BB approximation is a first degree 'bad' approximation because of the lack of radiative transfer especially at the
101   -;photosphere of stars.
102   -;REMARKS: BECAUSE WE STILL HAVEN'T SET THE DEFAULT VALUE FOR THE MAJORIY OF THE STELLAR POULATIONS,THE CORRESPONDING LINES ARE COMMENTED INSTEAD OF SETTING ARBITRARY VALUES.
103   -;CONSIDERED STARS: SPEC_TYPE(N=5) = OBAFG, LUM_CLASS(N=10) = IA+,IA,IAB,IB,II,III,IV,V,VI,VII
104   -;KM spectral types are not included because their UV part was not that important to excite the dust.
105   -;This will help ease and shorten the fitting procedure
106   -;If the user wants to use them without having to read the entirety of this plugin please contact the DustEmWrap team.
107   -
108 101 IF keyword_set(key) THEN BEGIN
109 102 ;stop
110 103 FOR i=0L,n_elements(key)-1 DO BEGIN
... ...