1 | include(dom.inc)
   2 | 
   3 |        SUBROUTINE PARAM(KCO,KC,KH,DCO,DC,DH,pathspec,ntbands,nb_base)
   4 | 
   5 |        include 'dom_constants.h'
   6 | 
   7 | !      IN  
   8 |        CHARACTER*80 :: pathspec,nb_base
   9 |        DOM_INT      :: ntbands
  10 | 
  11 | !      LOCAL
  12 |        DOM_INT      :: I, J
  13 | 
  14 | !      OUT
  15 |        DOM_REAL,DIMENSION(14,ntbands)::KCO,DCO,KC,DC,KH,DH
  16 | 
  17 |         IF (nb_base.eq.'SNB') THEN
  18 | 
  19 |           PRINT*, "  -> Using narrow band database SNB from EM2C"
  20 | 
  21 |           OPEN(UNIT=1010,FILE=trim(pathspec)//'/SNBH2O_K')
  22 |           OPEN(UNIT=1020,FILE=trim(pathspec)//'/SNBH2O_Phi')
  23 |           OPEN(UNIT=1030,FILE=trim(pathspec)//'/SNBCO2_K')
  24 |           OPEN(UNIT=1040,FILE=trim(pathspec)//'/SNBCO2_Phi')
  25 |           OPEN(UNIT=1050,FILE=trim(pathspec)//'/SNBCO_K')
  26 |           OPEN(UNIT=1060,FILE=trim(pathspec)//'/SNBCO_Phi')
  27 | 
  28 |           DO I=1,48
  29 |             READ(1050,*) (KCO(J,I),J=1,14)
  30 |             READ(1060,*) (DCO(J,I),J=1,14)
  31 |           END DO
  32 | 
  33 |           DO I=1,96
  34 |             READ(1030,*) (KC(J,I),J=1,14)
  35 |             READ(1040,*) (DC(J,I),J=1,14)
  36 |           END DO
  37 | 
  38 |           DO I=1,367
  39 |             READ(1010,*) (KH(J,I),J=1,14)
  40 |             READ(1020,*) (DH(J,I),J=1,14)
  41 |           END DO
  42 | 
  43 |           CLOSE(1010)
  44 |           CLOSE(1020)
  45 |           CLOSE(1030)
  46 |           CLOSE(1040)
  47 |           CLOSE(1050)
  48 |           CLOSE(1060)
  49 | 
  50 |         ELSEIF(nb_base.eq.'NBKM') THEN
  51 | 
  52 |           PRINT*, "  -> Using narrow band database NBKM from CETHIL"
  53 | 
  54 |           OPEN(UNIT=1010,FILE=trim(pathspec)//'/NBKMH2O')
  55 |           OPEN(UNIT=1020,FILE=trim(pathspec)//'/NBKMCO2')
  56 |           OPEN(UNIT=1030,FILE=trim(pathspec)//'/NBKMCO')
  57 | 
  58 |           DO I=1,48
  59 |             READ(1030,*) (KCO(J,I),J=1,14)
  60 |           END DO
  61 | 
  62 |           DO I=1,48
  63 |             READ(1030,*) (DCO(J,I),J=1,14)
  64 |           END DO
  65 | 
  66 |           DO I=1,96
  67 |             READ(1020,*) (KC(J,I),J=1,14)
  68 |           END DO
  69 | 
  70 |           DO I=1,96
  71 |             READ(1020,*) (DC(J,I),J=1,14)
  72 |           END DO
  73 | 
  74 |           DO I=1,367
  75 |             READ(1010,*) (KH(J,I),J=1,14)
  76 |           END DO
  77 | 
  78 |           DO I=1,367
  79 |             READ(1010,*) (DH(J,I),J=1,14)
  80 |           END DO
  81 | 
  82 |         ELSE
  83 |           PRINT*
  84 |           PRINT*, "ERROR : Unknown narrow band database"
  85 |           STOP
  86 |         ENDIF
  87 | 
  88 |        END SUBROUTINE PARAM