Function dustem_plugin_stellar_population, key=key, val=val, scope=scope, paramtag=paramtag,help=help ;+ ; NAME: ; dustem_plugin_stellar_population ; PURPOSE: ; replaces the default DUSTEM ISRF with a composite stellar spectrum ; CATEGORY: ; DUSTEM Wrapper ; CALLING SEQUENCE: ; dustem_plugin_stellar_population(key=key,val=val) ; INPUTS: ; None ; OPTIONAL INPUT PARAMETERS: ; key = input parameter number ; val = input parameter value ; OUTPUTS: ; None ; OPTIONAL OUTPUT PARAMETERS: ; None ; ACCEPTED KEY-WORDS: ; help = if set, print this help ; COMMON BLOCKS: ; None ; SIDE EFFECTS: ; None ; RESTRICTIONS: ; The dustem fortran code must be installed ; The dustem idl wrapper must be installed ; PROCEDURE: ; This is a dustem plugin ;- ;- HISTORY: ADDED STELLAR LUMINOSITY CLASSES IF keyword_set(help) THEN BEGIN doc_library,'dustem_plugin_stellar_population' goto,the_end ENDIF IF keyword_set(scope) THEN BEGIN out=0 goto, the_scope ENDIF IF keyword_set(paramtag) THEN BEGIN out=0 goto, the_paramtag ENDIF out=0;Do we need to initialize the output? ;===========Initializing the defaut stellar population structure=========== one_pop={radius:0.,temperature:0.,distance:0.,nstars:0.,popid:''} ;===========Counters=========== popnumber = 0. counto = 0. countb = 0. counta = 0. countf = 0. countg = 0. ;============================== ;===========Constants (in cgs)========== (except for the stellar population distance (in pc)) rsun2cm = 6.957e10 ;================NOTA BENE========================================================================== ;DATA IS NOW RETRIEVED FROM THIS TEXT FILE: "A Modern Mean Dwarf Stellar Color and Effective Temperature Sequence" ;NB: intermediate spectral classes '.5' are not taken into account - this will probably have to change ;NB: Also contact we need to contact the reasercher who wrote the text file because he said so in it. ;Other luminosity classes Should be included (only MS so far). ;BB approximation is a first degree 'bad' approximation because of the lack of radiative transfer especially at the ;photosphere of stars. ;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. ;CONSIDERED STARS: SPEC_TYPE(N=5) = OBAFG, LUM_CLASS(N=10) = IA+,IA,IAB,IB,II,III,IV,V,VI,VII ;KM spectral types are not included because their UV part was not that important to excite the dust. ;This will help ease and shorten the fitting procedure ;If the user wants to use them without having to read the entirety of this plugin please contact the DustEmWrap team. IF keyword_set(key) THEN BEGIN file = !dustem_wrap_soft_dir+'src/idl/EEM_dwarf_UBVIJHK_colors_Teff.txt' frmt = '(A,A,A,A,A,A,A,A,A,A,A,A,A,A,A,A,A,A,A,A,A,A,A,A,A,A,A,A,A,A,A,A,A)' readcol, file, V1,V2,V3,V4,V5,V6,V7,V8,V9,V10,V11,V12,V13,V14,V15,V16,V17,V18,V20,V21,V22,V23,V24,V25,V26,V27,V28,V29,V30,V31,V32, format=frmt, delim=' ', comment='#',/nan ;Paramters that are needed for the BB approximation (horrible approximation even for MS stars) SpecType_MS = V1 T_eff_MS = V2 R_star_MS = V7 lumclass = (strmid(key,2)).remove(-1) ; extracting the luminosity class strings params = strmid(key,0,/reverse_offset) ; extracting the indices (in the 'key' array) of the stellar parameters to fit ; Locating the O3 stellar population parameters - this new code structure accounts for the presence of different luminosity classes. testo3 = strupcase(strmid(key,0,2)) EQ 'O3' ;NB: test is like this so that user can use lowercase for spectral class ind = where(testo3,countindo3) ; retrieveing the defaut values from the txt database tsto3 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'O3' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tsto3,ctdato3) if ctdato3 ne 0 then begin T_eff_MS_o3 = double(T_eff_MS(indat)) R_star_MS_o3 = double(R_star_MS(indat)) endif IF countindo3 NE 0 THEN BEGIN lumclasso3 = lumclass(ind) nlmo3 = n_elements(lumclasso3) indo3_1 = fltarr(nlmo3) indo3_2 = indo3_1 indo3_3 = indo3_1 indo3_4 = indo3_1 countindo3_1 = indo3_1 countindo3_2 = indo3_1 countindo3_3 = indo3_1 countindo3_4 = indo3_1 FOR i=0L,nlmo3-1 DO BEGIN popnumber+=1 counto+=1 indo3_1[i] = where(testo3 and params EQ 1, countindo3_1x); index of radius of this stellar population in the val array countindo3_1[i] = countindo3_1x indo3_2[i] = where(testo3 and params EQ 2, countindo3_2x); index of temperature of this stellar population in the val array countindo3_2[i] = countindo3_2x indo3_3[i] = where(testo3 and params EQ 3, countindo3_3x); index of distance of this stellar population in the val array countindo3_3[i] = countindo3_3x indo3_4[i] = where(testo3 and params EQ 4, countindo3_4x); index of the nstars of this stellar population in the val array (Q: does it make sense to have an nstars>1.? If not upper limiting it in the general fitting procedure might be handy) countindo3_4[i] = countindo3_4x ENDFOR ENDIF ; Locating the O4 stellar population parameters testo4 = strupcase(strmid(key,0,2)) EQ 'O4' ind = where(testo4,countindo4) ; retrieveing the defaut values from the txt database tsto4 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'O4' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tsto4,ctdato4) if ctdato4 ne 0 then begin T_eff_MS_o4 = double(T_eff_MS(indat)) R_star_MS_o4 = double(R_star_MS(indat)) endif IF countindo4 NE 0 THEN BEGIN lumclasso4 = lumclass(ind) nlmo4 = n_elements(lumclasso4) indo4_1 = fltarr(nlmo4) indo4_2 = indo4_1 indo4_3 = indo4_1 indo4_4 = indo4_1 countindo4_1 = indo4_1 countindo4_2 = indo4_1 countindo4_3 = indo4_1 countindo4_4 = indo4_1 FOR i=0L,nlmo4-1 DO BEGIN popnumber+=1 counto+=1 indo4_1[i] = where(testo4 and params EQ 1, countindo4_1x); index of radius of this stellar population in the val array countindo4_1[i] = countindo4_1x indo4_2[i] = where(testo4 and params EQ 2, countindo4_2x); index of temperature of this stellar population in the val array countindo4_2[i] = countindo4_2x indo4_3[i] = where(testo4 and params EQ 3, countindo4_3x); index of distance of this stellar population in the val array countindo4_3[i] = countindo4_3x indo4_4[i] = where(testo4 and params EQ 4, countindo4_4x); index of the nstars of this stellar population in the val array countindo4_4[i] = countindo4_4x ENDFOR ENDIF ; Locating the O5 stellar population parameters testo5 = strupcase(strmid(key,0,2)) EQ 'O5' ind = where(testo5,countindo5) ; retrieveing the defaut values from the txt database tsto5 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'O5' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tsto5,ctdato5) if ctdato5 ne 0 then begin T_eff_MS_o5 = double(T_eff_MS(indat)) R_star_MS_o5 = double(R_star_MS(indat)) endif IF countindo5 NE 0 THEN BEGIN lumclasso5 = lumclass(ind) nlmo5 = n_elements(lumclasso5) indo5_1 = fltarr(nlmo5) indo5_2 = indo5_1 indo5_3 = indo5_1 indo5_4 = indo5_1 countindo5_1 = indo5_1 countindo5_2 = indo5_1 countindo5_3 = indo5_1 countindo5_4 = indo5_1 FOR i=0L,nlmo5-1 DO BEGIN popnumber+=1 counto+=1 indo5_1[i] = where(testo5 and params EQ 1, countindo5_1x); index of radius of this stellar population in the val array countindo5_1[i] = countindo5_1x indo5_2[i] = where(testo5 and params EQ 2, countindo5_2x); index of temperature of this stellar population in the val array countindo5_2[i] = countindo5_2x indo5_3[i] = where(testo5 and params EQ 3, countindo5_3x); index of distance of this stellar population in the val array countindo5_3[i] = countindo5_3x indo5_4[i] = where(testo5 and params EQ 4, countindo5_4x); index of the nstars of this stellar population in the val array countindo5_4[i] = countindo5_4x ENDFOR ENDIF ; Locating the O6 stellar population parameters testo6 = strupcase(strmid(key,0,2)) EQ 'O6' ind = where(testo6,countindo6) ; retrieveing the defaut values from the txt database tsto6 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'O6' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tsto6,ctdato6) if ctdato6 ne 0 then begin T_eff_MS_o6 = double(T_eff_MS(indat)) R_star_MS_o6 = double(R_star_MS(indat)) endif IF countindo6 NE 0 THEN BEGIN lumclasso6 = lumclass(ind) nlmo6 = long(n_elements(lumclasso6)) indo6_1 = fltarr(nlmo6) indo6_2 = indo6_1 indo6_3 = indo6_1 indo6_4 = indo6_1 countindo6_1 = indo6_1 countindo6_2 = indo6_1 countindo6_3 = indo6_1 countindo6_4 = indo6_1 FOR i=0L,nlmo6-1 DO BEGIN popnumber+=1 counto+=1 indo6_1[i] = where(testo6 and params EQ 1, countindo6_1x); index of radius of this stellar population in the val array countindo6_1[i] = countindo6_1x indo6_2[i] = where(testo6 and params EQ 2, countindo6_2x); index of temperature of this stellar population in the val array countindo6_2[i] = countindo6_2x indo6_3[i] = where(testo6 and params EQ 3, countindo6_3x); index of distance of this stellar population in the val array countindo6_3[i] = countindo6_3x indo6_4[i] = where(testo6 and params EQ 4, countindo6_4x); index of the nstars of this stellar population in the val array countindo6_4[i] = countindo6_4x ENDFOR ENDIF ; Locating the O7 stellar population parameters testo7 = strupcase(strmid(key,0,2)) EQ 'O7' ind = where(testo7,countindo7) ; retrieveing the defaut values from the txt database tsto7 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'O7' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tsto7,ctdato7) if ctdato7 ne 0 then begin T_eff_MS_o7 = double(T_eff_MS(indat)) R_star_MS_o7 = double(R_star_MS(indat)) endif IF countindo7 NE 0 THEN BEGIN lumclasso7 = lumclass(ind) nlmo7 = n_elements(lumclasso7) indo7_1 = fltarr(nlmo7) indo7_2 = indo7_1 indo7_3 = indo7_1 indo7_4 = indo7_1 countindo7_1 = indo7_1 countindo7_2 = indo7_1 countindo7_3 = indo7_1 countindo7_4 = indo7_1 FOR i=0L,nlmo7-1 DO BEGIN popnumber+=1 counto+=1 indo7_1[i] = where(testo7 and params EQ 1, countindo7_1x); index of radius of this stellar population in the val array countindo7_1[i] = countindo7_1x indo7_2[i] = where(testo7 and params EQ 2, countindo7_2x); index of temperature of this stellar population in the val array countindo7_2[i] = countindo7_2x indo7_3[i] = where(testo7 and params EQ 3, countindo7_3x); index of distance of this stellar population in the val array countindo7_3[i] = countindo7_3x indo7_4[i] = where(testo7 and params EQ 4, countindo7_4x); index of the nstars of this stellar population in the val array countindo7_4[i] = countindo7_4x ENDFOR ENDIF ; Locating the O8 stellar population parameters testo8 = strupcase(strmid(key,0,2)) EQ 'O8' ind = where(testo8,countindo8) ; retrieveing the defaut values from the txt database tsto8 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'O8' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tsto8,ctdato8) if ctdato8 ne 0 then begin T_eff_MS_o8 = double(T_eff_MS(indat)) R_star_MS_o8 = double(R_star_MS(indat)) endif IF countindo8 NE 0 THEN BEGIN lumclasso8 = lumclass(ind) nlmo8 = n_elements(lumclasso8) indo8_1 = fltarr(nlmo8) indo8_2 = indo8_1 indo8_3 = indo8_1 indo8_4 = indo8_1 countindo8_1 = indo8_1 countindo8_2 = indo8_1 countindo8_3 = indo8_1 countindo8_4 = indo8_1 FOR i=0L,nlmo8-1 DO BEGIN popnumber+=1 counto+=1 indo8_1[i] = where(testo8 and params EQ 1, countindo8_1x); index of radius of this stellar population in the val array countindo8_1[i] = countindo8_1x indo8_2[i] = where(testo8 and params EQ 2, countindo8_2x); index of temperature of this stellar population in the val array countindo8_2[i] = countindo8_2x indo8_3[i] = where(testo8 and params EQ 3, countindo8_3x); index of distance of this stellar population in the val array countindo8_3[i] = countindo8_3x indo8_4[i] = where(testo8 and params EQ 4, countindo8_4x); index of the nstars of this stellar population in the val array countindo8_4[i] = countindo8_4x ENDFOR ENDIF ; Locating the O9 stellar population parameters testo9 = strupcase(strmid(key,0,2)) EQ 'O9' ind = where(testo9,countindo9) ; retrieveing the defaut values from the txt database tsto9 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'O9' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tsto9,ctdato9) if ctdato9 ne 0 then begin T_eff_MS_o9 = double(T_eff_MS(indat)) R_star_MS_o9 = double(R_star_MS(indat)) endif IF countindo9 NE 0 THEN BEGIN lumclasso9 = strupcase(lumclass(ind)) nlmo9 = long(n_elements(lumclasso9)) indo9_1 = fltarr(nlmo9) indo9_2 = indo9_1 indo9_3 = indo9_1 indo9_4 = indo9_1 countindo9_1 = indo9_1 countindo9_2 = indo9_1 countindo9_3 = indo9_1 countindo9_4 = indo9_1 FOR i=0L,nlmo9-1 DO BEGIN popnumber+=1 counto+=1 indo9_1[i] = where(testo9 and params EQ 1, countindo9_1x); index of radius of this stellar population in the val array countindo9_1[i] = countindo9_1x indo9_2[i] = where(testo9 and params EQ 2, countindo9_2x); index of temperature of this stellar population in the val array countindo9_2[i] = countindo9_2x indo9_3[i] = where(testo9 and params EQ 3, countindo9_3x); index of distance of this stellar population in the val array countindo9_3[i] = countindo9_3x indo9_4[i] = where(testo9 and params EQ 4, countindo9_4x); index of the nstars of this stellar population in the val array countindo9_4[i] = countindo9_4x ENDFOR ENDIF ;====================================================== ; Locating the B0 stellar population parameters testb0 = strupcase(strmid(key,0,2)) EQ 'B0' ind = where(testb0,countindb0) ; retrieveing the defaut values from the txt database tstb0 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'B0' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstb0,ctdatb0) if ctdatb0 ne 0 then begin T_eff_MS_b0 = double(T_eff_MS(indat)) R_star_MS_b0 = double(R_star_MS(indat)) endif IF countindb0 NE 0 THEN BEGIN lumclassb0 = uniq(lumclass(ind)) nlmb0 = n_elements(lumclassb0) indb0_1 = fltarr(nlmb0) indb0_2 = indb0_1 indb0_3 = indb0_1 indb0_4 = indb0_1 countindb0_1 = indb0_1 countindb0_2 = indb0_1 countindb0_3 = indb0_1 countindb0_4 = indb0_1 FOR i=0L,nlmb0-1 DO BEGIN popnumber+=1 countb+=1 indb0_1[i] = where(testb0 and params EQ 1, countindb0_1x); index of radius of this stellar population in the val array countindb0_1[i] = countindb0_1x indb0_2[i] = where(testb0 and params EQ 2, countindb0_2x); index of temperature of this stellar population in the val array countindb0_2[i] = countindb0_2x indb0_3[i] = where(testb0 and params EQ 3, countindb0_3x); index of distance of this stellar population in the val array countindb0_3[i] = countindb0_3x indb0_4[i] = where(testb0 and params EQ 4, countindb0_4x); index of the nstars of this stellar population in the val array countindb0_4[i] = countindb0_4x ENDFOR ENDIF ; Locating the B1 stellar population parameters testb1 = strupcase(strmid(key,0,2)) EQ 'B1' ind = where(testb1,countindb1) ; retrieveing the defaut values from the txt database tstb1 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'B1' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstb1,ctdatb1) if ctdatb1 ne 0 then begin T_eff_MS_b1 = double(T_eff_MS(indat)) R_star_MS_b1 = double(R_star_MS(indat)) endif IF countindb1 NE 0 THEN BEGIN lumclassb1 = strupcase(lumclass(ind)) nlmb1 = long(n_elements(lumclassb1)) indb1_1 = fltarr(nlmb1) indb1_2 = indb1_1 indb1_3 = indb1_1 indb1_4 = indb1_1 countindb1_1 = indb1_1 countindb1_2 = indb1_1 countindb1_3 = indb1_1 countindb1_4 = indb1_1 FOR i=0L,nlmb1-1 DO BEGIN popnumber+=1 countb+=1 indb1_1[i] = where(testb1 and params EQ 1, countindb1_1x); index of radius of this stellar population in the val array countindb1_1[i] = countindb1_1x indb1_2[i] = where(testb1 and params EQ 2, countindb1_2x); index of temperature of this stellar population in the val array countindb1_2[i] = countindb1_2x indb1_3[i] = where(testb1 and params EQ 3, countindb1_3x); index of distance of this stellar population in the val array countindb1_3[i] = countindb1_3x indb1_4[i] = where(testb1 and params EQ 4, countindb1_4x); index of the nstars of this stellar population in the val array countindb1_4[i] = countindb1_4x ENDFOR ENDIF ; Locating the B2 stellar population parameters testb2 = strupcase(strmid(key,0,2)) EQ 'B2' ind = where(testb2,countindb2) ; retrieveing the defaut values from the txt database tstb2 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'B2' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstb2,ctdatb2) if ctdatb2 ne 0 then begin T_eff_MS_b2 = double(T_eff_MS(indat)) R_star_MS_b2 = double(R_star_MS(indat)) endif IF countindb2 NE 0 THEN BEGIN lumclassb2 = uniq(lumclass(ind)) nlmb2 = n_elements(lumclassb2) indb2_1 = fltarr(nlmb2) indb2_2 = indb2_1 indb2_3 = indb2_1 indb2_4 = indb2_1 countindb2_1 = indb2_1 countindb2_2 = indb2_1 countindb2_3 = indb2_1 countindb2_4 = indb2_1 FOR i=0L,nlmb2-1 DO BEGIN popnumber+=1 countb+=1 indb2_1[i] = where(testb2 and params EQ 1, countindb2_1x); index of radius of this stellar population in the val array countindb2_1[i] = countindb2_1x indb2_2[i] = where(testb2 and params EQ 2, countindb2_2x); index of temperature of this stellar population in the val array countindb2_2[i] = countindb2_2x indb2_3[i] = where(testb2 and params EQ 3, countindb2_3x); index of distance of this stellar population in the val array countindb2_3[i] = countindb2_3x indb2_4[i] = where(testb2 and params EQ 4, countindb2_4x); index of the nstars of this stellar population in the val array countindb2_4[i] = countindb2_4x ENDFOR ENDIF ; Locating the B3 stellar population parameters testb3 = strupcase(strmid(key,0,2)) EQ 'B3' ind = where(testb3,countindb3) ; retrieveing the defaut values from the txt database tstb3 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'B3' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstb3,ctdatb3) if ctdatb3 ne 0 then begin T_eff_MS_b3 = double(T_eff_MS(indat)) R_star_MS_b3 = double(R_star_MS(indat)) endif IF countindb3 NE 0 THEN BEGIN lumclassb3 = uniq(lumclass(ind)) nlmb3 = n_elements(lumclassb3) indb3_1 = fltarr(nlmb3) indb3_2 = indb3_1 indb3_3 = indb3_1 indb3_4 = indb3_1 countindb3_1 = indb3_1 countindb3_2 = indb3_1 countindb3_3 = indb3_1 countindb3_4 = indb3_1 FOR i=0L,nlmb3-1 DO BEGIN popnumber+=1 countb+=1 indb3_1[i] = where(testb3 and params EQ 1, countindb3_1x); index of radius of this stellar population in the val array countindb3_1[i] = countindb3_1x indb3_2[i] = where(testb3 and params EQ 2, countindb3_2x); index of temperature of this stellar population in the val array countindb3_2[i] = countindb3_2x indb3_3[i] = where(testb3 and params EQ 3, countindb3_3x); index of distance of this stellar population in the val array countindb3_3[i] = countindb3_3x indb3_4[i] = where(testb3 and params EQ 4, countindb3_4x); index of the nstars of this stellar population in the val array countindb3_4[i] = countindb3_4x ENDFOR ENDIF ; Locating the B4 stellar population parameters testb4 = strupcase(strmid(key,0,2)) EQ 'B4' ind = where(testb4,countindb4) ; retrieveing the defaut values from the txt database tstb4 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'B4' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstb4,ctdatb4) if ctdatb4 ne 0 then begin T_eff_MS_b4 = double(T_eff_MS(indat)) R_star_MS_b4 = double(R_star_MS(indat)) endif IF countindb4 NE 0 THEN BEGIN lumclassb4 = uniq(lumclass(ind)) nlmb4 = n_elements(lumclassb4) indb4_1 = fltarr(nlmb4) indb4_2 = indb4_1 indb4_3 = indb4_1 indb4_4 = indb4_1 countindb4_1 = indb4_1 countindb4_2 = indb4_1 countindb4_3 = indb4_1 countindb4_4 = indb4_1 FOR i=0L,nlmb4-1 DO BEGIN popnumber+=1 countb+=1 indb4_1[i] = where(testb4 and params EQ 1, countindb4_1x); index of radius of this stellar population in the val array countindb4_1[i] = countindb4_1x indb4_2[i] = where(testb4 and params EQ 2, countindb4_2x); index of temperature of this stellar population in the val array countindb4_2[i] = countindb4_2x indb4_3[i] = where(testb4 and params EQ 3, countindb4_3x); index of distance of this stellar population in the val array countindb4_3[i] = countindb4_3x indb4_4[i] = where(testb4 and params EQ 4, countindb4_4x); index of the nstars of this stellar population in the val array countindb4_4[i] = countindb4_4x ENDFOR ENDIF ; Locating the B5 stellar population parameters testb5 = strupcase(strmid(key,0,2)) EQ 'B5' ind = where(testb5,countindb5) ; retrieveing the defaut values from the txt database tstb5 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'B5' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstb5,ctdatb5) if ctdatb5 ne 0 then begin T_eff_MS_b5 = double(T_eff_MS(indat)) R_star_MS_b5 = double(R_star_MS(indat)) endif IF countindb5 NE 0 THEN BEGIN lumclassb5 = uniq(lumclass(ind)) nlmb5 = n_elements(lumclassb5) indb5_1 = fltarr(nlmb5) indb5_2 = indb5_1 indb5_3 = indb5_1 indb5_4 = indb5_1 countindb5_1 = indb5_1 countindb5_2 = indb5_1 countindb5_3 = indb5_1 countindb5_4 = indb5_1 FOR i=0L,nlmb5-1 DO BEGIN popnumber+=1 countb+=1 indb5_1[i] = where(testb5 and params EQ 1, countindb5_1x); index of radius of this stellar population in the val array countindb5_1[i] = countindb5_1x indb5_2[i] = where(testb5 and params EQ 2, countindb5_2x); index of temperature of this stellar population in the val array countindb5_2[i] = countindb5_2x indb5_3[i] = where(testb5 and params EQ 3, countindb5_3x); index of distance of this stellar population in the val array countindb5_3[i] = countindb5_3x indb5_4[i] = where(testb5 and params EQ 4, countindb5_4x); index of the nstars of this stellar population in the val array countindb5_4[i] = countindb5_4x ENDFOR ENDIF ; Locating the B6 stellar population parameters testb6 = strupcase(strmid(key,0,2)) EQ 'B6' ind = where(testb6,countindb6) ; retrieveing the defaut values from the txt database tstb6 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'B6' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstb6,ctdatb6) if ctdatb6 ne 0 then begin T_eff_MS_b6 = double(T_eff_MS(indat)) R_star_MS_b6 = double(R_star_MS(indat)) endif IF countindb6 NE 0 THEN BEGIN lumclassb6 = uniq(lumclass(ind)) nlmb6 = n_elements(lumclassb6) indb6_1 = fltarr(nlmb6) indb6_2 = indb6_1 indb6_3 = indb6_1 indb6_4 = indb6_1 countindb6_1 = indb6_1 countindb6_2 = indb6_1 countindb6_3 = indb6_1 countindb6_4 = indb6_1 FOR i=0L,nlmb6-1 DO BEGIN popnumber+=1 countb+=1 indb6_1[i] = where(testb6 and params EQ 1, countindb6_1x); index of radius of this stellar population in the val array countindb6_1[i] = countindb6_1x indb6_2[i] = where(testb6 and params EQ 2, countindb6_2x); index of temperature of this stellar population in the val array countindb6_2[i] = countindb6_2x indb6_3[i] = where(testb6 and params EQ 3, countindb6_3x); index of distance of this stellar population in the val array countindb6_3[i] = countindb6_3x indb6_4[i] = where(testb6 and params EQ 4, countindb6_4x); index of the nstars of this stellar population in the val array countindb6_4[i] = countindb6_4x ENDFOR ENDIF ; Locating the B7 stellar population parameters testb7 = strupcase(strmid(key,0,2)) EQ 'B7' ind = where(testb7,countindb7) ; retrieveing the defaut values from the txt database tstb7 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'B7' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstb7,ctdatb7) if ctdatb7 ne 0 then begin T_eff_MS_b7 = double(T_eff_MS(indat)) R_star_MS_b7 = double(R_star_MS(indat)) endif IF countindb7 NE 0 THEN BEGIN lumclassb7 = uniq(lumclass(ind)) nlmb7 = n_elements(lumclassb7) indb7_1 = fltarr(nlmb7) indb7_2 = indb7_1 indb7_3 = indb7_1 indb7_4 = indb7_1 countindb7_1 = indb7_1 countindb7_2 = indb7_1 countindb7_3 = indb7_1 countindb7_4 = indb7_1 FOR i=0L,nlmb7-1 DO BEGIN popnumber+=1 countb+=1 indb7_1[i] = where(testb7 and params EQ 1, countindb7_1x); index of radius of this stellar population in the val array countindb7_1[i] = countindb7_1x indb7_2[i] = where(testb7 and params EQ 2, countindb7_2x); index of temperature of this stellar population in the val array countindb7_2[i] = countindb7_2x indb7_3[i] = where(testb7 and params EQ 3, countindb7_3x); index of distance of this stellar population in the val array countindb7_3[i] = countindb7_3x indb7_4[i] = where(testb7 and params EQ 4, countindb7_4x); index of the nstars of this stellar population in the val array countindb7_4[i] = countindb7_4x ENDFOR ENDIF ; Locating the B8 stellar population parameters testb8 = strupcase(strmid(key,0,2)) EQ 'B8' ind = where(testb8,countindb8) ; retrieveing the defaut values from the txt database tstb8 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'B8' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstb8,ctdatb8) if ctdatb8 ne 0 then begin T_eff_MS_b8 = double(T_eff_MS(indat)) R_star_MS_b8 = double(R_star_MS(indat)) endif IF countindb8 NE 0 THEN BEGIN lumclassb8 = uniq(lumclass(ind)) nlmb8 = n_elements(lumclassb8) indb8_1 = fltarr(nlmb8) indb8_2 = indb8_1 indb8_3 = indb8_1 indb8_4 = indb8_1 countindb8_1 = indb8_1 countindb8_2 = indb8_1 countindb8_3 = indb8_1 countindb8_4 = indb8_1 FOR i=0L,nlmb8-1 DO BEGIN popnumber+=1 countb+=1 indb8_1[i] = where(testb8 and params EQ 1, countindb8_1x); index of radius of this stellar population in the val array countindb8_1[i] = countindb8_1x indb8_2[i] = where(testb8 and params EQ 2, countindb8_2x); index of temperature of this stellar population in the val array countindb8_2[i] = countindb8_2x indb8_3[i] = where(testb8 and params EQ 3, countindb8_3x); index of distance of this stellar population in the val array countindb8_3[i] = countindb8_3x indb8_4[i] = where(testb8 and params EQ 4, countindb8_4x); index of the nstars of this stellar population in the val array countindb8_4[i] = countindb8_4x ENDFOR ENDIF ; Locating the B9 stellar population parameters testb9 = strupcase(strmid(key,0,2)) EQ 'B9' ind = where(testb9,countindb9) ; retrieveing the defaut values from the txt database tstb9 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'B9' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstb9,ctdatb9) if ctdatb9 ne 0 then begin T_eff_MS_b9 = double(T_eff_MS(indat)) R_star_MS_b9 = double(R_star_MS(indat)) endif IF countindb9 NE 0 THEN BEGIN lumclassb9 = uniq(lumclass(ind)) nlmb9 = n_elements(lumclassb9) indb9_1 = fltarr(nlmb9) indb9_2 = indb9_1 indb9_3 = indb9_1 indb9_4 = indb9_1 countindb9_1 = indb9_1 countindb9_2 = indb9_1 countindb9_3 = indb9_1 countindb9_4 = indb9_1 FOR i=0L,nlmb9-1 DO BEGIN popnumber+=1 countb+=1 indb9_1[i] = where(testb9 and params EQ 1, countindb9_1x); index of radius of this stellar population in the val array countindb9_1[i] = countindb9_1x indb9_2[i] = where(testb9 and params EQ 2, countindb9_2x); index of temperature of this stellar population in the val array countindb9_2[i] = countindb9_2x indb9_3[i] = where(testb9 and params EQ 3, countindb9_3x); index of distance of this stellar population in the val array countindb9_3[i] = countindb9_3x indb9_4[i] = where(testb9 and params EQ 4, countindb9_4x); index of the nstars of this stellar population in the val array countindb9_4[i] = countindb9_4x ENDFOR ENDIF ;====================================================== ; Locating the A0 stellar population parameters testa0 = strupcase(strmid(key,0,2)) EQ 'A0' ind = where(testa0,countinda0) ; retrieveing the defaut values from the txt database tsta0 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'A0' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tsta0,ctdata0) if ctdata0 ne 0 then begin T_eff_MS_a0 = double(T_eff_MS(indat)) R_star_MS_a0 = double(R_star_MS(indat)) endif IF countinda0 NE 0 THEN BEGIN lumclassA0 = uniq(lumclass(ind)) nlmA0 = n_elements(lumclassA0) indA0_1 = fltarr(nlmA0) indA0_2 = indA0_1 indA0_3 = indA0_1 indA0_4 = indA0_1 countindA0_1 = indA0_1 countindA0_2 = indA0_1 countindA0_3 = indA0_1 countindA0_4 = indA0_1 FOR i=0L,nlmA0-1 DO BEGIN popnumber+=1 counta+=1 inda0_1[i] = where(testa0 and params EQ 1, countinda0_1x); index of radius of this stellar population in the val array countinda0_1[i] = countinda0_1x inda0_2[i] = where(testa0 and params EQ 2, countinda0_2x); index of temperature of this stellar population in the val array countinda0_2[i] = countinda0_2x inda0_3[i] = where(testa0 and params EQ 3, countinda0_3x); index of distance of this stellar population in the val array countinda0_3[i] = countinda0_3x inda0_4[i] = where(testa0 and params EQ 4, countinda0_4x); index of the nstars of this stellar population in the val array countinda0_4[i] = countinda0_4x ENDFOR ENDIF ; Locating the A1 stellar population parameters testa1 = strupcase(strmid(key,0,2)) EQ 'A1' ind = where(testa1,countinda1) ; retrieveing the defaut values from the txt database tsta1 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'A1' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tsta1,ctdata1) if ctdata1 ne 0 then begin T_eff_MS_a1 = double(T_eff_MS(indat)) R_star_MS_a1 = double(R_star_MS(indat)) endif IF countinda1 NE 0 THEN BEGIN lumclassA1 = lumclass(ind) nlmA1 = n_elements(lumclassA1) indA1_1 = fltarr(nlmA1) indA1_2 = indA1_1 indA1_3 = indA1_1 indA1_4 = indA1_1 countindA1_1 = indA1_1 countindA1_2 = indA1_1 countindA1_3 = indA1_1 countindA1_4 = indA1_1 FOR i=0L,nlmA1-1 DO BEGIN popnumber+=1 counta+=1 inda1_1[i] = where(testa1 and params EQ 1, countinda1_1x); index of radius of this stellar population in the val array countinda1_1[i] = countinda1_1x inda1_2[i] = where(testa1 and params EQ 2, countinda1_2x); index of temperature of this stellar population in the val array countinda1_2[i] = countinda1_2x inda1_3[i] = where(testa1 and params EQ 3, countinda1_3x); index of distance of this stellar population in the val array countinda1_3[i] = countinda1_3x inda1_4[i] = where(testa1 and params EQ 4, countinda1_4x); index of the nstars of this stellar population in the val array countinda1_4[i] = countinda1_4x ENDFOR ENDIF ; Locating the A2 stellar population parameters testa2 = strupcase(strmid(key,0,2)) EQ 'A2' ind = where(testa2,countinda2) ; retrieveing the defaut values from the txt database tsta2 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'A2' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tsta2,ctdata2) if ctdata2 ne 0 then begin T_eff_MS_a2 = double(T_eff_MS(indat)) R_star_MS_a2 = double(R_star_MS(indat)) endif IF countinda2 NE 0 THEN BEGIN lumclassA2 = lumclass(ind) nlmA2 = n_elements(lumclassA2) indA2_1 = fltarr(nlmA2) indA2_2 = indA2_1 indA2_3 = indA2_1 indA2_4 = indA2_1 countindA2_1 = indA2_1 countindA2_2 = indA2_1 countindA2_3 = indA2_1 countindA2_4 = indA2_1 FOR i=0L,nlmA2-1 DO BEGIN popnumber+=1 counta+=1 inda2_1[i] = where(testa2 and params EQ 1, countinda2_1x); index of radius of this stellar population in the val array countinda2_1[i] = countinda2_1x inda2_2[i] = where(testa2 and params EQ 2, countinda2_2x); index of temperature of this stellar population in the val array countinda2_2[i] = countinda2_2x inda2_3[i] = where(testa2 and params EQ 3, countinda2_3x); index of distance of this stellar population in the val array countinda2_3[i] = countinda2_3x inda2_4[i] = where(testa2 and params EQ 4, countinda2_4x); index of the nstars of this stellar population in the val array countinda2_4[i] = countinda2_4x ENDFOR ENDIF ; Locating the A3 stellar population parameters testa3 = strupcase(strmid(key,0,2)) EQ 'A3' ind = where(testa3,countinda3) ; retrieveing the defaut values from the txt database tsta3 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'A3' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tsta3,ctdata3) if ctdata3 ne 0 then begin T_eff_MS_a3 = double(T_eff_MS(indat)) R_star_MS_a3 = double(R_star_MS(indat)) endif IF countinda3 NE 0 THEN BEGIN lumclassA3 = lumclass(ind) nlmA3 = n_elements(lumclassA3) indA3_1 = fltarr(nlmA3) indA3_2 = indA3_1 indA3_3 = indA3_1 indA3_4 = indA3_1 countindA3_1 = indA3_1 countindA3_2 = indA3_1 countindA3_3 = indA3_1 countindA3_4 = indA3_1 FOR i=0L,nlmA3-1 DO BEGIN popnumber+=1 counta+=1 inda3_1[i] = where(testa3 and params EQ 1, countinda3_1x); index of radius of this stellar population in the val array countinda3_1[i] = countinda3_1x inda3_2[i] = where(testa3 and params EQ 2, countinda3_2x); index of temperature of this stellar population in the val array countinda3_2[i] = countinda3_2x inda3_3[i] = where(testa3 and params EQ 3, countinda3_3x); index of distance of this stellar population in the val array countinda3_3[i] = countinda3_3x inda3_4[i] = where(testa3 and params EQ 4, countinda3_4x); index of the nstars of this stellar population in the val array countinda3_4[i] = countinda3_4x ENDFOR ENDIF ; Locating the A4 stellar population parameters testa4 = strupcase(strmid(key,0,2)) EQ 'A4' ind = where(testa4,countinda4) ; retrieveing the defaut values from the txt database tsta4 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'A4' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tsta4,ctdata4) if ctdata4 ne 0 then begin T_eff_MS_a4 = double(T_eff_MS(indat)) R_star_MS_a4 = double(R_star_MS(indat)) endif IF countinda4 NE 0 THEN BEGIN lumclassA4 = lumclass(ind) nlmA4 = n_elements(lumclassA4) indA4_1 = fltarr(nlmA4) indA4_2 = indA4_1 indA4_3 = indA4_1 indA4_4 = indA4_1 countindA4_1 = indA4_1 countindA4_2 = indA4_1 countindA4_3 = indA4_1 countindA4_4 = indA4_1 FOR i=0L,nlmA4-1 DO BEGIN popnumber+=1 counta+=1 inda4_1[i] = where(testa4 and params EQ 1, countinda4_1x); index of radius of this stellar population in the val array countinda4_1[i] = countinda4_1x inda4_2[i] = where(testa4 and params EQ 2, countinda4_2x); index of temperature of this stellar population in the val array countinda4_2[i] = countinda4_2x inda4_3[i] = where(testa4 and params EQ 3, countinda4_3x); index of distance of this stellar population in the val array countinda4_3[i] = countinda4_3x inda4_4[i] = where(testa4 and params EQ 4, countinda4_4x); index of the nstars of this stellar population in the val array countinda4_4[i] = countinda4_4x ENDFOR ENDIF ; Locating the A5 stellar population parameters testa5 = strupcase(strmid(key,0,2)) EQ 'A5' ind = where(testa5,countinda5) ; retrieveing the defaut values from the txt database tsta5 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'A5' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tsta5,ctdata5) if ctdata5 ne 0 then begin T_eff_MS_a5 = double(T_eff_MS(indat)) R_star_MS_a5 = double(R_star_MS(indat)) endif IF countinda5 NE 0 THEN BEGIN lumclassA5 = lumclass(ind) nlmA5 = n_elements(lumclassA5) indA5_1 = fltarr(nlmA5) indA5_2 = indA5_1 indA5_3 = indA5_1 indA5_4 = indA5_1 countindA5_1 = indA5_1 countindA5_2 = indA5_1 countindA5_3 = indA5_1 countindA5_4 = indA5_1 FOR i=0L,nlmA5-1 DO BEGIN popnumber+=1 counta+=1 inda5_1[i] = where(testa5 and params EQ 1, countinda5_1x); index of radius of this stellar population in the val array countinda5_1[i] = countinda5_1x inda5_2[i] = where(testa5 and params EQ 2, countinda5_2x); index of temperature of this stellar population in the val array countinda5_2[i] = countinda5_2x inda5_3[i] = where(testa5 and params EQ 3, countinda5_3x); index of distance of this stellar population in the val array countinda5_3[i] = countinda5_3x inda5_4[i] = where(testa5 and params EQ 4, countinda5_4x); index of the nstars of this stellar population in the val array countinda5_4[i] = countinda5_4x ENDFOR ENDIF ; Locating the A6 stellar population parameters testa6 = strupcase(strmid(key,0,2)) EQ 'A6' ind = where(testa6,countinda6) ; retrieveing the defaut values from the txt database tsta6 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'A6' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tsta6,ctdata6) if ctdata6 ne 0 then begin T_eff_MS_a6 = double(T_eff_MS(indat)) R_star_MS_a6 = double(R_star_MS(indat)) endif IF countinda6 NE 0 THEN BEGIN lumclassA6 = lumclass(ind) nlmA6 = n_elements(lumclassA6) indA6_1 = fltarr(nlmA6) indA6_2 = indA6_1 indA6_3 = indA6_1 indA6_4 = indA6_1 countindA6_1 = indA6_1 countindA6_2 = indA6_1 countindA6_3 = indA6_1 countindA6_4 = indA6_1 FOR i=0L,nlmA6-1 DO BEGIN popnumber+=1 counta+=1 inda6_1[i] = where(testa6 and params EQ 1, countinda6_1x); index of radius of this stellar population in the val array countinda6_1[i] = countinda6_1x inda6_2[i] = where(testa6 and params EQ 2, countinda6_2x); index of temperature of this stellar population in the val array countinda6_2[i] = countinda6_2x inda6_3[i] = where(testa6 and params EQ 3, countinda6_3x); index of distance of this stellar population in the val array countinda6_3[i] = countinda6_3x inda6_4[i] = where(testa6 and params EQ 4, countinda6_4x); index of the nstars of this stellar population in the val array countinda6_4[i] = countinda6_4x ENDFOR ENDIF ; Locating the A7 stellar population parameters testa7 = strupcase(strmid(key,0,2)) EQ 'A7' ind = where(testa7,countinda7) ; retrieveing the defaut values from the txt database tsta7 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'A7' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tsta7,ctdata7) if ctdata7 ne 0 then begin T_eff_MS_a7 = double(T_eff_MS(indat)) R_star_MS_a7 = double(R_star_MS(indat)) endif IF countinda7 NE 0 THEN BEGIN lumclassA7 = lumclass(ind) nlmA7 = n_elements(lumclassA7) indA7_1 = fltarr(nlmA7) indA7_2 = indA7_1 indA7_3 = indA7_1 indA7_4 = indA7_1 countindA7_1 = indA7_1 countindA7_2 = indA7_1 countindA7_3 = indA7_1 countindA7_4 = indA7_1 FOR i=0L,nlmA7-1 DO BEGIN popnumber+=1 counta+=1 inda7_1[i] = where(testa7 and params EQ 1, countinda7_1x); index of radius of this stellar population in the val array countinda7_1[i] = countinda7_1x inda7_2[i] = where(testa7 and params EQ 2, countinda7_2x); index of temperature of this stellar population in the val array countinda7_2[i] = countinda7_2x inda7_3[i] = where(testa7 and params EQ 3, countinda7_3x); index of distance of this stellar population in the val array countinda7_3[i] = countinda7_3x inda7_4[i] = where(testa7 and params EQ 4, countinda7_4x); index of the nstars of this stellar population in the val array countinda7_4[i] = countinda7_4x ENDFOR ENDIF ; Locating the A8 stellar population parameters testa8 = strupcase(strmid(key,0,2)) EQ 'A8' ind = where(testa8,countinda8) ; retrieveing the defaut values from the txt database tsta8 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'A8' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tsta8,ctdata8) if ctdata8 ne 0 then begin T_eff_MS_a8 = double(T_eff_MS(indat)) R_star_MS_a8 = double(R_star_MS(indat)) endif IF countinda8 NE 0 THEN BEGIN lumclassA8 = lumclass(ind) nlmA8 = n_elements(lumclassA8) indA8_1 = fltarr(nlmA8) indA8_2 = indA8_1 indA8_3 = indA8_1 indA8_4 = indA8_1 countindA8_1 = indA8_1 countindA8_2 = indA8_1 countindA8_3 = indA8_1 countindA8_4 = indA8_1 FOR i=0L,nlmA8-1 DO BEGIN popnumber+=1 counta+=1 inda8_1[i] = where(testa8 and params EQ 1, countinda8_1x); index of radius of this stellar population in the val array countinda8_1[i] = countinda8_1x inda8_2[i] = where(testa8 and params EQ 2, countinda8_2x); index of temperature of this stellar population in the val array countinda8_2[i] = countinda8_2x inda8_3[i] = where(testa8 and params EQ 3, countinda8_3x); index of distance of this stellar population in the val array countinda8_3[i] = countinda8_3x inda8_4[i] = where(testa8 and params EQ 4, countinda8_4x); index of the nstars of this stellar population in the val array countinda8_4[i] = countinda8_4x ENDFOR ENDIF ; Locating the A9 stellar population parameters testa9 = strupcase(strmid(key,0,2)) EQ 'A9' ind = where(testa9,countinda9) ; retrieveing the defaut values from the txt database tsta9 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'A9' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tsta9,ctdata9) if ctdata9 ne 0 then begin T_eff_MS_a9 = double(T_eff_MS(indat)) R_star_MS_a9 = double(R_star_MS(indat)) endif IF countinda9 NE 0 THEN BEGIN lumclassA9 = lumclass(ind) nlmA9 = n_elements(lumclassA9) indA9_1 = fltarr(nlmA9) indA9_2 = indA9_1 indA9_3 = indA9_1 indA9_4 = indA9_1 countindA9_1 = indA9_1 countindA9_2 = indA9_1 countindA9_3 = indA9_1 countindA9_4 = indA9_1 FOR i=0L,nlmA9-1 DO BEGIN popnumber+=1 counta+=1 inda9_1[i] = where(testa9 and params EQ 1, countinda9_1x); index of radius of this stellar population in the val array countinda9_1[i] = countinda9_1x inda9_2[i] = where(testa9 and params EQ 2, countinda9_2x); index of temperature of this stellar population in the val array countinda9_2[i] = countinda9_2x inda9_3[i] = where(testa9 and params EQ 3, countinda9_3x); index of distance of this stellar population in the val array countinda9_3[i] = countinda9_3x inda9_4[i] = where(testa9 and params EQ 4, countinda9_4x); index of the nstars of this stellar population in the val array countinda9_4[i] = countinda9_4x ENDFOR ENDIF ;====================================================== ; Locating the F0 stellar population parameters testf0 = strupcase(strmid(key,0,2)) EQ 'F0' ind = where(testf0,countindf0) ; retrieveing the defaut values from the txt database tstf0 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'F0' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstf0,ctdatf0) if ctdatf0 ne 0 then begin T_eff_MS_f0 = double(T_eff_MS(indat)) R_star_MS_f0 = double(R_star_MS(indat)) endif IF countindf0 NE 0 THEN BEGIN lumclassF0 = lumclass(ind) nlmF0 = n_elements(lumclassF0) indF0_1 = fltarr(nlmF0) indF0_2 = indF0_1 indF0_3 = indF0_1 indF0_4 = indF0_1 countindF0_1 = indF0_1 countindF0_2 = indF0_1 countindF0_3 = indF0_1 countindF0_4 = indF0_1 FOR i=0L,nlmF0-1 DO BEGIN popnumber+=1 countf+=1 indf0_1[i] = where(testf0 and params EQ 1, countindf0_1x); index of radius of this stellar population in the val array countindf0_1[i] = countindf0_1x indf0_2[i] = where(testf0 and params EQ 2, countindf0_2x); index of temperature of this stellar population in the val array countindf0_2[i] = countindf0_2x indf0_3[i] = where(testf0 and params EQ 3, countindf0_3x); index of distance of this stellar population in the val array countindf0_3[i] = countindf0_3x indf0_4[i] = where(testf0 and params EQ 4, countindf0_4x); index of the nstars of this stellar population in the val array countindf0_4[i] = countindf0_4x ENDFOR ENDIF ; Locating the F1 stellar population parameters testf1 = strupcase(strmid(key,0,2)) EQ 'F1' ind = where(testf1,countindf1) ; retrieveing the defaut values from the txt database tstf1 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'F1' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstf1,ctdatf1) if ctdatf1 ne 0 then begin T_eff_MS_f1 = double(T_eff_MS(indat)) R_star_MS_f1 = double(R_star_MS(indat)) endif IF countindf1 NE 0 THEN BEGIN lumclassF1 = lumclass(ind) nlmF1 = n_elements(lumclassF1) indF1_1 = fltarr(nlmF1) indF1_2 = indF1_1 indF1_3 = indF1_1 indF1_4 = indF1_1 countindF1_1 = indF1_1 countindF1_2 = indF1_1 countindF1_3 = indF1_1 countindF1_4 = indF1_1 FOR i=0L,nlmF1-1 DO BEGIN popnumber+=1 countf+=1 indf1_1[i] = where(testf1 and params EQ 1, countindf1_1x); index of radius of this stellar population in the val array countindf1_1[i] = countindf1_1x indf1_2[i] = where(testf1 and params EQ 2, countindf1_2x); index of temperature of this stellar population in the val array countindf1_2[i] = countindf1_2x indf1_3[i] = where(testf1 and params EQ 3, countindf1_3x); index of distance of this stellar population in the val array countindf1_3[i] = countindf1_3x indf1_4[i] = where(testf1 and params EQ 4, countindf1_4x); index of the nstars of this stellar population in the val array countindf1_4[i] = countindf1_4x ENDFOR ENDIF ; Locating the F2 stellar population parameters testf2 = strupcase(strmid(key,0,2)) EQ 'F2' ind = where(testf2,countindf2) ; retrieveing the defaut values from the txt database tstf2 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'F2' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstf2,ctdatf2) if ctdatf2 ne 0 then begin T_eff_MS_f2 = double(T_eff_MS(indat)) R_star_MS_f2 = double(R_star_MS(indat)) endif IF countindf2 NE 0 THEN BEGIN lumclassF2 = lumclass(ind) nlmF2 = n_elements(lumclassF2) indF2_1 = fltarr(nlmF2) indF2_2 = testf2 indF2_3 = testf2 indF2_4 = testf2 countindF2_1 = testf2 countindF2_2 = testf2 countindF2_3 = testf2 countindF2_4 = testf2 FOR i=0L,nlmF2-1 DO BEGIN popnumber+=1 countf+=1 indf2_1[i] = where(testf2 and params EQ 1, countindf2_1x); index of radius of this stellar population in the val array countindf2_1[i] = countindf2_1x indf2_2[i] = where(testf2 and params EQ 2, countindf2_2x); index of temperature of this stellar population in the val array countindf2_2[i] = countindf2_2x indf2_3[i] = where(testf2 and params EQ 3, countindf2_3x); index of distance of this stellar population in the val array countindf2_3[i] = countindf2_3x indf2_4[i] = where(testf2 and params EQ 4, countindf2_4x); index of the nstars of this stellar population in the val array countindf2_4[i] = countindf2_4x ENDFOR ENDIF ; Locating the F3 stellar population parameters testf3 = strupcase(strmid(key,0,2)) EQ 'F3' ind = where(testf3,countindf3) ; retrieveing the defaut values from the txt database tstf3 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'F3' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstf3,ctdatf3) if ctdatf3 ne 0 then begin T_eff_MS_f3 = double(T_eff_MS(indat)) R_star_MS_f3 = double(R_star_MS(indat)) endif IF countindf3 NE 0 THEN BEGIN lumclassF3 = lumclass(ind) nlmF3 = n_elements(lumclassF3) indF3_1 = fltarr(nlmF3) indF3_2 = indF3_1 indF3_3 = indF3_1 indF3_4 = indF3_1 countindF3_1 = indF3_1 countindF3_2 = indF3_1 countindF3_3 = indF3_1 countindF3_4 = indF3_1 FOR i=0L,nlmF3-1 DO BEGIN popnumber+=1 countf+=1 indf3_1[i] = where(testf3 and params EQ 1, countindf3_1x); index of radius of this stellar population in the val array countindf3_1[i] = countindf3_1x indf3_2[i] = where(testf3 and params EQ 2, countindf3_2x); index of temperature of this stellar population in the val array countindf3_2[i] = countindf3_2x indf3_3[i] = where(testf3 and params EQ 3, countindf3_3x); index of distance of this stellar population in the val array countindf3_3[i] = countindf3_3x indf3_4[i] = where(testf3 and params EQ 4, countindf3_4x); index of the nstars of this stellar population in the val array countindf3_4[i] = countindf3_4x ENDFOR ENDIF ; Locating the F4 stellar population parameters testf4 = strupcase(strmid(key,0,2)) EQ 'F4' ind = where(testf4,countindf4) ; retrieveing the defaut values from the txt database tstf4 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'F4' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstf4,ctdatf4) if ctdatf4 ne 0 then begin T_eff_MS_f4 = double(T_eff_MS(indat)) R_star_MS_f4 = double(R_star_MS(indat)) endif IF countindf4 NE 0 THEN BEGIN lumclassF4 = lumclass(ind) nlmF4 = n_elements(lumclassF4) indF4_1 = fltarr(nlmF4) indF4_2 = indF4_1 indF4_3 = indF4_1 indF4_4 = indF4_1 countindF4_1 = indF4_1 countindF4_2 = indF4_1 countindF4_3 = indF4_1 countindF4_4 = indF4_1 FOR i=0L,nlmF4-1 DO BEGIN popnumber+=1 countf+=1 indf4_1[i] = where(testf4 and params EQ 1, countindf4_1x); index of radius of this stellar population in the val array countindf4_1[i] = countindf4_1x indf4_2[i] = where(testf4 and params EQ 2, countindf4_2x); index of temperature of this stellar population in the val array countindf4_2[i] = countindf4_2x indf4_3[i] = where(testf4 and params EQ 3, countindf4_3x); index of distance of this stellar population in the val array countindf4_3[i] = countindf4_3x indf4_4[i] = where(testf4 and params EQ 4, countindf4_4x); index of the nstars of this stellar population in the val array countindf4_4[i] = countindf4_4x ENDFOR ENDIF ; Locating the F5 stellar population parameters testf5 = strupcase(strmid(key,0,2)) EQ 'F5' ind = where(testf5,countindf5) ; retrieveing the defaut values from the txt database tstf5 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'F5' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstf5,ctdatf5) if ctdatf5 ne 0 then begin T_eff_MS_f5 = double(T_eff_MS(indat)) R_star_MS_f5 = double(R_star_MS(indat)) endif IF countindf5 NE 0 THEN BEGIN lumclassF5 = lumclass(ind) nlmF5 = n_elements(lumclassF5) indF5_1 = fltarr(nlmF5) indF5_2 = indF5_1 indF5_3 = indF5_1 indF5_4 = indF5_1 countindF5_1 = indF5_1 countindF5_2 = indF5_1 countindF5_3 = indF5_1 countindF5_4 = indF5_1 FOR i=0L,nlmF5-1 DO BEGIN popnumber+=1 countf+=1 indf5_1[i] = where(testf5 and params EQ 1, countindf5_1x); index of radius of this stellar population in the val array countindf5_1[i] = countindf5_1x indf5_2[i] = where(testf5 and params EQ 2, countindf5_2x); index of temperature of this stellar population in the val array countindf5_2[i] = countindf5_2x indf5_3[i] = where(testf5 and params EQ 3, countindf5_3x); index of distance of this stellar population in the val array countindf5_3[i] = countindf5_3x indf5_4[i] = where(testf5 and params EQ 4, countindf5_4x); index of the nstars of this stellar population in the val array countindf5_4[i] = countindf5_4x ENDFOR ENDIF ; Locating the F6 stellar population parameters testf6 = strupcase(strmid(key,0,2)) EQ 'F6' ind = where(testf6,countindf6) ; retrieveing the defaut values from the txt database tstf6 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'F6' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstf6,ctdatf6) if ctdatf6 ne 0 then begin T_eff_MS_f6 = double(T_eff_MS(indat)) R_star_MS_f6 = double(R_star_MS(indat)) endif IF countindf6 NE 0 THEN BEGIN lumclassF6 = lumclass(ind) nlmF6 = n_elements(lumclassF6) indF6_1 = fltarr(nlmF6) indF6_2 = indF6_1 indF6_3 = indF6_1 indF6_4 = indF6_1 countindF6_1 = indF6_1 countindF6_2 = indF6_1 countindF6_3 = indF6_1 countindF6_4 = indF6_1 FOR i=0L,nlmF6-1 DO BEGIN popnumber+=1 countf+=1 indf6_1[i] = where(testf6 and params EQ 1, countindf6_1x); index of radius of this stellar population in the val array countindf6_1[i] = countindf6_1x indf6_2[i] = where(testf6 and params EQ 2, countindf6_2x); index of temperature of this stellar population in the val array countindf6_2[i] = countindf6_2x indf6_3[i] = where(testf6 and params EQ 3, countindf6_3x); index of distance of this stellar population in the val array countindf6_3[i] = countindf6_3x indf6_4[i] = where(testf6 and params EQ 4, countindf6_4x); index of the nstars of this stellar population in the val array countindf6_4[i] = countindf6_4x ENDFOR ENDIF ; Locating the F7 stellar population parameters testf7 = strupcase(strmid(key,0,2)) EQ 'F7' ind = where(testf7,countindf7) ; retrieveing the defaut values from the txt database tstf7 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'F7' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstf7,ctdatf7) if ctdatf7 ne 0 then begin T_eff_MS_f7 = double(T_eff_MS(indat)) R_star_MS_f7 = double(R_star_MS(indat)) endif IF countindf7 NE 0 THEN BEGIN umclassF7 = lumclass(ind) nlmF7 = n_elements(lumclassF7) indF7_1 = fltarr(nlmF7) indF7_2 = indF7_1 indF7_3 = indF7_1 indF7_4 = indF7_1 countindF7_1 = indF7_1 countindF7_2 = indF7_1 countindF7_3 = indF7_1 countindF7_4 = indF7_1 FOR i=0L,nlmF7-1 DO BEGIN popnumber+=1 countf+=1 indf7_1[i] = where(testf7 and params EQ 1, countindf7_1x); index of radius of this stellar population in the val array countindf7_1[i] = countindf7_1x indf7_2[i] = where(testf7 and params EQ 2, countindf7_2x); index of temperature of this stellar population in the val array countindf7_2[i] = countindf7_2x indf7_3[i] = where(testf7 and params EQ 3, countindf7_3x); index of distance of this stellar population in the val array countindf7_3[i] = countindf7_3x indf7_4[i] = where(testf7 and params EQ 4, countindf7_4x); index of the nstars of this stellar population in the val array countindf7_4[i] = countindf7_4x ENDFOR ENDIF ; Locating the F8 stellar population parameters testf8 = strupcase(strmid(key,0,2)) EQ 'F8' ind = where(testf8,countindf8) ; retrieveing the defaut values from the txt database tstf8 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'F8' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstf8,ctdatf8) if ctdatf8 ne 0 then begin T_eff_MS_f8 = double(T_eff_MS(indat)) R_star_MS_f8 = double(R_star_MS(indat)) endif IF countindf8 NE 0 THEN BEGIN lumclassF8 = lumclass(ind) nlmF8 = n_elements(lumclassF8) indF8_1 = fltarr(nlmF8) indF8_2 = indF8_1 indF8_3 = indF8_1 indF8_4 = indF8_1 countindF8_1 = indF8_1 countindF8_2 = indF8_1 countindF8_3 = indF8_1 countindF8_4 = indF8_1 FOR i=0L,nlmF8-1 DO BEGIN popnumber+=1 countf+=1 indf8_1[i] = where(testf8 and params EQ 1, countindf8_1x); index of radius of this stellar population in the val array countindf8_1[i] = countindf8_1x indf8_2[i] = where(testf8 and params EQ 2, countindf8_2x); index of temperature of this stellar population in the val array countindf8_2[i] = countindf8_2x indf8_3[i] = where(testf8 and params EQ 3, countindf8_3x); index of distance of this stellar population in the val array countindf8_3[i] = countindf8_3x indf8_4[i] = where(testf8 and params EQ 4, countindf8_4x); index of the nstars of this stellar population in the val array countindf8_4[i] = countindf8_4x ENDFOR ENDIF ; Locating the F9 stellar population parameters testf9 = strupcase(strmid(key,0,2)) EQ 'F9' ind = where(testf9,countindf9) ; retrieveing the defaut values from the txt database tstf9 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'F9' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstf9,ctdatf9) if ctdatf9 ne 0 then begin T_eff_MS_f9 = double(T_eff_MS(indat)) R_star_MS_f9 = double(R_star_MS(indat)) endif IF countindf9 NE 0 THEN BEGIN lumclassF9 = lumclass(ind) nlmF9 = n_elements(lumclassF9) indF9_1 = fltarr(nlmF9) indF9_2 = indF9_1 indF9_3 = indF9_1 indF9_4 = indF9_1 countindF9_1 = indF9_1 countindF9_2 = indF9_1 countindF9_3 = indF9_1 countindF9_4 = indF9_1 FOR i=0L,nlmF9-1 DO BEGIN popnumber+=1 countf+=1 indf9_1[i] = where(testf9 and params EQ 1, countindf9_1x); index of radius of this stellar population in the val array countindf9_1[i] = countindf9_1x indf9_2[i] = where(testf9 and params EQ 2, countindf9_2x); index of temperature of this stellar population in the val array countindf9_2[i] = countindf9_2x indf9_3[i] = where(testf9 and params EQ 3, countindf9_3x); index of distance of this stellar population in the val array countindf9_3[i] = countindf9_3x indf9_4[i] = where(testf9 and params EQ 4, countindf9_4x); index of the nstars of this stellar population in the val array countindf9_4[i] = countindf9_4x ENDFOR ENDIF ;====================================================== ; Locating the G0 stellar population parameters testg0 = strupcase(strmid(key,0,2)) EQ 'G0' ind = where(testg0,countindg0) ; retrieveing the defaut values from the txt database tstg0 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'G0' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstg0,ctdatg0) if ctdatg0 ne 0 then begin T_eff_MS_g0 = double(T_eff_MS(indat)) R_star_MS_g0 = double(R_star_MS(indat)) endif IF countindg0 NE 0 THEN BEGIN lumclassG0 = lumclass(ind) nlmG0 = n_elements(lumclassG0) indG0_1 = fltarr(nlmG0) indG0_2 = indG0_1 indG0_3 = indG0_1 indG0_4 = indG0_1 countindG0_1 = indG0_1 countindG0_2 = indG0_1 countindG0_3 = indG0_1 countindG0_4 = indG0_1 FOR i=0L,nlmG0-1 DO BEGIN popnumber+=1 countg+=1 indg0_1[i] = where(testg0 and params EQ 1, countindg0_1x); index of radius of this stellar population in the val array countindg0_1[i] = countindg0_1x indg0_2[i] = where(testg0 and params EQ 2, countindg0_2x); index of temperature of this stellar population in the val array countindg0_2[i] = countindg0_2x indg0_3[i] = where(testg0 and params EQ 3, countindg0_3x); index of distance of this stellar population in the val array countindg0_3[i] = countindg0_3x indg0_4[i] = where(testg0 and params EQ 4, countindg0_4x); index of the nstars of this stellar population in the val array countindg0_4[i] = countindg0_4x ENDFOR ENDIF ; Locating the G1 stellar population parameters testg1 = strupcase(strmid(key,0,2)) EQ 'G1' ind = where(testg1,countindg1) ; retrieveing the defaut values from the txt database tstg1 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'G1' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstg1,ctdatg1) if ctdatg1 ne 0 then begin T_eff_MS_g1 = double(T_eff_MS(indat)) R_star_MS_g1 = double(R_star_MS(indat)) endif IF countindg1 NE 0 THEN BEGIN lumclassG1 = lumclass(ind) nlmG1 = n_elements(lumclassG1) indG1_1 = fltarr(nlmG1) indG1_2 = indG1_1 indG1_3 = indG1_1 indG1_4 = indG1_1 countindG1_1 = indG1_1 countindG1_2 = indG1_1 countindG1_3 = indG1_1 countindG1_4 = indG1_1 FOR i=0L,nlmG1-1 DO BEGIN popnumber+=1 countg+=1 indg1_1[i] = where(testg1 and params EQ 1, countindg1_1x); index of radius of this stellar population in the val array countindg1_1[i] = countindg1_1x indg1_2[i] = where(testg1 and params EQ 2, countindg1_2x); index of temperature of this stellar population in the val array countindg1_2[i] = countindg1_2x indg1_3[i] = where(testg1 and params EQ 3, countindg1_3x); index of distance of this stellar population in the val array countindg1_3[i] = countindg1_3x indg1_4[i] = where(testg1 and params EQ 4, countindg1_4x); index of the nstars of this stellar population in the val array countindg1_4[i] = countindg1_4x ENDFOR ENDIF ; Locating the G2 stellar population parameters testg2 = strupcase(strmid(key,0,2)) EQ 'G2' ind = where(testg2,countindg2) ;sun-like star ; retrieveing the defaut values from the txt database tstg2 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'G2' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstg2,ctdatg2) if ctdatg2 ne 0 then begin T_eff_MS_g2 = double(T_eff_MS(indat)) R_star_MS_g2 = double(R_star_MS(indat)) endif IF countindg2 NE 0 THEN BEGIN lumclassG2 = lumclass(ind) nlmG2 = n_elements(lumclassG2) indG2_1 = fltarr(nlmG2) indG2_2 = indG2_1 indG2_3 = indG2_1 indG2_4 = indG2_1 countindG2_1 = indG2_1 countindG2_2 = indG2_1 countindG2_3 = indG2_1 countindG2_4 = indG2_1 FOR i=0L,nlmG2-1 DO BEGIN popnumber+=1 countg+=1 indg2_1[i] = where(testg2 and params EQ 1, countindg2_1x); index of radius of this stellar population in the val array countindg2_1[i] = countindg2_1x indg2_2[i] = where(testg2 and params EQ 2, countindg2_2x); index of temperature of this stellar population in the val array countindg2_2[i] = countindg2_2x indg2_3[i] = where(testg2 and params EQ 3, countindg2_3x); index of distance of this stellar population in the val array countindg2_3[i] = countindg2_3x indg2_4[i] = where(testg2 and params EQ 4, countindg2_4x); index of the nstars of this stellar population in the val array countindg2_4[i] = countindg2_4x ENDFOR ENDIF ; Locating the G3 stellar population parameters testg3 = strupcase(strmid(key,0,2)) EQ 'G3' ind = where(testg3,countindg3) ; retrieveing the defaut values from the txt database tstg3 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'G3' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstg3,ctdatg3) if ctdatg3 ne 0 then begin T_eff_MS_g3 = double(T_eff_MS(indat)) R_star_MS_g3 = double(R_star_MS(indat)) endif IF countindg3 NE 0 THEN BEGIN lumclassG3 = lumclass(ind) nlmG3 = n_elements(lumclassG3) indG3_1 = fltarr(nlmG3) indG3_2 = indG3_1 indG3_3 = indG3_1 indG3_4 = indG3_1 countindG3_1 = indG3_1 countindG3_2 = indG3_1 countindG3_3 = indG3_1 countindG3_4 = indG3_1 FOR i=0L,nlmG3-1 DO BEGIN popnumber+=1 countg+=1 indg3_1[i] = where(testg3 and params EQ 1, countindg3_1x); index of radius of this stellar population in the val array countindg3_1[i] = countindg3_1x indg3_2[i] = where(testg3 and params EQ 2, countindg3_2x); index of temperature of this stellar population in the val array countindg3_2[i] = countindg3_2x indg3_3[i] = where(testg3 and params EQ 3, countindg3_3x); index of distance of this stellar population in the val array countindg3_3[i] = countindg3_3x indg3_4[i] = where(testg3 and params EQ 4, countindg3_4x); index of the nstars of this stellar population in the val array countindg3_4[i] = countindg3_4x ENDFOR ENDIF ; Locating the G4 stellar population parameters testg3 = strupcase(strmid(key,0,2)) EQ 'G4' ind = where(testg3,countindg4) ; retrieveing the defaut values from the txt database tstg4 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'G4' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstg4,ctdatg4) if ctdatg4 ne 0 then begin T_eff_MS_g4 = double(T_eff_MS(indat)) R_star_MS_g4 = double(R_star_MS(indat)) endif IF countindg4 NE 0 THEN BEGIN lumclassG4 = lumclass(ind) nlmG4 = n_elements(lumclassG4) indG4_1 = fltarr(nlmG4) indG4_2 = indG4_1 indG4_3 = indG4_1 indG4_4 = indG4_1 countindG4_1 = indG4_1 countindG4_2 = indG4_1 countindG4_3 = indG4_1 countindG4_4 = indG4_1 FOR i=0L,nlmG4-1 DO BEGIN popnumber+=1 countg+=1 indg4_1[i] = where(testg3 and params EQ 1, countindg4_1x); index of radius of this stellar population in the val array countindg4_1[i] = countindg4_1x indg4_2[i] = where(testg3 and params EQ 2, countindg4_2x); index of temperature of this stellar population in the val array countindg4_2[i] = countindg4_2x indg4_3[i] = where(testg3 and params EQ 3, countindg4_3x); index of distance of this stellar population in the val array countindg4_3[i] = countindg4_3x indg4_4[i] = where(testg3 and params EQ 4, countindg4_4x); index of the nstars of this stellar population in the val array countindg4_4[i] = countindg4_4x ENDFOR ENDIF ; Locating the G5 stellar population parameters testg5 = strupcase(strmid(key,0,2)) EQ 'G5' ind = where(testg5,countindg5) ; retrieveing the defaut values from the txt database tstg5 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'G5' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstg5,ctdatg5) if ctdatg5 ne 0 then begin T_eff_MS_g5 = double(T_eff_MS(indat)) R_star_MS_g5 = double(R_star_MS(indat)) endif IF countindg5 NE 0 THEN BEGIN lumclassG5 = lumclass(ind) nlmG5 = n_elements(lumclassG5) indG5_1 = fltarr(nlmG5) indG5_2 = indG5_1 indG5_3 = indG5_1 indG5_4 = indG5_1 countindG5_1 = indG5_1 countindG5_2 = indG5_1 countindG5_3 = indG5_1 countindG5_4 = indG5_1 FOR i=0L,nlmG5-1 DO BEGIN popnumber+=1 countg+=1 indg5_1[i] = where(testg5 and params EQ 1, countindg5_1x); index of radius of this stellar population in the val array countindg5_1[i] = countindg5_1x indg5_2[i] = where(testg5 and params EQ 2, countindg5_2x); index of temperature of this stellar population in the val array countindg5_2[i] = countindg5_2x indg5_3[i] = where(testg5 and params EQ 3, countindg5_3x); index of distance of this stellar population in the val array countindg5_3[i] = countindg5_3x indg5_4[i] = where(testg5 and params EQ 4, countindg5_4x); index of the nstars of this stellar population in the val array countindg5_4[i] = countindg5_4x ENDFOR ENDIF ; Locating the G6 stellar population parameters testg6 = strupcase(strmid(key,0,2)) EQ 'G6' ind = where(testg6,countindg6) ; retrieveing the defaut values from the txt database tstg6 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'G6' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstg6,ctdatg6) if ctdatg6 ne 0 then begin T_eff_MS_g6 = double(T_eff_MS(indat)) R_star_MS_g6 = double(R_star_MS(indat)) endif IF countindg6 NE 0 THEN BEGIN lumclassG6 = lumclass(ind) nlmG6 = n_elements(lumclassG6) indG6_1 = fltarr(nlmG6) indG6_2 = indG6_1 indG6_3 = indG6_1 indG6_4 = indG6_1 countindG6_1 = indG6_1 countindG6_2 = indG6_1 countindG6_3 = indG6_1 countindG6_4 = indG6_1 FOR i=0L,nlmG6-1 DO BEGIN popnumber+=1 countg+=1 indg6_1[i] = where(testg6 and params EQ 1, countindg6_1x); index of radius of this stellar population in the val array countindg6_1[i] = countindg6_1x indg6_2[i] = where(testg6 and params EQ 2, countindg6_2x); index of temperature of this stellar population in the val array countindg6_2[i] = countindg6_2x indg6_3[i] = where(testg6 and params EQ 3, countindg6_3x); index of distance of this stellar population in the val array countindg6_3[i] = countindg6_3x indg6_4[i] = where(testg6 and params EQ 4, countindg6_4x); index of the nstars of this stellar population in the val array countindg6_4[i] = countindg6_4x ENDFOR ENDIF ; Locating the G7 stellar population parameters testg7 = strupcase(strmid(key,0,2)) EQ 'G7' ind = where(testg7,countindg7) ; retrieveing the defaut values from the txt database tstg7 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'G7' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstg7,ctdatg7) if ctdatg7 ne 0 then begin T_eff_MS_g7 = double(T_eff_MS(indat)) R_star_MS_g7 = double(R_star_MS(indat)) endif IF countindg7 NE 0 THEN BEGIN lumclassG7 = lumclass(ind) nlmG7 = n_elements(lumclassG7) indG7_1 = fltarr(nlmG7) indG7_2 = indG7_1 indG7_3 = indG7_1 indG7_4 = indG7_1 countindG7_1 = indG7_1 countindG7_2 = indG7_1 countindG7_3 = indG7_1 countindG7_4 = indG7_1 FOR i=0L,nlmG7-1 DO BEGIN popnumber+=1 countg+=1 indg7_1[i] = where(testg7 and params EQ 1, countindg7_1x); index of radius of this stellar population in the val array countindg7_1[i] = countindg7_1x indg7_2[i] = where(testg7 and params EQ 2, countindg7_2x); index of temperature of this stellar population in the val array countindg7_2[i] = countindg7_2x indg7_3[i] = where(testg7 and params EQ 3, countindg7_3x); index of distance of this stellar population in the val array countindg7_3[i] = countindg7_3x indg7_4[i] = where(testg7 and params EQ 4, countindg7_4x); index of the nstars of this stellar population in the val array countindg7_4[i] = countindg7_4x ENDFOR ENDIF ; Locating the G8 stellar population parameters testg8 = strupcase(strmid(key,0,2)) EQ 'G8' ind = where(testg8,countindg8) ; retrieveing the defaut values from the txt database tstg8 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'G8' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstg8,ctdatg8) if ctdatg8 ne 0 then begin T_eff_MS_g8 = double(T_eff_MS(indat)) R_star_MS_g8 = double(R_star_MS(indat)) endif IF countindg8 NE 0 THEN BEGIN lumclassG8 = lumclass(ind) nlmG8 = n_elements(lumclassG8) indG8_1 = fltarr(nlmG8) indG8_2 = indG8_1 indG8_3 = indG8_1 indG8_4 = indG8_1 countindG8_1 = indG8_1 countindG8_2 = indG8_1 countindG8_3 = indG8_1 countindG8_4 = indG8_1 FOR i=0L,nlmG8-1 DO BEGIN popnumber+=1 countg+=1 indg8_1[i] = where(testg8 and params EQ 1, countindg8_1x); index of radius of this stellar population in the val array countindg8_1[i] = countindg8_1x indg8_2[i] = where(testg8 and params EQ 2, countindg8_2x); index of temperature of this stellar population in the val array countindg8_2[i] = countindg8_2x indg8_3[i] = where(testg8 and params EQ 3, countindg8_3x); index of distance of this stellar population in the val array countindg8_3[i] = countindg8_3x indg8_4[i] = where(testg8 and params EQ 4, countindg8_4x); index of the nstars of this stellar population in the val array countindg8_4[i] = countindg8_4x ENDFOR ENDIF ; Locating the G9 stellar population parameters testg9 = strupcase(strmid(key,0,2)) EQ 'G9' ind = where(testg9,countindg9) ; retrieveing the defaut values from the txt database tstg9 = strtrim(strmid(SpecType_MS,0,2),2) EQ 'G9' and strtrim(strmid(SpecType_MS,2,1),2) NE '.' indat = where(tstg9,ctdatg9) if ctdatg9 ne 0 then begin T_eff_MS_g9 = double(T_eff_MS(indat)) R_star_MS_g9 = double(R_star_MS(indat)) endif IF countindg9 NE 0 THEN BEGIN lumclassG9 = lumclass(ind) nlmG9 = n_elements(lumclassG9) indG9_1 = fltarr(nlmG9) indG9_2 = indG9_1 indG9_3 = indG9_1 indG9_4 = indG9_1 countindG9_1 = indG9_1 countindG9_2 = indG9_1 countindG9_3 = indG9_1 countindG9_4 = indG9_1 FOR i=0L,nlmG9-1 DO BEGIN popnumber+=1 countg+=1 indg9_1[i] = where(testg9 and params EQ 1, countindg9_1x); index of radius of this stellar population in the val array countindg9_1[i] = countindg9_1x indg9_2[i] = where(testg9 and params EQ 2, countindg9_2x); index of temperature of this stellar population in the val array countindg9_2[i] = countindg9_2x indg9_3[i] = where(testg9 and params EQ 3, countindg9_3x); index of distance of this stellar population in the val array countindg9_3[i] = countindg9_3x indg9_4[i] = where(testg9 and params EQ 4, countindg9_4x); index of the nstars of this stellar population in the val array countindg9_4[i] = countindg9_4x ENDFOR ENDIF comp_pop=replicate(one_pop,popnumber) ; Replication of the default and initialized stellar population structure to create the composite stellar population structure. All the tag values are by default set to te null value ;Filling the tags of the compsite stellar population structure using the indices defined above FOR i=0L,popnumber-1 DO BEGIN ; Looping over all the structure fields IF counto NE 0 THEN BEGIN ; Filling the structure of the O STELLAR POPULATION(S) FOR k=0L,counto-1 DO BEGIN IF countindo3 NE 0 THEN BEGIN FOR l=0L,nlmO3-1 DO BEGIN comp_pop(k+i+l).popid = 'O3'+lumclasso3(l)+'_stellar_population' IF countindo3_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indo3_1[l]))(0)*rsun2cm ELSE BEGIN ; IF lumclasso3(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ;Hypergiants or extremely luminous supergiants ; IF lumclasso3(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ;Luminos supergiants ; IF lumclasso3(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ;Intermediate-size luminous supergiants ; IF lumclasso3(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ;Less luminous supergiants ; IF lumclasso3(l) EQ 'II' THEN comp_pop(k+i+l).radius = ;Bright giants ; IF lumclasso3(l) EQ 'III' THEN comp_pop(k+i+l).radius = ;Normal giants ; IF lumclasso3(l) EQ 'IV' THEN comp_pop(k+i+l).radius = ;Subgiants IF lumclasso3(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_o3*rsun2cm;15.0*rsun2cm ;(cm,cgs);Dwarfs: Main-sequence stars ; IF lumclasso3(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ;Subdwarfs ; IF lumclasso3(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ;White dwarfs ENDELSE IF countindo3_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indo3_2[l]))(0) ELSE BEGIN ; IF lumclasso3(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasso3(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasso3(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasso3(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasso3(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasso3(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasso3(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasso3(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_o3;44.9e3 ;(K) ; IF lumclasso3(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasso3(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindo3_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indo3_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 ; should be lower limited in the general procedure because this is the average distance between two stars IF countindo3_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indo3_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmo3 ENDIF IF countindo4 NE 0 THEN BEGIN FOR l=0L,nlmO4-1 DO BEGIN comp_pop(k+i+l).popid = 'O4'+lumclasso4(l)+'_stellar_population' IF countindo4_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indo4_1[l]))(0)*rsun2cm ELSE BEGIN ; IF lumclasso4(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ;Hypergiants or extremely luminous supergiants ; IF lumclasso4(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ;Luminos supergiants ; IF lumclasso4(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ;Intermediate-size luminous supergiants ; IF lumclasso4(l) EQ 'IB' THEN comp_pop(k+i+l).radius ;Less luminous supergiants ; IF lumclasso4(l) EQ 'II' THEN comp_pop(k+i+l).radius = ;Bright giants ; IF lumclasso4(l) EQ 'III' THEN comp_pop(k+i+l).radius = ;Normal giants ; IF lumclasso4(l) EQ 'IV' THEN comp_pop(k+i+l).radius = ;Subgiants IF lumclasso4(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_o4*rsun2cm;13.43*rsun2cm ;(cm,cgs);Dwarfs: Main-sequence stars ; IF lumclasso4(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ;Subdwarfs ; IF lumclasso4(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ;White dwarfs ENDELSE IF countindo4_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indo4_2[l]))(0) ELSE BEGIN ; IF lumclasso4(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasso4(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasso4(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasso4(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasso4(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasso4(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasso4(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasso4(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_o4;42.9e3 ; IF lumclasso4(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasso4(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindo4_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indo4_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindo4_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indo4_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmo4 ENDIF IF countindo5 NE 0 THEN BEGIN FOR l=0L,nlmo5-1 DO BEGIN comp_pop(k+i+l).popid = 'O5'+lumclasso5(l)+'_stellar_population' IF countindo5_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indo5_1[l]))(0)*rsun2cm ELSE BEGIN ; IF lumclasso5(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ;Hypergiants or extremely luminous supergiants ; IF lumclasso5(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ;Luminos supergiants ; IF lumclasso5(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ;Intermediate-size luminous supergiants ; IF lumclasso5(l) EQ 'IB' THEN comp_pop(k+i+l).radius ;Less luminous supergiants ; IF lumclasso5(l) EQ 'II' THEN comp_pop(k+i+l).radius = ;Bright giants ; IF lumclasso5(l) EQ 'III' THEN comp_pop(k+i+l).radius = ;Normal giants ; IF lumclasso5(l) EQ 'IV' THEN comp_pop(k+i+l).radius = ;Subgiants IF lumclasso5(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_o5*rsun2cm;12*rsun2cm ;(cm,cgs);Dwarfs: Main-sequence stars ; IF lumclasso5(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ;Subdwarfs ; IF lumclasso5(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ;White dwarfs ENDELSE IF countindo5_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indo5_2[l]))(0) ELSE BEGIN ; IF lumclasso5(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasso5(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasso5(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasso5(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasso5(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasso5(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasso5(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasso5(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_o5;41.4e3 ; IF lumclasso5(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasso5(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindo5_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indo5_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindo5_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indo5_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmo5 ENDIF IF countindo6 NE 0 THEN BEGIN FOR l=0L,nlmo6-1 DO BEGIN comp_pop(k+i+l).popid = 'O6'+lumclasso6(l)+'_stellar_population' IF countindo6_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indo6_1[l]))(0)*rsun2cm ELSE BEGIN ; IF lumclasso6(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclasso6(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclasso6(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclasso6(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclasso6(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclasso6(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclasso6(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclasso6(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_o6*rsun2cm;10.71*rsun2cm ; IF lumclasso6(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclasso6(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindo6_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indo6_2[l]))(0) ELSE BEGIN ; IF lumclasso6(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasso6(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasso6(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasso6(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasso6(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasso6(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasso6(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasso6(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_o6;39.5e3 ; IF lumclasso6(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasso6(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindo6_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indo6_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindo6_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indo6_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmo6 ENDIF IF countindo7 NE 0 THEN BEGIN FOR l=0L,nlmo7-1 DO BEGIN comp_pop(k+l+i).popid = 'O7'+lumclasso7(l)+'_stellar_population' IF countindo7_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indo7_1[l]))(0) ELSE BEGIN ; IF lumclasso7(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclasso7(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclasso7(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclasso7(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclasso7(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclasso7(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclasso7(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclasso7(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_o7*rsun2cm;9.52*rsun2cm ; IF lumclasso7(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclasso7(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindo7_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indo7_2[l]))(0) ELSE BEGIN ; IF lumclasso7(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasso7(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasso7(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasso7(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasso7(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasso7(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasso7(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasso7(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_o7;37.1e3 ; IF lumclasso7(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasso7(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindo7_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indo7_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindo7_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indo7_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmo7 ENDIF IF countindo8 NE 0 THEN BEGIN FOR l=0L,nlmo8-1 DO BEGIN comp_pop(k+i+l).popid = 'O8'+lumclasso8(l)+'_stellar_population' IF countindo8_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indo8_1[l]))(0) ELSE BEGIN ; IF lumclasso8(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclasso8(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclasso8(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclasso8(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclasso8(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclasso8(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclasso8(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclasso8(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_o8*rsun2cm;8.5*rsun2cm ; IF lumclasso8(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclasso8(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindo8_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indo8_2[l]))(0) ELSE BEGIN ; IF lumclasso8(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasso8(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasso8(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasso8(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasso8(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasso8(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasso8(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasso8(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_o8*rsun2cm;35.1e3 ; IF lumclasso8(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasso8(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindo8_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indo8_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindo8_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indo8_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmo8 ENDIF IF countindo9 NE 0 THEN BEGIN FOR l=0L,nlmo9-1 DO BEGIN comp_pop(k+i+l).popid = 'O9'+lumclasso9(l)+'_stellar_population' IF countindo9_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indo9_1[l]))(0) ELSE BEGIN ; IF lumclasso9(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclasso9(l) EQ 'IA' THEN comp_pop(k+i+l).radius = IF lumclasso9(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = 20.0*rsun2cm ;taken from Alnitak Aa's data. Even though its spec type is O9.5 (not O9) ; IF lumclasso9(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclasso9(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclasso9(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclasso9(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclasso9(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_o9*rsun2cm;7.51*rsun2cm ; IF lumclasso9(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclasso9(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindo9_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indo9_2[l]))(0) ELSE BEGIN ; IF lumclasso9(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasso9(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = IF lumclasso9(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = 29.5e3 ;taken from Alnitak Aa's data. Even though its spec type is O9.5 (not O9) ; IF lumclasso9(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasso9(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasso9(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasso9(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasso9(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_o9*rsun2cm;33.3e3 ; IF lumclasso9(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasso9(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindo9_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indo9_3[l]))(0) ELSE comp_pop(k+i+l).distance = 1.0E+00 IF countindo9_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indo9_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmo9 ENDIF ENDFOR counto=0L i+=k-1 ENDIF ;============================================================================================== IF countb NE 0 THEN BEGIN ; Filling the structure of the B STELLAR POPULATION(S) FOR k=0L,countb-1 DO BEGIN IF countindb0 NE 0 THEN BEGIN FOR l=0L,nlmb0-1 DO BEGIN comp_pop(k+i+l).popid = 'B0'+lumclassb0(l)+'_stellar_population' IF countindb0_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indb0_1[l]))(0) ELSE BEGIN ; IF lumclassb0(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassb0(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassb0(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassb0(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassb0(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassb0(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassb0(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassb0(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_b0*rsun2cm;7.16*rsun2cm ; IF lumclassb0(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassb0(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindb0_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indb0_2[l]))(0) ELSE BEGIN ; IF lumclasb0(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasb0(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasb0(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasb0(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasb0(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasb0(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasb0(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasb0(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_b0;31.4e3 ; IF lumclasb0(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasb0(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindb0_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indb0_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindb0_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indb0_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmb0 ENDIF IF countindb1 NE 0 THEN BEGIN FOR l=0L,nlmb1-1 DO BEGIN comp_pop(k+i+l).popid = 'B1'+lumclassb1(l)+'_stellar_population' IF countindb1_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indb1_1[l]))(0) ELSE BEGIN ; IF lumclassb1(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassb1(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassb1(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassb1(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassb1(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassb1(l) EQ 'III' THEN comp_pop(k+i+l).radius = IF lumclassb1(l) EQ 'IV' THEN comp_pop(k+i+l).radius = 7.3*rsun2cm ;taken from Alnitak Ab's data. IF lumclassb1(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_b1*rsun2cm;5.71*rsun2cm ; IF lumclassb1(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassb1(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindb1_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indb1_2[l]))(0) ELSE BEGIN ; IF lumclassb1(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclassb1(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasb1(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclassb1(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclassb1(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclassb1(l) EQ 'III' THEN comp_pop(k+i+l).temperature = IF lumclassb1(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = 29.0e3 ;taken from Alnitak Ab's data. IF lumclassb1(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_b1*rsun2cm;26.0e3 ; IF lumclassb1(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclassb1(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindb1_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indb1_3[l]))(0) ELSE comp_pop(k+i+l).distance = 1.0E+00 IF countindb1_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indb1_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmb1 ENDIF IF countindb2 NE 0 THEN BEGIN FOR l=0L,nlmb2-1 DO BEGIN comp_pop(k+i+l).popid = 'B2'+lumclassb2(l)+'_stellar_population' IF countindb2_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indb2_1[l]))(0) ELSE BEGIN ; IF lumclassb2(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassb2(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassb2(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassb2(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassb2(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassb2(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassb2(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassb2(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_b2*rsun2cm;4.06*rsun2cm ; IF lumclassb2(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassb2(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindb2_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indb2_2[l]))(0) ELSE BEGIN ; IF lumclasb2(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasb2(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasb2(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasb2(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasb2(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasb2(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasb2(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasb2(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_b2;20.6e3 ; IF lumclasb2(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasb2(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindb2_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indb2_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindb2_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indb2_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmb2 ENDIF IF countindb3 NE 0 THEN BEGIN FOR l=0L,nlmb3-1 DO BEGIN comp_pop(k+i+l).popid = 'B3'+lumclassb3(l)+'_stellar_population' IF countindb3_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indb3_1[l]))(0) ELSE BEGIN ; IF lumclassb3(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassb3(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassb3(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassb3(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassb3(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassb3(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassb3(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassb3(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_b3*rsun2cm;3.61*rsun2cm ; IF lumclassb3(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassb3(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindb3_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indb3_2[l]))(0) ELSE BEGIN ; IF lumclasb3(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasb3(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasb3(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasb3(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasb3(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasb3(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasb3(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasb3(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_b3;17.0e3 ; IF lumclasb3(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasb3(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindb3_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indb3_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindb3_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indb3_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmb3 ENDIF IF countindb4 NE 0 THEN BEGIN FOR l=0L,nlmb4-1 DO BEGIN comp_pop(k+i+l).popid = 'B4'+lumclassb4(l)+'_stellar_population' IF countindb4_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indb4_1[l]))(0) ELSE BEGIN ; IF lumclassb4(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassb4(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassb4(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassb4(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassb4(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassb4(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassb4(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassb4(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_b4*rsun2cm;3.46*rsun2cm ; IF lumclassb4(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassb4(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindb4_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indb4_2[l]))(0) ELSE BEGIN ; IF lumclasb4(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasb4(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasb4(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasb4(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasb4(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasb4(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasb4(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_b4;16.4e3 ; IF lumclasb4(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasb4(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindb4_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indb4_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindb4_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indb4_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmb4 ENDIF IF countindb5 NE 0 THEN BEGIN FOR l=0L,nlmb5-1 DO BEGIN comp_pop(k+i+l).popid = 'B5'+lumclassb5(l)+'_stellar_population' IF countindb5_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indb5_1[l]))(0) ELSE BEGIN ; IF lumclassb5(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassb5(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassb5(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassb5(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassb5(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassb5(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassb5(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassb5(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_b5*rsun2cm;3.36*rsun2cm ; IF lumclassb5(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassb5(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindb5_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indb5_2[l]))(0) ELSE BEGIN ; IF lumclasb5(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasb5(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasb5(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasb5(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasb5(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasb5(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasb5(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasb5(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_b5;15.7e3 ; IF lumclasb5(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasb5(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindb5_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indb5_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindb5_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indb5_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmb5 ENDIF IF countindb6 NE 0 THEN BEGIN FOR l=0L,nlmb6-1 DO BEGIN comp_pop(k+i+l).popid = 'B6'+lumclassb6(l)+'_stellar_population' IF countindb6_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indb6_1[l]))(0) ELSE BEGIN ; IF lumclassb6(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassb6(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassb6(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassb6(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassb6(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassb6(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassb6(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassb6(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_b6*rsun2cm;3.27*rsun2cm ; IF lumclassb6(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassb6(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindb6_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indb6_2[l]))(0) ELSE BEGIN ; IF lumclasb6(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasb6(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasb6(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasb6(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasb6(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasb6(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasb6(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasb6(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_b6;14.5e3 ; IF lumclasb6(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasb6(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindb6_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indb6_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindb6_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indb6_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmb6 ENDIF IF countindb7 NE 0 THEN BEGIN FOR l=0L,nlmb7-1 DO BEGIN comp_pop(k+i+l).popid = 'B7'+lumclassb7(l)+'_stellar_population' IF countindb7_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indb7_1[l]))(0) ELSE BEGIN ; IF lumclassb7(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassb7(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassb7(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassb7(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassb7(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassb7(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassb7(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassb7(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_b7*rsun2cm;2.94*rsun2cm ; IF lumclassb7(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassb7(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindb7_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indb7_2[l]))(0) ELSE BEGIN ; IF lumclasb7(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasb7(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasb7(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasb7(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasb7(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasb7(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasb7(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasb7(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_b7;14.0e3 ; IF lumclasb7(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasb7(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindb7_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indb7_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindb7_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indb7_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmb7 ENDIF IF countindb8 NE 0 THEN BEGIN FOR l=0L,nlmb7-1 DO BEGIN comp_pop(k+i+l).popid = 'B8'+lumclassb8(l)+'_stellar_population' IF countindb8_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indb8_1[l]))(0) ELSE BEGIN ; IF lumclassb8(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassb9(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassb8(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassb8(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassb8(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassb8(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassb8(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassb8(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_b8*rsun2cm;2.86*rsun2cm ; IF lumclassb8(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassb8(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindb8_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indb8_2[l]))(0) ELSE BEGIN ; IF lumclasb8(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasb8(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasb8(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasb8(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasb8(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasb8(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasb8(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasb8(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_b8;12.3e3 ; IF lumclasb8(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasb8(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindb8_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indb8_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindb8_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indb8_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmb8 ENDIF IF countindb9 NE 0 THEN BEGIN FOR l=0L,nlmb9-1 DO BEGIN comp_pop(k+i+l).popid = 'B9'+lumclassb9(l)+'_stellar_population' IF countindb9_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indb9_1[l]))(0) ELSE BEGIN ; IF lumclassb9(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassb9(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassb9(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassb9(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassb9(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassb9(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassb9(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassb9(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_b9*rsun2cm;2.49*rsun2cm ; IF lumclassb9(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassb9(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindb9_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indb9_2[l]))(0) ELSE BEGIN ; IF lumclasb9(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasb9(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasb9(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasb9(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasb9(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasb9(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasb9(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasb9(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_b9;10.7e3 ; IF lumclasb9(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasb9(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindb9_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indb9_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindb9_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indb9_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmb9 ENDIF ENDFOR countb=0. i+=k-1 ENDIF ;============================================================================================== IF counta NE 0 THEN BEGIN ; Filling the structure of the A STELLAR POPULATION(S) FOR k=0L+i,counta-1+i DO BEGIN IF countinda0 NE 0 THEN BEGIN FOR l=0L,nlma0-1 DO BEGIN comp_pop(k+i+l).popid = 'A0'+lumclassa0(l)+'_stellar_population' IF countinda0_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(inda0_1[l]))(0) ELSE BEGIN ; IF lumclassa0(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassa0(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassa0(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassa0(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassa0(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassa0(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassa0(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassa0(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_a0*rsun2cm;2.193*rsun2cm ; IF lumclassa0(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassa0(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countinda0_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(inda0_2[l]))(0) ELSE BEGIN ; IF lumclasa0(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasa0(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasa0(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasa0(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasa0(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasa0(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasa0(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasa0(l) EQ 'V' THEN comp_pop(k+i+l).temperature =T_eff_MS_a0; 9.7e3 ; IF lumclasa0(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasa0(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countinda0_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(inda0_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countinda0_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(inda0_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlma0 ENDIF IF countinda1 NE 0 THEN BEGIN FOR l=0L,nlma1-1 DO BEGIN comp_pop(k+i+l).popid = 'A1'+lumclassa1(l)+'_stellar_population' IF countinda1_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(inda1_1[l]))(0) ELSE BEGIN ; IF lumclassa1(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassa1(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassa1(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassa1(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassa1(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassa1(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassa1(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassa1(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_a1*rsun2cm;2.136*rsun2cm ; IF lumclassa1(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassa1(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countinda1_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(inda1_2[l]))(0) ELSE BEGIN ; IF lumclasa1(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasa1(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasa1(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasa1(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasa1(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasa1(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasa1(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasa1(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_a1;9.3e3 ; IF lumclasa1(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasa1(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countinda1_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(inda1_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countinda1_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(inda1_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ; K+=1 ENDFOR k+=nlma1 ENDIF IF countinda2 NE 0 THEN BEGIN FOR l=0L,nlma2-1 DO BEGIN comp_pop(k+i+l).popid = 'A2'+lumclassa2(l)+'_stellar_population' IF countinda2_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(inda2_1[l]))(0) ELSE BEGIN ; IF lumclassa2(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassa2(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassa2(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassa2(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassa2(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassa2(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassa2(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassa2(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_a2*rsun2cm;2.117*rsun2cm ; IF lumclassa2(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassa2(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countinda2_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(inda2_2[l]))(0) ELSE BEGIN ; IF lumclasa2(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasa2(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasa2(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasa2(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasa2(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasa2(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasa2(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasa2(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_a2;8.8e3 ; IF lumclasa2(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasa2(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countinda2_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(inda2_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countinda2_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(inda2_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlma2 ENDIF IF countinda3 NE 0 THEN BEGIN FOR l=0L,nlma3-1 DO BEGIN comp_pop(k+i+l).popid = 'A3'+lumclassa3(l)+'_stellar_population' IF countinda3_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(inda3_1[l]))(0) ELSE BEGIN ; IF lumclassa3(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassa3(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassa3(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassa3(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassa3(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassa3(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassa3(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassa3(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_a3*rsun2cm;1.861*rsun2cm ; IF lumclassa3(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassa3(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countinda3_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(inda3_2[l]))(0) ELSE BEGIN ; IF lumclasa3(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasa3(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasa3(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasa3(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasa3(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasa3(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasa3(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasa3(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_a3;8.6e3 ; IF lumclasa3(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasa3(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countinda3_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(inda3_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countinda3_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(inda3_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlma3 ENDIF IF countinda4 NE 0 THEN BEGIN FOR l=0L,nlma4-1 DO BEGIN comp_pop(k+i+l).popid = 'A4'+lumclassa4(l)+'_stellar_population' IF countinda4_1 NE 0 THEN comp_pop(k+i+l).radius = (val(inda4_1))(0) ELSE BEGIN ; IF lumclassa4(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassa4(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassa4(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassa4(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassa4(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassa4(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassa4(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassa4(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_a4*rsun2cm;1.794*rsun2cm ; IF lumclassa4(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassa4(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countinda4_2 NE 0 THEN comp_pop(k+i+l).temperature = (val(inda4_2))(0) ELSE BEGIN ; IF lumclasa4(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasa4(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasa4(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasa4(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasa4(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasa4(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasa4(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasa4(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_a4;8.25e3 ; IF lumclasa4(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasa4(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countinda4_3 NE 0 THEN comp_pop(k+i+l).distance = (val(inda4_3))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countinda4_4 NE 0 THEN comp_pop(k+i+l).nstars = (val(inda4_4))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlma4 ENDIF IF countinda5 NE 0 THEN BEGIN FOR l=0L,nlma5-1 DO BEGIN comp_pop(k+i+l).popid = 'A5'+lumclassa5(l)+'_stellar_population' IF countinda5_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(inda5_1[l]))(0) ELSE BEGIN ; IF lumclassa5(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassa5(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassa5(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassa5(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassa5(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassa5(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassa5(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassa5(l) EQ 'V' THEN omp_pop(k+i+l).radius = R_star_MS_a5*rsun2cm;1.785*rsun2cm ; IF lumclassa5(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassa5(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countinda5_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(inda5_2[l]))(0) ELSE BEGIN ; IF lumclasa5(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasa5(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasa5(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasa5(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasa5(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasa5(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasa5(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasa5(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_a5;8.1e3 ; IF lumclasa5(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasa5(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countinda5_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(inda5_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countinda5_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(inda5_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlma5 ENDIF IF countinda6 NE 0 THEN BEGIN FOR l=0L,nlma6-1 DO BEGIN comp_pop(k+i+l).popid = 'A6'+lumclassa6(l)+'_stellar_population' IF countinda6_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(inda6_1[l]))(0) ELSE BEGIN ; IF lumclassa6(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassa6(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassa6(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassa6(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassa6(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassa6(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassa6(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassa6(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_a6*rsun2cm;1.775*rsun2cm ; IF lumclassa6(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassa6(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countinda6_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(inda6_2[l]))(0) ELSE BEGIN ; IF lumclasa6(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasa6(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasa6(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasa6(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasa6(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasa6(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasa6(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasa6(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_a6;7.91e3 ; IF lumclasa6(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasa6(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countinda6_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(inda6_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countinda6_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(inda6_5[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlma6 ENDIF IF countinda7 NE 0 THEN BEGIN FOR l=0L,nlma7-1 DO BEGIN comp_pop(k+i+l).popid = 'A7'+lumclassa7(l)+'_stellar_population' IF countinda7_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(inda7_1[l]))(0) ELSE BEGIN ; IF lumclassa7(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassa7(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassa7(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassa7(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassa7(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassa7(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassa7(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassa7(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_a7*rsun2cm;1.75*rsun2cm ; IF lumclassa7(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassa7(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countinda7_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(inda7_2[l]))(0) ELSE BEGIN ; IF lumclasa7(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasa7(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasa7(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasa7(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasa7(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasa7(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasa7(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_a7;7.76e3 ; IF lumclasa7(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasa7(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countinda7_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(inda7_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countinda7_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(inda7_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlma7 ENDIF IF countinda8 NE 0 THEN BEGIN FOR l=0L,nlma8-1 DO BEGIN comp_pop(k+i+l).popid = 'A8'+lumclassa8(l)+'_stellar_population' IF countinda8_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(inda8_1[l]))(0) ELSE BEGIN ; IF lumclassa8(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassa8(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassa8(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassa8(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassa8(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassa8(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassa8(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassa8(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_a8*rsun2cm;1.747*rsun2cm ; IF lumclassa8(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassa8(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countinda8_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(inda8_2[l]))(0) ELSE BEGIN ; IF lumclasa8(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasa8(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasa8(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasa8(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasa8(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasa8(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasa8(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasa8(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_a8;7.590e3 ; IF lumclasa8(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasa8(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countinda8_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(inda8_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countinda8_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(inda8_4[l]))(0) ELSE comp_pop(k+i+l).nstars =1; stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlma8 ENDIF IF countinda9 NE 0 THEN BEGIN FOR l=0L,nlma9-1 DO BEGIN comp_pop(k+i+l).popid = 'A9'+lumclassa9(l)+'_stellar_population' IF countinda9_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(inda9_1[l]))(0) ELSE BEGIN ; IF lumclassa9(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassa9(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassa9(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassa9(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassa9(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassa9(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassa9(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassa9(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_a9*rsun2cm;1.747*rsun2cm ; IF lumclassa9(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassa9(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countinda9_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(inda9_2[l]))(0) ELSE begin ; IF lumclasa9(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclasa9(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclasa9(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclasa9(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclasa9(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclasa9(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclasa9(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclasa9(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_a9;7.4e3 ; IF lumclasa9(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclasa9(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countinda9_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(inda9_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countinda9_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(inda9_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlma9 ENDIF ENDFOR counta=0. i+=k-1 ENDIF ;============================================================================================== IF countf NE 0 THEN BEGIN ; Filling the structure of the F STELLAR POPULATION(S) FOR k=0L+i,countf-1+i DO BEGIN IF countindf0 NE 0 THEN BEGIN FOR l=0L,nlmaf0-1 DO BEGIN comp_pop(k+i+l).popid = 'F0'+lumclassf0(l)+'_stellar_population' IF countindf0_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indf0_1[l]))(0) ELSE BEGIN ; IF lumclassf0(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassf0(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassf0(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassf0(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassf0(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassf0(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassf0(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassf0(l) EQ 'V' THEN omp_pop(k+i+l).radius = R_star_MS_f0*rsun2cm;1.728*rsun2cm ; IF lumclassf0(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassf0(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindf0_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indf0_2[l]))(0) ELSE BEGIN ; IF lumclassf0(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclassf0(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclassf0(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclassf0(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclassf0(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclassf0(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclassf0(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclassf0(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_f0;7.22e3 ; IF lumclassf0(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclassf0(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindf0_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indf0_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindf0_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indf0_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmf0 ENDIF IF countindf1 NE 0 THEN BEGIN FOR l=0L,nlmaf1-1 DO BEGIN comp_pop(k+i+l).popid = 'F1'+lumclassf1(l)+'_stellar_population' IF countindf1_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indf1_1[l]))(0) ELSE BEGIN ; IF lumclassf1(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassf1(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassf1(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassf1(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassf1(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassf1(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassf1(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassf1(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_f1*rsun2cm;1.679*rsun2cm ; IF lumclassf1(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassf1(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindf1_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indf1_2[l]))(0) ELSE BEGIN ; IF lumclassf1(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclassf1(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclassf1(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclassf1(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclassf1(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclassf1(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclassf1(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclassf1(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_f1;7.02e3 ; IF lumclassf1(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclassf1(l) EQ 'VIII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindf1_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indf1_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindf1_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indf1_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmf1 ENDIF IF countindf2 NE 0 THEN BEGIN FOR l=0L,nlmaf2-1 DO BEGIN comp_pop(k+i+l).popid = 'F2'+lumclassf2(l)+'_stellar_population' IF countindf2_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indf2_1[l]))(0) ELSE BEGIN ; IF lumclassf2(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassf2(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassf2(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassf2(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassf2(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassf2(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassf2(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassf2(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_f2*rsun2cm;1.622*rsun2cm ; IF lumclassf2(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassf2(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindf2_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indf2_2[l]))(0) ELSE BEGIN ; IF lumclassf2(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclassf2(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclassf2(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclassf2(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclassf2(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclassf2(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclassf2(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclassf2(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_f2;6.820e3 ; IF lumclassf2(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclassf2(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindf2_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indf2_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindf2_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indf2_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmf2 ENDIF IF countindf3 NE 0 THEN BEGIN FOR l=0L,nlmaf3-1 DO BEGIN comp_pop(k+i+l).popid = 'F3'+lumclassf3(l)+'_stellar_population' IF countindf3_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indf3_1[l]))(0) ELSE BEGIN ; IF lumclassf3(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassf3(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassf3(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassf3(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassf3(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassf3(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassf3(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassf3(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_f3*rsun2cm;1.578*rsun2cm ; IF lumclassf3(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassf3(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindf3_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indf3_2[l]))(0) ELSE BEGIN ; IF lumclassf3(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclassf3(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclassf3(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclassf3(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclassf3(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclassf3(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclassf3(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclassf3(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_f3;6.75e3 ; IF lumclassf3(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclassf3(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindf3_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indf3_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindf3_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indf3_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmf3 ENDIF IF countindf4 NE 0 THEN BEGIN FOR l=0L,nlmaf4-1 DO BEGIN comp_pop(k+i+l).popid = 'F4'+lumclassf4(l)+'_stellar_population' IF countindf4_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indf4_1[l]))(0) ELSE BEGIN ; IF lumclassf4(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassf4(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassf4(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassf4(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassf4(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassf4(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassf4(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassf4(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_f4*rsun2cm;1.533*rsun2cm ; IF lumclassf4(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassf4(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindf4_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indf4_2[l]))(0) ELSE BEGIN ; IF lumclassf4(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclassf4(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclassf4(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclassf4(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclassf4(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclassf4(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclassf4(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclassf4(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_f4;6.67e3 ; IF lumclassf4(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclassf4(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindf4_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indf4_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindf4_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indf4_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmf4 ENDIF IF countindf5 NE 0 THEN BEGIN FOR l=0L,nlmaf5-1 DO BEGIN comp_pop(k+i+l).popid = 'F5'+lumclassf5(l)+'_stellar_population' IF countindf5_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indf5_1[l]))(0) ELSE BEGIN ; IF lumclassf5(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassf5(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassf5(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassf5(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassf5(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassf5(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassf5(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassf5(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_f5*rsun2cm;1.473*rsun2cm ; IF lumclassf5(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassf5(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindf5_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indf5_2[l]))(0) ELSE BEGIN ; IF lumclassf5(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclassf5(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclassf5(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclassf5(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclassf5(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclassf5(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclassf5(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclassf5(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_f5;6.55e3 ; IF lumclassf5(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclassf5(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindf5_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indf5_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindf5_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indf5_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmf5 ENDIF IF countindf6 NE 0 THEN BEGIN FOR l=0L,nlmaf6-1 DO BEGIN comp_pop(k+i+l).popid = 'F6'+lumclassf6(l)+'_stellar_population' IF countindf6_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indf6_1[l]))(0) ELSE BEGIN ; IF lumclassf6(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassf6(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassf6(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassf6(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassf6(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassf6(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassf6(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassf6(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_f6*rsun2cm;1.359*rsun2cm ; IF lumclassf6(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassf6(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindf6_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indf6_2[l]))(0) ELSE BEGIN ; IF lumclassf6(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclassf6(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclassf6(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclassf6(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclassf6(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclassf6(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclassf6(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclassf6(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_f6;6.35e3 ; IF lumclassf6(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclassf6(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindf6_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indf6_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindf6_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indf6_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmf6 ENDIF IF countindf7 NE 0 THEN BEGIN FOR l=0L,nlmaf7-1 DO BEGIN comp_pop(k+i+l).popid = 'F7'+lumclassf7(l)+'_stellar_population' IF countindf7_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indf7_1[l]))(0) ELSE BEGIN ; IF lumclassf7(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassf7(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassf7(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassf7(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassf7(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassf7(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassf7(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassf7(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_f7*rsun2cm;1.324*rsun2cm ; IF lumclassf7(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassf7(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindf7_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indf7_2[l]))(0) ELSE BEGIN ; IF lumclassf7(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclassf7(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclassf7(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclassf7(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclassf7(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclassf7(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclassf7(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclassf7(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_f7;6.28e3 ; IF lumclassf7(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclassf7(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindf7_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indf7_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindf7_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indf7_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmf7 ENDIF IF countindf8 NE 0 THEN BEGIN FOR l=0L,nlmaf8-1 DO BEGIN comp_pop(k+i+l).popid = 'F8'+lumclassf8(l)+'_stellar_population' IF countindf8_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indf8_1[l]))(0) ELSE BEGIN ; IF lumclassf8(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassf8(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassf8(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassf8(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassf8(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassf8(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassf8(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassf8(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_f8*rsun2cm;1.221*rsun2cm ; IF lumclassf8(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassf8(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindf8_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indf8_2[l]))(0) ELSE BEGIN ; IF lumclassf8(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclassf8(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclassf8(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclassf8(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclassf8(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclassf8(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclassf8(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclassf8(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_f8;6.18e3 ; IF lumclassf8(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclassf8(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindf8_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indf8_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindf8_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indf8_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmf8 ENDIF IF countindf9 NE 0 THEN BEGIN FOR l=0L,nlmaf9-1 DO BEGIN comp_pop(k+i+l).popid = 'F9'+lumclassf9(l)+'_stellar_population' IF countindf9_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indf9_1[l]))(0) ELSE BEGIN ; IF lumclassf9(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassf9(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassf9(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassf9(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassf9(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassf9(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassf9(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassf9(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_f9*rsun2cm;1.167*rsun2cm ; IF lumclassf9(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassf9(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindf9_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indf9_2[l]))(0) ELSE BEGIN ; IF lumclassf9(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclassf9(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclassf9(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclassf9(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclassf9(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclassf9(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclassf9(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclassf9(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_f9;6.05e3 ; IF lumclassf9(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclassf9(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindf9_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indf9_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindf9_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indf9_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmf9 ENDIF ENDFOR countf=0. i+=k-1 ENDIF ;============================================================================================== IF countg NE 0 THEN BEGIN ; Filling the structure of the G STELLAR POPULATION(S) FOR k=0L+i,countg-1+i DO BEGIN IF countindg0 NE 0 THEN BEGIN FOR l=0L,nlmg0-1 DO BEGIN comp_pop(k+i+l).popid = 'G0'+lumclassg0(l)+'_stellar_population' IF countindg0_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indg0_1[l]))(0) ELSE BEGIN ; IF lumclassg0(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassg0(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassg0(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassg0(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassg0(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassg0(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassg0(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassg0(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_g0*rsun2cm;1.1*rsun2cm ; IF lumclassg0(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassg0(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindg0_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indg0_2[l]))(0) ELSE BEGIN ; IF lumclassfg0(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg0(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg0(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg0(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg0(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg0(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg0(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclassfg0(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_g0;5.93e3 ; IF lumclassfg0(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg0(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindg0_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indg0_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindg0_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indg0_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmg0 ENDIF IF countindg1 NE 0 THEN BEGIN FOR l=0L,nlmg1-1 DO BEGIN comp_pop(k+i+l).popid = 'G1'+lumclassg1(l)+'_stellar_population' IF countindg1_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indg1_1[l]))(0) ELSE BEGIN ; IF lumclassg1(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassg1(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassg1(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassg1(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassg1(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassg1(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassg1(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassg1(l) EQ 'V' THEN omp_pop(k+i+l).radius = R_star_MS_g1*rsun2cm;1.06*rsun2cm ; IF lumclassg1(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassg1(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindg1_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indg1_2[l]))(0) ELSE begin ; IF lumclassfg1(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg1(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg1(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg1(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg1(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg1(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg1(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclassfg1(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_g1;5.86e3 ; IF lumclassfg1(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg1(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindg1_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indg1_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindg1_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indg1_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmg1 ENDIF IF countindg2 NE 0 THEN BEGIN FOR l=0L,nlmg2-1 DO BEGIN comp_pop(k+i+l).popid = 'G2'+lumclassg2(l)+'_stellar_population' IF countindg2_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indg2_1[l]))(0) ELSE BEGIN ; IF lumclassg2(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassg2(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassg2(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassg2(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassg2(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassg2(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassg2(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassg2(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_g2*rsun2cm;1.012*rsun2cm ; IF lumclassg2(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassg2(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindg2_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indg2_2[l]))(0) ELSE BEGIN ; IF lumclassfg2(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg2(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg2(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg2(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg2(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg2(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg2(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclassfg2(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_g2;5.77e3 ; IF lumclassfg2(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg2(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindg2_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indg2_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindg2_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indg2_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmg2 ENDIF IF countindg3 NE 0 THEN BEGIN FOR l=0L,nlmg3-1 DO BEGIN comp_pop(k+i+l).popid = 'G3'+lumclassg3(l)+'_stellar_population' IF countindg3_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indg3_1[l]))(0) ELSE BEGIN ; IF lumclassg3(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassg3(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassg3(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassg3(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassg3(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassg3(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassg3(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassg3(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_g3*rsun2cm;1.002*rsun2cm ; IF lumclassg3(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassg3(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindg3_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indg3_2[l]))(0) ELSE BEGIN ; IF lumclassfg3(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg3(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg3(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg3(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg3(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg3(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg3(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclassfg3(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_g3;5.72e3 ; IF lumclassfg3(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg3(l) EQ 'VII' THEN comp_pop(k+i+l).temperature ENDELSE IF countindg3_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indg3_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindg3_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indg3_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmg3 ENDIF IF countindg4 NE 0 THEN BEGIN FOR l=0L,nlmg4-1 DO BEGIN comp_pop(k+i+l).popid = 'G4'+lumclassg4(l)+'_stellar_population' IF countindg4_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indg4_1[l]))(0) ELSE BEGIN ; IF lumclassg4(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassg4(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassg4(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassg4(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassg4(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassg4(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassg4(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassg4(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_g4*rsun2cm;0.991*rsun2cm ; IF lumclassg4(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassg4(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindg4_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indg4_2[l]))(0) ELSE BEGIN ; IF lumclassfg4(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg4(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg4(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg4(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg4(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg4(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg4(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclassfg4(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_g4;5.68e3 ; IF lumclassfg4(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg4(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindg4_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indg4_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindg4_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indg4_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmg4 ENDIF IF countindg5 NE 0 THEN BEGIN FOR l=0L,nlmg5-1 DO BEGIN comp_pop(k+i+l).popid = 'G5'+lumclassg5(l)+'_stellar_population' IF countindg5_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indg5_1[l]))(0) ELSE BEGIN ; IF lumclassg5(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassg5(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassg5(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassg5(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassg5(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassg5(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassg5(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassg5(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_g5*rsun2cm;0.977*rsun2cm ; IF lumclassg5(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassg5(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindg5_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indg5_2[l]))(0) ELSE BEGIn ; IF lumclassfg5(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg5(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg5(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg5(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg5(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg5(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg5(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclassfg5(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_g5;5.66e3 ; IF lumclassfg5(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg5(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindg5_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indg5_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindg5_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indg5_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmg5 ENDIF IF countindg6 NE 0 THEN BEGIN FOR l=0L,nlmg6-1 DO BEGIN comp_pop(k+i+l).popid = 'G6'+lumclassg6(l)+'_stellar_population' IF countindg6_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indg6_1[l]))(0) ELSE BEGIN ; IF lumclassg6(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassg6(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassg6(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassg6(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassg6(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassg6(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassg6(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassg6(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_g6*rsun2cm;0.949*rsun2cm ; IF lumclassg6(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassg6(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindg6_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indg6_2[l]))(0) ELSE BEGIN ; IF lumclassfg6(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg6(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg6(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg6(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg6(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg6(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg6(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclassfg6(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_g6;5.6e3 ; IF lumclassfg6(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg6(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindg6_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indg6_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindg6_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indg6_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmg6 ENDIF IF countindg7 NE 0 THEN BEGIN FOR l=0L,nlmg7-1 DO BEGIN comp_pop(k+i+l).popid = 'G7'+lumclassg7(l)+'_stellar_population' IF countindg7_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indg7_1[l]))(0) ELSE BEGIN ; IF lumclassg7(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassg7(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassg7(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassg7(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassg7(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassg7(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassg7(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassg7(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_g7*rsun2cm;0.927*rsun2cm ; IF lumclassg7(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassg7(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindg7_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indg7_2[l]))(0) ELSE BEGIN ; IF lumclassfg7(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg7(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg7(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg7(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg7(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg7(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg7(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclassfg7(l) EQ 'V' THEN ccomp_pop(k+i+l).temperature = T_eff_MS_g7;5.55e3 ; IF lumclassfg7(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg7(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindg7_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indg7_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindg7_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indg7_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmg7 ENDIF IF countindg8 NE 0 THEN BEGIN FOR l=0L,nlmg8-1 DO BEGIN comp_pop(k+i+l).popid = 'G8'+lumclassg8(l)+'_stellar_population' IF countindg8_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indg8_1[l]))(0) ELSE BEGIN ; IF lumclassg8(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassg8(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassg8(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassg8(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassg8(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassg8(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassg8(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassg8(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_g8*rsun2cm;0.914*rsun2cm ; IF lumclassg8(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassg8(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindg8_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indg8_2[l]))(0) ELSE BEGIN ; IF lumclassfg8(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg8(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg8(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg8(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg8(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg8(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg8(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclassfg8(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_g8;5.48e3 ; IF lumclassfg8(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg8(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindg8_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indg8_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindg8_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indg8_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmg8 ENDIF IF countindg9 NE 0 THEN BEGIN FOR l=0L,nlmg9-1 DO BEGIN comp_pop(k+i+l).popid = 'G9'+lumclassg9(l)+'_stellar_population' IF countindg9_1[l] NE 0 THEN comp_pop(k+i+l).radius = (val(indg9_1[l]))(0) ELSE BEGIN ; IF lumclassg9(l) EQ 'IA+' THEN comp_pop(k+i+l).radius = ; IF lumclassg9(l) EQ 'IA' THEN comp_pop(k+i+l).radius = ; IF lumclassg9(l) EQ 'IAB' THEN comp_pop(k+i+l).radius = ; IF lumclassg9(l) EQ 'IB' THEN comp_pop(k+i+l).radius = ; IF lumclassg9(l) EQ 'II' THEN comp_pop(k+i+l).radius = ; IF lumclassg9(l) EQ 'III' THEN comp_pop(k+i+l).radius = ; IF lumclassg9(l) EQ 'IV' THEN comp_pop(k+i+l).radius = IF lumclassg9(l) EQ 'V' THEN comp_pop(k+i+l).radius = R_star_MS_g9*rsun2cm;0.853*rsun2cm ; IF lumclassg9(l) EQ 'VI' THEN comp_pop(k+i+l).radius = ; IF lumclassg9(l) EQ 'VII' THEN comp_pop(k+i+l).radius = ENDELSE IF countindg9_2[l] NE 0 THEN comp_pop(k+i+l).temperature = (val(indg9_2[l]))(0) ELSE BEGIN ; IF lumclassfg9(l) EQ 'IA+' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg9(l) EQ 'IA' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg9(l) EQ 'IAB' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg9(l) EQ 'IB' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg9(l) EQ 'II' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg9(l) EQ 'III' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg9(l) EQ 'IV' THEN comp_pop(k+i+l).temperature = IF lumclassfg9(l) EQ 'V' THEN comp_pop(k+i+l).temperature = T_eff_MS_g9;5.38e3 ; IF lumclassfg9(l) EQ 'VI' THEN comp_pop(k+i+l).temperature = ; IF lumclassfg9(l) EQ 'VII' THEN comp_pop(k+i+l).temperature = ENDELSE IF countindg9_3[l] NE 0 THEN comp_pop(k+i+l).distance = (val(indg9_3[l]))(0) ELSE comp_pop(k+i+l).distance = 10.0 IF countindg9_4[l] NE 0 THEN comp_pop(k+i+l).nstars = (val(indg9_4[l]))(0) ELSE comp_pop(k+i+l).nstars = 1;stellar_density*(4.*!pi/3)*comp_pop(k+i+l).distance ;K+=1 ENDFOR k+=nlmg9 ENDIF ENDFOR countg=0. i+=k-1 ENDIF ;============================================================================================== ENDFOR st=((*!dustem_params).isrf) c2a = 3e18 ;speed of light in ansgtroms/s (because of the Astron's PLANCK function) pc2cm = 3.086e18 ;cm (cgs) wave_angstrom = st.lambisrf*1.e4 ;mic to Angstrom (Astron Planck's function uses wavelengths in Angstroms) stellar_component=fltarr(n_elements(st)) ; array of zeros to contain the new ISRF values. FOR i=0L,n_elements(comp_pop.popid)-1 DO BEGIN ; Looping over all the stellar populations omega = ((comp_pop.radius)[i]/((comp_pop.distance)[i]*pc2cm))^2 ; Dilution factor of a stellar population Inu = planck(wave_angstrom,(comp_pop.temperature)[i])*(wave_angstrom)^2/c2a ; ergs/cm2/s/Hz/sr ; ;THIS PART IS ONLY FOR ORION. ; ;IT SHOULD NOT BE PROVIDED IN THE RELEASE ; if isa(!dustem_current) then begin ; ; idff = where(tag_names(*!dustem_plugin) EQ 'FREEFREE',ct_ff) ; ; if ct_ff ne 0 then begin ; ; ;look for freefree temperature in the pd vector (sysvar) ; ;we need an index ; ; id_fftmp = where(strupcase(strmid((*(*!dustem_fit).param_descs),14)) EQ 'FREEFREE_1',ct_fftmp) ; if ct_fftmp ne 0 then Te = ((*(*!dustem_fit).current_param_values)[id_fftmp])[0]*(((*(*!dustem_fit).param_init_values))[id_fftmp])[0] ; ; test if Te is correct ; ; Beta2 = 2.e-10*(Te)^(-3/4) ;cm^3/s-1 ; Rs_p3 = 1.1355d+57 ;Rs^3 for orion nebula ; nH_stromgren = sqrt(Rs_p3*4*!pi*Beta2/(3*(!const.sigma)*1.0d-04*((comp_pop.temperature)[i])^4*!pi)) ; ; ; Inu *= exp(((*!dustem_current).ext.ext_tot)*(-nh_stromgren)*4/3*Rs_p3^(1/3)/1.0d21) ; ; ; endif ; ; ENDIF If !dustem_dim and isa(!dustem_current) then Inu *= exp(((*!dustem_current).ext.ext_tot)*(-(*!dustem_HCD))/1.0d21) stellar_component=stellar_component+(comp_pop.nstars)[i]*omega*Inu ;no *4 because of Lambert's cosine law. ENDFOR out = stellar_component ;============This block creates a composite (9V) stellar population structure for when the function isn't used as a plugin=================== ENDIF ELSE BEGIN popnumber = 5. comp_pop = replicate(one_pop,popnumber) stellar_density = 0.14 ;stars/pc^3 ;###BY DEFAULT ALL STELLAR POPULATIONS ARE ON THE MS### ;BECAUSE LESS MASSIVE (COLDER) STARS ARE MORE ABUNDANT, THEY WILL BE CHOSEN TO BUILD THE COMPOSITE STELLAR STRUCTURE ;THIS MEANS THAT THE CHOSEN SPECTRAL TYPE IS 9V (MAIN SEQUENCE) ;CHOSEN STELLAR POPULATIONS ;O9V ;B9V ;A9V ;F9V ;G9V comp_pop(0).popid = 'O9V_stellar_population' comp_pop(0).radius = 7.51*rsun2cm comp_pop(0).temperature = 33.3e3 comp_pop(0).distance = 10.0 comp_pop(0).nstars = stellar_density*(4.*!pi/3)*comp_pop(0).distance comp_pop(1).popid = 'B9V_stellar_population' comp_pop(1).radius = 2.49*rsun2cm comp_pop(1).temperature = 10.7e3 comp_pop(1).distance = 10.0 comp_pop(1).nstars = stellar_density*(4.*!pi/3)*comp_pop(1).distance comp_pop(2).popid = 'A9V_stellar_population' comp_pop(2).radius = 1.747*rsun2cm comp_pop(2).temperature = 7.4e3 comp_pop(2).distance = 10.0 comp_pop(2).nstars = stellar_density*(4.*!pi/3)*comp_pop(2).distance comp_pop(3).popid = 'F9V_stellar_population' comp_pop(3).radius = 1.167*rsun2cm comp_pop(3).temperature = 6.05e3 comp_pop(3).distance = 10.0 comp_pop(3).nstars = stellar_density*(4.*!pi/3)*comp_pop(3).distance comp_pop(4).popid = 'G9V_stellar_population' comp_pop(4).radius = 0.853*rsun2cm comp_pop(4).temperature = 5.38e3 comp_pop(4).distance = 10.0 comp_pop(4).nstars = stellar_density*(4.*!pi/3)*comp_pop(4).distance st=((*!dustem_params).isrf) c2a = 3e18 ;speed of light in ansgtroms/s (because of the Astron's PLANCK function) pc2cm = 3.086e18 ;cm (cgs) wave_angstrom = st.lambisrf*1.e4 ;mic to Angstrom (Astron Planck's function uses wavelengths in Angstroms) stellar_component=fltarr(n_elements(st)) ; array of zeros to contain the new ISRF values. FOR i=0L,n_elements(comp_pop.popid)-1 DO BEGIN ; Looping over all the stellar populations omega = ((comp_pop.radius)[i]/((comp_pop.distance)[i]*pc2cm))^2 ; Dilution factor of a stellar population Inu = 4*planck(wave_angstrom,(comp_pop.temperature)[i])*(wave_angstrom)^2/c2a ; ergs/cm2/s/Hz/sr ;*4 and not using Lambert's cosine law stellar_component=stellar_component+(comp_pop.nstars)[i]*omega*Inu ENDFOR out = stellar_component ENDELSE the_scope: scope='STELLAR_POPULATION' the_paramtag: paramtag=['R_star (R_sol)','T_BB (K)','D (pc)','N_stars'] the_end: return, out end