pro initPartitionFunction, qtname, nmole, nspeci, isovec, q296, qcoef

  buf=''
  nspeci = 0
  nmole = 0
  d = 0.d
  d1 = 0.d
  d2 = 0.d
  d3 = 0.d
  d4 = 0.d
  
  OPENR, fptr, qtname, /get_lun, ERROR = err
  IF (err NE 0) then return
  
  readf, fptr, buf; 	/* read first text line */ 
  readf, fptr, buf

  reads,buf,nmole,nspeci,format='(I,I)'
  
  isovec = intarr(nmole)
  q296 = dblarr(nspeci)
  qcoef = dblarr(nspeci, 3, 4)

  readf, fptr, buf; 	/* read isonm title line */ 
  for i=1, nmole-1 do begin &$ ; 	/* read isotope vector */ 
	readf, fptr, buf &$
	reads,buf,j,format='(I)' &$
	isovec[i] = isovec[i-1] + j &$
  end      

  readf, fptr, buf;    /* read last isotope vector element */ 
  readf, fptr, buf;    /* read title line */ 

  for i=0, nspeci-1 do begin &$
	readf, fptr, buf &$
	reads,buf,d,format='(D)' &$
	q296[i] = d &$
  end
  
  for range = 0, 2 do begin &$
	readf, fptr, buf &$ ;     /* read title line */ 
     for i=0, nspeci-1 do begin &$
		readf, fptr, buf &$
		reads,buf,d1,d2,d3,d4 &$
		qcoef[i,range,*] = [d1,d2,d3,d4] &$
	end &$
  end 

  close, fptr

end