program create2

c        Reads in plot3d grid and solution files and generates the  
c        restart file for nparc3d. 
c        
c-----------------------------------------------------------------------

      parameter ( nimax =  162 )
      parameter ( njmax =  121 )
      parameter ( nkmax =   71 )

      dimension q  (nimax,njmax,nkmax,5)
      dimension x  (nimax,njmax,nkmax)
      dimension y  (nimax,njmax,nkmax)
      dimension z  (nimax,njmax,nkmax)

c-----------------------------------------------------------------------

c...Set some constants

      gamma  = 1.4

      gp1    = gamma + 1.0
      gm1    = gamma - 1.0

c...Read in the grid.

      open ( unit=7, file='run.x', form='unformatted' )
      read (7) nblk
      read (7) ni, nj, nk
      read (7) ((( x(i,j,k), i=1,ni), j=1,nj), k=1,nk),
     &         ((( y(i,j,k), i=1,ni), j=1,nj), k=1,nk),
     &         ((( z(i,j,k), i=1,ni), j=1,nj), k=1,nk)

c...Read in the solution.

      open ( unit=8, file='run.q', form='unformatted' )
      read (8) nblk
      read (8) ni, nj, nk
      read (8) t1, t2, t3, t4
      read (8) (((( q(i,j,k,m), i=1,ni), j=1,nj), k=1,nk), m=1,5)

c...Nondimensionalize x and q for output to nparc

      rlen =1.0

      do  i = 1, ni
        do  j = 1, nj
          do  k = 1, nk
            x(i,j,k)   = x(i,j,k) / rlen
            y(i,j,k)   = y(i,j,k) / rlen
            z(i,j,k)   = z(i,j,k) / rlen
            q(i,j,k,1) = q(i,j,k,1)
            q(i,j,k,2) = q(i,j,k,2)
            q(i,j,k,3) = q(i,j,k,3)
            q(i,j,k,4) = q(i,j,k,5) 
          enddo
        enddo
      enddo

c...Output the nparc3d restart file (unformatted, fort.2)

      open ( 2, file='fort.2', form='unformatted' )

      nc1 = 0

      write ( 2 ) nc1
      write ( 2 ) ni, nj
      write ( 2 ) ((( x(i,j,k), i=1,ni), j=1,nj), k=1,nk),
     &            ((( y(i,j,k), i=1,ni), j=1,nj), k=1,nk)
      write ( 2 ) (((( q(i,j,k,m), i=1,ni), j=1,nj), k=1,nk), m=1,4 )

c-----------------------------------------------------------------------

      stop
      end