c      ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c      cc     Initializing                                         ccc
c      cc                  init.f                                  ccc
c      ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
       subroutine init( temp, rhom )
        implicit real*8 ( a-h, o-z )
        integer nt, nr
        dimension temp( 1 ), rhom( 1 )
        common /param/ ndat, n
c
        open ( unit = 80, file = 'init.par' , status = 'unknown'
     &       , access = 'sequential', form = 'formatted'        )
        read(80,2000) tinit
        read(80,2000) tstep
        read(80,2010) nt
        read(80,2000) rinit
        read(80,2000) rstep
        read(80,2010) nr
        read(80,2010) n      
        close ( 80 )
 2000   format( 10x,f8.0)
 2010   format( 10x,i8)
c
        ndat = 0
        do 100    it=1, nt
           do 200 ir=1, nr
             ndat = ndat + 1
              temp ( ndat ) = 1.d1**( tinit + (it - 1)*tstep )
              rhom ( ndat ) = 1.d1**( rinit + (ir - 1)*rstep )
 200       continue
 100    continue
        return
c
        end