init_ck.f [SRC] [CPP] [JOB] [SCAN]
src



   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:
mcecile.f [src] - 415