Geometric information and definitions
| Type | Visibility | Attributes | Name | Initial | |||
|---|---|---|---|---|---|---|---|
| type(treelmesh_type), | public | :: | tree |
tree data type |
|||
| type(tem_BC_prop_type), | public | :: | boundary |
boundary information as stored on disk |
|||
| type(mus_geomIncrHead_type), | public, | allocatable | :: | geomIncr(:) |
The header type containing all the geometry increase information |
||
| logical, | public | :: | dynamicGeom | = | .false. |
Logical to define whether geometry increase is active or not |
|
| integer, | public, | allocatable | :: | posInBndID(:) |
Tree element position in the boundary_ID( nDir, nElems) in bc_prop_type it has a size of tree%nElems How to use: do iElem = 1, tree%nElems posInBndID = posInBndID( iElem ) ! current element has boundary only if posInBndID>0 ! else posInBndID = -1 if (posInBnd > 0 ) bnd_ID(1:nDir) = bc_prop%boundary_ID( 1:nDir, posInBndID ) end if end do |
||
| integer, | public, | allocatable | :: | posInQVal(:) |
Tree element position in the qVal( nDir, nElems) in bc_prop_type it has a size of tree%nElems How to use: do iElem = 1, tree%nElems posInQVal = posInQVal( iElem ) ! current element has qVal if posInQVal>0 else posInQVal = -1 if (posInQVal > 0 ) qVal(1:nDir) = bc_prop%qVal( 1:nDir, posInQVal ) end if end do |
||
| integer, | public, | allocatable | :: | levelPointer(:) |
tree element position in level descriptor total list it has a size of tree%nElems How to use: do iElem = 1, tree%nElems treeID = tree%treeID( iElem ) level = tem_levelOf( treeID ) posInTotal = levelPointer( iElem ) treeID = LevelDesc( iLevel )%total( posInTotal ) end do |
||
| integer, | public, | allocatable | :: | bcLevelPointer(:) |
Boundary element poisition in the levelwise globBC%elemLvl(:)%elem%val It has a size of geometry%boundary%property%nElems. It is used in tracking to extract value stored in boundary types. Hot to use this access normal direction of boundary element: do iElem = 1, tree%nElems level = tem_levelOf( treeID ) posInBndID = posInBndID(iElem) if (posInBndID > 0) then BCIDs = bc_prop%boundary_ID(:, posInBndID) minBCID = minval(BCIDs, BCIDs > 0) posInBcElem = bcLevelPointer(posInBndID) normal = globBC%elemLvl(iLevel)%normal(posInBcElem) end if end do |
||
| integer, | public, | allocatable | :: | minBcID(:) |
Minimum bcID for each boundary element. if a element has more than one boundary then use minBcID which depends on boundary order in seeder configuration. |
||
| type(mus_IBM_globType), | public | :: | globIBM |
immersed boundary data |
|||
| real(kind=rk), | public, | allocatable | :: | bndForce(:,:) |
Contains Forces on boundary elements computed using momentum exchange method. This will be used to derive_bndForce routine to compute force of certain boundaries. Forces are stored in level-independent fashion as geometry%boundaryID loaded from mesh files. Dim1: geometry%boundaryi%property%nElems Dim2: 3 |
||
| real(kind=rk), | public, | allocatable | :: | bndMoment(:,:) |
Contains Moments on boundary elements computed using momentum exchange method. This will be used to derive_bndMoment routine to compute moment of certain boundaries. Forces are stored in level-independent fashion as geometry%boundaryID loaded from mesh files. Dim1: geometry%boundaryi%property%nElems Dim2: 3 |