1 |       program testban
   2 | cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
   3 | c     module
   4 | 
   5 | c     variables
   6 |       implicit none
   7 |       include 'propradia.inc'
   8 |       include 'propradiabis.inc'
   9 | 
  10 |       include 'cecile.inc'
  11 |       include 'radiatif.inc'
  12 | 
  13 |       include 'entre.inc'
  14 | 
  15 |       integer it, i, ii
  16 |       double precision rt
  17 | c      double precision temperat, ptot, fmco, fmco2, fmh2o
  18 | 
  19 | 
  20 | ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
  21 | c     lecture de la structure en bande commune a tous les spectres
  22 | c     moleculaires
  23 |     
  24 | 
  25 |       open(unit=9, file='SNBWN')
  26 | c     lecture initial pour boucle conditionnelle     
  27 |       ibande =1 
  28 |       read(9,*) eta(ibande), delta_eta(ibande)
  29 | 
  30 | c     lecture boucle     
  31 |       ibande=1
  32 |       do while (eta(ibande).ge.0)
  33 |          read(9,*) eta(ibande+1), delta_eta(ibande+1)
  34 |          ibande=ibande+1
  35 |       enddo
  36 | c     si on decouvre le nombre de bande
  37 | c     nbande =ibande - 1
  38 |       close(9)
  39 | 
  40 | c     test ecriture
  41 | c      do ibande=1, nbande
  42 | c         print *, eta(ibande), delta_eta(ibande)
  43 | c      enddo
  44 | 
  45 | c------------------------------------------------------------------
  46 |       
  47 |       call parambgaz
  48 | c     correspondance index bande et index param de bande
  49 |       call paramind
  50 | 
  51 | c     test ecriture
  52 | c      do ibande=1 ,nbande
  53 |          
  54 | c      call genere_b (1,ibande)
  55 | c      print *, ibande, ICO(ibande), ICO2(ibande), IH2O(ibande)
  56 | c      print *, '   ',LICO(ibande), LICO2(ibande), LIH2O(ibande)
  57 | c      read(*,*)
  58 | 
  59 | c      enddo
  60 | 
  61 | c------------------------------------------------------------------
  62 |       temp(1)=2400.
  63 |       call modbgazinterp(temp(1),rt,it)
  64 | 
  65 | c     test ecriture
  66 | c       print *, rt ,it
  67 | c       read(*,*)
  68 | c-----------------------------------------------------------------
  69 |       do ii=1,1000
  70 | 
  71 | c      ibande =300
  72 |       call genere_b(1,ibande)
  73 | 
  74 | 
  75 |       ptot=1.
  76 |       fm(1,1)=0.
  77 |       fm(2,1)=0.6
  78 |       fm(3,1)=0.8
  79 | 
  80 | c      call modbgazdissoc(ptot,temperat,fmco,fmco2,fmh2o)
  81 |       call modbgaz(ptot,temp(1),fm(1,1),fm(2,1),fm(3,1))
  82 |       print *,''
  83 |       print *,'ibande',ibande
  84 |       print *, 'les parametres'
  85 |       print *,'--------------'
  86 |       print *, LICO(ibande)
  87 |       print *, LICO2(ibande)
  88 |       print *, LIH2O(ibande)
  89 | 
  90 | c      do i=1,3
  91 | c         
  92 | c      print *, kgb6p(i,ibande), dinv(i,ibande),
  93 | c     &         gamma(i,ibande), 
  94 | c     &         phi(i,ibande) 
  95 | c         
  96 | c      enddo
  97 | 
  98 | c      read(*,*)
  99 |       enddo
 100 | 
 101 | c      print *, dinv_piv(3,11,45)
 102 | c      print *, kgb_piv(3,11,45)
 103 | 
 104 | c----------------------------------------------------------------
 105 | c     melange de gaz
 106 |       in =5
 107 |       iin=1
 108 | c      call modbgaz
 109 | c     attention ibande est fixe ?
 110 |       do iin=1, n_mx
 111 | 
 112 |       fm(3,iin)=0.1*iin
 113 | 
 114 |       call modbgaz(ptot,temp(1),fm(1,1),fm(2,1),fm(3,1))
 115 | 
 116 |       print *,''
 117 |       print *, 'kgbar et phig'
 118 |       print *, kgbar(1,iin), phig(1,iin)
 119 | 
 120 |       enddo
 121 | 
 122 | ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
 123 |       end