pro get_ice2,first,last

openw, 2, '/home/jeneen/intensity/idl/out_1.5.ice' 

for i=first, last do begin
    ice = '/tmp80/jeneen/ice/ice.' + strtrim(i,2)
    dark = '/tmp80/jeneen/dark/dark.' + strtrim(i,2)
    if (file_exist(ice) and file_exist(dark)) then begin
      itmp = rfits(ice,head=ihd)
      dtmp = rfits(dark,head=dhd)

      t_ref = strtrim(sxpar(ihd, 'T_REF'),1)
      missvals = sxpar(ihd, 'MISSVALS')
      obs_dist = sxpar(ihd, 'OBS_DIST')
      imean = sxpar(ihd, 'DATAMEAN')
      dmean = sxpar(dhd, 'DATAMEAN')
      data = itmp - dtmp

      tot = float(0)
      for x = 465, 564 do begin
        for y = 465, 564 do begin
          tot = tot + data(x,y)
        end
      end
       
      itot = float(0)
      for x = 465, 564 do begin
        for y = 465, 564 do begin
          itot = itot + itmp(x,y)
        end
      end

      dtot = float(0)
      for x = 465, 564 do begin
        for y = 465, 564 do begin
          dtot = dtot + dtmp(x,y)
        end
      end

      avg = tot/10000
      iavg = itot/10000
      davg = dtot/10000

     print, strtrim(t_ref,2)+'  '+strtrim(avg,2)+'  '+strtrim(iavg,2) +'  '+strtrim(davg,2) +'  '+strtrim(missvals,2) +'  '+strtrim(imean,2)+'  '+strtrim(dmean,2)+'  '+strtrim(obs_dist,2)
     printf,2, strtrim(t_ref,2)+'  '+strtrim(avg,2)+'  '+strtrim(iavg,2) +'  '+strtrim(davg,2) +'  '+strtrim(missvals,2) +'  '+strtrim(imean,2)+'  '+strtrim(dmean,2)+'  '+strtrim(obs_dist,2)
  end
end

close, 2


end