subroutine h2opf(t,pf)
c   
c     partition function for H2Ofrom EXOMOILA data
c
      INCLUDE 'PARAMS.FOR'
      dimension ttab(10000),pftab(10000)
c   
      data init /1/
c
      if(init.eq.1) then
         open(67,file='./data/h2o_exomol.pf',status='old')
         do i=1,10000
            read(67,*) ttab(i),pftab(i)
         end do 
         close(67)
         init=0
      end if
c   
      itab=ifix(real(t))
      pf=pftab(itab)+(t-ttab(itab))*(pftab(itab+1)-pftab(itab))
      return
      end