1 |       SUBROUTINE TPARAM
   2 | C*****************************************************************
   3 | C                                                                *
   4 | C   THIS SUBROUTINE READS THE MODEL PARAMETERS                   *
   5 | c    prog taine param
   6 | c
   7 | c     in --> nothing
   8 | c     out--> common intaine.inc
   9 | c            k: coeff d'absorption moyen sur une bande etroite   
  10 | c               (eta et delta_eta et temperature fixes par la
  11 | c               structure de bande de Taine) par unite d'atmos
  12 | c               unites: cm-1 atm-1
  13 | c
  14 | c            dinv: moyennne de l'inverse de la 1/2 largeur de raie
  15 | c                  a mi-hauteur.
  16 | c                  unites: cm
  17 | C                                                                *
  18 | C*****************************************************************
  19 | 
  20 | 
  21 | 
  22 | c l'utilisation de intaine necessite tjrs l'appel a gazsuie
  23 | c pour nbande_mx et ngaz_mx
  24 |       include 'gazsuie.inc'
  25 |       include 'intaine.inc'
  26 | c      include 'gaz.inc'
  27 | c      REAL KCO,KC,KH
  28 | c      COMMON /SPCD/DCO(14,48),DC(14,367),DH(14,367)
  29 | c      COMMON /SPCK/KCO(14,48),KC(14,367),KH(14,367)
  30 | c      include gaz.inc avec ajout indice pour nature du gaz
  31 | c      double precision k(ngaz_mx,nbande_mx,ntemp_mx)
  32 | c      double precision dinv(ngaz_mx,nbande_mx,ntemp_mx)
  33 |       integer I,J 
  34 | 
  35 | 
  36 | 
  37 |       OPEN(UNIT=45,FILE='SNBCO')
  38 |       OPEN(UNIT=46,FILE='SNBCO2')
  39 |       OPEN(UNIT=47,FILE='SNBH2O')
  40 | C
  41 | C   READING THE PARAMETERS
  42 | C
  43 | 
  44 | c      lecture pour le co=1
  45 |       DO  I=1,48
  46 | c      READ(45,*) (KCO(1,J,I),J=1,12)
  47 |       READ(45,*) (k(1,I,J),J=1,12)
  48 |       enddo
  49 | 
  50 |       DO  I=1,48
  51 | c      READ(45,*) (DCO(1,J,I),J=1,12)
  52 |       READ(45,*) (dinv(1,I,J),J=1,12)
  53 |       enddo
  54 | 
  55 | c   lecture pour le co2=2
  56 |       DO 3 I=1,96
  57 | c3     READ(46,*) (KC(2,J,I),J=1,14)
  58 | 3     READ(46,*) (k(2,I,J),J=1,14)
  59 | 
  60 |       DO 4 I=1,96
  61 | c4     READ(46,*) (DC(2,J,I),J=1,14)
  62 | 4     READ(46,*) (dinv(2,I,J),J=1,14)
  63 | 
  64 | c   lecture pour le h20=3
  65 |       DO 5 I=1,367
  66 | c5     READ(47,*) (KH(3,J,I),J=1,14)
  67 | 5     READ(47,*) (k(3,I,J),J=1,14)
  68 | 
  69 |       DO 6 I=1,367
  70 | c6     READ(47,*) (DH(3,J,I),J=1,14)
  71 | 6     READ(47,*) (dinv(3,I,J),J=1,14)
  72 | 
  73 |       close(45)
  74 |       close(46)
  75 |       close(47)
  76 | 
  77 | 
  78 |       RETURN
  79 |       END
  80 | 
  81 |