1 |       SUBROUTINE trajet_fonc(w, bi, bj)
   2 |       IMPLICIT NONE
   3 | c
   4 | c&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
   5 | c expression de la fonstion psi a integrer
   6 | c     in -->
   7 | c           in, iin, w, bi , bj
   8 | c
   9 | c
  10 | c     out -->
  11 | c           psi
  12 | c           dpsi6d?  cf cecile.inc
  13 | c&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
  14 | c
  15 | c.......................................................................
  16 | c
  17 |       INCLUDE 'cecile.inc'
  18 |       include 'propradia.inc'
  19 | c      include 'radiatif.inc'
  20 | c      include 'entre.inc'
  21 | 
  22 |       double precision w
  23 |       double precision bi, bj
  24 | 
  25 | c.......................................................................
  26 | c      DOUBLE PRECISION tmp1,tmp2,tmp3,tmp4,tmp5
  27 | c      INTEGER itmp1,itmp2,itmp3,itmp4,itmp5
  28 | c.......................................................................
  29 | 
  30 | 
  31 | 
  32 | 
  33 | 
  34 |        
  35 | 
  36 | 
  37 | c -----{calcul de psi...}
  38 |      
  39 |         psi(in,iin)=psi(in,iin)+w*delta_eta(ibande)*1.d+2*(bi-bj)
  40 | 
  41 |         var_psi(in,iin)=var_psi(in,iin)+
  42 |      &       (w*delta_eta(ibande)*1.d+2*(bi-bj))**2
  43 | 
  44 | 
  45 |         !if ((in.eq.4).and. (iin.eq.5))then
  46 |         !if ((in.eq.iin))then
  47 |         !   print *, 'in, iin, psi(in, iin), w'
  48 |         !   print *, in, iin, psi(in, iin), w, (bi-bj), itir
  49 |         !   read *
  50 |         !endif
  51 |         if (in .eq. 60) then
  52 |            psi60bande(ibande)= psi60bande(ibande)+ psi(in,iin)
  53 |         endif
  54 | c
  55 | c.......................................................................
  56 | c
  57 |       RETURN
  58 |       END