Nullify normals for 2D or 1D domains.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(block_object), | intent(inout) | :: | self | Block. |
elemental subroutine nullify_normals(self)
!< Nullify normals for 2D or 1D domains.
class(block_object), intent(inout) :: self !< Block.
if (self%signature%is_null_x) then
self%face_i%normal = (self%face_i%normal.paral.ex) + (0._R8P * ey ) + (0._R8P * ez )
self%face_j%normal = (0._R8P * ex ) + (self%face_j%normal.paral.ey) + (self%face_j%normal.paral.ez)
self%face_k%normal = (0._R8P * ex ) + (self%face_k%normal.paral.ey) + (self%face_k%normal.paral.ez)
endif
if (self%signature%is_null_y) then
self%face_i%normal = (self%face_i%normal.paral.ex) + (0._R8P * ey ) + (self%face_i%normal.paral.ez)
self%face_j%normal = (0._R8P * ex ) + (self%face_j%normal.paral.ey) + (0._R8P * ez )
self%face_k%normal = (self%face_k%normal.paral.ex) + (0._R8P * ey ) + (self%face_k%normal.paral.ez)
endif
if (self%signature%is_null_z) then
self%face_i%normal = (self%face_i%normal.paral.ex) + (self%face_i%normal.paral.ey) + (0._R8P * ez )
self%face_j%normal = (self%face_j%normal.paral.ex) + (self%face_j%normal.paral.ey) + (0._R8P * ez )
self%face_k%normal = (0._R8P * ex ) + (0._R8P * ey ) + (self%face_k%normal.paral.ez)
endif
endsubroutine nullify_normals