This routine find the treeID on the opposite side neighbor of the periodic plane for current leaf node
write(dbgUnit(5),*) '1.1 qVal ', qVal
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=rk), | intent(in) | :: | elemBary(3) |
current element barycenter |
||
| integer, | intent(in) | :: | iDir |
Current boundary neighbor direction |
||
| integer(kind=long_k), | intent(inout) | :: | bc_ID |
treeiD of opposite neighbor with periodic plane |
||
| real(kind=rk), | intent(out) | :: | qVal |
distance from boundary for direction iDir Todo HK: maybe turn this optional, why is it relevant for periodic? |
||
| logical, | intent(inout) | :: | unKnownBnd |
Set to true if a neighbor with no property is encountered |
||
| integer, | intent(in) | :: | neighbor_pos |
position of neighbor treeID in the proto%node list |
||
| integer, | intent(in) | :: | neighbor_level |
level of the periodic boundary neighbor node |
||
| type(levelValues_type), | intent(in) | :: | leVal |
level value of parent node |
||
| type(sdr_protoTree_type), | intent(in) | :: | proto |
preliminary tree |
||
| type(sdr_geometry_type), | intent(in) | :: | geometry |
contains all geometrical objects |
||
| type(treelmesh_type), | intent(in) | :: | meshUniverse |
contains bounding cube information |