Rocstar  1.0
Rocstar multiphysics simulation application
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
TURB_floLesUniFiltCCJ.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 J-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_floLesUniFiltCCJ.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_flolesunifiltccj( 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_floLesUniFiltCCJ.F90,v $'
82 
83  CALL registerfunction( global,'TURB_FloLesUniFiltCCJ',&
84  'TURB_floLesUniFiltCCJ.F90' )
85 
86 ! integration over J-direction -----------------------------------------------
87 
88  IF (ndel(dirj)==filwidth_zero) THEN
89 
90 ! - no filtering in J-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(dirj)==filwidth_one) .OR. &
96  (ndel(dirj)==filwidth_two)) THEN
97 
98  DO k=kbeg,kend
99  DO i=ibeg,iend
100 
101 ! ----- perform integration
102 
103  DO j=jbeg+1,jend-1
104  ijkc = indijk(i ,j ,k ,icoff,ijcoff)
105  DO l=idbeg,idend
106  filtvar(l,ijkc)=fact1(ndel(dirj))*(fvar(l,ijkc-icoff) + &
107  fvar(l,ijkc+icoff))+ &
108  fact2(ndel(dirj))*fvar(l,ijkc)
109  END DO ! l
110  ENDDO ! j
111 
112 ! ----- perform integration at the edges
113 
114  ijkcb = indijk(i ,jbeg ,k ,icoff,ijcoff)
115  ijkce = indijk(i ,jend ,k ,icoff,ijcoff)
116 
117  DO l=idbeg,idend
118  tmpb(l,1)=2._rfreal*fvar(l,ijkcb)-fvar(l,ijkcb+icoff)
119  tmpe(l,1)=2._rfreal*fvar(l,ijkce)-fvar(l,ijkce-icoff)
120  filtvar(l,ijkcb)=fact1(ndel(dirj))*(tmpb(l,1)+fvar(l,ijkcb+icoff))+ &
121  fact2(ndel(dirj))*fvar(l,ijkcb)
122  filtvar(l,ijkce)=fact1(ndel(dirj))*(tmpe(l,1)+fvar(l,ijkce-icoff))+ &
123  fact2(ndel(dirj))*fvar(l,ijkce)
124  ENDDO ! l
125 
126  ENDDO ! i
127  ENDDO ! k
128 
129  ELSEIF (ndel(dirj)==filwidth_four) THEN
130 
131  DO k=kbeg,kend
132  DO i=ibeg,iend
133 
134 ! ----- perform integration
135 
136  DO j=jbeg+2,jend-2
137  ijkc = indijk(i ,j ,k ,icoff,ijcoff)
138  DO l=idbeg,idend
139  filtvar(l,ijkc)=fact1(ndel(dirj))*(fvar(l,ijkc-2*icoff) + &
140  fvar(l,ijkc+2*icoff))+ &
141  fact2(ndel(dirj))*(fvar(l,ijkc+icoff) + &
142  fvar(l,ijkc)+fvar(l,ijkc-icoff))
143  ENDDO ! l
144  ENDDO ! j
145 
146 ! ----- perform integration at the edges
147 
148  ijkcb = indijk(i ,jbeg ,k ,icoff,ijcoff)
149  ijkce = indijk(i ,jend ,k ,icoff,ijcoff)
150  ijkcb1= indijk(i ,jbeg+1 ,k ,icoff,ijcoff)
151  ijkce1= indijk(i ,jend-1 ,k ,icoff,ijcoff)
152 
153  DO l=idbeg,idend
154  tmpb(l,1)=2._rfreal*fvar(l,ijkcb)-fvar(l,ijkcb+icoff)
155  tmpe(l,1)=2._rfreal*fvar(l,ijkce)-fvar(l,ijkce-icoff)
156  filtvar(l,ijkcb1)=fact1(ndel(dirj))*(tmpb(l,1)+fvar(l,ijkcb1+2*icoff))+ &
157  fact2(ndel(dirj))*(fvar(l,ijkcb1+icoff)+fvar(l,ijkcb1)+ &
158  fvar(l,ijkcb1-icoff))
159  filtvar(l,ijkce1)=fact1(ndel(dirj))*(tmpe(l,1)+fvar(l,ijkce1-2*icoff))+ &
160  fact2(ndel(dirj))*(fvar(l,ijkce1+icoff)+fvar(l,ijkce1)+ &
161  fvar(l,ijkce1-icoff))
162  ENDDO ! l
163 
164  DO l=idbeg,idend
165  tmpb(l,2)=2._rfreal*fvar(l,ijkcb)-fvar(l,ijkcb+2*icoff)
166  tmpe(l,2)=2._rfreal*fvar(l,ijkce)-fvar(l,ijkce-2*icoff)
167  filtvar(l,ijkcb)=fact1(ndel(dirj))*(tmpb(l,2)+fvar(l,ijkcb+2*icoff))+ &
168  fact2(ndel(dirj))*(fvar(l,ijkcb+icoff)+fvar(l,ijkcb)+ &
169  tmpb(l,1))
170  filtvar(l,ijkce)=fact1(ndel(dirj))*(tmpe(l,2)+fvar(l,ijkce-2*icoff))+ &
171  fact2(ndel(dirj))*(fvar(l,ijkce-icoff)+fvar(l,ijkce)+ &
172  tmpe(l,1))
173  ENDDO ! l
174  ENDDO ! i
175  ENDDO ! k
176  ENDIF ! nDel
177 
178 ! finalize --------------------------------------------------------------------
179 
180  CALL deregisterfunction( global )
181 
182 END SUBROUTINE turb_flolesunifiltccj
183 
184 !******************************************************************************
185 !
186 ! RCS Revision history:
187 !
188 ! $Log: TURB_floLesUniFiltCCJ.F90,v $
189 ! Revision 1.5 2008/12/06 08:44:43 mtcampbe
190 ! Updated license.
191 !
192 ! Revision 1.4 2008/11/19 22:17:55 mtcampbe
193 ! Added Illinois Open Source License/Copyright
194 !
195 ! Revision 1.3 2004/06/03 02:12:49 wasistho
196 ! expand CC-filtering to all dummy layers
197 !
198 ! Revision 1.2 2004/03/12 02:55:36 wasistho
199 ! changed rocturb routine names
200 !
201 ! Revision 1.1 2004/03/08 23:35:45 wasistho
202 ! changed turb nomenclature
203 !
204 ! Revision 1.3 2003/05/16 05:44:12 wasistho
205 ! modified array range of CC-filtered
206 !
207 ! Revision 1.2 2003/05/15 02:57:06 jblazek
208 ! Inlined index function.
209 !
210 ! Revision 1.1 2002/10/14 23:55:30 wasistho
211 ! Install Rocturb
212 !
213 !******************************************************************************
214 
215 
216 
217 
218 
219 
220 
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
turb_flolesunifiltccj
subroutine turb_flolesunifiltccj(global, nDum, ibeg, iend, jbeg, jend, kbeg, kend, iCOff, ijCOff, nDel, idBeg, idEnd, fact1, fact2, fVar, filtVar)
Definition: TURB_floLesUniFiltCCJ.F90:52
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