1 |       SUBROUTINE parambgaz
   2 | C*****************************************************************
   3 | C                                                                *
   4 | C   THIS SUBROUTINE READS THE MODEL PARAMETERS                   *
   5 | c   lecture des parametres du modele de bande pour les gaz
   6 | c   c'est la subroutine param du prog de Taine
   7 | c
   8 | c     in --> nothing
   9 | c     out--> common propradia.inc
  10 | c            kgb_piv: coeff d'absorption moyen sur une bande etroite   
  11 | c               (eta et delta_eta et temperature fixes par la
  12 | c               structure de bande de Taine) par unite d'atmos
  13 | c               unites: cm-1 atm-1
  14 | c
  15 | c            dinv_piv: moyennne de l'inverse de la 1/2 largeur de 
  16 | c                  raie a mi-hauteur.
  17 | c                  unites: cm
  18 | C                                                                *
  19 | C*****************************************************************
  20 | 
  21 | 
  22 |       implicit none
  23 |       include 'propradia.inc'
  24 |  
  25 | c      include 'gaz.inc'
  26 | c      REAL KCO,KC,KH
  27 | c      COMMON /SPCD/DCO(14,48),DC(14,367),DH(14,367)
  28 | c      COMMON /SPCK/KCO(14,48),KC(14,367),KH(14,367)
  29 | c      include gaz.inc avec ajout indice pour nature du gaz
  30 | c      double precision k(ngaz_mx,nbande_mx,ntemp_mx)
  31 | c      double precision dinv(ngaz_mx,nbande_mx,ntemp_mx)
  32 | c     i est l'indexe de bande tronque
  33 | c     j est l'indexe de la temperature
  34 |       integer I,J 
  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,*) (kgb_piv(1,J,I),J=1,14)
  48 |        enddo
  49 | 
  50 |       DO  I=1,48
  51 | c      READ(45,*) (DCO(1,J,I),J=1,12)
  52 |       READ(45,*) (dinv_piv(1,J,I),J=1,14)
  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,*) (kgb_piv(2,J,I),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_piv(2,J,I),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,*) (kgb_piv(3,J,I),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_piv(3,J,I),J=1,14)
  72 | 
  73 |       close(45)
  74 |       close(46)
  75 |       close(47)
  76 | 
  77 | 
  78 |       RETURN
  79 |       END