Rocstar  1.0
Rocstar multiphysics simulation application
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
RFLO_CentralDissipation.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 artificial numerical dissipation (JST type).
26 !
27 ! Description: the dissipation consists of a blend of 2nd- and 4th-order
28 ! differences (undivided).
29 !
30 ! Input: region = data of current region.
31 !
32 ! Output: region%levels%mixt%diss = dissipative fluxes.
33 !
34 ! Notes: the TVD type of pressure switch is taken from:
35 ! Turkel, E.; Swanson, R.C.; Vatsa, V.N.; White, J.A.: Multigrid
36 ! for Hypersonic Viscous Two- and Three-Dimensional Flows. AIAA
37 ! Paper 91-1572, 1991.
38 !
39 !******************************************************************************
40 !
41 ! $Id: RFLO_CentralDissipation.F90,v 1.3 2008/12/06 08:44:26 mtcampbe Exp $
42 !
43 ! Copyright: (c) 2001 by the University of Illinois
44 !
45 !******************************************************************************
46 
47 SUBROUTINE rflo_centraldissipation( region )
48 
49  USE moddatatypes
50  USE moddatastruct, ONLY : t_region
51  USE modinterfaces, ONLY : rflo_getdimensphys, rflo_getcelloffset
52  USE moderror
53  USE modparameters
54  IMPLICIT NONE
55 
56 #include "Indexing.h"
57 
58 ! ... parameters
59  TYPE(t_region) :: region
60 
61 ! ... loop variables
62  INTEGER :: i, j, k, ii, jj, kk
63 
64 ! ... local variables
65  INTEGER :: ipcbeg, ipcend, jpcbeg, jpcend, kpcbeg, kpcend
66  INTEGER :: ilev, icoff, ijcoff, ijkc0, ijkcm1, ijkcp1, ijkcp2, pswitchtype
67 
68  REAL(RFREAL) :: beta, vis2, vis4, eval, pmax, eps2, eps4, fd(5)
69  REAL(RFREAL) :: pswitchomega, ptvd, psum
70  REAL(RFREAL), POINTER :: cv(:,:), dv(:,:), diss(:,:), srad(:,:), dp(:)
71 
72 !******************************************************************************
73 
74  CALL registerfunction( region%global,'RFLO_CentralDissipation',&
75  'RFLO_CentralDissipation.F90' )
76 
77 ! get dimensions and pointers -------------------------------------------------
78 
79  ilev = region%currLevel
80 
81  CALL rflo_getdimensphys( region,ilev,ipcbeg,ipcend, &
82  jpcbeg,jpcend,kpcbeg,kpcend )
83  CALL rflo_getcelloffset( region,ilev,icoff,ijcoff )
84 
85  cv => region%levels(ilev)%mixt%cv
86  dv => region%levels(ilev)%mixt%dv
87  diss => region%levels(ilev)%mixt%diss
88  srad => region%levels(ilev)%mixt%srad
89  dp => region%work1D
90 
91 ! get coefficients and switch type --------------------------------------------
92 
93  beta = region%mixtInput%betrk(region%irkStep)
94  vis2 = beta*region%mixtInput%vis2
95  vis4 = beta*region%mixtInput%vis4
96  pswitchtype = region%mixtInput%pSwitchType
97  pswitchomega = region%mixtInput%pSwitchOmega
98 
99 ! dissipation in i-direction --------------------------------------------------
100 
101  DO k=kpcbeg,kpcend
102  DO j=jpcbeg,jpcend
103 
104 ! --- pressure switch
105 
106  ii = 0
107  IF (pswitchtype == pswitch_std) THEN
108  DO i=ipcbeg-1,ipcend+1
109  ii = ii + 1
110  ijkc0 = indijk(i ,j,k,icoff,ijcoff)
111  ijkcm1 = indijk(i-1,j,k,icoff,ijcoff)
112  ijkcp1 = indijk(i+1,j,k,icoff,ijcoff)
113  dp(ii) = abs(( dv(dv_mixt_pres,ijkcp1)- &
114  2._rfreal*dv(dv_mixt_pres,ijkc0 )+ &
115  dv(dv_mixt_pres,ijkcm1))/ &
116  ( dv(dv_mixt_pres,ijkcp1)+ &
117  2._rfreal*dv(dv_mixt_pres,ijkc0 )+ &
118  dv(dv_mixt_pres,ijkcm1)))
119  ENDDO
120  ELSE
121  DO i=ipcbeg-1,ipcend+1
122  ii = ii + 1
123  ijkc0 = indijk(i ,j,k,icoff,ijcoff)
124  ijkcm1 = indijk(i-1,j,k,icoff,ijcoff)
125  ijkcp1 = indijk(i+1,j,k,icoff,ijcoff)
126  psum = dv(dv_mixt_pres,ijkcp1)+ &
127  2._rfreal*dv(dv_mixt_pres,ijkc0 )+ &
128  dv(dv_mixt_pres,ijkcm1)
129  ptvd = abs(dv(dv_mixt_pres,ijkcp1)-dv(dv_mixt_pres,ijkc0 ))+ &
130  abs(dv(dv_mixt_pres,ijkc0 )-dv(dv_mixt_pres,ijkcm1))
131  dp(ii) = abs( dv(dv_mixt_pres,ijkcp1)- &
132  2._rfreal*dv(dv_mixt_pres,ijkc0 )+ &
133  dv(dv_mixt_pres,ijkcm1))/ &
134  ((1._rfreal-pswitchomega)*ptvd+pswitchomega*psum)
135  ENDDO
136  ENDIF
137 
138 ! --- dissipative fluxes at I+1/2
139 
140  ii = 0
141  DO i=ipcbeg-1,ipcend
142  ii = ii + 1
143  ijkc0 = indijk(i ,j,k,icoff,ijcoff)
144  ijkcm1 = indijk(i-1,j,k,icoff,ijcoff)
145  ijkcp1 = indijk(i+1,j,k,icoff,ijcoff)
146  ijkcp2 = indijk(i+2,j,k,icoff,ijcoff)
147  eval = 0.5_rfreal*(srad(icoord,ijkc0)+srad(icoord,ijkcp1) + &
148  max(srad(jcoord,ijkc0)+srad(jcoord,ijkcp1), &
149  srad(kcoord,ijkc0)+srad(kcoord,ijkcp1)))
150  pmax = max(dp(ii),dp(ii+1))
151  eps2 = eval*vis2*pmax
152  eps4 = eval*vis4
153  eps4 = dim(eps4,eps2)
154  fd(1) = eps2*(cv(cv_mixt_dens,ijkcp1)-cv(cv_mixt_dens,ijkc0)) - &
155  eps4*(cv(cv_mixt_dens,ijkcp2)- &
156  3._rfreal*(cv(cv_mixt_dens,ijkcp1)- &
157  cv(cv_mixt_dens,ijkc0 ))- &
158  cv(cv_mixt_dens,ijkcm1))
159  fd(2) = eps2*(cv(cv_mixt_xmom,ijkcp1)-cv(cv_mixt_xmom,ijkc0)) - &
160  eps4*(cv(cv_mixt_xmom,ijkcp2)- &
161  3._rfreal*(cv(cv_mixt_xmom,ijkcp1)- &
162  cv(cv_mixt_xmom,ijkc0 ))- &
163  cv(cv_mixt_xmom,ijkcm1))
164  fd(3) = eps2*(cv(cv_mixt_ymom,ijkcp1)-cv(cv_mixt_ymom,ijkc0)) - &
165  eps4*(cv(cv_mixt_ymom,ijkcp2)- &
166  3._rfreal*(cv(cv_mixt_ymom,ijkcp1)- &
167  cv(cv_mixt_ymom,ijkc0 ))- &
168  cv(cv_mixt_ymom,ijkcm1))
169  fd(4) = eps2*(cv(cv_mixt_zmom,ijkcp1)-cv(cv_mixt_zmom,ijkc0)) - &
170  eps4*(cv(cv_mixt_zmom,ijkcp2)- &
171  3._rfreal*(cv(cv_mixt_zmom,ijkcp1)- &
172  cv(cv_mixt_zmom,ijkc0 ))- &
173  cv(cv_mixt_zmom,ijkcm1))
174  fd(5) = eps2*(cv(cv_mixt_ener,ijkcp1)-cv(cv_mixt_ener,ijkc0)) - &
175  eps4*(cv(cv_mixt_ener,ijkcp2)- &
176  3._rfreal*(cv(cv_mixt_ener,ijkcp1)- &
177  cv(cv_mixt_ener,ijkc0 ))- &
178  cv(cv_mixt_ener,ijkcm1))
179 
180  diss(cv_mixt_dens,ijkc0 ) = diss(cv_mixt_dens,ijkc0 ) + fd(1)
181  diss(cv_mixt_xmom,ijkc0 ) = diss(cv_mixt_xmom,ijkc0 ) + fd(2)
182  diss(cv_mixt_ymom,ijkc0 ) = diss(cv_mixt_ymom,ijkc0 ) + fd(3)
183  diss(cv_mixt_zmom,ijkc0 ) = diss(cv_mixt_zmom,ijkc0 ) + fd(4)
184  diss(cv_mixt_ener,ijkc0 ) = diss(cv_mixt_ener,ijkc0 ) + fd(5)
185 
186  diss(cv_mixt_dens,ijkcp1) = diss(cv_mixt_dens,ijkcp1) - fd(1)
187  diss(cv_mixt_xmom,ijkcp1) = diss(cv_mixt_xmom,ijkcp1) - fd(2)
188  diss(cv_mixt_ymom,ijkcp1) = diss(cv_mixt_ymom,ijkcp1) - fd(3)
189  diss(cv_mixt_zmom,ijkcp1) = diss(cv_mixt_zmom,ijkcp1) - fd(4)
190  diss(cv_mixt_ener,ijkcp1) = diss(cv_mixt_ener,ijkcp1) - fd(5)
191  ENDDO ! i
192 
193  ENDDO ! j
194  ENDDO ! k
195 
196 ! dissipation in j-direction --------------------------------------------------
197 
198  DO k=kpcbeg,kpcend
199  DO i=ipcbeg,ipcend
200 
201 ! --- pressure switch
202 
203  jj = 0
204  IF (pswitchtype == pswitch_std) THEN
205  DO j=jpcbeg-1,jpcend+1
206  jj = jj + 1
207  ijkc0 = indijk(i,j ,k,icoff,ijcoff)
208  ijkcm1 = indijk(i,j-1,k,icoff,ijcoff)
209  ijkcp1 = indijk(i,j+1,k,icoff,ijcoff)
210  dp(jj) = abs(( dv(dv_mixt_pres,ijkcp1)- &
211  2._rfreal*dv(dv_mixt_pres,ijkc0 )+ &
212  dv(dv_mixt_pres,ijkcm1))/ &
213  ( dv(dv_mixt_pres,ijkcp1)+ &
214  2._rfreal*dv(dv_mixt_pres,ijkc0 )+ &
215  dv(dv_mixt_pres,ijkcm1)))
216  ENDDO
217  ELSE
218  DO j=jpcbeg-1,jpcend+1
219  jj = jj + 1
220  ijkc0 = indijk(i,j ,k,icoff,ijcoff)
221  ijkcm1 = indijk(i,j-1,k,icoff,ijcoff)
222  ijkcp1 = indijk(i,j+1,k,icoff,ijcoff)
223  psum = dv(dv_mixt_pres,ijkcp1)+ &
224  2._rfreal*dv(dv_mixt_pres,ijkc0 )+ &
225  dv(dv_mixt_pres,ijkcm1)
226  ptvd = abs(dv(dv_mixt_pres,ijkcp1)-dv(dv_mixt_pres,ijkc0 ))+ &
227  abs(dv(dv_mixt_pres,ijkc0 )-dv(dv_mixt_pres,ijkcm1))
228  dp(jj) = abs( dv(dv_mixt_pres,ijkcp1)- &
229  2._rfreal*dv(dv_mixt_pres,ijkc0 )+ &
230  dv(dv_mixt_pres,ijkcm1))/ &
231  ((1._rfreal-pswitchomega)*ptvd+pswitchomega*psum)
232  ENDDO
233  ENDIF
234 
235 ! --- dissipative fluxes at J+1/2
236 
237  jj = 0
238  DO j=jpcbeg-1,jpcend
239  jj = jj + 1
240  ijkc0 = indijk(i,j ,k,icoff,ijcoff)
241  ijkcm1 = indijk(i,j-1,k,icoff,ijcoff)
242  ijkcp1 = indijk(i,j+1,k,icoff,ijcoff)
243  ijkcp2 = indijk(i,j+2,k,icoff,ijcoff)
244  eval = 0.5_rfreal*(srad(jcoord,ijkc0)+srad(jcoord,ijkcp1) + &
245  max(srad(icoord,ijkc0)+srad(icoord,ijkcp1), &
246  srad(kcoord,ijkc0)+srad(kcoord,ijkcp1)))
247  pmax = max(dp(jj),dp(jj+1))
248  eps2 = eval*vis2*pmax
249  eps4 = eval*vis4
250  eps4 = dim(eps4,eps2)
251  fd(1) = eps2*(cv(cv_mixt_dens,ijkcp1)-cv(cv_mixt_dens,ijkc0)) - &
252  eps4*(cv(cv_mixt_dens,ijkcp2)- &
253  3._rfreal*(cv(cv_mixt_dens,ijkcp1)- &
254  cv(cv_mixt_dens,ijkc0 ))- &
255  cv(cv_mixt_dens,ijkcm1))
256  fd(2) = eps2*(cv(cv_mixt_xmom,ijkcp1)-cv(cv_mixt_xmom,ijkc0)) - &
257  eps4*(cv(cv_mixt_xmom,ijkcp2)- &
258  3._rfreal*(cv(cv_mixt_xmom,ijkcp1)- &
259  cv(cv_mixt_xmom,ijkc0 ))- &
260  cv(cv_mixt_xmom,ijkcm1))
261  fd(3) = eps2*(cv(cv_mixt_ymom,ijkcp1)-cv(cv_mixt_ymom,ijkc0)) - &
262  eps4*(cv(cv_mixt_ymom,ijkcp2)- &
263  3._rfreal*(cv(cv_mixt_ymom,ijkcp1)- &
264  cv(cv_mixt_ymom,ijkc0 ))- &
265  cv(cv_mixt_ymom,ijkcm1))
266  fd(4) = eps2*(cv(cv_mixt_zmom,ijkcp1)-cv(cv_mixt_zmom,ijkc0)) - &
267  eps4*(cv(cv_mixt_zmom,ijkcp2)- &
268  3._rfreal*(cv(cv_mixt_zmom,ijkcp1)- &
269  cv(cv_mixt_zmom,ijkc0 ))- &
270  cv(cv_mixt_zmom,ijkcm1))
271  fd(5) = eps2*(cv(cv_mixt_ener,ijkcp1)-cv(cv_mixt_ener,ijkc0)) - &
272  eps4*(cv(cv_mixt_ener,ijkcp2)- &
273  3._rfreal*(cv(cv_mixt_ener,ijkcp1)- &
274  cv(cv_mixt_ener,ijkc0 ))- &
275  cv(cv_mixt_ener,ijkcm1))
276 
277  diss(cv_mixt_dens,ijkc0 ) = diss(cv_mixt_dens,ijkc0 ) + fd(1)
278  diss(cv_mixt_xmom,ijkc0 ) = diss(cv_mixt_xmom,ijkc0 ) + fd(2)
279  diss(cv_mixt_ymom,ijkc0 ) = diss(cv_mixt_ymom,ijkc0 ) + fd(3)
280  diss(cv_mixt_zmom,ijkc0 ) = diss(cv_mixt_zmom,ijkc0 ) + fd(4)
281  diss(cv_mixt_ener,ijkc0 ) = diss(cv_mixt_ener,ijkc0 ) + fd(5)
282 
283  diss(cv_mixt_dens,ijkcp1) = diss(cv_mixt_dens,ijkcp1) - fd(1)
284  diss(cv_mixt_xmom,ijkcp1) = diss(cv_mixt_xmom,ijkcp1) - fd(2)
285  diss(cv_mixt_ymom,ijkcp1) = diss(cv_mixt_ymom,ijkcp1) - fd(3)
286  diss(cv_mixt_zmom,ijkcp1) = diss(cv_mixt_zmom,ijkcp1) - fd(4)
287  diss(cv_mixt_ener,ijkcp1) = diss(cv_mixt_ener,ijkcp1) - fd(5)
288  ENDDO ! j
289 
290  ENDDO ! i
291  ENDDO ! k
292 
293 ! dissipation in k-direction --------------------------------------------------
294 
295  DO j=jpcbeg,jpcend
296  DO i=ipcbeg,ipcend
297 
298 ! --- pressure switch
299 
300  kk = 0
301  IF (pswitchtype == pswitch_std) THEN
302  DO k=kpcbeg-1,kpcend+1
303  kk = kk + 1
304  ijkc0 = indijk(i,j,k ,icoff,ijcoff)
305  ijkcm1 = indijk(i,j,k-1,icoff,ijcoff)
306  ijkcp1 = indijk(i,j,k+1,icoff,ijcoff)
307  dp(kk) = abs(( dv(dv_mixt_pres,ijkcp1)- &
308  2._rfreal*dv(dv_mixt_pres,ijkc0 )+ &
309  dv(dv_mixt_pres,ijkcm1))/ &
310  ( dv(dv_mixt_pres,ijkcp1)+ &
311  2._rfreal*dv(dv_mixt_pres,ijkc0 )+ &
312  dv(dv_mixt_pres,ijkcm1)))
313  ENDDO
314  ELSE
315  DO k=kpcbeg-1,kpcend+1
316  kk = kk + 1
317  ijkc0 = indijk(i,j,k ,icoff,ijcoff)
318  ijkcm1 = indijk(i,j,k-1,icoff,ijcoff)
319  ijkcp1 = indijk(i,j,k+1,icoff,ijcoff)
320  psum = dv(dv_mixt_pres,ijkcp1)+ &
321  2._rfreal*dv(dv_mixt_pres,ijkc0 )+ &
322  dv(dv_mixt_pres,ijkcm1)
323  ptvd = abs(dv(dv_mixt_pres,ijkcp1)-dv(dv_mixt_pres,ijkc0 ))+ &
324  abs(dv(dv_mixt_pres,ijkc0 )-dv(dv_mixt_pres,ijkcm1))
325  dp(kk) = abs( dv(dv_mixt_pres,ijkcp1)- &
326  2._rfreal*dv(dv_mixt_pres,ijkc0 )+ &
327  dv(dv_mixt_pres,ijkcm1))/ &
328  ((1._rfreal-pswitchomega)*ptvd+pswitchomega*psum)
329  ENDDO
330  ENDIF
331 
332 ! --- dissipative fluxes at K+1/2
333 
334  kk = 0
335  DO k=kpcbeg-1,kpcend
336  kk = kk + 1
337  ijkc0 = indijk(i,j,k ,icoff,ijcoff)
338  ijkcm1 = indijk(i,j,k-1,icoff,ijcoff)
339  ijkcp1 = indijk(i,j,k+1,icoff,ijcoff)
340  ijkcp2 = indijk(i,j,k+2,icoff,ijcoff)
341  eval = 0.5_rfreal*(srad(kcoord,ijkc0)+srad(kcoord,ijkcp1) + &
342  max(srad(icoord,ijkc0)+srad(icoord,ijkcp1), &
343  srad(jcoord,ijkc0)+srad(jcoord,ijkcp1)))
344  pmax = max(dp(kk),dp(kk+1))
345  eps2 = eval*vis2*pmax
346  eps4 = eval*vis4
347  eps4 = dim(eps4,eps2)
348  fd(1) = eps2*(cv(cv_mixt_dens,ijkcp1)-cv(cv_mixt_dens,ijkc0)) - &
349  eps4*(cv(cv_mixt_dens,ijkcp2)- &
350  3._rfreal*(cv(cv_mixt_dens,ijkcp1)- &
351  cv(cv_mixt_dens,ijkc0 ))- &
352  cv(cv_mixt_dens,ijkcm1))
353  fd(2) = eps2*(cv(cv_mixt_xmom,ijkcp1)-cv(cv_mixt_xmom,ijkc0)) - &
354  eps4*(cv(cv_mixt_xmom,ijkcp2)- &
355  3._rfreal*(cv(cv_mixt_xmom,ijkcp1)- &
356  cv(cv_mixt_xmom,ijkc0 ))- &
357  cv(cv_mixt_xmom,ijkcm1))
358  fd(3) = eps2*(cv(cv_mixt_ymom,ijkcp1)-cv(cv_mixt_ymom,ijkc0)) - &
359  eps4*(cv(cv_mixt_ymom,ijkcp2)- &
360  3._rfreal*(cv(cv_mixt_ymom,ijkcp1)- &
361  cv(cv_mixt_ymom,ijkc0 ))- &
362  cv(cv_mixt_ymom,ijkcm1))
363  fd(4) = eps2*(cv(cv_mixt_zmom,ijkcp1)-cv(cv_mixt_zmom,ijkc0)) - &
364  eps4*(cv(cv_mixt_zmom,ijkcp2)- &
365  3._rfreal*(cv(cv_mixt_zmom,ijkcp1)- &
366  cv(cv_mixt_zmom,ijkc0 ))- &
367  cv(cv_mixt_zmom,ijkcm1))
368  fd(5) = eps2*(cv(cv_mixt_ener,ijkcp1)-cv(cv_mixt_ener,ijkc0)) - &
369  eps4*(cv(cv_mixt_ener,ijkcp2)- &
370  3._rfreal*(cv(cv_mixt_ener,ijkcp1)- &
371  cv(cv_mixt_ener,ijkc0 ))- &
372  cv(cv_mixt_ener,ijkcm1))
373 
374  diss(cv_mixt_dens,ijkc0 ) = diss(cv_mixt_dens,ijkc0 ) + fd(1)
375  diss(cv_mixt_xmom,ijkc0 ) = diss(cv_mixt_xmom,ijkc0 ) + fd(2)
376  diss(cv_mixt_ymom,ijkc0 ) = diss(cv_mixt_ymom,ijkc0 ) + fd(3)
377  diss(cv_mixt_zmom,ijkc0 ) = diss(cv_mixt_zmom,ijkc0 ) + fd(4)
378  diss(cv_mixt_ener,ijkc0 ) = diss(cv_mixt_ener,ijkc0 ) + fd(5)
379 
380  diss(cv_mixt_dens,ijkcp1) = diss(cv_mixt_dens,ijkcp1) - fd(1)
381  diss(cv_mixt_xmom,ijkcp1) = diss(cv_mixt_xmom,ijkcp1) - fd(2)
382  diss(cv_mixt_ymom,ijkcp1) = diss(cv_mixt_ymom,ijkcp1) - fd(3)
383  diss(cv_mixt_zmom,ijkcp1) = diss(cv_mixt_zmom,ijkcp1) - fd(4)
384  diss(cv_mixt_ener,ijkcp1) = diss(cv_mixt_ener,ijkcp1) - fd(5)
385  ENDDO ! k
386 
387  ENDDO ! i
388  ENDDO ! j
389 
390 ! finalize --------------------------------------------------------------------
391 
392  CALL deregisterfunction( region%global )
393 
394 END SUBROUTINE rflo_centraldissipation
395 
396 !******************************************************************************
397 !
398 ! RCS Revision history:
399 !
400 ! $Log: RFLO_CentralDissipation.F90,v $
401 ! Revision 1.3 2008/12/06 08:44:26 mtcampbe
402 ! Updated license.
403 !
404 ! Revision 1.2 2008/11/19 22:17:37 mtcampbe
405 ! Added Illinois Open Source License/Copyright
406 !
407 ! Revision 1.1 2004/11/29 20:51:38 wasistho
408 ! lower to upper case
409 !
410 ! Revision 1.10 2003/11/20 16:40:38 mdbrandy
411 ! Backing out RocfluidMP changes from 11-17-03
412 !
413 ! Revision 1.6 2003/05/15 02:57:03 jblazek
414 ! Inlined index function.
415 !
416 ! Revision 1.5 2002/09/05 17:40:21 jblazek
417 ! Variable global moved into regions().
418 !
419 ! Revision 1.4 2002/07/25 00:36:48 jblazek
420 ! Option for TVD type pressure switch.
421 !
422 ! Revision 1.3 2002/02/21 23:25:05 jblazek
423 ! Blocks renamed as regions.
424 !
425 ! Revision 1.2 2002/01/28 23:55:22 jblazek
426 ! Added flux computation (central scheme).
427 !
428 ! Revision 1.1 2002/01/23 03:51:25 jblazek
429 ! Added low-level time-stepping routines.
430 !
431 !******************************************************************************
432 
433 
434 
435 
436 
437 
438 
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
max
Vector_n max(const Array_n_const &v1, const Array_n_const &v2)
Definition: Vector_n.h:354
rflo_centraldissipation
subroutine rflo_centraldissipation(region)
Definition: RFLO_CentralDissipation.F90:47
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
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
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
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
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