1 | include(dom.inc)
2 |
3 | SUBROUTINE readtab(infile, kabs_all, nYH, nYC, nYCo, nT, nquad, &
4 | & model, dataref)
5 |
6 | IMPLICIT NONE
7 |
8 | CHARACTER*64 :: infile,model
9 | DOM_INT :: nYH,nYC,nYCO
10 | DOM_INT :: nT,nquad
11 | DOM_INT :: i,j,k,iT,iq
12 | DOM_REAL,DIMENSION(8) :: dataref
13 | DOM_REAL,DIMENSION(nYH,nYC,nYCo,nT,nquad) :: kabs_all
14 |
15 | ! ------------------!
16 | ! READ GLOBAL DATAS !
17 | ! ------------------!
18 |
19 | OPEN(12,FILE=infile,FORM='FORMATTED')
20 | READ(12,*)
21 | READ(12,*)
22 | READ(12,*)
23 | READ(12,*)
24 | READ(12,*)
25 |
26 | ! -----------------!
27 | ! READ TABLE DATA !
28 | ! -----------------!
29 |
30 | DO i=1,nYH
31 | DO j=1,nYC
32 | DO k=1,nYCO
33 | DO iT=1,nT
34 |
35 | READ(12,*) (kabs_all(i,j,k,iT,iq),iq=1,nquad)
36 |
37 | ENDDO
38 | ENDDO
39 | ENDDO
40 | ENDDO
41 |
42 | IF (model.eq.'TFSCK') READ(12,*) (dataref(iq),iq=1,8)
43 |
44 | CLOSE(12)
45 |
46 | END SUBROUTINE readtab
47 |
readtab.F could be called by: