0d0 enddo dir 1 buffid 2 dir bufflen 0 if axis eq

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: uff(i,4) = buff(i,3) buff(i,2) = buff(i,1) enddo subroutine take3( axis, dir, u, n1, n2, n3 ) use caf_intrinsics integer axis, dir, n1, n2, n3 double precision u( n1, n2, n3 ) dir = -1 buff_id = 2 + dir buff_len = 0 implicit none integer i3, i2, i1, buff_len,buff_id integer i, kk, indx buff_id = 3 + dir indx = 0 if( axis .eq. 1 )then if( dir .eq. -1 )then include 'cafnpb.h' include 'globals.h' dir = -1 if( axis .eq. 1 )then do i3=2,n3-1 do i2=2,n2-1 indx = indx + 1 u(n1,i2,i3) = buff(indx, buff_id ) enddo enddo endif integer i3, i2, i1, buff_len,buff_id do i3=2,n3-1 do i2=2,n2-1 buff_len = buff_len + 1 buff(buff_len,buff_id ) = u( 2, enddo enddo > integer axis, dir, n1, n2, n3 double precision u( n1, n2, n3 ) i2,i3) integer i3, i2, i1 buff(1:buff_len,buff_id+1)[nbr(axis,dir,k)] = buff(1:buff_len,buff_id) else if( dir .eq. +1 ) then do i3=2,n3-1 do i2=2,n2-1 buff_len = buff_len + 1 buff(buff_len, buff_id ) = u( n1-1, i2,i3) enddo enddo > buff(1:buff_len,buff_id+1)[nbr(axis,dir,k)] = buff(1:buff_len,buff_i...
View Full Document

Ask a homework question - tutors are online