1 | SUBROUTINE init_ck
2 | IMPLICIT NONE
3 |
4 | DOUBLE PRECISION phi
5 |
6 | INCLUDE 'tab_ck.inc'
7 |
8 | OPEN(10,FILE='tab_ck.fic')
9 | DO i_phitab=0,n_phitab
10 | atab(i_phitab,0)=0.D+0
11 | DO i_gtab=1,n_gtab
12 | READ(10,*) atab(i_phitab,i_gtab)
13 | ENDDO
14 | c READ(10,'(9(1x,e12.5))') phi,lphitab(i_phitab)
15 | c & ,atab(i_phitab,0)
16 | c & ,atab(i_phitab,1)
17 | c & ,atab(i_phitab,2)
18 | c & ,atab(i_phitab,3)
19 | c & ,atab(i_phitab,4)
20 | c & ,atab(i_phitab,5)
21 | c & ,atab(i_phitab,6)
22 | ENDDO
23 | CLOSE(10)
24 |
25 | OPEN(10,FILE='g_ck.fic')
26 | DO i_gtab=0,n_gtab
27 | READ(10,*) gtab(i_gtab),fgtab(i_gtab)
28 | ENDDO
29 | CLOSE(10)
30 |
31 | OPEN(10,FILE='lphi_ck.fic')
32 | DO i_phitab=0,n_phitab
33 | READ(10,*) phitab(i_phitab),lphitab(i_phitab)
34 | ENDDO
35 | CLOSE(10)
36 |
37 | RETURN
38 | END
init_ck.f could be called by: