+++ /dev/null
-
-c---------------------------------------------------------------------
-c---------------------------------------------------------------------
-
- subroutine setcoeff
-
-c---------------------------------------------------------------------
-c---------------------------------------------------------------------
-
- implicit none
-
- include 'applu.incl'
-
-c---------------------------------------------------------------------
-c local variables
-c---------------------------------------------------------------------
-
-
-c---------------------------------------------------------------------
-c set up coefficients
-c---------------------------------------------------------------------
- dxi = 1.0d+00 / ( nx0 - 1 )
- deta = 1.0d+00 / ( ny0 - 1 )
- dzeta = 1.0d+00 / ( nz0 - 1 )
-
- tx1 = 1.0d+00 / ( dxi * dxi )
- tx2 = 1.0d+00 / ( 2.0d+00 * dxi )
- tx3 = 1.0d+00 / dxi
-
- ty1 = 1.0d+00 / ( deta * deta )
- ty2 = 1.0d+00 / ( 2.0d+00 * deta )
- ty3 = 1.0d+00 / deta
-
- tz1 = 1.0d+00 / ( dzeta * dzeta )
- tz2 = 1.0d+00 / ( 2.0d+00 * dzeta )
- tz3 = 1.0d+00 / dzeta
-
- ii1 = 2
- ii2 = nx0 - 1
- ji1 = 2
- ji2 = ny0 - 2
- ki1 = 3
- ki2 = nz0 - 1
-
-c---------------------------------------------------------------------
-c diffusion coefficients
-c---------------------------------------------------------------------
- dx1 = 0.75d+00
- dx2 = dx1
- dx3 = dx1
- dx4 = dx1
- dx5 = dx1
-
- dy1 = 0.75d+00
- dy2 = dy1
- dy3 = dy1
- dy4 = dy1
- dy5 = dy1
-
- dz1 = 1.00d+00
- dz2 = dz1
- dz3 = dz1
- dz4 = dz1
- dz5 = dz1
-
-c---------------------------------------------------------------------
-c fourth difference dissipation
-c---------------------------------------------------------------------
- dssp = ( max (dx1, dy1, dz1 ) ) / 4.0d+00
-
-c---------------------------------------------------------------------
-c coefficients of the exact solution to the first pde
-c---------------------------------------------------------------------
- ce(1,1) = 2.0d+00
- ce(1,2) = 0.0d+00
- ce(1,3) = 0.0d+00
- ce(1,4) = 4.0d+00
- ce(1,5) = 5.0d+00
- ce(1,6) = 3.0d+00
- ce(1,7) = 5.0d-01
- ce(1,8) = 2.0d-02
- ce(1,9) = 1.0d-02
- ce(1,10) = 3.0d-02
- ce(1,11) = 5.0d-01
- ce(1,12) = 4.0d-01
- ce(1,13) = 3.0d-01
-
-c---------------------------------------------------------------------
-c coefficients of the exact solution to the second pde
-c---------------------------------------------------------------------
- ce(2,1) = 1.0d+00
- ce(2,2) = 0.0d+00
- ce(2,3) = 0.0d+00
- ce(2,4) = 0.0d+00
- ce(2,5) = 1.0d+00
- ce(2,6) = 2.0d+00
- ce(2,7) = 3.0d+00
- ce(2,8) = 1.0d-02
- ce(2,9) = 3.0d-02
- ce(2,10) = 2.0d-02
- ce(2,11) = 4.0d-01
- ce(2,12) = 3.0d-01
- ce(2,13) = 5.0d-01
-
-c---------------------------------------------------------------------
-c coefficients of the exact solution to the third pde
-c---------------------------------------------------------------------
- ce(3,1) = 2.0d+00
- ce(3,2) = 2.0d+00
- ce(3,3) = 0.0d+00
- ce(3,4) = 0.0d+00
- ce(3,5) = 0.0d+00
- ce(3,6) = 2.0d+00
- ce(3,7) = 3.0d+00
- ce(3,8) = 4.0d-02
- ce(3,9) = 3.0d-02
- ce(3,10) = 5.0d-02
- ce(3,11) = 3.0d-01
- ce(3,12) = 5.0d-01
- ce(3,13) = 4.0d-01
-
-c---------------------------------------------------------------------
-c coefficients of the exact solution to the fourth pde
-c---------------------------------------------------------------------
- ce(4,1) = 2.0d+00
- ce(4,2) = 2.0d+00
- ce(4,3) = 0.0d+00
- ce(4,4) = 0.0d+00
- ce(4,5) = 0.0d+00
- ce(4,6) = 2.0d+00
- ce(4,7) = 3.0d+00
- ce(4,8) = 3.0d-02
- ce(4,9) = 5.0d-02
- ce(4,10) = 4.0d-02
- ce(4,11) = 2.0d-01
- ce(4,12) = 1.0d-01
- ce(4,13) = 3.0d-01
-
-c---------------------------------------------------------------------
-c coefficients of the exact solution to the fifth pde
-c---------------------------------------------------------------------
- ce(5,1) = 5.0d+00
- ce(5,2) = 4.0d+00
- ce(5,3) = 3.0d+00
- ce(5,4) = 2.0d+00
- ce(5,5) = 1.0d-01
- ce(5,6) = 4.0d-01
- ce(5,7) = 3.0d-01
- ce(5,8) = 5.0d-02
- ce(5,9) = 4.0d-02
- ce(5,10) = 3.0d-02
- ce(5,11) = 1.0d-01
- ce(5,12) = 3.0d-01
- ce(5,13) = 2.0d-01
-
- return
- end
-
-
