This subroutine creates the required stencil.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| type(tem_stencilHeader_type), | intent(out) | :: | stencil |
stencil type to be defined |
||
| character(len=*) | :: | stencilKind |
stencil kind to decide which create function to call |
subroutine tem_create_stencil( stencil, stencilKind ) ! -------------------------------------------------------------------- ! !> stencil type to be defined type( tem_stencilHeader_type), intent(out) :: stencil !> stencil kind to decide which create function to call character(len=*) :: stencilKind ! -------------------------------------------------------------------- ! write(logUnit(1),*) 'Creating stencil kind: '//trim(stencilKind) select case (trim(stencilKind)) case ('d3q19') call init( me = stencil, & & QQN = 18, & & QQ = 19, & & useAll = .true., & & nDims = 3, & & label = 'd3q19', & & cxDir = d3q19_cxDir ) case ('d3q13') call init( me = stencil, & & QQN = 12, & & QQ = 13, & & useAll = .true., & & nDims = 3, & & label = 'd3q13', & & cxDir = d3q13_cxDir ) case ('d3q27') call init( me = stencil, & & QQN = 26, & & QQ = 27, & & useAll = .true., & & nDims = 3, & & label = 'd3q27', & & cxDir = d3q27_cxDir ) case ('d3q7') call init( me = stencil, & & QQN = 6, & & QQ = 7, & & useAll = .true., & & nDims = 3, & & label = 'd3q7', & & cxDir = d3q7_cxDir ) case ('d3q6', 'flekkoy') call init( me = stencil, & & QQN = 6, & & QQ = 6, & & useAll = .false., & & nDims = 3, & & label = 'd3q6', & & cxDir = d3q6_cxDir ) case ('d2q9') call init( me = stencil, & & QQN = 8, & & QQ = 9, & & useAll = .true., & & nDims = 2, & & label = 'd2q9', & & cxDir = d2q9_cxDir ) case ('d2q5') call init( me = stencil, & & QQN = 4, & & QQ = 5, & & useAll = .true., & & nDims = 2, & & label = 'd2q5', & & cxDir = d2q5_cxDir ) case ('d1q3') call init( me = stencil, & & QQN = 2, & & QQ = 3, & & useAll = .true., & & nDims = 1, & & label = 'd1q3', & & cxDir = d1q3_cxDir ) case ('default') write(logUnit(1),*) "Stencil kind Not defined yet!!" call tem_abort() end select end subroutine tem_create_stencil