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