Rocstar  1.0
Rocstar multiphysics simulation application
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
RFLO_RoeDissipFirst.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 numerical dissipation based on 1st-order
26 ! Roe`s upwind scheme.
27 !
28 ! Description: none.
29 !
30 ! Input: region = data of current region.
31 !
32 ! Output: region%levels%mixt%diss = dissipative fluxes.
33 !
34 ! Notes: none.
35 !
36 !******************************************************************************
37 !
38 ! $Id: RFLO_RoeDissipFirst.F90,v 1.3 2008/12/06 08:44:27 mtcampbe Exp $
39 !
40 ! Copyright: (c) 2003 by the University of Illinois
41 !
42 !******************************************************************************
43 
44 SUBROUTINE rflo_roedissipfirst( region )
45 
46  USE moddatatypes
47  USE moddatastruct, ONLY : t_region
48  USE modinterfaces, ONLY : rflo_getdimensphys, rflo_getcelloffset, &
49  rflo_getnodeoffset
50  USE moderror
51  USE modparameters
52  IMPLICIT NONE
53 
54 #include "Indexing.h"
55 
56 ! ... parameters
57  TYPE(t_region) :: region
58 
59 ! ... loop variables
60  INTEGER :: i, j, k
61 
62 ! ... local variables
63  INTEGER :: ipcbeg, ipcend, jpcbeg, jpcend, kpcbeg, kpcend
64  INTEGER :: ilev, icoff, ijcoff, inoff, ijnoff, ijkc0, ijkcm1, ijkn
65  INTEGER :: indcp, indmol, indsvel
66 
67  REAL(RFREAL) :: beta5, ds, nx, ny, nz, svel, rgas, gaml, gamr, ggm1, gam1, &
68  rl, ul, vl, wl, pl, hl, rr, ur, vr, wr, pr, hr, rav, dd, &
69  dd1, uav, vav, wav, hav, q2a, c2a, cav, uvw, du, eabs1, &
70  eabs2, eabs5, h1, h2, h3, h4, h5, epsentr, delta, fd(5)
71  REAL(RFREAL), POINTER :: cv(:,:), dv(:,:), gv(:,:), diss(:,:), &
72  si(:,:), sj(:,:), sk(:,:), &
73  sivel(:), sjvel(:), skvel(:)
74 
75 !******************************************************************************
76 
77  CALL registerfunction( region%global,'RFLO_RoeDissipFirst',&
78  'RFLO_RoeDissipFirst.F90' )
79 
80 ! get dimensions and pointers -------------------------------------------------
81 
82  ilev = region%currLevel
83 
84  CALL rflo_getdimensphys( region,ilev,ipcbeg,ipcend, &
85  jpcbeg,jpcend,kpcbeg,kpcend )
86  CALL rflo_getcelloffset( region,ilev,icoff,ijcoff )
87  CALL rflo_getnodeoffset( region,ilev,inoff,ijnoff )
88 
89  cv => region%levels(ilev)%mixt%cv
90  dv => region%levels(ilev)%mixt%dv
91  gv => region%levels(ilev)%mixt%gv
92  si => region%levels(ilev)%grid%si
93  sj => region%levels(ilev)%grid%sj
94  sk => region%levels(ilev)%grid%sk
95  sivel => region%levels(ilev)%grid%siVel
96  sjvel => region%levels(ilev)%grid%sjVel
97  skvel => region%levels(ilev)%grid%skVel
98  diss => region%levels(ilev)%mixt%diss
99 
100 ! get coefficients ------------------------------------------------------------
101 
102  beta5 = 0.5_rfreal*region%mixtInput%betrk(region%irkStep)
103  epsentr = region%mixtInput%epsEntr
104  indsvel = region%levels(ilev)%grid%indSvel
105  indcp = region%levels(ilev)%mixt%indCp
106  indmol = region%levels(ilev)%mixt%indMol
107 
108 ! dissipation in i-direction --------------------------------------------------
109 
110  DO k=kpcbeg,kpcend
111  DO j=jpcbeg,jpcend
112  DO i=ipcbeg,ipcend+1
113  ijkc0 = indijk(i ,j,k,icoff,ijcoff)
114  ijkcm1 = indijk(i-1,j,k,icoff,ijcoff)
115  ijkn = indijk(i ,j,k,inoff,ijnoff)
116  ds = sqrt(si(xcoord,ijkn)*si(xcoord,ijkn)+ &
117  si(ycoord,ijkn)*si(ycoord,ijkn)+ &
118  si(zcoord,ijkn)*si(zcoord,ijkn))
119  nx = si(xcoord,ijkn)/ds
120  ny = si(ycoord,ijkn)/ds
121  nz = si(zcoord,ijkn)/ds
122  svel = sivel(ijkn*indsvel)/ds
123  ds = ds*beta5
124 
125 ! ----- left and right states
126 
127  rl = cv(cv_mixt_dens,ijkcm1)
128  ul = dv(dv_mixt_uvel,ijkcm1)
129  vl = dv(dv_mixt_vvel,ijkcm1)
130  wl = dv(dv_mixt_wvel,ijkcm1)
131  pl = dv(dv_mixt_pres,ijkcm1)
132  rgas = 8314.3_rfreal/gv(gv_mixt_mol,ijkcm1*indmol)
133  gaml = gv(gv_mixt_cp,ijkcm1*indcp)/(gv(gv_mixt_cp,ijkcm1*indcp)-rgas)
134  ggm1 = gaml/(gaml-1._rfreal)
135  hl = ggm1*pl/rl + 0.5_rfreal*(ul*ul+vl*vl+wl*wl)
136 
137  rr = cv(cv_mixt_dens,ijkc0)
138  ur = dv(dv_mixt_uvel,ijkc0)
139  vr = dv(dv_mixt_vvel,ijkc0)
140  wr = dv(dv_mixt_wvel,ijkc0)
141  pr = dv(dv_mixt_pres,ijkc0)
142  rgas = 8314.3_rfreal/gv(gv_mixt_mol,ijkc0*indmol)
143  gamr = gv(gv_mixt_cp,ijkc0*indcp)/(gv(gv_mixt_cp,ijkc0*indcp)-rgas)
144  ggm1 = gamr/(gamr-1._rfreal)
145  hr = ggm1*pr/rr + 0.5_rfreal*(ur*ur+vr*vr+wr*wr)
146 
147 ! ----- Roe`s average
148 
149  rav = sqrt(rl*rr)
150  dd = rav/rl
151  dd1 = 1._rfreal/(1._rfreal+dd)
152  uav = (ul+dd*ur)*dd1
153  vav = (vl+dd*vr)*dd1
154  wav = (wl+dd*wr)*dd1
155  hav = (hl+dd*hr)*dd1
156  q2a = 0.5_rfreal*(uav*uav+vav*vav+wav*wav)
157  gam1 = 0.5_rfreal*(gaml+gamr) - 1._rfreal
158  c2a = gam1*(hav-q2a)
159  cav = sqrt(c2a)
160  uvw = uav*nx + vav*ny + wav*nz - svel
161  du = (ur-ul)*nx + (vr-vl)*ny + (wr-wl)*nz
162 
163 ! ----- eigenvalues
164 
165  h1 = abs(uvw - cav)
166  h2 = abs(uvw)
167  eabs5 = abs(uvw + cav)
168  delta = epsentr*eabs5
169  eabs1 = entropy_corr1( h1,delta )
170  eabs2 = entropy_corr1( h2,delta )
171 
172 ! ----- upwind fluxes
173 
174  h1 = rav*cav*du
175  h2 = eabs1*(pr-pl - h1)/(2._rfreal*c2a)
176  h3 = eabs2*(rr-rl - (pr-pl)/c2a)
177  h4 = eabs2*rav
178  h5 = eabs5*(pr-pl + h1)/(2._rfreal*c2a)
179 
180  fd(1) = h2 + h3 + h5
181  fd(2) = h2*(uav-cav*nx) + h3*uav + h4*(ur-ul-du*nx) + &
182  h5*(uav+cav*nx)
183  fd(3) = h2*(vav-cav*ny) + h3*vav + h4*(vr-vl-du*ny) + &
184  h5*(vav+cav*ny)
185  fd(4) = h2*(wav-cav*nz) + h3*wav + h4*(wr-wl-du*nz) + &
186  h5*(wav+cav*nz)
187  fd(5) = h2*(hav-cav*uvw) + h3*q2a + h4*(uav*(ur-ul)+ &
188  vav*(vr-vl)+wav*(wr-wl)-uvw*du) + &
189  h5*(hav+cav*uvw)
190 
191  fd(1) = fd(1)*ds
192  fd(2) = fd(2)*ds
193  fd(3) = fd(3)*ds
194  fd(4) = fd(4)*ds
195  fd(5) = fd(5)*ds
196 
197  diss(cv_mixt_dens,ijkc0 ) = diss(cv_mixt_dens,ijkc0 ) - fd(1)
198  diss(cv_mixt_xmom,ijkc0 ) = diss(cv_mixt_xmom,ijkc0 ) - fd(2)
199  diss(cv_mixt_ymom,ijkc0 ) = diss(cv_mixt_ymom,ijkc0 ) - fd(3)
200  diss(cv_mixt_zmom,ijkc0 ) = diss(cv_mixt_zmom,ijkc0 ) - fd(4)
201  diss(cv_mixt_ener,ijkc0 ) = diss(cv_mixt_ener,ijkc0 ) - fd(5)
202 
203  diss(cv_mixt_dens,ijkcm1) = diss(cv_mixt_dens,ijkcm1) + fd(1)
204  diss(cv_mixt_xmom,ijkcm1) = diss(cv_mixt_xmom,ijkcm1) + fd(2)
205  diss(cv_mixt_ymom,ijkcm1) = diss(cv_mixt_ymom,ijkcm1) + fd(3)
206  diss(cv_mixt_zmom,ijkcm1) = diss(cv_mixt_zmom,ijkcm1) + fd(4)
207  diss(cv_mixt_ener,ijkcm1) = diss(cv_mixt_ener,ijkcm1) + fd(5)
208  ENDDO ! i
209  ENDDO ! j
210  ENDDO ! k
211 
212 ! dissipation in j-direction --------------------------------------------------
213 
214  DO k=kpcbeg,kpcend
215  DO i=ipcbeg,ipcend
216  DO j=jpcbeg,jpcend+1
217  ijkc0 = indijk(i,j ,k,icoff,ijcoff)
218  ijkcm1 = indijk(i,j-1,k,icoff,ijcoff)
219  ijkn = indijk(i,j ,k,inoff,ijnoff)
220  ds = sqrt(sj(xcoord,ijkn)*sj(xcoord,ijkn)+ &
221  sj(ycoord,ijkn)*sj(ycoord,ijkn)+ &
222  sj(zcoord,ijkn)*sj(zcoord,ijkn))
223  nx = sj(xcoord,ijkn)/ds
224  ny = sj(ycoord,ijkn)/ds
225  nz = sj(zcoord,ijkn)/ds
226  svel = sjvel(ijkn*indsvel)/ds
227  ds = ds*beta5
228 
229 ! ----- left and right states
230 
231  rl = cv(cv_mixt_dens,ijkcm1)
232  ul = dv(dv_mixt_uvel,ijkcm1)
233  vl = dv(dv_mixt_vvel,ijkcm1)
234  wl = dv(dv_mixt_wvel,ijkcm1)
235  pl = dv(dv_mixt_pres,ijkcm1)
236  rgas = 8314.3_rfreal/gv(gv_mixt_mol,ijkcm1*indmol)
237  gaml = gv(gv_mixt_cp,ijkcm1*indcp)/(gv(gv_mixt_cp,ijkcm1*indcp)-rgas)
238  ggm1 = gaml/(gaml-1._rfreal)
239  hl = ggm1*pl/rl + 0.5_rfreal*(ul*ul+vl*vl+wl*wl)
240 
241  rr = cv(cv_mixt_dens,ijkc0)
242  ur = dv(dv_mixt_uvel,ijkc0)
243  vr = dv(dv_mixt_vvel,ijkc0)
244  wr = dv(dv_mixt_wvel,ijkc0)
245  pr = dv(dv_mixt_pres,ijkc0)
246  rgas = 8314.3_rfreal/gv(gv_mixt_mol,ijkc0*indmol)
247  gamr = gv(gv_mixt_cp,ijkc0*indcp)/(gv(gv_mixt_cp,ijkc0*indcp)-rgas)
248  ggm1 = gamr/(gamr-1._rfreal)
249  hr = ggm1*pr/rr + 0.5_rfreal*(ur*ur+vr*vr+wr*wr)
250 
251 ! ----- Roe`s average
252 
253  rav = sqrt(rl*rr)
254  dd = rav/rl
255  dd1 = 1._rfreal/(1._rfreal+dd)
256  uav = (ul+dd*ur)*dd1
257  vav = (vl+dd*vr)*dd1
258  wav = (wl+dd*wr)*dd1
259  hav = (hl+dd*hr)*dd1
260  q2a = 0.5_rfreal*(uav*uav+vav*vav+wav*wav)
261  gam1 = 0.5_rfreal*(gaml+gamr) - 1._rfreal
262  c2a = gam1*(hav-q2a)
263  cav = sqrt(c2a)
264  uvw = uav*nx + vav*ny + wav*nz - svel
265  du = (ur-ul)*nx + (vr-vl)*ny + (wr-wl)*nz
266 
267 ! ----- eigenvalues
268 
269  h1 = abs(uvw - cav)
270  h2 = abs(uvw)
271  eabs5 = abs(uvw + cav)
272  delta = epsentr*eabs5
273  eabs1 = entropy_corr1( h1,delta )
274  eabs2 = entropy_corr1( h2,delta )
275 
276 ! ----- upwind fluxes
277 
278  h1 = rav*cav*du
279  h2 = eabs1*(pr-pl - h1)/(2._rfreal*c2a)
280  h3 = eabs2*(rr-rl - (pr-pl)/c2a)
281  h4 = eabs2*rav
282  h5 = eabs5*(pr-pl + h1)/(2._rfreal*c2a)
283 
284  fd(1) = h2 + h3 + h5
285  fd(2) = h2*(uav-cav*nx) + h3*uav + h4*(ur-ul-du*nx) + &
286  h5*(uav+cav*nx)
287  fd(3) = h2*(vav-cav*ny) + h3*vav + h4*(vr-vl-du*ny) + &
288  h5*(vav+cav*ny)
289  fd(4) = h2*(wav-cav*nz) + h3*wav + h4*(wr-wl-du*nz) + &
290  h5*(wav+cav*nz)
291  fd(5) = h2*(hav-cav*uvw) + h3*q2a + h4*(uav*(ur-ul)+ &
292  vav*(vr-vl)+wav*(wr-wl)-uvw*du) + &
293  h5*(hav+cav*uvw)
294 
295  fd(1) = fd(1)*ds
296  fd(2) = fd(2)*ds
297  fd(3) = fd(3)*ds
298  fd(4) = fd(4)*ds
299  fd(5) = fd(5)*ds
300 
301  diss(cv_mixt_dens,ijkc0 ) = diss(cv_mixt_dens,ijkc0 ) - fd(1)
302  diss(cv_mixt_xmom,ijkc0 ) = diss(cv_mixt_xmom,ijkc0 ) - fd(2)
303  diss(cv_mixt_ymom,ijkc0 ) = diss(cv_mixt_ymom,ijkc0 ) - fd(3)
304  diss(cv_mixt_zmom,ijkc0 ) = diss(cv_mixt_zmom,ijkc0 ) - fd(4)
305  diss(cv_mixt_ener,ijkc0 ) = diss(cv_mixt_ener,ijkc0 ) - fd(5)
306 
307  diss(cv_mixt_dens,ijkcm1) = diss(cv_mixt_dens,ijkcm1) + fd(1)
308  diss(cv_mixt_xmom,ijkcm1) = diss(cv_mixt_xmom,ijkcm1) + fd(2)
309  diss(cv_mixt_ymom,ijkcm1) = diss(cv_mixt_ymom,ijkcm1) + fd(3)
310  diss(cv_mixt_zmom,ijkcm1) = diss(cv_mixt_zmom,ijkcm1) + fd(4)
311  diss(cv_mixt_ener,ijkcm1) = diss(cv_mixt_ener,ijkcm1) + fd(5)
312  ENDDO ! j
313  ENDDO ! i
314  ENDDO ! k
315 
316 ! dissipation in k-direction --------------------------------------------------
317 
318  DO j=jpcbeg,jpcend
319  DO i=ipcbeg,ipcend
320  DO k=kpcbeg,kpcend+1
321  ijkc0 = indijk(i,j,k ,icoff,ijcoff)
322  ijkcm1 = indijk(i,j,k-1,icoff,ijcoff)
323  ijkn = indijk(i,j,k ,inoff,ijnoff)
324  ds = sqrt(sk(xcoord,ijkn)*sk(xcoord,ijkn)+ &
325  sk(ycoord,ijkn)*sk(ycoord,ijkn)+ &
326  sk(zcoord,ijkn)*sk(zcoord,ijkn))
327  nx = sk(xcoord,ijkn)/ds
328  ny = sk(ycoord,ijkn)/ds
329  nz = sk(zcoord,ijkn)/ds
330  svel = skvel(ijkn*indsvel)/ds
331  ds = ds*beta5
332 
333 ! ----- left and right states
334 
335  rl = cv(cv_mixt_dens,ijkcm1)
336  ul = dv(dv_mixt_uvel,ijkcm1)
337  vl = dv(dv_mixt_vvel,ijkcm1)
338  wl = dv(dv_mixt_wvel,ijkcm1)
339  pl = dv(dv_mixt_pres,ijkcm1)
340  rgas = 8314.3_rfreal/gv(gv_mixt_mol,ijkcm1*indmol)
341  gaml = gv(gv_mixt_cp,ijkcm1*indcp)/(gv(gv_mixt_cp,ijkcm1*indcp)-rgas)
342  ggm1 = gaml/(gaml-1._rfreal)
343  hl = ggm1*pl/rl + 0.5_rfreal*(ul*ul+vl*vl+wl*wl)
344 
345  rr = cv(cv_mixt_dens,ijkc0)
346  ur = dv(dv_mixt_uvel,ijkc0)
347  vr = dv(dv_mixt_vvel,ijkc0)
348  wr = dv(dv_mixt_wvel,ijkc0)
349  pr = dv(dv_mixt_pres,ijkc0)
350  rgas = 8314.3_rfreal/gv(gv_mixt_mol,ijkc0*indmol)
351  gamr = gv(gv_mixt_cp,ijkc0*indcp)/(gv(gv_mixt_cp,ijkc0*indcp)-rgas)
352  ggm1 = gamr/(gamr-1._rfreal)
353  hr = ggm1*pr/rr + 0.5_rfreal*(ur*ur+vr*vr+wr*wr)
354 
355 ! ----- Roe`s average
356 
357  rav = sqrt(rl*rr)
358  dd = rav/rl
359  dd1 = 1._rfreal/(1._rfreal+dd)
360  uav = (ul+dd*ur)*dd1
361  vav = (vl+dd*vr)*dd1
362  wav = (wl+dd*wr)*dd1
363  hav = (hl+dd*hr)*dd1
364  q2a = 0.5_rfreal*(uav*uav+vav*vav+wav*wav)
365  gam1 = 0.5_rfreal*(gaml+gamr) - 1._rfreal
366  c2a = gam1*(hav-q2a)
367  cav = sqrt(c2a)
368  uvw = uav*nx + vav*ny + wav*nz - svel
369  du = (ur-ul)*nx + (vr-vl)*ny + (wr-wl)*nz
370 
371 ! ----- eigenvalues
372 
373  h1 = abs(uvw - cav)
374  h2 = abs(uvw)
375  eabs5 = abs(uvw + cav)
376  delta = epsentr*eabs5
377  eabs1 = entropy_corr1( h1,delta )
378  eabs2 = entropy_corr1( h2,delta )
379 
380 ! ----- upwind fluxes
381 
382  h1 = rav*cav*du
383  h2 = eabs1*(pr-pl - h1)/(2._rfreal*c2a)
384  h3 = eabs2*(rr-rl - (pr-pl)/c2a)
385  h4 = eabs2*rav
386  h5 = eabs5*(pr-pl + h1)/(2._rfreal*c2a)
387 
388  fd(1) = h2 + h3 + h5
389  fd(2) = h2*(uav-cav*nx) + h3*uav + h4*(ur-ul-du*nx) + &
390  h5*(uav+cav*nx)
391  fd(3) = h2*(vav-cav*ny) + h3*vav + h4*(vr-vl-du*ny) + &
392  h5*(vav+cav*ny)
393  fd(4) = h2*(wav-cav*nz) + h3*wav + h4*(wr-wl-du*nz) + &
394  h5*(wav+cav*nz)
395  fd(5) = h2*(hav-cav*uvw) + h3*q2a + h4*(uav*(ur-ul)+ &
396  vav*(vr-vl)+wav*(wr-wl)-uvw*du) + &
397  h5*(hav+cav*uvw)
398 
399  fd(1) = fd(1)*ds
400  fd(2) = fd(2)*ds
401  fd(3) = fd(3)*ds
402  fd(4) = fd(4)*ds
403  fd(5) = fd(5)*ds
404 
405  diss(cv_mixt_dens,ijkc0 ) = diss(cv_mixt_dens,ijkc0 ) - fd(1)
406  diss(cv_mixt_xmom,ijkc0 ) = diss(cv_mixt_xmom,ijkc0 ) - fd(2)
407  diss(cv_mixt_ymom,ijkc0 ) = diss(cv_mixt_ymom,ijkc0 ) - fd(3)
408  diss(cv_mixt_zmom,ijkc0 ) = diss(cv_mixt_zmom,ijkc0 ) - fd(4)
409  diss(cv_mixt_ener,ijkc0 ) = diss(cv_mixt_ener,ijkc0 ) - fd(5)
410 
411  diss(cv_mixt_dens,ijkcm1) = diss(cv_mixt_dens,ijkcm1) + fd(1)
412  diss(cv_mixt_xmom,ijkcm1) = diss(cv_mixt_xmom,ijkcm1) + fd(2)
413  diss(cv_mixt_ymom,ijkcm1) = diss(cv_mixt_ymom,ijkcm1) + fd(3)
414  diss(cv_mixt_zmom,ijkcm1) = diss(cv_mixt_zmom,ijkcm1) + fd(4)
415  diss(cv_mixt_ener,ijkcm1) = diss(cv_mixt_ener,ijkcm1) + fd(5)
416  ENDDO ! k
417  ENDDO ! i
418  ENDDO ! j
419 
420 ! finalize --------------------------------------------------------------------
421 
422  CALL deregisterfunction( region%global )
423 
424 ! #############################################################################
425 
426  CONTAINS
427 
428  REAL(RFREAL) FUNCTION entropy_corr1( z,d )
429 
430  REAL(RFREAL) :: z, d
431 
432  IF (z > d) THEN
433  entropy_corr1 = z
434  ELSE
435  entropy_corr1 = 0.5_rfreal*(z*z+d*d)/d
436  ENDIF
437 
438  END FUNCTION entropy_corr1
439 
440 END SUBROUTINE rflo_roedissipfirst
441 
442 !******************************************************************************
443 !
444 ! RCS Revision history:
445 !
446 ! $Log: RFLO_RoeDissipFirst.F90,v $
447 ! Revision 1.3 2008/12/06 08:44:27 mtcampbe
448 ! Updated license.
449 !
450 ! Revision 1.2 2008/11/19 22:17:38 mtcampbe
451 ! Added Illinois Open Source License/Copyright
452 !
453 ! Revision 1.1 2004/11/29 20:51:40 wasistho
454 ! lower to upper case
455 !
456 ! Revision 1.6 2003/11/20 16:40:40 mdbrandy
457 ! Backing out RocfluidMP changes from 11-17-03
458 !
459 ! Revision 1.2 2003/10/01 23:52:10 jblazek
460 ! Corrected bug in moving noslip wall BC and grid speeds.
461 !
462 ! Revision 1.1 2003/05/29 17:28:43 jblazek
463 ! Implemented Roe scheme.
464 !
465 !******************************************************************************
466 
467 
468 
469 
470 
471 
472 
d
const NT & d
Definition: rational_rotation.h:73
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
sqrt
double sqrt(double d)
Definition: double.h:73
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
rflo_getnodeoffset
subroutine rflo_getnodeoffset(region, iLev, iNodeOffset, ijNodeOffset)
Definition: RFLO_GetDimensions.F90:267
entropy_corr1
REAL(RFREAL) function entropy_corr1(z, d)
Definition: RFLO_RoeDissipFirst.F90:428
z
void int int int REAL REAL REAL * z
Definition: write.cpp:76
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
rflo_roedissipfirst
subroutine rflo_roedissipfirst(region)
Definition: RFLO_RoeDissipFirst.F90:44
i
blockLoc i
Definition: read.cpp:79
rflo_getcelloffset
subroutine rflo_getcelloffset(region, iLev, iCellOffset, ijCellOffset)
Definition: RFLO_GetDimensions.F90:235
Direction_2::delta
RT delta(int i) const
Definition: Direction_2.h:159
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