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