Rocstar  1.0
Rocstar multiphysics simulation application
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
TURB_floRansSARoe2ndFlux.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: compute SA convective fluxes based on 2nd-order Roe upwind scheme.
26 !
27 ! Description: none.
28 !
29 ! Input: region = data of current region.
30 !
31 ! Output: region%levels%turb%rhs = convective fluxes added to the residual.
32 !
33 ! Notes: uses MUSCL scheme with kappa=1/3.
34 !
35 !******************************************************************************
36 !
37 ! $Id: TURB_floRansSARoe2ndFlux.F90,v 1.4 2008/12/06 08:44:44 mtcampbe Exp $
38 !
39 ! Copyright: (c) 2003 by the University of Illinois
40 !
41 !******************************************************************************
42 
43 SUBROUTINE turb_floranssaroe2ndflux( region )
44 
45  USE moddatatypes
46  USE moddatastruct, ONLY : t_region
47  USE modinterfaces, ONLY : rflo_getdimensphys, rflo_getcelloffset, &
48  rflo_getnodeoffset
49  USE turb_modinterfaces, ONLY : turb_floranssaroefluxpatch
50  USE moderror
51  USE modparameters
52  USE turb_modparameters
53  IMPLICIT NONE
54 
55 #include "Indexing.h"
56 
57 ! ... parameters
58  TYPE(t_region) :: region
59 
60 ! ... loop variables
61  INTEGER :: i, j, k, ipatch
62 
63 ! ... local variables
64  INTEGER :: ilev, ipcbeg, ipcend, jpcbeg, jpcend, kpcbeg, kpcend
65  INTEGER :: icoff, ijcoff, inoff, ijnoff, ijkc0, ijkc1, ijkc2, ijkc3, ijkn
66  INTEGER :: indsvel
67 
68  REAL(RFREAL) :: limfac, limfac3, rvolref, vola, eps2(3), dvar(4), dvarm(4), &
69  dvarp(4), qsl, qsr, deltl(4), deltr(4), fc, eps2n, &
70  rl, rr, ul, ur, vl, vr, wl, wr, qsrl, qsrr
71  REAL(RFREAL), POINTER :: dv(:,:), si(:,:), sj(:,:), sk(:,:), vol(:)
72  REAL(RFREAL), POINTER :: tcv(:,:), trhs(:,:), sivel(:), sjvel(:), skvel(:)
73 
74 ! ... functions
75  REAL(RFREAL) :: muscl3, af, bf, eps
76 
77 !******************************************************************************
78 ! limiter function
79 
80  muscl3(af,bf,eps) = (bf*(2._rfreal*af*af+eps)+af*(bf*bf+2._rfreal*eps))/ &
81  (2._rfreal*af*af+2._rfreal*bf*bf-af*bf+ &
82  3._rfreal*eps+1.e-30_rfreal)
83 
84  CALL registerfunction( region%global,'TURB_FloRansSARoe2ndFlux',&
85  'TURB_floRansSARoe2ndFlux.F90' )
86 
87 ! get dimensions and pointers -------------------------------------------------
88 
89  ilev = region%currLevel
90 
91  CALL rflo_getdimensphys( region,ilev,ipcbeg,ipcend, &
92  jpcbeg,jpcend,kpcbeg,kpcend )
93  CALL rflo_getcelloffset( region,ilev,icoff,ijcoff )
94  CALL rflo_getnodeoffset( region,ilev,inoff,ijnoff )
95 
96  dv => region%levels(ilev)%mixt%dv
97  tcv => region%levels(ilev)%turb%cv
98  trhs => region%levels(ilev)%turb%rhs
99  vol => region%levels(ilev)%grid%vol
100  si => region%levels(ilev)%grid%si
101  sj => region%levels(ilev)%grid%sj
102  sk => region%levels(ilev)%grid%sk
103  sivel => region%levels(ilev)%grid%siVel
104  sjvel => region%levels(ilev)%grid%sjVel
105  skvel => region%levels(ilev)%grid%skVel
106  indsvel = region%levels(ilev)%grid%indSvel
107  limfac = region%mixtInput%limFac
108 
109 ! normalise epsilon^2 for all limited variables (rho, u, v, w, p) -------------
110 
111  limfac3 = limfac*limfac*limfac
112  rvolref = 1._rfreal/region%global%limVolRef**1.5_rfreal
113  eps2(1) = limfac3*region%global%limRef(1)*region%global%limRef(1)*rvolref
114  eps2(2) = limfac3*region%global%limRef(2)*region%global%limRef(2)*rvolref
115  eps2(3) = limfac3*region%global%limRef(3)*region%global%limRef(3)*rvolref
116 
117 ! flux in i-direction (except through boundary) -------------------------------
118 
119  DO k=kpcbeg,kpcend
120  DO j=jpcbeg,jpcend
121  DO i=ipcbeg+1,ipcend
122  ijkc0 = indijk(i ,j,k,icoff,ijcoff)
123  ijkc1 = indijk(i-1,j,k,icoff,ijcoff)
124  ijkc2 = indijk(i-2,j,k,icoff,ijcoff)
125  ijkc3 = indijk(i+1,j,k,icoff,ijcoff)
126  ijkn = indijk(i ,j,k,inoff,ijnoff)
127 
128  vola = (0.5_rfreal*(vol(ijkc0)+vol(ijkc1)))**1.5_rfreal
129 
130 ! ----- limited extrapolations
131 
132  dvar(1) = tcv(cv_sa_nutil,ijkc0) - tcv(cv_sa_nutil,ijkc1)
133  dvar(2) = dv(dv_mixt_uvel,ijkc0) - dv(dv_mixt_uvel,ijkc1)
134  dvar(3) = dv(dv_mixt_vvel,ijkc0) - dv(dv_mixt_vvel,ijkc1)
135  dvar(4) = dv(dv_mixt_wvel,ijkc0) - dv(dv_mixt_wvel,ijkc1)
136 
137  dvarm(1) = tcv(cv_sa_nutil,ijkc1) - tcv(cv_sa_nutil,ijkc2)
138  dvarm(2) = dv(dv_mixt_uvel,ijkc1) - dv(dv_mixt_uvel,ijkc2)
139  dvarm(3) = dv(dv_mixt_vvel,ijkc1) - dv(dv_mixt_vvel,ijkc2)
140  dvarm(4) = dv(dv_mixt_wvel,ijkc1) - dv(dv_mixt_wvel,ijkc2)
141 
142  dvarp(1) = tcv(cv_sa_nutil,ijkc3) - tcv(cv_sa_nutil,ijkc0)
143  dvarp(2) = dv(dv_mixt_uvel,ijkc3) - dv(dv_mixt_uvel,ijkc0)
144  dvarp(3) = dv(dv_mixt_vvel,ijkc3) - dv(dv_mixt_vvel,ijkc0)
145  dvarp(4) = dv(dv_mixt_wvel,ijkc3) - dv(dv_mixt_wvel,ijkc0)
146 
147  eps2n = eps2(1)*vola
148  deltl(1) = 0.5_rfreal*muscl3(dvar(1) ,dvarm(1),eps2n)
149  deltr(1) = 0.5_rfreal*muscl3(dvarp(1),dvar(1) ,eps2n)
150  eps2n = eps2(2)*vola
151  deltl(2) = 0.5_rfreal*muscl3(dvar(2) ,dvarm(2),eps2n)
152  deltr(2) = 0.5_rfreal*muscl3(dvarp(2),dvar(2) ,eps2n)
153  deltl(3) = 0.5_rfreal*muscl3(dvar(3) ,dvarm(3),eps2n)
154  deltr(3) = 0.5_rfreal*muscl3(dvarp(3),dvar(3) ,eps2n)
155  deltl(4) = 0.5_rfreal*muscl3(dvar(4) ,dvarm(4),eps2n)
156  deltr(4) = 0.5_rfreal*muscl3(dvarp(4),dvar(4) ,eps2n)
157 
158 ! ----- left and right states
159 
160  rl = tcv(cv_sa_nutil,ijkc1) + deltl(1)
161  ul = dv(dv_mixt_uvel,ijkc1) + deltl(2)
162  vl = dv(dv_mixt_vvel,ijkc1) + deltl(3)
163  wl = dv(dv_mixt_wvel,ijkc1) + deltl(4)
164  qsrl = (ul*si(xcoord,ijkn)+vl*si(ycoord,ijkn)+ &
165  wl*si(zcoord,ijkn)-sivel(ijkn*indsvel))*rl
166 
167  rr = tcv(cv_sa_nutil,ijkc0) - deltr(1)
168  ur = dv(dv_mixt_uvel,ijkc0) - deltr(2)
169  vr = dv(dv_mixt_vvel,ijkc0) - deltr(3)
170  wr = dv(dv_mixt_wvel,ijkc0) - deltr(4)
171  qsrr = (ur*si(xcoord,ijkn)+vr*si(ycoord,ijkn)+ &
172  wr*si(zcoord,ijkn)-sivel(ijkn*indsvel))*rr
173 
174 ! ----- fluxes
175 
176  fc = 0.5_rfreal*(qsrl +qsrr )
177  trhs(cv_sa_nutil,ijkc0) = trhs(cv_sa_nutil,ijkc0) + fc
178  trhs(cv_sa_nutil,ijkc1) = trhs(cv_sa_nutil,ijkc1) - fc
179  ENDDO ! i
180  ENDDO ! j
181  ENDDO ! k
182 
183 ! flux in j-direction (except through boundary) -------------------------------
184 
185  DO k=kpcbeg,kpcend
186  DO j=jpcbeg+1,jpcend
187  DO i=ipcbeg,ipcend
188  ijkc0 = indijk(i,j ,k,icoff,ijcoff)
189  ijkc1 = indijk(i,j-1,k,icoff,ijcoff)
190  ijkc2 = indijk(i,j-2,k,icoff,ijcoff)
191  ijkc3 = indijk(i,j+1,k,icoff,ijcoff)
192  ijkn = indijk(i,j ,k,inoff,ijnoff)
193 
194  vola = (0.5_rfreal*(vol(ijkc0)+vol(ijkc1)))**1.5_rfreal
195 
196 ! ----- limited extrapolations
197 
198  dvar(1) = tcv(cv_sa_nutil,ijkc0) - tcv(cv_sa_nutil,ijkc1)
199  dvar(2) = dv(dv_mixt_uvel,ijkc0) - dv(dv_mixt_uvel,ijkc1)
200  dvar(3) = dv(dv_mixt_vvel,ijkc0) - dv(dv_mixt_vvel,ijkc1)
201  dvar(4) = dv(dv_mixt_wvel,ijkc0) - dv(dv_mixt_wvel,ijkc1)
202 
203  dvarm(1) = tcv(cv_sa_nutil,ijkc1) - tcv(cv_sa_nutil,ijkc2)
204  dvarm(2) = dv(dv_mixt_uvel,ijkc1) - dv(dv_mixt_uvel,ijkc2)
205  dvarm(3) = dv(dv_mixt_vvel,ijkc1) - dv(dv_mixt_vvel,ijkc2)
206  dvarm(4) = dv(dv_mixt_wvel,ijkc1) - dv(dv_mixt_wvel,ijkc2)
207 
208  dvarp(1) = tcv(cv_sa_nutil,ijkc3) - tcv(cv_sa_nutil,ijkc0)
209  dvarp(2) = dv(dv_mixt_uvel,ijkc3) - dv(dv_mixt_uvel,ijkc0)
210  dvarp(3) = dv(dv_mixt_vvel,ijkc3) - dv(dv_mixt_vvel,ijkc0)
211  dvarp(4) = dv(dv_mixt_wvel,ijkc3) - dv(dv_mixt_wvel,ijkc0)
212 
213  eps2n = eps2(1)*vola
214  deltl(1) = 0.5_rfreal*muscl3(dvar(1) ,dvarm(1),eps2n)
215  deltr(1) = 0.5_rfreal*muscl3(dvarp(1),dvar(1) ,eps2n)
216  eps2n = eps2(2)*vola
217  deltl(2) = 0.5_rfreal*muscl3(dvar(2) ,dvarm(2),eps2n)
218  deltr(2) = 0.5_rfreal*muscl3(dvarp(2),dvar(2) ,eps2n)
219  deltl(3) = 0.5_rfreal*muscl3(dvar(3) ,dvarm(3),eps2n)
220  deltr(3) = 0.5_rfreal*muscl3(dvarp(3),dvar(3) ,eps2n)
221  deltl(4) = 0.5_rfreal*muscl3(dvar(4) ,dvarm(4),eps2n)
222  deltr(4) = 0.5_rfreal*muscl3(dvarp(4),dvar(4) ,eps2n)
223 
224 ! ----- left and right states
225 
226  rl = tcv(cv_sa_nutil,ijkc1) + deltl(1)
227  ul = dv(dv_mixt_uvel,ijkc1) + deltl(2)
228  vl = dv(dv_mixt_vvel,ijkc1) + deltl(3)
229  wl = dv(dv_mixt_wvel,ijkc1) + deltl(4)
230  qsrl = (ul*sj(xcoord,ijkn)+vl*sj(ycoord,ijkn)+ &
231  wl*sj(zcoord,ijkn)-sjvel(ijkn*indsvel))*rl
232 
233  rr = tcv(cv_sa_nutil,ijkc0) - deltr(1)
234  ur = dv(dv_mixt_uvel,ijkc0) - deltr(2)
235  vr = dv(dv_mixt_vvel,ijkc0) - deltr(3)
236  wr = dv(dv_mixt_wvel,ijkc0) - deltr(4)
237  qsrr = (ur*sj(xcoord,ijkn)+vr*sj(ycoord,ijkn)+ &
238  wr*sj(zcoord,ijkn)-sjvel(ijkn*indsvel))*rr
239 
240 ! ----- fluxes
241 
242  fc = 0.5_rfreal*(qsrl +qsrr )
243  trhs(cv_sa_nutil,ijkc0) = trhs(cv_sa_nutil,ijkc0) + fc
244  trhs(cv_sa_nutil,ijkc1) = trhs(cv_sa_nutil,ijkc1) - fc
245  ENDDO ! i
246  ENDDO ! j
247  ENDDO ! k
248 
249 ! flux in k-direction (except through boundary) -------------------------------
250 
251  DO k=kpcbeg+1,kpcend
252  DO j=jpcbeg,jpcend
253  DO i=ipcbeg,ipcend
254  ijkc0 = indijk(i,j,k ,icoff,ijcoff)
255  ijkc1 = indijk(i,j,k-1,icoff,ijcoff)
256  ijkc2 = indijk(i,j,k-2,icoff,ijcoff)
257  ijkc3 = indijk(i,j,k+1,icoff,ijcoff)
258  ijkn = indijk(i,j,k ,inoff,ijnoff)
259 
260  vola = (0.5_rfreal*(vol(ijkc0)+vol(ijkc1)))**1.5_rfreal
261 
262 ! ----- limited extrapolations
263 
264  dvar(1) = tcv(cv_sa_nutil,ijkc0) - tcv(cv_sa_nutil,ijkc1)
265  dvar(2) = dv(dv_mixt_uvel,ijkc0) - dv(dv_mixt_uvel,ijkc1)
266  dvar(3) = dv(dv_mixt_vvel,ijkc0) - dv(dv_mixt_vvel,ijkc1)
267  dvar(4) = dv(dv_mixt_wvel,ijkc0) - dv(dv_mixt_wvel,ijkc1)
268 
269  dvarm(1) = tcv(cv_sa_nutil,ijkc1) - tcv(cv_sa_nutil,ijkc2)
270  dvarm(2) = dv(dv_mixt_uvel,ijkc1) - dv(dv_mixt_uvel,ijkc2)
271  dvarm(3) = dv(dv_mixt_vvel,ijkc1) - dv(dv_mixt_vvel,ijkc2)
272  dvarm(4) = dv(dv_mixt_wvel,ijkc1) - dv(dv_mixt_wvel,ijkc2)
273 
274  dvarp(1) = tcv(cv_sa_nutil,ijkc3) - tcv(cv_sa_nutil,ijkc0)
275  dvarp(2) = dv(dv_mixt_uvel,ijkc3) - dv(dv_mixt_uvel,ijkc0)
276  dvarp(3) = dv(dv_mixt_vvel,ijkc3) - dv(dv_mixt_vvel,ijkc0)
277  dvarp(4) = dv(dv_mixt_wvel,ijkc3) - dv(dv_mixt_wvel,ijkc0)
278 
279  eps2n = eps2(1)*vola
280  deltl(1) = 0.5_rfreal*muscl3(dvar(1) ,dvarm(1),eps2n)
281  deltr(1) = 0.5_rfreal*muscl3(dvarp(1),dvar(1) ,eps2n)
282  eps2n = eps2(2)*vola
283  deltl(2) = 0.5_rfreal*muscl3(dvar(2) ,dvarm(2),eps2n)
284  deltr(2) = 0.5_rfreal*muscl3(dvarp(2),dvar(2) ,eps2n)
285  deltl(3) = 0.5_rfreal*muscl3(dvar(3) ,dvarm(3),eps2n)
286  deltr(3) = 0.5_rfreal*muscl3(dvarp(3),dvar(3) ,eps2n)
287  deltl(4) = 0.5_rfreal*muscl3(dvar(4) ,dvarm(4),eps2n)
288  deltr(4) = 0.5_rfreal*muscl3(dvarp(4),dvar(4) ,eps2n)
289 
290 ! ----- left and right states
291 
292  rl = tcv(cv_sa_nutil,ijkc1) + deltl(1)
293  ul = dv(dv_mixt_uvel,ijkc1) + deltl(2)
294  vl = dv(dv_mixt_vvel,ijkc1) + deltl(3)
295  wl = dv(dv_mixt_wvel,ijkc1) + deltl(4)
296  qsrl = (ul*sk(xcoord,ijkn)+vl*sk(ycoord,ijkn)+ &
297  wl*sk(zcoord,ijkn)-skvel(ijkn*indsvel))*rl
298 
299  rr = tcv(cv_sa_nutil,ijkc0) - deltr(1)
300  ur = dv(dv_mixt_uvel,ijkc0) - deltr(2)
301  vr = dv(dv_mixt_vvel,ijkc0) - deltr(3)
302  wr = dv(dv_mixt_wvel,ijkc0) - deltr(4)
303  qsrr = (ur*sk(xcoord,ijkn)+vr*sk(ycoord,ijkn)+ &
304  wr*sk(zcoord,ijkn)-skvel(ijkn*indsvel))*rr
305 
306 ! ----- fluxes
307 
308  fc = 0.5_rfreal*(qsrl +qsrr )
309  trhs(cv_sa_nutil,ijkc0) = trhs(cv_sa_nutil,ijkc0) + fc
310  trhs(cv_sa_nutil,ijkc1) = trhs(cv_sa_nutil,ijkc1) - fc
311  ENDDO ! i
312  ENDDO ! j
313  ENDDO ! k
314 
315 ! fluxes through boundaries ---------------------------------------------------
316 
317  DO ipatch=1,region%nPatches
318  CALL turb_floranssaroefluxpatch( region, &
319  region%levels(ilev)%patches(ipatch) )
320  ENDDO
321 
322 ! finalize --------------------------------------------------------------------
323 
324  CALL deregisterfunction( region%global )
325 
326 END SUBROUTINE turb_floranssaroe2ndflux
327 
328 !******************************************************************************
329 !
330 ! RCS Revision history:
331 !
332 ! $Log: TURB_floRansSARoe2ndFlux.F90,v $
333 ! Revision 1.4 2008/12/06 08:44:44 mtcampbe
334 ! Updated license.
335 !
336 ! Revision 1.3 2008/11/19 22:17:56 mtcampbe
337 ! Added Illinois Open Source License/Copyright
338 !
339 ! Revision 1.2 2004/03/12 02:55:36 wasistho
340 ! changed rocturb routine names
341 !
342 ! Revision 1.1 2004/03/08 23:35:46 wasistho
343 ! changed turb nomenclature
344 !
345 ! Revision 1.2 2003/10/27 23:10:26 wasistho
346 ! established RaNS 2nd order upwind scheme
347 !
348 ! Revision 1.1 2003/10/27 04:54:38 wasistho
349 ! added RaNS upwind schemes
350 !
351 !
352 !******************************************************************************
353 
354 
355 
356 
357 
358 
359 
k
j indices k indices k
Definition: Indexing.h:6
kpcbeg
**********************************************************************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 kpcbeg
Definition: RFLO_GetDimensPhys.h:27
moderror::registerfunction
subroutine registerfunction(global, funName, fileName)
Definition: ModError.F90:449
jpcbeg
**********************************************************************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 jpcbeg
Definition: RFLO_GetDimensPhys.h:27
turb_modinterfaces
Definition: TURB_ModInterfaces.F90:39
turb_modparameters
Definition: TURB_ModParameters.F90:39
ipcend
**********************************************************************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 ipcend
Definition: RFLO_GetDimensPhys.h:27
rflo_getnodeoffset
subroutine rflo_getnodeoffset(region, iLev, iNodeOffset, ijNodeOffset)
Definition: RFLO_GetDimensions.F90:267
ipcbeg
**********************************************************************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 ipcbeg
Definition: RFLO_GetDimensPhys.h:27
i
blockLoc i
Definition: read.cpp:79
rflo_getcelloffset
subroutine rflo_getcelloffset(region, iLev, iCellOffset, ijCellOffset)
Definition: RFLO_GetDimensions.F90:235
turb_floranssaroefluxpatch
subroutine turb_floranssaroefluxpatch(region, patch)
Definition: TURB_floRansSARoeFluxPatch.F90:45
j
j indices j
Definition: Indexing.h:6
jpcend
**********************************************************************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 jpcend
Definition: RFLO_GetDimensPhys.h:27
turb_floranssaroe2ndflux
subroutine turb_floranssaroe2ndflux(region)
Definition: TURB_floRansSARoe2ndFlux.F90:43
moderror::deregisterfunction
subroutine deregisterfunction(global)
Definition: ModError.F90:469
rflo_getdimensphys
subroutine rflo_getdimensphys(region, iLev, ipcbeg, ipcend, jpcbeg, jpcend, kpcbeg, kpcend)
Definition: RFLO_GetDimensions.F90:109