diff options
Diffstat (limited to 'src/Dissipation_4_2.F90')
| -rw-r--r-- | src/Dissipation_4_2.F90 | 18 |
1 files changed, 18 insertions, 0 deletions
diff --git a/src/Dissipation_4_2.F90 b/src/Dissipation_4_2.F90 index 5b419c9..84746e7 100644 --- a/src/Dissipation_4_2.F90 +++ b/src/Dissipation_4_2.F90 @@ -35,6 +35,7 @@ subroutine dissipation_4_2 (var, ni, nj, nk, bb, gsize, delta, epsilon, rhs) if ( bb(1) == 0 ) then il = 1 + gsize(1) else +!$omp parallel workshare rhs(1,:,:) = rhs(1,:,:) + & ( a(1,1) * var(1,:,:) + a(2,1) * var(2,:,:) + & a(3,1) * var(3,:,:) ) * idel @@ -49,12 +50,14 @@ subroutine dissipation_4_2 (var, ni, nj, nk, bb, gsize, delta, epsilon, rhs) ( a(2,4) * var(2,:,:) + a(3,4) * var(3,:,:) + & a(4,4) * var(4,:,:) + a(5,4) * var(5,:,:) + & a(6,4) * var(6,:,:) ) * idel +!$omp end parallel workshare il = 5 end if if ( bb(2) == 0 ) then ir = ni - gsize(1) else +!$omp parallel workshare rhs(ni-3,:,:) = rhs(ni-3,:,:) + & ( a(2,4) * var(ni-1,:,:) + a(3,4) * var(ni-2,:,:) + & a(4,4) * var(ni-3,:,:) + a(5,4) * var(ni-4,:,:) + & @@ -69,15 +72,18 @@ subroutine dissipation_4_2 (var, ni, nj, nk, bb, gsize, delta, epsilon, rhs) rhs(ni,:,:) = rhs(ni,:,:) + & ( a(1,1) * var(ni,:,:) + a(2,1) * var(ni-1,:,:) + & a(3,1) * var(ni-2,:,:) ) * idel +!$omp end parallel workshare ir = ni - 4 end if +!$omp parallel workshare rhs(il:ir,:,:) = rhs(il:ir,:,:) + & ( -6.0_wp * var(il:ir,:,:) + & 4.0_wp * ( var(il-1:ir-1,:,:) + & var(il+1:ir+1,:,:) ) - & ( var(il-2:ir-2,:,:) + & var(il+2:ir+2,:,:) ) ) * idel +!$omp end parallel workshare if ( zero_derivs_y == 0 ) then call set_coeff ( a ) @@ -91,6 +97,7 @@ subroutine dissipation_4_2 (var, ni, nj, nk, bb, gsize, delta, epsilon, rhs) if ( bb(3) == 0 ) then jl = 1 + gsize(2) else +!$omp parallel workshare rhs(:,1,:) = rhs(:,1,:) + & ( a(1,1) * var(:,1,:) + a(2,1) * var(:,2,:) + & a(3,1) * var(:,3,:) ) * idel @@ -105,12 +112,14 @@ subroutine dissipation_4_2 (var, ni, nj, nk, bb, gsize, delta, epsilon, rhs) ( a(2,4) * var(:,2,:) + a(3,4) * var(:,3,:) + & a(4,4) * var(:,4,:) + a(5,4) * var(:,5,:) + & a(6,4) * var(:,6,:) ) * idel +!$omp end parallel workshare jl = 5 end if if ( bb(4) == 0 ) then jr = nj - gsize(2) else +!$omp parallel workshare rhs(:,nj-3,:) = rhs(:,nj-3,:) + & ( a(2,4) * var(:,nj-1,:) + a(3,4) * var(:,nj-2,:) + & a(4,4) * var(:,nj-3,:) + a(5,4) * var(:,nj-4,:) + & @@ -125,15 +134,18 @@ subroutine dissipation_4_2 (var, ni, nj, nk, bb, gsize, delta, epsilon, rhs) rhs(:,nj,:) = rhs(:,nj,:) + & ( a(1,1) * var(:,nj,:) + a(2,1) * var(:,nj-1,:) + & a(3,1) * var(:,nj-2,:) ) * idel +!$omp end parallel workshare jr = nj - 4 end if +!$omp parallel workshare rhs(:,jl:jr,:) = rhs(:,jl:jr,:) + & ( -6.0_wp * var(:,jl:jr,:) + & 4.0_wp * ( var(:,jl-1:jr-1,:) + & var(:,jl+1:jr+1,:) ) - & ( var(:,jl-2:jr-2,:) + & var(:,jl+2:jr+2,:) ) ) * idel +!$omp end parallel workshare end if if ( zero_derivs_z == 0 ) then @@ -148,6 +160,7 @@ subroutine dissipation_4_2 (var, ni, nj, nk, bb, gsize, delta, epsilon, rhs) if ( bb(5) == 0 ) then kl = 1 + gsize(3) else +!$omp parallel workshare rhs(:,:,1) = rhs(:,:,1) + & ( a(1,1) * var(:,:,1) + a(2,1) * var(:,:,2) + & a(3,1) * var(:,:,3) ) * idel @@ -162,12 +175,14 @@ subroutine dissipation_4_2 (var, ni, nj, nk, bb, gsize, delta, epsilon, rhs) ( a(2,4) * var(:,:,2) + a(3,4) * var(:,:,3) + & a(4,4) * var(:,:,4) + a(5,4) * var(:,:,5) + & a(6,4) * var(:,:,6) ) * idel +!$omp end parallel workshare kl = 5 end if if ( bb(6) == 0 ) then kr = nk - gsize(3) else +!$omp parallel workshare rhs(:,:,nk-3) = rhs(:,:,nk-3) + & ( a(2,4) * var(:,:,nk-1) + a(3,4) * var(:,:,nk-2) + & a(4,4) * var(:,:,nk-3) + a(5,4) * var(:,:,nk-4) + & @@ -182,15 +197,18 @@ subroutine dissipation_4_2 (var, ni, nj, nk, bb, gsize, delta, epsilon, rhs) rhs(:,:,nk) = rhs(:,:,nk) + & ( a(1,1) * var(:,:,nk) + a(2,1) * var(:,:,nk-1) + & a(3,1) * var(:,:,nk-2) ) * idel +!$omp end parallel workshare kr = nk - 4 end if +!$omp parallel workshare rhs(:,:,kl:kr) = rhs(:,:,kl:kr) + & ( -6.0_wp * var(:,:,kl:kr) + & 4.0_wp * ( var(:,:,kl-1:kr-1) + & var(:,:,kl+1:kr+1) ) - & ( var(:,:,kl-2:kr-2) + & var(:,:,kl+2:kr+2) ) ) * idel +!$omp end parallel workshare end if contains |