examples/smpi/NAS/LU/setbv.f

   1
   2 c---------------------------------------------------------------------
   3 c---------------------------------------------------------------------
   4
   5       subroutine setbv
   6
   7 c---------------------------------------------------------------------
   8 c---------------------------------------------------------------------
   9
  10 c---------------------------------------------------------------------
  11 c   set the boundary values of dependent variables
  12 c---------------------------------------------------------------------
  13
  14       implicit none
  15
  16       include 'applu.incl'
  17
  18 c---------------------------------------------------------------------
  19 c   local variables
  20 c---------------------------------------------------------------------
  21       integer i, j, k
  22       integer iglob, jglob
  23
  24 c---------------------------------------------------------------------
  25 c   set the dependent variable values along the top and bottom faces
  26 c---------------------------------------------------------------------
  27       do j = 1, ny
  28          jglob = jpt + j
  29          do i = 1, nx
  30            iglob = ipt + i
  31             call exact( iglob, jglob, 1, u( 1, i, j, 1 ) )
  32             call exact( iglob, jglob, nz, u( 1, i, j, nz ) )
  33          end do
  34       end do
  35
  36 c---------------------------------------------------------------------
  37 c   set the dependent variable values along north and south faces
  38 c---------------------------------------------------------------------
  39       IF (west.eq.-1) then
  40          do k = 1, nz
  41             do i = 1, nx
  42                iglob = ipt + i
  43                call exact( iglob, 1, k, u( 1, i, 1, k ) )
  44             end do
  45          end do
  46       END IF
  47
  48       IF (east.eq.-1) then
  49           do k = 1, nz
  50              do i = 1, nx
  51                 iglob = ipt + i
  52                 call exact( iglob, ny0, k, u( 1, i, ny, k ) )
  53              end do
  54           end do
  55       END IF
  56
  57 c---------------------------------------------------------------------
  58 c   set the dependent variable values along east and west faces
  59 c---------------------------------------------------------------------
  60       IF (north.eq.-1) then
  61          do k = 1, nz
  62             do j = 1, ny
  63                jglob = jpt + j
  64                call exact( 1, jglob, k, u( 1, 1, j, k ) )
  65             end do
  66          end do
  67       END IF
  68
  69       IF (south.eq.-1) then
  70          do k = 1, nz
  71             do j = 1, ny
  72                   jglob = jpt + j
  73             call exact( nx0, jglob, k, u( 1, nx, j, k ) )
  74             end do
  75          end do
  76       END IF
  77
  78       return
  79       end