program hitrandriver

c **********************************************************************
c *                                                                    *
c *                                                                    *
c *                                                                    *
c *                                                    Leslie Young    *
c *                                                    02 July 1999    *
c *                                                                    *
c **********************************************************************

        implicit double precision (o-z, a-h)
        include "parameters.dat"

        dimension 

c       ---------------------------------------------------------
c            Arrays : Model atmosphere
c       ---------------------------------------------------------
     &   ifund(5)

c       ---------------------------------------------------------
c            Arrays : line listing
c       ---------------------------------------------------------
     &  , wav(nd4),str(nd4),wid(nd4),eng(nd4),nlines(5),nl(5)
     &  , stot(5),wmean(5),wstart(5),wend(5)
c     &  , wstep(5)

c       ---------------------------------------------------------
c                           Data for Titan 
c       ---------------------------------------------------------
        data ifund/1,5,7,8,9/

c       ----------------------------------------------------------
c             Read in line listing for fundamentals
c       ----------------------------------------------------------
        n = 1

        do imol = 1, 5
           nl(imol) = n
           jb = ifund(imol)
c           selmin = 2809.
c           selmax = 3210.
           call sread(jb,wav(n),str(n),wid(n),eng(n),nlines(imol)
     &      ,stot(imol),wmean(imol),wstart(imol),wend(imol))
c           wdopref(imol)=dsqrt(two*rkb*tref/rmass(imol)/amu)
c     &      *wmean(imol)/cs
c           opp0(imol) = stot(imol)/wdopref(imol) 
c           wstep(imol) = wdopref(imol)
           
           write(*,*), imol, wmean(imol), stot(imol), nlines(imol)

           n = n + nlines(imol)

        end do



        end