;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;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. ; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;The following code will read in the data into idl ;(after they are uncompressed with gunzip) ncol=[200,200,1.,1,.05,.001,fltarr(10)]*1e18 ; a guess at the columndensities td='taufiledirectory/' ;Where you store the tauh2's bb=1 ;the b you want - integer only bbs=strtrim(string(fix(bb)),2) openr,1,td+'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) tauh2=lj0*0 for j=0,15 do tauh2=tauh2+v0(j+1)*ncol(j)/1e21 h2a=exp(tauh2) close,1 end