flood_periphery_diagonal Subroutine

private subroutine flood_periphery_diagonal(proto, node_pos, treeID, qVal, iColor)

This routine checks for qVal of the periphery and floods if qVal < 0 and the node in that direction is fluid and not intersected by boundary

Arguments

Type IntentOptional Attributes Name
type(sdr_protoTree_type), intent(inout) :: proto

The proto tree description with all the data enabling the flooding.

integer, intent(in) :: node_pos

node position in protoTree

integer(kind=long_k), intent(in) :: treeID

treeID of current node

real(kind=rk) :: qVal(:)

qVal for all 26 neighbor directions, should be calculated already

integer, intent(in) :: iColor

Color to do the flooding in