diff options
Diffstat (limited to 'src/include/centered_second2.h')
| -rw-r--r-- | src/include/centered_second2.h | 56 |
1 files changed, 28 insertions, 28 deletions
diff --git a/src/include/centered_second2.h b/src/include/centered_second2.h index c29fd41..e7aca4b 100644 --- a/src/include/centered_second2.h +++ b/src/include/centered_second2.h @@ -1,39 +1,39 @@ ! Calculation of centered differences. ! $Header$ -if ( eh_mask(i,j,k) .ge. 0 ) then - if ( ( .not. btest(eh_mask(i,j,k),0) ) .and. & - ( .not. btest(eh_mask(i,j,k),1) ) ) then - dfx(i,j,k) = idx * ( f(i+1,j,k) - f(i-1,j,k) ) - else if ( btest(eh_mask(i,j,k),0) ) then - dfx(i,j,k) = idx * ( -three * f(i,j,k) + & - four * f(i+1,j,k) - f(i+2,j,k) ) +if ( eh_mask(i,j,k,l) .ge. 0 ) then + if ( ( .not. btest(eh_mask(i,j,k,l),0) ) .and. & + ( .not. btest(eh_mask(i,j,k,l),1) ) ) then + dfx(i,j,k,l) = idx * ( f(i+1,j,k,l) - f(i-1,j,k,l) ) + else if ( btest(eh_mask(i,j,k,l),0) ) then + dfx(i,j,k,l) = idx * ( -three * f(i,j,k,l) + & + four * f(i+1,j,k,l) - f(i+2,j,k,l) ) else - dfx(i,j,k) = idx * ( three * f(i,j,k) - & - four * f(i-1,j,k) + f(i-2,j,k) ) + dfx(i,j,k,l) = idx * ( three * f(i,j,k,l) - & + four * f(i-1,j,k,l) + f(i-2,j,k,l) ) end if - if ( ( .not. btest(eh_mask(i,j,k),2) ) .and. & - ( .not. btest(eh_mask(i,j,k),3) ) ) then - dfy(i,j,k) = idy * ( f(i,j+1,k) - f(i,j-1,k) ) - else if ( btest(eh_mask(i,j,k),2) ) then - dfy(i,j,k) = idy * ( -three * f(i,j,k) + & - four * f(i,j+1,k) - f(i,j+2,k) ) + if ( ( .not. btest(eh_mask(i,j,k,l),2) ) .and. & + ( .not. btest(eh_mask(i,j,k,l),3) ) ) then + dfy(i,j,k,l) = idy * ( f(i,j+1,k,l) - f(i,j-1,k,l) ) + else if ( btest(eh_mask(i,j,k,l),2) ) then + dfy(i,j,k,l) = idy * ( -three * f(i,j,k,l) + & + four * f(i,j+1,k,l) - f(i,j+2,k,l) ) else - dfy(i,j,k) = idy * ( three * f(i,j,k) - & - four * f(i,j-1,k) + f(i,j-2,k) ) + dfy(i,j,k,l) = idy * ( three * f(i,j,k,l) - & + four * f(i,j-1,k,l) + f(i,j-2,k,l) ) end if - if ( ( .not. btest(eh_mask(i,j,k),4) ) .and. & - ( .not. btest(eh_mask(i,j,k),5) ) ) then - dfz(i,j,k) = idz * ( f(i,j,k+1) - f(i,j,k-1) ) - else if ( btest(eh_mask(i,j,k),4) ) then - dfz(i,j,k) = idz * ( -three * f(i,j,k) + & - four * f(i,j,k+1) - f(i,j,k+2) ) + if ( ( .not. btest(eh_mask(i,j,k,l),4) ) .and. & + ( .not. btest(eh_mask(i,j,k,l),5) ) ) then + dfz(i,j,k,l) = idz * ( f(i,j,k+1,l) - f(i,j,k-1,l) ) + else if ( btest(eh_mask(i,j,k,l),4) ) then + dfz(i,j,k,l) = idz * ( -three * f(i,j,k,l) + & + four * f(i,j,k+1,l) - f(i,j,k+2,l) ) else - dfz(i,j,k) = idz * ( three * f(i,j,k) - & - four * f(i,j,k-1) + f(i,j,k-2) ) + dfz(i,j,k,l) = idz * ( three * f(i,j,k,l) - & + four * f(i,j,k-1,l) + f(i,j,k-2,l) ) end if else - dfx(i,j,k) = zero - dfy(i,j,k) = zero - dfz(i,j,k) = zero + dfx(i,j,k,l) = zero + dfy(i,j,k,l) = zero + dfz(i,j,k,l) = zero end if |