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

      subroutine setcoeff

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

      implicit none

      include 'header.h'

c---------------------------------------------------------------------
c  local variables
c---------------------------------------------------------------------


c---------------------------------------------------------------------
c   set up coefficients
c---------------------------------------------------------------------

c     The following three settings are based on average cell size,
c     not actual cell size

      dxi   = 1.0d0 / ( dble(gx_size)/dble(x_zones) - 1 )
      deta  = 1.0d0 / ( dble(gy_size)/dble(y_zones) - 1 )
      dzeta = 1.0d0 / ( dble(gz_size - 1) )

      tx1 = 1.0d0 / ( dxi * dxi )
      tx2 = 1.0d0 / ( 2.0d0 * dxi )
      tx3 = 1.0d0 / dxi

      ty1 = 1.0d0 / ( deta * deta )
      ty2 = 1.0d0 / ( 2.0d0 * deta )
      ty3 = 1.0d0 / deta

      tz1 = 1.0d0 / ( dzeta * dzeta )
      tz2 = 1.0d0 / ( 2.0d0 * dzeta )
      tz3 = 1.0d0 / dzeta

c---------------------------------------------------------------------
c   diffusion coefficients
c---------------------------------------------------------------------
      dx1 = 0.75d0
      dx2 = dx1
      dx3 = dx1
      dx4 = dx1
      dx5 = dx1

      dy1 = 0.75d0
      dy2 = dy1
      dy3 = dy1
      dy4 = dy1
      dy5 = dy1

      dz1 = 1.00d0
      dz2 = dz1
      dz3 = dz1
      dz4 = dz1
      dz5 = dz1

c---------------------------------------------------------------------
c   fourth difference dissipation
c---------------------------------------------------------------------
      dssp = ( max (dx1, dy1, dz1 ) ) / 4.0d0

c---------------------------------------------------------------------
c   coefficients of the exact solution to the first pde
c---------------------------------------------------------------------
      ce(1,1) = 2.0d0
      ce(1,2) = 0.0d0
      ce(1,3) = 0.0d0
      ce(1,4) = 4.0d0
      ce(1,5) = 5.0d0
      ce(1,6) = 3.0d0
      ce(1,7) = 5.0d-1
      ce(1,8) = 2.0d-2
      ce(1,9) = 1.0d-2
      ce(1,10) = 3.0d-2
      ce(1,11) = 5.0d-1
      ce(1,12) = 4.0d-1
      ce(1,13) = 3.0d-1

c---------------------------------------------------------------------
c   coefficients of the exact solution to the second pde
c---------------------------------------------------------------------
      ce(2,1) = 1.0d0
      ce(2,2) = 0.0d0
      ce(2,3) = 0.0d0
      ce(2,4) = 0.0d0
      ce(2,5) = 1.0d0
      ce(2,6) = 2.0d0
      ce(2,7) = 3.0d0
      ce(2,8) = 1.0d-2
      ce(2,9) = 3.0d-2
      ce(2,10) = 2.0d-2
      ce(2,11) = 4.0d-1
      ce(2,12) = 3.0d-1
      ce(2,13) = 5.0d-1

c---------------------------------------------------------------------
c   coefficients of the exact solution to the third pde
c---------------------------------------------------------------------
      ce(3,1) = 2.0d0
      ce(3,2) = 2.0d0
      ce(3,3) = 0.0d0
      ce(3,4) = 0.0d0
      ce(3,5) = 0.0d0
      ce(3,6) = 2.0d0
      ce(3,7) = 3.0d0
      ce(3,8) = 4.0d-2
      ce(3,9) = 3.0d-2
      ce(3,10) = 5.0d-2
      ce(3,11) = 3.0d-1
      ce(3,12) = 5.0d-1
      ce(3,13) = 4.0d-1

c---------------------------------------------------------------------
c   coefficients of the exact solution to the fourth pde
c---------------------------------------------------------------------
      ce(4,1) = 2.0d0
      ce(4,2) = 2.0d0
      ce(4,3) = 0.0d0
      ce(4,4) = 0.0d0
      ce(4,5) = 0.0d0
      ce(4,6) = 2.0d0
      ce(4,7) = 3.0d0
      ce(4,8) = 3.0d-2
      ce(4,9) = 5.0d-2
      ce(4,10) = 4.0d-2
      ce(4,11) = 2.0d-1
      ce(4,12) = 1.0d-1
      ce(4,13) = 3.0d-1

c---------------------------------------------------------------------
c   coefficients of the exact solution to the fifth pde
c---------------------------------------------------------------------
      ce(5,1) = 5.0d0
      ce(5,2) = 4.0d0
      ce(5,3) = 3.0d0
      ce(5,4) = 2.0d0
      ce(5,5) = 1.0d-1
      ce(5,6) = 4.0d-1
      ce(5,7) = 3.0d-1
      ce(5,8) = 5.0d-2
      ce(5,9) = 4.0d-2
      ce(5,10) = 3.0d-2
      ce(5,11) = 1.0d-1
      ce(5,12) = 3.0d-1
      ce(5,13) = 2.0d-1

      return
      end