Rocstar  1.0
Rocstar multiphysics simulation application
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
TURB_floLesUniFiltCCI.F90
Go to the documentation of this file.
1 ! *********************************************************************
2 ! * Rocstar Simulation Suite *
3 ! * Copyright@2015, Illinois Rocstar LLC. All rights reserved. *
4 ! * *
5 ! * Illinois Rocstar LLC *
6 ! * Champaign, IL *
7 ! * www.illinoisrocstar.com *
8 ! * sales@illinoisrocstar.com *
9 ! * *
10 ! * License: See LICENSE file in top level of distribution package or *
11 ! * http://opensource.org/licenses/NCSA *
12 ! *********************************************************************
13 ! *********************************************************************
14 ! * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, *
15 ! * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES *
16 ! * OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND *
17 ! * NONINFRINGEMENT. IN NO EVENT SHALL THE CONTRIBUTORS OR *
18 ! * COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER *
19 ! * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, *
20 ! * Arising FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE *
21 ! * USE OR OTHER DEALINGS WITH THE SOFTWARE. *
22 ! *********************************************************************
23 !******************************************************************************
24 !
25 ! Purpose: Perform cell to cell uniform filtering in I-direction.
26 !
27 ! Description: The filtering is performed using fixed filter coefficients.
28 ! If filterwidth is zero in this direction the cell-variable is
29 ! simply copied to the filtered cell-variable.
30 !
31 ! Input: global = global data, only used in register function
32 ! nDum = number of dummy cells
33 ! i,j,kbeg, i,j,kend = cell index boundaries
34 ! iCOff, ijCOff = cell Offset
35 ! nDel = three components filter width parameter
36 ! idBeg = begin variable index to be filtered
37 ! idEnd = end variable index to be filtered
38 ! fVar = cell variable to be filtered
39 !
40 ! Output: filtVar = filtered cell variable
41 !
42 ! Notes: Mother routine = LesUniFiltCC.
43 !
44 !******************************************************************************
45 !
46 ! $Id: TURB_floLesUniFiltCCI.F90,v 1.5 2008/12/06 08:44:43 mtcampbe Exp $
47 !
48 ! Copyright: (c) 2001 by the University of Illinois
49 !
50 !******************************************************************************
51 
52 SUBROUTINE turb_flolesunifiltcci( global,nDum,ibeg,iend,jbeg,jend,kbeg,kend, &
53  icoff,ijcoff,ndel,idbeg,idend,fact1,fact2,fvar, &
54  filtvar )
55 
56  USE moddatatypes
57  USE modglobal, ONLY : t_global
58  USE modturbulence
59  USE moderror
60  USE turb_modparameters
61  IMPLICIT NONE
62 
63 #include "Indexing.h"
64 
65 ! ... parameters
66  TYPE(t_global), POINTER :: global
67  INTEGER :: ndum,ibeg,iend,jbeg,jend,kbeg,kend,icoff,ijcoff
68  INTEGER :: ndel(diri:dirk),idbeg,idend
69  REAL(RFREAL) :: fact1(filwidth_four),fact2(filwidth_four)
70  REAL(RFREAL), POINTER :: fvar(:,:),filtvar(:,:)
71 
72 ! ... loop variables
73  INTEGER :: i, j, k, l, ijkc, ijkcb, ijkce, ijkcb1, ijkce1
74 
75 ! ... local variables
76  CHARACTER(CHRLEN) :: rcsidentstring
77  REAL(RFREAL) :: tmpb(idbeg:idend,2),tmpe(idbeg:idend,2)
78 
79 !******************************************************************************
80 
81  rcsidentstring = '$RCSfile: TURB_floLesUniFiltCCI.F90,v $'
82 
83  CALL registerfunction( global,'TURB_FloLesUniFiltCCI',&
84  'TURB_floLesUniFiltCCI.F90' )
85 
86 ! integration over I-direction -----------------------------------------------
87 
88  IF (ndel(diri)==filwidth_zero) THEN
89 
90 ! - no filtering in I-direction is needed and we copy fVar into filtVar
91 
92  filtvar(idbeg:idend,:)=fvar(idbeg:idend,:) ! no simple copy filtvar=fvar
93  ! array sizes can be different
94 
95  ELSEIF ((ndel(diri)==filwidth_one) .OR. &
96  (ndel(diri)==filwidth_two)) THEN
97 
98  DO k=kbeg,kend
99  DO j=jbeg,jend
100 
101 ! ----- perform integration
102 
103  DO i=ibeg+1,iend-1
104  ijkc = indijk(i ,j ,k ,icoff,ijcoff)
105  DO l=idbeg,idend
106  filtvar(l,ijkc)=fact1(ndel(diri))*(fvar(l,ijkc-1)+fvar(l,ijkc+1))+ &
107  fact2(ndel(diri))*fvar(l,ijkc)
108  END DO ! l
109  ENDDO ! i
110 
111 ! ----- perform integration at the edges
112 
113  ijkcb = indijk(ibeg ,j ,k ,icoff,ijcoff)
114  ijkce = indijk(iend ,j ,k ,icoff,ijcoff)
115 
116  DO l=idbeg,idend
117  tmpb(l,1)=2._rfreal*fvar(l,ijkcb)-fvar(l,ijkcb+1)
118  tmpe(l,1)=2._rfreal*fvar(l,ijkce)-fvar(l,ijkce-1)
119  filtvar(l,ijkcb)=fact1(ndel(diri))*(tmpb(l,1)+fvar(l,ijkcb+1))+ &
120  fact2(ndel(diri))*fvar(l,ijkcb)
121  filtvar(l,ijkce)=fact1(ndel(diri))*(tmpe(l,1)+fvar(l,ijkce-1))+ &
122  fact2(ndel(diri))*fvar(l,ijkce)
123  ENDDO ! l
124  ENDDO ! j
125  ENDDO ! k
126 
127  ELSEIF (ndel(diri)==filwidth_four) THEN
128 
129  DO k=kbeg,kend
130  DO j=jbeg,jend
131 
132 ! ----- perform integration
133 
134  DO i=ibeg+2,iend-2
135  ijkc = indijk(i ,j ,k ,icoff,ijcoff)
136  DO l=idbeg,idend
137  filtvar(l,ijkc)=fact1(ndel(diri))*(fvar(l,ijkc-2)+fvar(l,ijkc+2))+ &
138  fact2(ndel(diri))*(fvar(l,ijkc+1)+fvar(l,ijkc)+ &
139  fvar(l,ijkc-1))
140  ENDDO ! l
141  ENDDO ! i
142 
143 ! ----- perform integration at the edges
144 
145  ijkcb = indijk(ibeg ,j ,k ,icoff,ijcoff)
146  ijkce = indijk(iend ,j ,k ,icoff,ijcoff)
147  ijkcb1= indijk(ibeg+1 ,j ,k ,icoff,ijcoff)
148  ijkce1= indijk(iend-1 ,j ,k ,icoff,ijcoff)
149 
150  DO l=idbeg,idend
151  tmpb(l,1)=2._rfreal*fvar(l,ijkcb)-fvar(l,ijkcb+1)
152  tmpe(l,1)=2._rfreal*fvar(l,ijkce)-fvar(l,ijkce-1)
153  filtvar(l,ijkcb1)=fact1(ndel(diri))*(tmpb(l,1)+fvar(l,ijkcb1+2))+ &
154  fact2(ndel(diri))*(fvar(l,ijkcb1+1)+fvar(l,ijkcb1)+ &
155  fvar(l,ijkcb1-1))
156  filtvar(l,ijkce1)=fact1(ndel(diri))*(tmpe(l,1)+fvar(l,ijkce1-2))+ &
157  fact2(ndel(diri))*(fvar(l,ijkce1+1)+fvar(l,ijkce1)+ &
158  fvar(l,ijkce1-1))
159  ENDDO ! l
160 
161  DO l=idbeg,idend
162  tmpb(l,2)=2._rfreal*fvar(l,ijkcb)-fvar(l,ijkcb+2)
163  tmpe(l,2)=2._rfreal*fvar(l,ijkce)-fvar(l,ijkce-2)
164  filtvar(l,ijkcb)=fact1(ndel(diri))*(tmpb(l,2)+fvar(l,ijkcb+2))+ &
165  fact2(ndel(diri))*(fvar(l,ijkcb+1)+fvar(l,ijkcb)+ &
166  tmpb(l,1))
167  filtvar(l,ijkce)=fact1(ndel(diri))*(tmpe(l,2)+fvar(l,ijkce-2))+ &
168  fact2(ndel(diri))*(fvar(l,ijkce-1)+fvar(l,ijkce)+ &
169  tmpe(l,1))
170  ENDDO ! l
171  ENDDO ! j
172  ENDDO ! k
173 
174  ENDIF ! nDel
175 
176 ! finalize --------------------------------------------------------------------
177 
178  CALL deregisterfunction( global )
179 
180 END SUBROUTINE turb_flolesunifiltcci
181 
182 !******************************************************************************
183 !
184 ! RCS Revision history:
185 !
186 ! $Log: TURB_floLesUniFiltCCI.F90,v $
187 ! Revision 1.5 2008/12/06 08:44:43 mtcampbe
188 ! Updated license.
189 !
190 ! Revision 1.4 2008/11/19 22:17:55 mtcampbe
191 ! Added Illinois Open Source License/Copyright
192 !
193 ! Revision 1.3 2004/06/03 02:12:45 wasistho
194 ! expand CC-filtering to all dummy layers
195 !
196 ! Revision 1.2 2004/03/12 02:55:36 wasistho
197 ! changed rocturb routine names
198 !
199 ! Revision 1.1 2004/03/08 23:35:45 wasistho
200 ! changed turb nomenclature
201 !
202 ! Revision 1.3 2003/05/16 05:44:03 wasistho
203 ! modified array range of CC-filtered
204 !
205 ! Revision 1.2 2003/05/15 02:57:06 jblazek
206 ! Inlined index function.
207 !
208 ! Revision 1.1 2002/10/14 23:55:30 wasistho
209 ! Install Rocturb
210 !
211 !******************************************************************************
212 
213 
214 
215 
216 
217 
218 
ibeg
**********************************************************************Rocstar Simulation Suite Illinois Rocstar LLC All rights reserved ****Illinois Rocstar LLC IL **www illinoisrocstar com **sales illinoisrocstar com WITHOUT WARRANTY OF ANY **EXPRESS OR INCLUDING BUT NOT LIMITED TO THE WARRANTIES **OF FITNESS FOR A PARTICULAR PURPOSE AND **NONINFRINGEMENT IN NO EVENT SHALL THE CONTRIBUTORS OR **COPYRIGHT HOLDERS BE LIABLE FOR ANY DAMAGES OR OTHER WHETHER IN AN ACTION OF TORT OR **Arising OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE **USE OR OTHER DEALINGS WITH THE SOFTWARE **********************************************************************INTERFACE SUBROUTINE ibeg
Definition: RFLO_GetPatchIndices.h:27
k
j indices k indices k
Definition: Indexing.h:6
moderror::registerfunction
subroutine registerfunction(global, funName, fileName)
Definition: ModError.F90:449
turb_modparameters
Definition: TURB_ModParameters.F90:39
iend
**********************************************************************Rocstar Simulation Suite Illinois Rocstar LLC All rights reserved ****Illinois Rocstar LLC IL **www illinoisrocstar com **sales illinoisrocstar com WITHOUT WARRANTY OF ANY **EXPRESS OR INCLUDING BUT NOT LIMITED TO THE WARRANTIES **OF FITNESS FOR A PARTICULAR PURPOSE AND **NONINFRINGEMENT IN NO EVENT SHALL THE CONTRIBUTORS OR **COPYRIGHT HOLDERS BE LIABLE FOR ANY DAMAGES OR OTHER WHETHER IN AN ACTION OF TORT OR **Arising OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE **USE OR OTHER DEALINGS WITH THE SOFTWARE **********************************************************************INTERFACE SUBROUTINE knode iend
Definition: RFLO_CalcGradDummySymm.h:32
i
blockLoc i
Definition: read.cpp:79
j
j indices j
Definition: Indexing.h:6
jend
**********************************************************************Rocstar Simulation Suite Illinois Rocstar LLC All rights reserved ****Illinois Rocstar LLC IL **www illinoisrocstar com **sales illinoisrocstar com WITHOUT WARRANTY OF ANY **EXPRESS OR INCLUDING BUT NOT LIMITED TO THE WARRANTIES **OF FITNESS FOR A PARTICULAR PURPOSE AND **NONINFRINGEMENT IN NO EVENT SHALL THE CONTRIBUTORS OR **COPYRIGHT HOLDERS BE LIABLE FOR ANY DAMAGES OR OTHER WHETHER IN AN ACTION OF TORT OR **Arising OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE **USE OR OTHER DEALINGS WITH THE SOFTWARE **********************************************************************INTERFACE SUBROUTINE knode jend
Definition: RFLO_CalcGradDummySymm.h:32
modglobal::t_global
Definition: ModGlobal.F90:50
jbeg
**********************************************************************Rocstar Simulation Suite Illinois Rocstar LLC All rights reserved ****Illinois Rocstar LLC IL **www illinoisrocstar com **sales illinoisrocstar com WITHOUT WARRANTY OF ANY **EXPRESS OR INCLUDING BUT NOT LIMITED TO THE WARRANTIES **OF FITNESS FOR A PARTICULAR PURPOSE AND **NONINFRINGEMENT IN NO EVENT SHALL THE CONTRIBUTORS OR **COPYRIGHT HOLDERS BE LIABLE FOR ANY DAMAGES OR OTHER WHETHER IN AN ACTION OF TORT OR **Arising OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE **USE OR OTHER DEALINGS WITH THE SOFTWARE **********************************************************************INTERFACE SUBROUTINE knode jbeg
Definition: RFLO_CalcGradDummySymm.h:32
kbeg
**********************************************************************Rocstar Simulation Suite Illinois Rocstar LLC All rights reserved ****Illinois Rocstar LLC IL **www illinoisrocstar com **sales illinoisrocstar com WITHOUT WARRANTY OF ANY **EXPRESS OR INCLUDING BUT NOT LIMITED TO THE WARRANTIES **OF FITNESS FOR A PARTICULAR PURPOSE AND **NONINFRINGEMENT IN NO EVENT SHALL THE CONTRIBUTORS OR **COPYRIGHT HOLDERS BE LIABLE FOR ANY DAMAGES OR OTHER WHETHER IN AN ACTION OF TORT OR **Arising OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE **USE OR OTHER DEALINGS WITH THE SOFTWARE **********************************************************************INTERFACE SUBROUTINE knode kbeg
Definition: RFLO_CalcGradDummySymm.h:32
moderror::deregisterfunction
subroutine deregisterfunction(global)
Definition: ModError.F90:469
turb_flolesunifiltcci
subroutine turb_flolesunifiltcci(global, nDum, ibeg, iend, jbeg, jend, kbeg, kend, iCOff, ijCOff, nDel, idBeg, idEnd, fact1, fact2, fVar, filtVar)
Definition: TURB_floLesUniFiltCCI.F90:52