;This program plots quicklook data from the raw FUSE science pipline. ;It was started on 8-13-99 by SRM @h1abs @convolute @newjoin @h2label @h2data @atlas @rescol @fm @make_kernel gd='n' read,'Get data?',gd ;This is the current ism (ie h2) offset. ismoff=.15 ;A few constants that may or may not be used c=2.9979e5 sig=5.66956e-5;erg/cm^2/sec/k^4 col=255 cs=1.25 ang=string("305b) ;Angstrom symbol in idl for plots fst='ergs cm!e-2!n s!e-1!n '+ang+'!e-1!n' ;a flux string for plots ;The wavelength grid dlin=.01 linearize=900.+dlin*findgen(290/dlin) lminmax=minmax(linearize) ;The stellar model distpc=1./4.3*1000 ;hipparcos is (4.3+-.81m'') t=20000 g=4 atm=atlas(t,g) l=transpose(atm(0,*)*10) ; To convert cloudy atlas to flux in erg/(cm^2 s Ang) f=transpose(atm(1,*)*((c*1e5)*1e8/(atm(0,*)*10)^2)*!pi*4) inrange=where(l lt 3200) l=l(inrange) f=f(inrange) rv=3.2 mbmv=.13 ibmv=-.2 ;infered for B3V from fitz&massa ebv=mbmv-ibmv ;ccmpar=[4.603,0.991,-0.143,.7,3.536,.5],x0,gam,c1,c2,c3,c4 ccmpar=[4.543,1.219,-0.966,0.321,6.64,.533] ;from EBB's fm_fit modifyccm=fm(l,ccmpar(0),ccmpar(1),ccmpar(2),ccmpar(3),ccmpar(4),ccmpar(5)) dust=modifyccm ;alam=(dust/rv+1.)*rv*ebv alam=dust*ebv tdust=10.^(-0.4*alam) ftdust=f*tdust gg=10d0^g ;g is log GM/R^2 cgs bigg=6.673d-8 ;G in cgs masssol=1.989d33 ;solar mass in grams rsol=6.96e10 ;solar radius in cm radius=sqrt(bigg*masssol/gg) print,radius/rsol pc=3.085e18 ;parsec in cm flam=ftdust*(radius/(distpc*pc))^2 fl=interpol(flam,l,linearize) ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;The h2 optical depth models stored in the files called ;tauh2n21b*j0-15v0.dat contain negative optical depths as ;a function of wavelength (.01 Ang grid) computed for an ;maximum column density of 10^21 (cm^-2). ; ;To get another column you must scale downward ;by multiplying by the ratio of ncolh2/1e21. ; ;The file code for tauh2n21b1j0-15v0.dat ;tauh2 >>>> (negative) optical depth for h2 ;n21 >>>> computed for column density of 1e21 ;b1 >>>> the doppler parameter in units of km/s ;j0-15 >>>> rotational states from 0-15 ;v0 >>>> vibrational state is 0 ; ;The data are stored as 17 double precision vectors 59000 long. ;The first vector is the wavelength grid. ;The 1-16 vectors contain the individual rotational state templates. ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; bb=5 ;the b for the templates bbs=strtrim(string(fix(bb)),2) openr,1,'/home/stephan/idlpro/h2/h2abspvppn/tauh2n21b'+bbs+'j0-15v0.dat' length=59000l ;the l means a long interger (it has no sign) v0=assoc(1,dblarr(length)) lj0=v0(0) ncol=[100,100,1.,1,.05,.001,fltarr(10)]*1e18 ; a guess tauh2=lj0*0 for j=0,15 do tauh2=tauh2+v0(j+1)*ncol(j)/1e21 h2a=exp(tauh2) close,1 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;HI model is generated by the function h1abs.pro ;which is used more like a procedure. ;h1 is the returned absorption vector ;lall is the returned wavelength vector ;hcol is the input column (log N (cm^-2)) ;hvkms is the input doppler shift (km/s) ;hb is the input doppler parameter (km/s) ;xphi are the individual cross-sections ;lama are the individual wavelength vectors they are not shifted ;tall is the combined optical depth as a function of wavelength on the lall grid ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; yn='' read,'redo h1',yn if yn eq 'y' then begin hcol=[22.] hvkms=[0] hb=[30] h1=h1abs(lall,hcol,hvkms,hb,xphi,lama,tall) ;h1c=interpol(h1,lall,linearize) endif tauh1=interpol(tall,lall,linearize) ;combine h1 and h2 and allow for velocity offset ;tauh1h2=tauh1+tauh2 tauh1h2=tauh1+interpol(tauh2(0:n_elements(linearize)-1),(1.+ismoff/1100.)*linearize,linearize) ;convolve to resolution of instrument and attenuate stellar model ;minmod=6e-13*h2a respower=12000. fwhm=1000./respower ;in angstroms kk=make_kernel(grid=linearize,fwhm=fwhm,/gauss) h=convol(exp(tauh1h2),kk) minmod=h*fl ;The channels averaged over all orbits as given by B-G if gd eq 'y' then begin yr=[0,1e-10] dr='/home/stephan/FUSE/CAMMdata/' files=findfile(dr+'*0606*') nraw=n_elements(files) print,[transpose(string(indgen(nraw))),transpose(files)] segs=fltarr(2^14,nraw,6) for fn=0,n_elements(files)-1 do begin file=files(fn) im=mrdfits(file,1) print,'Reading file',file segs(0,fn,0)=im.wave segs(0,fn,1)=im.flux segs(0,fn,2)=im.error segs(0,fn,3)=im.quality segs(0,fn,4)=im.counts segs(0,fn,5)=im.cntserr endfor window,0,xsize=500,ysize=900 !p.multi=[0,1,8] for fn=0,n_elements(files)-1 do begin lam=segs(*,fn,0) xr=minmax(lam) filetit=strmid(files(fn),40,17)+' '+string(fn,'(I3)') plot,lam,smooth(segs(*,fn,1),6),psym=10,ticklen=-.025,charsize=cs,title=filetit,xtitle='Wavelength ('+ang+')',ytit='Flux ('+fst+')',color=col,ys=1,xs=1,yr=yr,xr=xr,ymargin=[0,0] oplot,l,flam endfor ;This is section prepares the data for joining ;Read in the data and rebin down by factor of 2 ;and define wavelength, flux and count vectors as noted below for i=0,7 do print,i,' ',strmid(files(i),strpos(files(i),'000')+3,5) ;a rebin down parameter rpow=14.-1. ;need to rescale the sic channels that didn't have all the star is the slit ;all the time rscls1=1. rscls2=1. ;select data -- cl** are the true counts for weighting ll1a=rebin(segs(*,0,0),2^rpow) ;wavelength fl1a=rebin(segs(*,0,1),2^rpow) ;flux nl1a=(rebin(segs(*,0,2),2^rpow))/sqrt(2.^(14.-rpow));fluxnoisedropsbyroot2perreb el1a=segs(*,0,5) ;sqrt of total counts? cl1a=rebin(el1a^2,2^rpow)*2^(14.-rpow) ;total counts? to conserve cnts mult 2 ; ls1a=rebin(segs(*,1,0),2^rpow) fs1a=rebin(segs(*,1,1),2^rpow)/rscls1 ns1a=(rebin(segs(*,1,2),2^rpow)/rscls1)/sqrt(2.^(14.-rpow)) es1a=segs(*,1,5) cs1a=rebin(es1a^2,2^rpow)*2^(14.-rpow) ; ll1b=rebin(segs(*,2,0),2^rpow) fl1b=rebin(segs(*,2,1),2^rpow) nl1b=(rebin(segs(*,2,2),2^rpow))/sqrt(2.^(14.-rpow)) el1b=segs(*,2,5) cl1b=rebin(el1b^2,2^rpow)*2^(14.-rpow) ; ls1b=rebin(segs(*,3,0),2^rpow) fs1b=rebin(segs(*,3,1),2^rpow)/rscls1 ns1b=(rebin(segs(*,3,2),2^rpow)/rscls1)/sqrt(2.^(14.-rpow)) es1b=segs(*,3,5) cs1b=rebin(es1b^2,2^rpow)*2^(14.-rpow) ; ll2a=rebin(segs(*,4,0),2^rpow) fl2a=rebin(segs(*,4,1),2^rpow) nl2a=(rebin(segs(*,4,2),2^rpow))/sqrt(2.^(14.-rpow)) el2a=segs(*,4,5) cl2a=rebin(el2a^2,2^rpow)*2^(14.-rpow) ; ls2a=rebin(segs(*,5,0),2^rpow) fs2a=rebin(segs(*,5,1),2^rpow)/rscls2 ns2a=(rebin(segs(*,5,2),2^rpow)/rscls2)/sqrt(2.^(14.-rpow)) es2a=segs(*,5,5) cs2a=rebin(es2a^2,2^rpow)*2^(14.-rpow) ; ll2b=rebin(segs(*,6,0),2^rpow) fl2b=rebin(segs(*,6,1),2^rpow) nl2b=(rebin(segs(*,6,2),2^rpow))/sqrt(2.^(14.-rpow)) el2b=segs(*,6,5) cl2b=rebin(el2b^2,2^rpow)*2^(14.-rpow) ; ls2b=rebin(segs(*,7,0),2^rpow) fs2b=rebin(segs(*,7,1),2^rpow)/rscls2 ns2b=(rebin(segs(*,7,2),2^rpow)/rscls2)/sqrt(2.^(14.-rpow)) es2b=segs(*,7,5) cs2b=rebin(es2b^2,2^rpow)*2^(14.-rpow) ; ;Begin the join program. Will linearize to linearize (here .01 900 to 1190) newjoin,linearize,ll1a,fl1a,nl1a,cl1a,ls1a,fs1a,ns1a,cs1a,ll1b,fl1b,nl1b,cl1b,ls1b,fs1b,ns1b,cs1b,ll2a,fl2a,nl2a,cl2a,ls2a,fs2a,ns2a,cs2a,ll2b,fl2b,nl2b,cl2b,ls2b,fs2b,ns2b,cs2b,ftot,ntot,ctot,shifts print,[transpose(string(shifts.shift)),transpose(shifts.label)] endif ;This ends the read in and join segments program window,1,xsize=1200,ysize=500 !p.multi=0 xr=lminmax yr=[0,.8e-10] plot,linearize,ftot,psym=10,ticklen=-.025,charsize=cs,title=filetit,xtitle='Wavength ('+ang+')',ytit='Flux ('+fst+')',ys=1,yr=yr,xs=1,xr=xr ;plot,linearize(1500:*),smooth(ftot(1500:*),3),psym=10,ticklen=-.025,charsize=cs,title=filetit,xtitle='Wavength ('+ang+')',ytit='Flux ('+fst+')',ys=1,xs=1,xr=[-10,10]+1180 oiread,oidat ;stru {wav:0.0d0,el:0.0,eu:0.0,gl:0,gu:0,a:0.0,gam:0.0,f:0.0,lgf:0.0,llf:0.0} oi0=where(oidat.el eq 0.) oi0=oidat(oi0).wav noi=n_elements(oi0) ariread,aridat ;stru {wav:0.0d0,el:0.0,eu:0.0,gl:0,gu:0,a:0.0,gam:0.0,f:0.0,lgf:0.0,llf:0.0} ari0=where(aridat.el eq 0.) ari0=aridat(ari0).wav nari=n_elements(ari0) niread,nidat ;stru {wav:0.0d0,el:0.0,eu:0.0,gl:0,gu:0,a:0.0,gam:0.0,f:0.0,lgf:0.0,llf:0.0} ni0=where(nidat.el eq 0.) ni0=nidat(ni0).wav nni=n_elements(ni0) niiread,niidat ;stru {wav:0.0d0,el:0.0,eu:0.0,gl:0,gu:0,a:0.0,gam:0.0,f:0.0,lgf:0.0,llf:0.0} ;nii0=where(niidat.el eq 0.) ;nii0=niidat(nii0).wav nii0=niidat.wav nnii=n_elements(nii0) piiread,piidat ;stru {wav:0.0d0,el:0.0,eu:0.0,gl:0,gu:0,a:0.0,gam:0.0,f:0.0,lgf:0.0,llf:0.0} ;pii0=where(piidat.el eq 0.) ;pii0=piidat(pii0).wav pii0=piidat.wav npii=n_elements(pii0) ciread,cidat ;stru {wav:0.0d0,el:0.0,eu:0.0,gl:0,gu:0,a:0.0,gam:0.0,f:0.0,lgf:0.0,llf:0.0} ci0=where(cidat.el eq 0.) ci0=cidat(ci0).wav ;ci0=cidat.wav nci=n_elements(ci0) ciiread,ciidat ;stru {wav:0.0d0,el:0.0,eu:0.0,gl:0,gu:0,a:0.0,gam:0.0,f:0.0,lgf:0.0,llf:0.0} cii0=where(ciidat.el eq 0.) cii0=ciidat(cii0).wav ncii=n_elements(cii0) si_iiread,si_iidat ;stru {wav:0.0d0,el:0.0,eu:0.0,gl:0,gu:0,a:0.0,gam:0.0,f:0.0,lgf:0.0,llf:0.0} si_ii0=where(si_iidat.el eq 0.) si_ii0=si_iidat(si_ii0).wav nsi_ii=n_elements(si_ii0) siiread,siidat ;stru {wav:0.0d0,el:0.0,eu:0.0,gl:0,gu:0,a:0.0,gam:0.0,f:0.0,lgf:0.0,llf:0.0} sii0=where(siidat.el eq 0.) sii0=siidat(sii0).wav nsii=n_elements(sii0) feiiread,feiidat ;stru {wav:0.0d0,el:0.0,eu:0.0,gl:0,gu:0,a:0.0,gam:0.0,f:0.0,lgf:0.0,llf:0.0} feii0=where(feiidat.el eq 0.) feii0=feiidat(feii0).wav ;feii0=feiidat.wav nfeii=n_elements(feii0) feiiiread,feiiidat ;stru {wav:0.0d0,el:0.0,eu:0.0,gl:0,gu:0,a:0.0,gam:0.0,f:0.0,lgf:0.0,llf:0.0} ;feiii0=where(feiiidat.el eq 0.) ;feiii0=feiiidat(feiii0).wav feiii0=feiiidat.wav nfeiii=n_elements(feiii0) cs=1.25 hc='' read,'big hc',hc if hc eq 'y' then begin set_plot,'ps' device,filename='CAMM0606combspec.ps',xsize=32,ysize=32,/inches,/color,/times,/isolatin1,xoff=2,yoff=1,portrait=1 !p.font=0 rescol yr=[0,.8];[0,1.25] renorm=1e-10 !p.multi=[0,1,10] for i=900,1170,30 do begin plot,linearize,ftot/renorm,psym=10,xr=[i,i+30],xtitle='Wavelength ('+ang+')',ytitle='Flux (x10!U-10!N ergs cm!U-2!N s!U-1!N '+ang+'!U-1!N)',charsize=cs,yr=yr,color=9,ys=1 oplot,linearize+ismoff,fl/renorm,psym=10,color=5 oplot,linearize,minmod/renorm,psym=10,color=2 oplot,linearize,ftot/(10*renorm),linestyle=2,color=9 xyouts,i+29,.125,'1/10',charsize='2',color=9 ytr=[.1,.9]*!cymax ;LABELS ;h2label,jpp,vpp,lminmax,lamoff,cs,indata for j=0,5 do h2label,j,0,[i,i+30],ismoff/1065,.6,ncol(j)*1e21,bb,indata,2 atomicy=.9*!cymax glam=where(oi0 ge i and oi0 le i+30) if glam(0) ne -1 then begin for j=0,noi-1 do plots,[1,1]*oi0(j)+ismoff,ytr,color=5,linestyle=2,noclip=0,psym=0 xyouts,oi0(glam)+ismoff,atomicy*(0*oi0(glam)+1),'OI '+string(oi0(glam),'(f7.2)'),orient=90,color=5,charsize=.7 endif glam=where(ari0 ge i and ari0 le i+30) if glam(0) ne -1 then begin for j=0,nari-1 do plots,[1,1]*ari0(j)+ismoff,ytr,color=5,linestyle=2,noclip=0,psym=0 xyouts,ari0(glam)+ismoff,atomicy*(0*ari0(glam)+1),'ArI '+string(ari0(glam),'(f7.2)'),orient=90,color=5,charsize=.7 endif glam=where(ni0 ge i and ni0 le i+30) if glam(0) ne -1 then begin for j=0,nni-1 do plots,[1,1]*ni0(j)+ismoff,ytr,color=6,linestyle=2,noclip=0,psym=0 xyouts,ni0(glam)+ismoff,atomicy*(0*ni0(glam)+1),'NI '+string(ni0(glam),'(f7.2)'),orient=90,color=6,charsize=.7 endif glam=where(nii0 ge i and nii0 le i+30) if glam(0) ne -1 then begin for j=0,nnii-1 do plots,[1,1]*nii0(j)+ismoff,ytr,color=6,linestyle=2,noclip=0,psym=0 xyouts,nii0(glam)+ismoff,atomicy*(0*nii0(glam)+1),'NII '+string(nii0(glam),'(f7.2)'),orient=90,color=6,charsize=.7 endif glam=where(pii0 ge i and pii0 le i+30) if glam(0) ne -1 then begin for j=0,npii-1 do plots,[1,1]*pii0(j)+ismoff,ytr,color=7,linestyle=2,noclip=0,psym=0 xyouts,pii0(glam)+ismoff,atomicy*(0*pii0(glam)+1),'PII '+string(pii0(glam),'(f7.2)'),orient=90,color=7,charsize=.7 endif glam=where(ci0 ge i and ci0 le i+30) if glam(0) ne -1 then begin for j=0,nci-1 do plots,[1,1]*ci0(j)+ismoff,ytr,color=7,linestyle=2,noclip=0,psym=0 xyouts,ci0(glam)+ismoff,atomicy*(0*ci0(glam)+1),'CI '+string(ci0(glam),'(f7.2)'),orient=90,color=7,charsize=.7 endif glam=where(cii0 ge i and cii0 le i+30) if glam(0) ne -1 then begin for j=0,ncii-1 do plots,[1,1]*cii0(j)+ismoff,ytr,color=7,linestyle=2,noclip=0,psym=0 xyouts,cii0(glam)+ismoff,atomicy*(0*cii0(glam)+1),'CII '+string(cii0(glam),'(f7.2)'),orient=90,color=7,charsize=.7 endif glam=where(si_ii0 ge i and si_ii0 le i+30) if glam(0) ne -1 then begin for j=0,nsi_ii-1 do plots,[1,1]*si_ii0(j)+ismoff,ytr,color=8,linestyle=2,noclip=0,psym=0 xyouts,si_ii0(glam)+ismoff,atomicy*(0*si_ii0(glam)+1),'SiII '+string(si_ii0(glam),'(f7.2)'),orient=90,color=8,charsize=.7 endif glam=where(sii0 ge i and sii0 le i+30) if glam(0) ne -1 then begin for j=0,nsii-1 do plots,[1,1]*sii0(j)+ismoff,ytr,color=2,linestyle=2,noclip=0,psym=0 xyouts,sii0(glam)+ismoff,atomicy*(0*sii0(glam)+1),'SII '+string(sii0(glam),'(f7.2)'),orient=90,color=2,charsize=.7 endif glam=where(feii0 ge i and feii0 le i+30) if glam(0) ne -1 then begin for j=0,nfeii-1 do plots,[1,1]*feii0(j)+ismoff,ytr,color=3,linestyle=2,noclip=0,psym=0 xyouts,feii0(glam)+ismoff,atomicy*(0*feii0(glam)+1),'FeII '+string(feii0(glam),'(f7.2)'),orient=90,color=3,charsize=.7 endif glam=where(feiii0 ge i and feiii0 le i+30) if glam(0) ne -1 then begin for j=0,nfeiii-1 do plots,[1,1]*feiii0(j)+ismoff,ytr,color=4,linestyle=2,noclip=0,psym=0 xyouts,feiii0(glam)+ismoff,atomicy*(0*feiii0(glam)+1),'FeIII '+string(feiii0(glam),'(f7.2)'),orient=90,color=4,charsize=.7 endif ;Some Emission lines goto,skip_emission esi_iv=[1066.636,1066.650,1122.485,1128.325,1128.325,1128.34] nesi_iv=n_elements(esi_iv) esiv=[1062.664,1072.973,1073.519,1098.930,1099.481] nesiv=n_elements(esiv) eciii=[977.020,1174.933,1175.263,1175.590,1175.711,1175.987,1176.370] neciii=n_elements(eciii) epv=[1117.977,1128.008] nepv=n_elements(epv) eniii=[989.799 ,991.511,991.577] neniii=n_elements(eniii) elams=[esi_iv,esiv,eciii,epv,eniii] nelams=n_elements(elams) emstlab=[strarr(nesi_iv)+'SiIV',strarr(nesiv)+'SIV',strarr(neciii)+'CIII',strarr(nepv)+'PV',strarr(neniii)+'NIII'] ecolor=[intarr(nesi_iv)+2,intarr(nesiv)+3,intarr(neciii)+4,intarr(nepv)+5,intarr(neniii)+6] glam=where(elams ge i and elams le i+30) if glam(0) ne -1 then begin for j=0,nelams-1 do plots,[1,1]*elams(j)-.33,[.3,.7]*!cymax,color=4,linestyle=2,noclip=0 xyouts,elams(glam)-.33,!cymax*.3*(0*elams(glam)+1),emstlab(glam)+' '+string(elams(glam),'(f7.2)'),orient=-90,color=ecolor(glam),charsize=.7 endif skip_emission:dum=0 endfor ;Footnotes --- cmm2=' (cm!e-2!n) ' kms=' (km s!e-1!n) ' xyouts,1190,yr(1)*.9,'Kurucz Atlas'+'!c!c*'+string(t,'(I6)')+'!c!c*'+string(g,'(F4.1)')+'!c!c* E(B-V)='+string(ebv,'(F5.2)')+'Extinction Curve '+string(ccmpar,'(6f7.3)')+'!c!c* Dist='+string(distpc,'(f6.0)')+'pc'+'!c!c* '+string(n_elements(hcol),'(I2)')+' Component '+string('HI','(A3,1x)')+' Absorption Model'+'!c!c* Columns '+cmm2+'= '+string(10.^hcol,'(2(g7.2,1x))')+'!c!c* Velocities '+kms+'= '+string(hvkms,'(2(f6.1,1x))')+'!c!c* 1 Component '+string('H2','(A3,1x)')+' Absorption Model'+'!c!c* Jn Columns '+cmm2+'= '+string(ncol(where(ncol ne 0)),'(6(g7.2,1x))')+'!c!c* Velocities '+kms+'= ?'+'!c!c* b '+kms+'='+bbs+'!c!c* Resolution='+string(respower,'(I6)')+'!c!cStephan R. McCandliss - JHU!c!c'+systime(0),color=9,font=0,charsize=.6 device,/close !p.font=-1 pm0 set_plot,'x' endif oplot,lj0,minmod,color=255 ;oplot,lj1+.05,h2av1*4e-13,color=255 ; oplot,lj0-.2,minmod*h1c ;oplot,lj1+.15,h2av1*5.5e-13,color=255 ;oplot,lj0+.34,minmod,color=255 ;oplot,lj0+.34,fl h2s=[1062.8821,1064.6061,1063.4603,1066.9004,1064.9944,1070.1411,1067.4786,1074.3125,1070.8997,1079.4004,1049.3674,1051.0332,1049.9598,1053.2839,1051.4980,1056.4716,1053.9760,1060.5806,1057.3800,1065.5964,1061.6968,1071.4967,1066.9054,1078.2660,1072.9918] close,1 cii=[1036.337,1037.018] ovi=[1031.926,1037.617] siv=[1062.664,1072.973,1073.519] ciiiex=[1174.933,1175.263,1175.590,1175.711,1175.987,1176.370] svi=[933.378,944.523] pv=[1117.977,1128.008] ciii=[977.020] niii=[989.799 ,991.511,991.577] vline,ciiiex c=2.998e5 wset,0 kmspan=500 xti='Velocity (km s!e-1!n)' !p.multi=[0,1,2] plot,(linearize-cii(0))/cii(0)*c,ftot,xr=[-1,1]*kmspan,psym=10,xs=9,xtitle=xti,charsize=cs,ytitle='Flux ('+fst+')' oplot,(linearize-cii(0)+ismoff)/cii(0)*c,minmod,color=255 axis,0,!cymax,xax=1,xr=([-1,1]*kmspan/c+1.)*cii(0),xs=1,xtitle='CII ' +string(cii(0),'(f8.3)'),charsize=cs*1.1 plot,(linearize-cii(1))/cii(1)*c,ftot,xr=[-1,1]*kmspan,psym=10,xs=9,xtitle=xti,charsize=cs,ytitle='Flux ('+fst+')' oplot,(linearize-cii(1))/cii(1)*c,minmod,color=255 axis,0,!cymax,xax=1,xr=([-1,1]*kmspan/c+1.)*cii(1),xs=1,xtitle='CII ' +string(cii(1),'(f8.3)'),charsize=cs*1.1 set_plot,'x' !p.font=-1 !p.multi=[0,1,2] set_plot,'ps' device,file='fluxandsignaltonoise.ps',/isolatin1,/times,/landscape !p.font=0 plot,linearize,sqrt(ctot),psym=10,xtitle='Wavelength ('+ang+')',ytitle='(Counts)!e1/2!n (0.01'+ang+')!e-1!n',title=strmid(files(0),32,4)+'-- from combined spectral total',color=9 plot,linearize,ftot,psym=10,xtitle='Wavelength ('+ang+')',ytitle='Flux ('+fst+')',title=strmid(files(0),32,4)+'-- from combined spectral total',color=9 oplot,linearize,minmod,color=2 device,/close set_plot,'x' !p.font=-1 !p.multi=0 goto,skipprint rlin=900+.1*findgen(290/.1) rftot=rebin(ftot,2900) rntot=rebin(ntot,2900) rminmod=rebin(float(minmod),2900) rfl=rebin(float(fl),2900) openw,1,'EG01rebinnedcombined09082000.dat' headout=head(5)+' '+head(6)+'* E(B-V)='+string(ebv,'(F5.2)')+'* CCM Extinction Curve'+'* Dist='+string(distpc,'(f6.0)')+'pc'+'* '+string(n_elements(hcol),'(I2)')+' Component '+string('HI','(A3,1x)')+' Absorption Model'+' Columns '+cmm2+'= '+string(10.^hcol,'(2(g7.2,1x))')+' Velocities '+kms+'= '+string(hvkms,'(2(f6.1,1x))')+'* 1 Component '+string('H2','(A3,1x)')+' Absorption Model'+' Jn Columns '+cmm2+'= '+string(ncol(where(ncol ne 0)),'(6(g7.2,1x))')+'* Velocities '+kms+'= ?'+'* b '+kms+'='+bbs+'* Resolution='+string(respower,'(I6)') printf,1,headout printf,1,'lam favg noise model+abs model' printf,1,[transpose(rlin(50:2899-12)),transpose(rftot(50:2899-12)),transpose(rntot(50:2899-12)),transpose(rminmod(50:2899-12)),transpose(rfl(50:2899-12))] plot,rlin,rftot oplot,rlin,rminmod oplot,rlin,rfl close,1 skipprint:dum=0 end