c      cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c      cc          init2 routine                                  ccc
c      cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
       subroutine init2 ( rho , tem ,nn )
       implicit real*8 (a-h,o-z)
       real*8 logt , lrho
       dimension tem(nn), rho(nn)
       common /mika2/ tinit , tstep , rinit , rstep , nt , nr
       common /nika2/ a , z , ic 
c
        ic = nt*nr
            n = 0
       do 100 it = 1, nt 
                 ct = float( it )
               logt = tinit + ( ct - 1 )*tstep
          do 120  ir = 1, nr
                 cr = float( ir )
               lrho = rinit + ( cr - 1 )*rstep
                  n = n + 1
              if ( n .eq. nn ) call errcall(2)
             tem( n ) = 1.d1**logt
             rho( n ) = 1.d1**lrho
c            write(*,*) 'tem,rho (i) = ',tem(n),rho(n),n
 120      continue
 100   continue
        return
        end