Skip to content

Commit 80ddc85

Browse files
ChrisZYJChrisZYJ
committed
remove beta scaling
1 parent 17d6308 commit 80ddc85

File tree

1 file changed

+5
-4
lines changed

1 file changed

+5
-4
lines changed

‎src/simulation/m_weno.fpp‎

Lines changed: 5 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -841,6 +841,7 @@ contains
841841

842842
if (.not. teno) then
843843
! (Balsara & Shu, 2000) Page 11 Section III.a
844+
! Note: paratheses are needed to group the terms before '+ weno_eps' to avoid unintended floating point errors
844845
beta(0) = ( v_rs_ws_${XYZ}$ (j - 3, k, l, i)*( 547d0*v_rs_ws_${XYZ}$ (j - 3, k, l, i) & !&
845846
- 3882d0*v_rs_ws_${XYZ}$ (j - 2, k, l, i) & !&
846847
+ 4642d0*v_rs_ws_${XYZ}$ (j - 1, k, l, i) & !&
@@ -850,7 +851,7 @@ contains
850851
+ 7042d0*v_rs_ws_${XYZ}$ (j , k, l, i) ) & !&
851852
+ v_rs_ws_${XYZ}$ (j - 1, k, l, i)*( 11003d0*v_rs_ws_${XYZ}$ (j - 1, k, l, i) & !&
852853
- 9402d0*v_rs_ws_${XYZ}$ (j , k, l, i) ) & !&
853-
+ 2107d0*v_rs_ws_${XYZ}$ (j, k, l, i)**2d0 ) / 240d0 & !&
854+
+ 2107d0*v_rs_ws_${XYZ}$ (j, k, l, i)**2d0 ) & !&
854855
+ weno_eps !&
855856

856857
beta(1) = ( v_rs_ws_${XYZ}$ (j - 2, k, l, i)*( 267d0*v_rs_ws_${XYZ}$ (j - 2, k, l, i) & !&
@@ -862,7 +863,7 @@ contains
862863
+ 1922d0*v_rs_ws_${XYZ}$ (j + 1, k, l, i) ) & !&
863864
+ v_rs_ws_${XYZ}$ (j , k, l, i)*( 3443d0*v_rs_ws_${XYZ}$ (j , k, l, i) & !&
864865
- 2522d0*v_rs_ws_${XYZ}$ (j + 1, k, l, i) ) & !&
865-
+ 547d0*v_rs_ws_${XYZ}$ (j + 1, k, l, i)**2d0 ) / 240d0 & !&
866+
+ 547d0*v_rs_ws_${XYZ}$ (j + 1, k, l, i)**2d0 ) & !&
866867
+ weno_eps !&
867868

868869
beta(2) = ( v_rs_ws_${XYZ}$ (j - 1, k, l, i)*( 547d0*v_rs_ws_${XYZ}$ (j - 1, k, l, i) & !&
@@ -874,7 +875,7 @@ contains
874875
+ 1602d0*v_rs_ws_${XYZ}$ (j + 2, k, l, i) ) & !&
875876
+ v_rs_ws_${XYZ}$ (j + 1, k, l, i)*( 2843d0*v_rs_ws_${XYZ}$ (j + 1, k, l, i) & !&
876877
- 1642d0*v_rs_ws_${XYZ}$ (j + 2, k, l, i) ) & !&
877-
+ 267d0*v_rs_ws_${XYZ}$ (j + 2, k, l, i)**2d0 ) / 240d0 & !&
878+
+ 267d0*v_rs_ws_${XYZ}$ (j + 2, k, l, i)**2d0 ) & !&
878879
+ weno_eps !&
879880

880881
beta(3) = ( v_rs_ws_${XYZ}$ (j , k, l, i)*( 2107d0*v_rs_ws_${XYZ}$ (j , k, l, i) & !&
@@ -886,7 +887,7 @@ contains
886887
+ 4642d0*v_rs_ws_${XYZ}$ (j + 3, k, l, i) ) & !&
887888
+ v_rs_ws_${XYZ}$ (j + 2, k, l, i)*( 7043d0*v_rs_ws_${XYZ}$ (j + 2, k, l, i) & !&
888889
- 3882d0*v_rs_ws_${XYZ}$ (j + 3, k, l, i) ) & !&
889-
+ 547d0*v_rs_ws_${XYZ}$ (j + 3, k, l, i)**2d0 ) / 240d0 & !&
890+
+ 547d0*v_rs_ws_${XYZ}$ (j + 3, k, l, i)**2d0 ) & !&
890891
+ weno_eps !&
891892

892893
else ! TENO

0 commit comments

Comments
 (0)