subroutine chemsoot(fvsoot, fmoletyl, P, temperature, dt)
c*********************************************************************
c ce programme calcul la formation des suies
c c'est un post-traitement sur un solution de combustion
c c'est justifie car 
c     1/ on ne fabrique pas d'especes a reinjecter (faux!!!)
c     2/ la constante de temps est immense ==> la perte
c        d'espece C2H2 qu'on introduit est tres faible...
c
c     in: fmol c2h2 fraction molaire de C2h2 produite
c                   pendant DELTA t
c         la pression en pascal (1,013 e+5=atmos)
c         la temperturre en Kelvin
c         l'ecart de temps dt
c
c     out: fv C(s)  fraction volumique d'atome de Carbone
c                   (1er etape (pure chimie)/ 4: Modele de TIEN etc)
c                   (A introduire dans le fichier.scheme de CHEMKIN)
c*********************************************************************
c autre Q: gerer le transport des suies dans code de fokion 
c autre Q: gerer le radiatif ses bons 
      
c     module
c     neant

      implicit none
      double precision  dt, fvsoot,fmoletyl, P, temperature
      double precision vmol_carb, R, somfmol

      parameter (vmol_carb=0.00615) !volume molaire du carbone
      parameter (R= 8.315)          !cte des gaz parfaits
      parameter (somfmol = 1.)      !somme locale des frac molaire




c---------------------------------------------------------------------
c (1er etape (pure chimie)/ 4: Modele de TIEN etc)
c (A introduire dans le fichier.scheme de CHEMKIN)

      if ( dt .eq. 0.) then
         ! pour une seconde
         fvsoot = (1.e+4*dexp(-21100/temperature)*fmoletyl*vmol_carb)
     & /((R*temperature/P)*somfmol)






      else
         print *, 'CHEMSOOT: WARNING: instationnary solution not 
     &             yet implemented'
      endif

c----------------------------------------------------------------------
! the tree other steps are not yet impanted







c----------------------------------------------------------------------
      return
      end