1 c---------------------------------------------------------------------
2 c---------------------------------------------------------------------
4 subroutine exact_solution(xi,eta,zeta,dtemp)
6 c---------------------------------------------------------------------
7 c---------------------------------------------------------------------
9 c---------------------------------------------------------------------
10 c this function returns the exact solution at point xi, eta, zeta
11 c---------------------------------------------------------------------
15 double precision xi, eta, zeta, dtemp(5)
20 > xi*(ce(m,2) + xi*(ce(m,5) + xi*(ce(m,8) + xi*ce(m,11)))) +
21 > eta*(ce(m,3) + eta*(ce(m,6) + eta*(ce(m,9) + eta*ce(m,12))))+
22 > zeta*(ce(m,4) + zeta*(ce(m,7) + zeta*(ce(m,10) +