Liste de commandes a ajouter dans CHIMERE pour effectuer des etudes de sensibilite.
v. Decembre 98.
Pour retrouver les insertions dans les fortrans -> grep lme *.f
-------------------------------------------------------------------------
1. Creer le fichier de common 'sensib.h', contenant:

      character*15 charsensib
      character*6 datecar
      common /sensib1/ charsensib
      common /sensib3/ datecar
      
      double precision spoll
      common /sensib2/ spoll
      
-------------------------------------------------------------------------
2. Lire la chaine de caractere donnant le nom de l'espece a etudier.
Dans CHIMV5 -> routine Binipoll.f
Dans CHIMV6 -> routine Binisens.f

Ajouter:

c-lme - Ajout de l'appel  'sensib.h'     
      include 'sensib.h'
c-lme - Affectation du nom de la 1ere variable etudiee dans charsensib  
	   if (np.eq.1)charsensib=charspec
c-lme -Ouverture des fichiers de resultats pour awrateloc haute frequence
      open(86,file=datecar//'_sens_rates_Surf_'//charsensib)
      open(87,file=datecar//'_sens_rates_Mixd_'//charsensib)

-------------------------------------------------------------------------
3. Ajouter les fonctionnalites de creation de fichiers XMGR

en ajoutant la routine Boutlme.f
et en l'appelant dans:
- Binteghour.f:

if(mod(nostep,nprint).eq.0) call outprint
if(mod(nostep,nprint).eq.0) call outlme

- Bintegrun.f:
call outprint
call outlme
-------------------------------------------------------------------------
4. Verifier que l'adjoint des rates est bien initialise dans Binichem.f

      do nr=1,nreac
         do nh=0,nhourrunmax
	    do nv=1,nverti
	       wrate(nr,nh,nv) = dun
	    enddo
	 enddo
      enddo

-------------------------------------------------------------------------
5. Verifier que le poids des rates est bien interpole dans Blocvalues.f

c Sensitivity to reactions rates  
      do nb=1,nbox
         do nr=1,nreac
	    wrateloc(nr,nb) = w2*wrate(nr,ihourrun,nb)
     $                      + w1*wrate(nr,ihourrun+1,nb)	    
	 enddo
      enddo

et adjointise

c  Interpolation of reactions rates (lme)

      do nv=1,nverti
         do nr=1,nreac
            awrate(nr,ihourrun,nv) = awrate(nr,ihourrun,nv) + 
     $                               awrateloc(nr,nv)*w2
            awrate(nr,ihourrun+1,nv) = awrate(nr,ihourrun+1,nv) + 
     $                               awrateloc(nr,nv)*w1
         enddo
      enddo

-------------------------------------------------------------------------
6. Ouverture des fichiers de resultats dans Biniread.f

c-lme - Appel de 'sensib.h'      
      include 'sensib.h'
c-lme  Opening files for XMGR results

      rewind(33)
      read(33,'8a')datecar
      rewind(33)
      
      open(70,file=datecar//'_1opspec_'//fnout)
      open(71,file=datecar//'_2opspec_'//fnout)
      open(72,file=datecar//'_3opspec_'//fnout)
      open(73,file=datecar//'_4opspec_'//fnout)
      open(74,file=datecar//'_5opspec_'//fnout)
      open(75,file=datecar//'_6opspec_'//fnout)

      open(80,file=datecar//'_1opspec_24h_'//fnout)
      open(81,file=datecar//'_2opspec_24h_'//fnout)
      open(82,file=datecar//'_3opspec_24h_'//fnout)
      open(83,file=datecar//'_4opspec_24h_'//fnout)
      open(84,file=datecar//'_5opspec_24h_'//fnout)
      open(85,file=datecar//'_6opspec_24h_'//fnout)

