1 |         SUBROUTINE modbsuie(fvsuie)
   2 | C**********************************************************************
   3 | C
   4 | C   THIS SUBROUTINE calcule les trois parametres du modele de gaz de
   5 | c                   bande etroite pour le melange de gaz
   6 | c                   pour la maille ---IIN----------- ou on se trouve
   7 | c                  
   8 | c     in -->  profils 
   9 | c             in et tous les iin par cecile.inc
  10 | c               
  11 | c
  12 | c     out --> kgbar et phi maille d'emission
  13 | c             kgbar pour toutes les autres 
  14 | c
  15 | c
  16 | c
  17 | c     ATTENTION: Toute modif du modele doit etre accompagne des modif
  18 | c     identique dans trajet.f lors des calculs des sensibilites
  19 | c
  20 | C*********************************************************************
  21 |      
  22 | 
  23 | c     declarations 
  24 |       implicit none
  25 |       include 'cecile.inc'
  26 | c     necessaire pour radiatif
  27 |       
  28 |       include 'propradia.inc'
  29 | c     car il faut ngaz_mx
  30 |       include 'propradiabis.inc'
  31 | 
  32 |       include 'radiatif.inc'
  33 | c     pour kgbar et phi
  34 |       include 'entre.inc'
  35 | c     pour profils temp ptot et concentrations
  36 |       double precision fvsuie !fv suie locale ou je suis en iin
  37 | 
  38 | 
  39 | 
  40 | ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
  41 | c       lee et Tien pour sft98
  42 | c       sur une bande etroite ks = ksbar = ks pour eta(ibande) par exemple
  43 | c       fois cent  car eta et cm-1 et va en m-1
  44 | 
  45 |         !ksbar(iin) = 5.5D+0 * eta(ibande)* 1.d+2 * fv(iin)
  46 |         ksbar(iin) = 5.5D+0 * eta(ibande)* 1.d+2 * fvsuie
  47 | 
  48 |       RETURN
  49 |       END