gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

REFINEMENT: Take absolute values in the error calculation.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2013-12-19 14:47:12 -02:00
parent 072f234e0f
commit 85b92c8090
1 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
src/refinement.F90
View File

@ -283,16 +283,16 @@ module refinement
! !
fr = pdata%q(iqt,ip1,j,k) - pdata%q(iqt,i ,j,k) fr = pdata%q(iqt,ip1,j,k) - pdata%q(iqt,i ,j,k)
fl = pdata%q(iqt,im1,j,k) - pdata%q(iqt,i ,j,k) fl = pdata%q(iqt,im1,j,k) - pdata%q(iqt,i ,j,k)
fc = pdata%q(iqt,ip1,j,k) + pdata%q(iqt,im1,j,k) & fc = abs(pdata%q(iqt,ip1,j,k)) + abs(pdata%q(iqt,im1,j,k)) &
+ 2.0d+00 * pdata%q(iqt,i,j,k) + 2.0d+00 * abs(pdata%q(iqt,i,j,k))
fx = (fr + fl) / (abs(fr) + abs(fl) + epsref * fc) fx = (fr + fl) / (abs(fr) + abs(fl) + epsref * fc)
! calculate the second derivative error along the Y direction ! calculate the second derivative error along the Y direction
! !
fr = pdata%q(iqt,i,jp1,k) - pdata%q(iqt,i,j ,k) fr = pdata%q(iqt,i,jp1,k) - pdata%q(iqt,i,j ,k)
fl = pdata%q(iqt,i,jm1,k) - pdata%q(iqt,i,j ,k) fl = pdata%q(iqt,i,jm1,k) - pdata%q(iqt,i,j ,k)
fc = pdata%q(iqt,i,jp1,k) + pdata%q(iqt,i,jm1,k) & fc = abs(pdata%q(iqt,i,jp1,k)) + abs(pdata%q(iqt,i,jm1,k)) &
+ 2.0d+00 * pdata%q(iqt,i,j,k) + 2.0d+00 * abs(pdata%q(iqt,i,j,k))
fy = (fr + fl) / (abs(fr) + abs(fl) + epsref * fc) fy = (fr + fl) / (abs(fr) + abs(fl) + epsref * fc)
#if NDIMS == 3 #if NDIMS == 3
@ -300,8 +300,8 @@ module refinement
! !
fr = pdata%q(iqt,i,j,kp1) - pdata%q(iqt,i,j,k ) fr = pdata%q(iqt,i,j,kp1) - pdata%q(iqt,i,j,k )
fl = pdata%q(iqt,i,j,km1) - pdata%q(iqt,i,j,k ) fl = pdata%q(iqt,i,j,km1) - pdata%q(iqt,i,j,k )
fc = pdata%q(iqt,i,j,kp1) + pdata%q(iqt,i,j,km1) & fc = abs(pdata%q(iqt,i,j,kp1)) + abs(pdata%q(iqt,i,j,km1)) &
+ 2.0d+00 * pdata%q(iqt,i,j,k) + 2.0d+00 * abs(pdata%q(iqt,i,j,k))
fz = (fr + fl) / (abs(fr) + abs(fl) + epsref * fc) fz = (fr + fl) / (abs(fr) + abs(fl) + epsref * fc)
#endif /* NDIMS == 3 */ #endif /* NDIMS == 3 */