program titplot
        implicit double precision (o-z, a-h)
        include "parameters.dat"
        character fname*64

c       ---------------------------------------------------------
c            Arrays : Model atmosphere
c       ---------------------------------------------------------
        dimension zalt(nd3),pres(nd3),temp(nd3),dens(nd3,0:4)
     &   ,col(nd3,0:4)


        write(unit=*, fmt='(19h Enter file name > ,$)') 
        read (unit=*, fmt=900) fname
        call atmread2(fname,1.d3,natm,zalt,pres,temp,dens,col)

        open(unit=60,file='titan.plot')
          do n = 1, natm
             dn0 = dens(n,0)
             f1 = dens(n,1)/dn0
             f2 = dens(n,2)/dn0
             f3 = dens(n,3)/dn0
             f4 = dens(n,4)/dn0
             write(unit=60,fmt=901) 1.e-5*zalt(n),pres(n),temp(n)
     &        ,f1,f2,f3,f4
          end do
        close(unit=60)
 900    format(a)
 901    format(1x,f8.2,e12.4,1x,f8.2,4e12.4)
        stop
        end