pro h2data,vpp,lminmax,indata ;;;indata(numlines,8) ;;;; v' v'' j' j'' lam f g'' atot ;;;; '=upperlevel and ''=lowerlevel dir='/home/stephan/idlpro/h2/' if vpp eq 0 then begin openr,1,dir+'highjsh2vpnvpp0.dat' numlines=1580 endif if vpp eq 1 then begin openr,1,dir+'highjsh2vpnvpp1.dat' numlines=1580 endif if vpp eq 2 then begin numlines =1580 openr,1,dir+'highjsh2vpnvpp2.dat' endif if vpp eq 3 then begin numlines =1580 openr,1,dir+'highjsh2vpnvpp3.dat' endif if vpp eq 4 then begin numlines =1580 openr,1,dir+'highjsh2vpnvpp4.dat' endif if vpp eq 5 then begin numlines =1580 openr,1,dir+'highjsh2vpnvpp5.dat' endif if vpp eq 6 then begin numlines =1580 openr,1,dir+'highjsh2vpnvpp6.dat' endif if vpp eq 7 then begin numlines =1580 openr,1,dir+'highjsh2vpnvpp7.dat' endif indatas=dblarr(8,numlines) readf,1,indatas close,1 ;select j'' ;j0=where(indatas(3,*) eq jpp) ;indata=indatas(*,j0) indata=indatas ;numline=n_elements(j0) bandpass=where(indata(4,*) ge lminmax(0) and indata(4,*) le lminmax(1)) if bandpass(0) ne -1 then indata=indata(*,bandpass) else goto,leave ;order of indatas is v' v'' j' j'' lam f g gamma'' ps=where(indata(2,*) lt indata(3,*)) qs=where(indata(2,*) eq indata(3,*)) rs=where(indata(2,*) gt indata(3,*)) leave:print,'leaving h2data' return end