1 | include(dom.inc)
   2 | 
   3 |        SUBROUTINE PARAM(KCO,KC,KH,DCO,DC,DH,pathspec,ntbands)
   4 | 
   5 |        include 'dom_constants.h'
   6 | 
   7 | !      IN  
   8 |        CHARACTER*80 :: pathspec
   9 |        DOM_INT      :: ntbands
  10 | 
  11 | !      LOCAL
  12 |        CHARACTER*80 :: kh2ofile, kco2file, kcofile
  13 |        CHARACTER*80 :: ph2ofile, pco2file, pcofile
  14 |        DOM_INT      :: I, J
  15 | 
  16 | !      OUT
  17 |        DOM_REAL,DIMENSION(14,ntbands)::KCO,DCO,KC,DC,KH,DH
  18 | 
  19 |        kh2ofile = pathspec(1:len_trim(pathspec))//'/SNBH2O_K'
  20 |        ph2ofile = pathspec(1:len_trim(pathspec))//'/SNBH2O_Phi'
  21 |        kco2file = pathspec(1:len_trim(pathspec))//'/SNBCO2_K'
  22 |        pco2file = pathspec(1:len_trim(pathspec))//'/SNBCO2_Phi'
  23 |        kcofile  = pathspec(1:len_trim(pathspec))//'/SNBCO_K'
  24 |        pcofile  = pathspec(1:len_trim(pathspec))//'/SNBCO_Phi'
  25 | 
  26 |        OPEN(UNIT=1010,FILE=kh2ofile)
  27 |        OPEN(UNIT=1020,FILE=ph2ofile)
  28 |        OPEN(UNIT=1030,FILE=kco2file)
  29 |        OPEN(UNIT=1040,FILE=pco2file)
  30 |        OPEN(UNIT=1050,FILE=kcofile)
  31 |        OPEN(UNIT=1060,FILE=pcofile)
  32 | 
  33 |        DO I=1,48
  34 |          READ(1050,*) (KCO(J,I),J=1,12)
  35 |          READ(1060,*) (DCO(J,I),J=1,12)
  36 |        END DO
  37 | 
  38 |        DO I=1,96
  39 |          READ(1030,*) (KC(J,I),J=1,14)
  40 |          READ(1040,*) (DC(J,I),J=1,14)
  41 |        END DO
  42 | 
  43 |        DO I=1,367
  44 |          READ(1010,*) (KH(J,I),J=1,14)
  45 |          READ(1020,*) (DH(J,I),J=1,14)
  46 |        END DO
  47 | 
  48 |        CLOSE(1010)
  49 |        CLOSE(1020)
  50 |        CLOSE(1030)
  51 |        CLOSE(1040)
  52 |        CLOSE(1050)
  53 |        CLOSE(1060)
  54 | 
  55 |        END SUBROUTINE PARAM