Rocstar  1.0
Rocstar multiphysics simulation application
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
RADI_DiffRadFluxPatch.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 diffusion approximation, FLDSRC or ROSS, radiative fluxes
26 ! through a patch by using an average of variables.
27 !
28 ! Description: none.
29 !
30 ! Input: region = data of current region
31 ! patch = current patch.
32 !
33 ! Output: region%levels%mixt%rhs = radiative fluxes added to the residual.
34 !
35 ! Notes: In addition to the energy source term due to radiation (div(qr)),
36 ! radiant fluxes normal to surface (cell faces) are stored in data
37 ! structure to be delivered later to combustion module, through Genx,
38 ! for burn-rate calculations. The radiant fluxes are stored always in
39 ! positive directions w.r.t lbound = 1, 3, and 5. Opposite sign
40 ! should be applied later for lbound = 2, 4, and 6 before transfering
41 ! the radiant fluxes through Genx.
42 !
43 !******************************************************************************
44 !
45 ! $Id: RADI_DiffRadFluxPatch.F90,v 1.4 2008/12/06 08:44:37 mtcampbe Exp $
46 !
47 ! Copyright: (c) 2001 by the University of Illinois
48 !
49 !******************************************************************************
50 
51 SUBROUTINE radi_diffradfluxpatch( region,patch )
52 
53  USE moddatatypes
54  USE modbndpatch, ONLY : t_patch
55  USE moddatastruct, ONLY : t_region
56  USE modinterfaces, ONLY : rflo_getpatchindices, rflo_getpatchdirection, &
57  rflo_getcelloffset, rflo_getnodeoffset
58  USE moderror
59  USE modparameters
60  USE radi_modparameters
61  IMPLICIT NONE
62 
63 #include "Indexing.h"
64 
65 ! ... parameters
66  TYPE(t_region) :: region
67  TYPE(t_patch) :: patch
68 
69 ! ... loop variables
70  INTEGER :: i, j, k, n1, n2
71 
72 ! ... local variables
73  INTEGER :: ilev, lbound, bctype, distrib, flowmodel, bcopt
74  INTEGER :: ibeg, iend, jbeg, jend, kbeg, kend, idir, jdir, kdir
75  INTEGER :: icoff, ijcoff, inoff, ijnoff, ijkcd, ijkcb0, ijkcb1, ijknb
76  INTEGER :: inode, jnode, knode
77  INTEGER :: idbeg, idend, jdbeg, jdend, kdbeg, kdend, i2d, noff
78 
79  REAL(RFREAL) :: sgn, gfaca, flima, coefa, tempa, tburn, twall, mrate
80  REAL(RFREAL) :: stboltz, rati, beta, qrx, qry, qrz, fr, modsf, sf(3)
81  REAL(RFREAL), POINTER :: rhs(:,:), sface(:,:), vals(:,:), grad(:,:)
82  REAL(RFREAL), POINTER :: rdv(:,:), coef(:,:), gofact(:), flim(:), qr(:)
83  REAL(RFREAL), POINTER :: wvint(:,:)
84 
85 !******************************************************************************
86 
87  CALL registerfunction( region%global,'RADI_DiffRadFluxPatch',&
88  'RADI_DiffRadFluxPatch.F90' )
89 
90 ! get dimensions and pointers -------------------------------------------------
91 
92  ilev = region%currLevel
93  lbound = patch%lbound
94 
95  CALL rflo_getpatchindices( region,patch,ilev,ibeg,iend,jbeg,jend,kbeg,kend )
96  CALL rflo_getpatchdirection( patch,idir,jdir,kdir )
97  CALL rflo_getcelloffset( region,ilev,icoff,ijcoff )
98  CALL rflo_getnodeoffset( region,ilev,inoff,ijnoff )
99 
100  bctype = patch%bcType
101  noff = abs(patch%l1end-patch%l1beg) + 1
102  distrib = patch%mixt%distrib
103  flowmodel = region%mixtInput%flowModel
104 
105  rhs => region%levels(ilev)%mixt%rhs
106  vals => patch%mixt%vals
107  rdv => region%levels(ilev)%radi%dv
108  flim => region%levels(ilev)%radi%fluxLim
109  coef => region%levels(ilev)%radi%radCoef
110  gofact=> region%levels(ilev)%radi%goFact
111  wvint => region%levels(ilev)%radi%wvInt
112 
113 ! take the right face vector and make it point outwards
114 
115  sgn = +1._rfreal
116  inode = 0
117  jnode = 0
118  knode = 0
119  IF (lbound==2 .OR. lbound==4 .OR. lbound==6) THEN
120  sgn = -1._rfreal
121  inode = -idir
122  jnode = -jdir
123  knode = -kdir
124  ENDIF
125 
126 ! get the appropriate face vector, gradients and radiant flux
127 
128  IF (lbound==1 .OR. lbound==2) THEN
129  sface => region%levels(ilev)%grid%si
130  grad => region%levels(ilev)%mixt%gradi
131  qr => region%levels(ilev)%radi%qri
132  ELSE IF (lbound==3 .OR. lbound==4) THEN
133  sface => region%levels(ilev)%grid%sj
134  grad => region%levels(ilev)%mixt%gradj
135  qr => region%levels(ilev)%radi%qrj
136  ELSE IF (lbound==5 .OR. lbound==6) THEN
137  sface => region%levels(ilev)%grid%sk
138  grad => region%levels(ilev)%mixt%gradk
139  qr => region%levels(ilev)%radi%qrk
140  ENDIF
141 
142 ! get radiation and rk-stage constants
143 
144  stboltz = region%radiInput%stBoltz
145  rati = -16._rfreal
146  beta = region%mixtInput%betrk(region%irkStep)
147 
148 ! stationary grid -------------------------------------------------------------
149 ! slip wall
150 
151  IF (bctype>=bc_slipwall .AND. bctype<=bc_slipwall+bc_range) THEN
152 
153  DO k=kbeg,kend
154  DO j=jbeg,jend
155  DO i=ibeg,iend
156  ijkcb0 = indijk(i ,j ,k ,icoff,ijcoff) ! boundary
157  ijkcb1 = indijk(i+idir ,j+jdir ,k+kdir ,icoff,ijcoff) ! interior
158  ijknb = indijk(i+inode,j+jnode,k+knode,inoff,ijnoff)
159 
160  gfaca = gofact(ijkcb0)
161  flima = 0.5_rfreal*(3._rfreal*flim(ijkcb0)-flim(ijkcb1))
162  tempa = 0.5_rfreal*(3._rfreal*rdv(dv_radi_teff,ijkcb0)- &
163  rdv(dv_radi_teff,ijkcb1))
164  coefa = 0.5_rfreal*(3._rfreal*coef(ijkcb0,radi_coeff_extinct)- &
165  coef(ijkcb1,radi_coeff_extinct))
166 
167  coefa = gfaca*rati*flima*stboltz*tempa**3/coefa
168 
169  qrx = coefa*grad(gr_mixt_tx,ijknb)
170  qry = coefa*grad(gr_mixt_ty,ijknb)
171  qrz = coefa*grad(gr_mixt_tz,ijknb)
172 
173  sf(1) = sgn*sface(xcoord,ijknb)
174  sf(2) = sgn*sface(ycoord,ijknb)
175  sf(3) = sgn*sface(zcoord,ijknb)
176 
177  fr = qrx*sf(1) + qry*sf(2) + qrz*sf(3)
178 
179  rhs(cv_mixt_ener,ijkcb0) = rhs(cv_mixt_ener,ijkcb0) - fr*beta
180 
181 ! ------- store rad. flux normal to surf. in positive direction (hence sgn* )
182  modsf = sqrt( sf(1)*sf(1) + sf(2)*sf(2) + sf(3)*sf(3) )
183  qr(ijknb) = -sgn*fr/modsf
184  ENDDO ! i
185  ENDDO ! j
186  ENDDO ! k
187 
188 ! noslip wall
189 
190  ELSE IF (bctype>=bc_noslipwall .AND. bctype<=bc_noslipwall+bc_range) THEN
191 
192  bcopt = patch%mixt%switches(bcswi_noslip_adiabat)
193 
194  DO k=kbeg,kend
195  DO j=jbeg,jend
196  DO i=ibeg,iend
197  ijkcb0 = indijk(i ,j ,k ,icoff,ijcoff) ! boundary
198  ijknb = indijk(i+inode,j+jnode,k+knode,inoff,ijnoff)
199 
200 ! ------- adiabatic wall
201  IF (bcopt == bcopt_adiabat) THEN
202  tempa = rdv(dv_radi_teff,ijkcb0)
203  ELSE
204 
205 ! ------- prescribed wall temperature
206  IF (lbound==1 .OR. lbound==2) THEN
207  n1 = j - jbeg
208  n2 = k - kbeg
209  ELSE IF (lbound==3 .OR. lbound==4) THEN
210  n1 = k - kbeg
211  n2 = i - ibeg
212  ELSE IF (lbound==5 .OR. lbound==6) THEN
213  n1 = i - ibeg
214  n2 = j - jbeg
215  ENDIF
216  i2d = distrib * indij(n1,n2,noff)
217  twall = vals(bcdat_noslip_twall,i2d)
218  tempa = twall
219  ENDIF
220 
221  gfaca = gofact(ijkcb0)
222  flima = flim(ijkcb0)
223  coefa = coef(ijkcb0,radi_coeff_extinct)
224  coefa = gfaca*rati*flima*stboltz*tempa**3/coefa
225 
226  qrx = coefa*grad(gr_mixt_tx,ijknb)
227  qry = coefa*grad(gr_mixt_ty,ijknb)
228  qrz = coefa*grad(gr_mixt_tz,ijknb)
229 
230  sf(1) = sgn*sface(xcoord,ijknb)
231  sf(2) = sgn*sface(ycoord,ijknb)
232  sf(3) = sgn*sface(zcoord,ijknb)
233 
234  fr = qrx*sf(1) + qry*sf(2) + qrz*sf(3)
235 
236  rhs(cv_mixt_ener,ijkcb0) = rhs(cv_mixt_ener,ijkcb0) - fr*beta
237 
238 ! ------- store rad. flux normal to surf. in positive direction (hence sgn* )
239  modsf = sqrt( sf(1)*sf(1) + sf(2)*sf(2) + sf(3)*sf(3) )
240  qr(ijknb) = -sgn*fr/modsf
241  ENDDO ! i
242  ENDDO ! j
243  ENDDO ! k
244 
245 ! injection boundary
246 
247  ELSE IF (bctype>=bc_injection .AND. bctype<=bc_injection+bc_range) THEN
248 
249 
250  DO k=kbeg,kend
251  DO j=jbeg,jend
252  DO i=ibeg,iend
253  ijkcb0 = indijk(i ,j ,k ,icoff,ijcoff) ! boundary
254  ijkcb1 = indijk(i+idir ,j+jdir ,k+kdir ,icoff,ijcoff) ! interior
255  ijknb = indijk(i+inode,j+jnode,k+knode,inoff,ijnoff)
256 
257  IF (lbound==1 .OR. lbound==2) THEN
258  n1 = j - jbeg
259  n2 = k - kbeg
260  ELSE IF (lbound==3 .OR. lbound==4) THEN
261  n1 = k - kbeg
262  n2 = i - ibeg
263  ELSE IF (lbound==5 .OR. lbound==6) THEN
264  n1 = i - ibeg
265  n2 = j - jbeg
266  ENDIF
267  i2d = distrib * indij(n1,n2,noff)
268  mrate = vals(bcdat_inject_mfrate,i2d)
269  tburn = vals(bcdat_inject_temp ,i2d)
270  gfaca = gofact(ijkcb0)
271 
272  IF (mrate > 0._rfreal) THEN ! surface burning
273  flima = flim(ijkcb0)
274  tempa = tburn
275  coefa = coef(ijkcb0,radi_coeff_extinct)
276 
277  ELSE ! not burning - slip/noslip wall
278  IF (flowmodel == flow_euler) THEN
279  flima = 0.5_rfreal*(3._rfreal*flim(ijkcb0)-flim(ijkcb1))
280  tempa = 0.5_rfreal*(3._rfreal*rdv(dv_radi_teff,ijkcb0)- &
281  rdv(dv_radi_teff,ijkcb1))
282  coefa = 0.5_rfreal*(3._rfreal*coef(ijkcb0,radi_coeff_extinct)- &
283  coef(ijkcb1,radi_coeff_extinct))
284 
285  ELSE ! treated as NS adiabatic
286  flima = flim(ijkcb0)
287  tempa = rdv(dv_radi_teff,ijkcb0)
288  coefa = coef(ijkcb0,radi_coeff_extinct)
289  ENDIF
290  ENDIF
291 
292  coefa = gfaca*rati*flima*stboltz*tempa**3/coefa
293 
294  qrx = coefa*grad(gr_mixt_tx,ijknb)
295  qry = coefa*grad(gr_mixt_ty,ijknb)
296  qrz = coefa*grad(gr_mixt_tz,ijknb)
297 
298  sf(1) = sgn*sface(xcoord,ijknb)
299  sf(2) = sgn*sface(ycoord,ijknb)
300  sf(3) = sgn*sface(zcoord,ijknb)
301 
302  fr = qrx*sf(1) + qry*sf(2) + qrz*sf(3)
303 
304  rhs(cv_mixt_ener,ijkcb0) = rhs(cv_mixt_ener,ijkcb0) - fr*beta
305 
306 ! ------- store rad. flux normal to surf. in positive direction (hence sgn* )
307  modsf = sqrt( sf(1)*sf(1) + sf(2)*sf(2) + sf(3)*sf(3) )
308  qr(ijknb) = -sgn*fr/modsf
309  ENDDO ! i
310  ENDDO ! j
311  ENDDO ! k
312 
313 ! non-conforming region interface
314 
315  ELSE IF (bctype>=bc_regionint .AND. bctype<=bc_regionint +bc_range) THEN
316 
317  ELSE IF (bctype>=bc_regnonconf .AND. bctype<=bc_regnonconf+bc_range) THEN
318 
319 ! everything else
320 
321  ELSE
322 
323  DO k=kbeg,kend
324  DO j=jbeg,jend
325  DO i=ibeg,iend
326  ijkcb0 = indijk(i ,j ,k ,icoff,ijcoff) ! boundary
327  ijkcd = indijk(i-idir ,j-jdir ,k-kdir ,icoff,ijcoff) ! dummy
328  ijknb = indijk(i+inode,j+jnode,k+knode,inoff,ijnoff)
329 
330  gfaca = 0.5_rfreal*(gofact(ijkcb0)+gofact(ijkcd))
331  flima = 0.5_rfreal*(flim(ijkcb0)+flim(ijkcd))
332  tempa = 0.5_rfreal*(rdv(dv_radi_teff,ijkcb0) + &
333  rdv(dv_radi_teff,ijkcd))
334  coefa = 0.5_rfreal*(coef(ijkcb0,radi_coeff_extinct) + &
335  coef(ijkcd ,radi_coeff_extinct))
336 
337  coefa = gfaca*rati*flima*stboltz*tempa**3/coefa
338 
339  qrx = coefa*grad(gr_mixt_tx,ijknb)
340  qry = coefa*grad(gr_mixt_ty,ijknb)
341  qrz = coefa*grad(gr_mixt_tz,ijknb)
342 
343  sf(1) = sgn*sface(xcoord,ijknb)
344  sf(2) = sgn*sface(ycoord,ijknb)
345  sf(3) = sgn*sface(zcoord,ijknb)
346 
347  fr = qrx*sf(1) + qry*sf(2) + qrz*sf(3)
348 
349  rhs(cv_mixt_ener,ijkcb0) = rhs(cv_mixt_ener,ijkcb0) - fr*beta
350 
351 ! ------- store rad. flux normal to surf. in positive direction (hence sgn* )
352  modsf = sqrt( sf(1)*sf(1) + sf(2)*sf(2) + sf(3)*sf(3) )
353  qr(ijknb) = -sgn*fr/modsf
354  ENDDO ! i
355  ENDDO ! j
356  ENDDO ! k
357 
358  ENDIF ! bcType
359 
360 ! add contribution of grad(T) from patch for intensity calculation later
361 
362  IF (lbound==1 .OR. lbound==2) THEN
363  idbeg = ibeg - idir*region%nDumCells
364  idend = iend - idir*region%nDumCells
365  jdbeg = jbeg - region%nDumCells
366  jdend = jend + region%nDumCells
367  kdbeg = kbeg - region%nDumCells
368  kdend = kend + region%nDumCells
369 
370  ELSE IF (lbound==3 .OR. lbound==4) THEN
371  jdbeg = jbeg - jdir*region%nDumCells
372  jdend = jend - jdir*region%nDumCells
373  idbeg = ibeg - region%nDumCells
374  idend = iend + region%nDumCells
375  kdbeg = kbeg - region%nDumCells
376  kdend = kend + region%nDumCells
377 
378  ELSE IF (lbound==5 .OR. lbound==6) THEN
379  kdbeg = kbeg - kdir*region%nDumCells
380  kdend = kend - kdir*region%nDumCells
381  jdbeg = jbeg - region%nDumCells
382  jdend = jend + region%nDumCells
383  idbeg = ibeg - region%nDumCells
384  idend = iend + region%nDumCells
385  ENDIF
386 
387  DO k=kdbeg,kdend
388  DO j=jdbeg,jdend
389  DO i=idbeg,idend
390  ijkcb0 = indijk(i ,j ,k ,icoff,ijcoff) ! boundary
391  ijknb = indijk(i+inode,j+jnode,k+knode,inoff,ijnoff)
392 
393  wvint(xcoord,ijkcb0) = wvint(xcoord,ijkcb0) + grad(gr_mixt_tx,ijknb)
394  wvint(ycoord,ijkcb0) = wvint(ycoord,ijkcb0) + grad(gr_mixt_ty,ijknb)
395  wvint(zcoord,ijkcb0) = wvint(zcoord,ijkcb0) + grad(gr_mixt_tz,ijknb)
396  ENDDO ! i
397  ENDDO ! j
398  ENDDO ! k
399 
400 ! finalize --------------------------------------------------------------------
401 
402  CALL deregisterfunction( region%global )
403 
404 END SUBROUTINE radi_diffradfluxpatch
405 
406 !******************************************************************************
407 !
408 ! RCS Revision history:
409 !
410 ! $Log: RADI_DiffRadFluxPatch.F90,v $
411 ! Revision 1.4 2008/12/06 08:44:37 mtcampbe
412 ! Updated license.
413 !
414 ! Revision 1.3 2008/11/19 22:17:49 mtcampbe
415 ! Added Illinois Open Source License/Copyright
416 !
417 ! Revision 1.2 2006/08/19 15:40:11 mparmar
418 ! Renamed patch variables
419 !
420 ! Revision 1.1 2004/09/22 02:35:49 wasistho
421 ! changed file nomenclature from lower to upper case
422 !
423 ! Revision 1.5 2004/09/22 01:32:10 wasistho
424 ! switch LFD to FLD for flux limited diffusion
425 !
426 ! Revision 1.4 2004/09/18 17:41:54 wasistho
427 ! install Limited Flux Diffusion radiation
428 !
429 ! Revision 1.3 2003/08/13 21:18:12 wasistho
430 ! corrected sign and rk-stage of radiation flux
431 !
432 ! Revision 1.2 2003/07/30 22:24:34 wasistho
433 ! enter part and smoke data into radiation
434 !
435 ! Revision 1.1 2003/07/17 01:16:59 wasistho
436 ! initial activation rocrad
437 !
438 !
439 !******************************************************************************
440 
441 
442 
443 
444 
445 
446 
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
rflo_getpatchdirection
subroutine rflo_getpatchdirection(patch, idir, jdir, kdir)
Definition: RFLO_GetPatchDimens.F90:215
inode
**********************************************************************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 inode
Definition: RFLO_CalcGradDummySymm.h:32
k
j indices k indices k
Definition: Indexing.h:6
rhs
NT rhs
Definition: rational_rotation.h:86
moderror::registerfunction
subroutine registerfunction(global, funName, fileName)
Definition: ModError.F90:449
rflo_getpatchindices
subroutine rflo_getpatchindices(region, patch, iLev, ibeg, iend, jbeg, jend, kbeg, kend)
Definition: RFLO_GetPatchDimens.F90:61
sqrt
double sqrt(double d)
Definition: double.h:73
rflo_getnodeoffset
subroutine rflo_getnodeoffset(region, iLev, iNodeOffset, ijNodeOffset)
Definition: RFLO_GetDimensions.F90:267
jdir
**********************************************************************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 jdir
Definition: RFLO_CalcGradDummyConn.h:30
patch
Definition: patch.h:74
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
idir
**********************************************************************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 idir
Definition: RFLO_CalcGradDummyConn.h:30
rflo_getcelloffset
subroutine rflo_getcelloffset(region, iLev, iCellOffset, ijCellOffset)
Definition: RFLO_GetDimensions.F90:235
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
kdir
**********************************************************************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 kdir
Definition: RFLO_CalcGradDummySymm.h:32
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
jnode
**********************************************************************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 jnode
Definition: RFLO_CalcGradDummySymm.h:32
modbndpatch::t_patch
Definition: ModBndPatch.F90:91
radi_diffradfluxpatch
subroutine radi_diffradfluxpatch(region, patch)
Definition: RADI_DiffRadFluxPatch.F90:51
radi_modparameters
Definition: RADI_ModParameters.F90:39