Blame view

src/idl/inspect_emiles_ssp.pro 661 Bytes
2d9acde0   Jean-Philippe Bernard   First commit
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
PRO inspect_emiles_ssp

loadct,13
;dir='/Users/jpb/Soft_Libraries/dustem-wrapper_idl/Data/SSPs/ssp_emiles-ch/'
dir='/Users/jpb/Soft_Libraries/dustem-wrapper_idl/Data/SSPs/eMILES_Ch-young-MetalPoorRemoved/'
files=file_search(dir+'*.fits')
Nfiles=n_elements(files)
yr=[1.e-5,1.e-2]
xtit='Wavelenght [mic]'
ytit='Intensity []'
colours=long(range_gen(Nfiles,[0,255]))
FOR i=0L,Nfiles-1 DO BEGIN
  message,'reading '+files[i],/continue
  d=readfits(files[i],h)
  make_axis,h,wav
  wav=wav/1.e10*1.e6
  IF i EQ 0 THEN BEGIN
    cgplot,wav,d,/xlog,/ylog,yrange=yr,xtit=xtit,ytit=ytit
  ENDIF ELSE BEGIN
    cgoplot,wav,d,color=colours[i]
    stop
  ENDELSE
ENDFOR

END