2 c---------------------------------------------------------------------
3 c---------------------------------------------------------------------
6 c---------------------------------------------------------------------
7 c---------------------------------------------------------------------
9 c---------------------------------------------------------------------
11 c set the initial values of independent variables based on tri-linear
12 c interpolation of boundary values in the computational space.
14 c---------------------------------------------------------------------
20 c---------------------------------------------------------------------
22 c---------------------------------------------------------------------
25 double precision xi, eta, zeta
26 double precision pxi, peta, pzeta
27 double precision ue_1jk(5),ue_nx0jk(5),ue_i1k(5),
28 > ue_iny0k(5),ue_ij1(5),ue_ijnz(5)
32 zeta = ( dble (k-1) ) / (nz-1)
35 IF (jglob.ne.1.and.jglob.ne.ny0) then
36 eta = ( dble (jglob-1) ) / (ny0-1)
39 IF (iglob.ne.1.and.iglob.ne.nx0) then
40 xi = ( dble (iglob-1) ) / (nx0-1)
41 call exact (1,jglob,k,ue_1jk)
42 call exact (nx0,jglob,k,ue_nx0jk)
43 call exact (iglob,1,k,ue_i1k)
44 call exact (iglob,ny0,k,ue_iny0k)
45 call exact (iglob,jglob,1,ue_ij1)
46 call exact (iglob,jglob,nz,ue_ijnz)
48 pxi = ( 1.0d+00 - xi ) * ue_1jk(m)
50 peta = ( 1.0d+00 - eta ) * ue_i1k(m)
52 pzeta = ( 1.0d+00 - zeta ) * ue_ij1(m)
55 u( m, i, j, k ) = pxi + peta + pzeta
56 > - pxi * peta - peta * pzeta - pzeta * pxi
57 > + pxi * peta * pzeta