@@ -1251,9 +1251,9 @@ contains
12511251 if (Re_size(i) > 0 ) Re_R(i) = 0._wp
12521252 $:GPU_LOOP(parallelism= ' [seq]' )
12531253 do q = 1 , Re_size(i)
1254- Re_L(i) = qL_prim_rs${XYZ}$_vf(j, k, l, E_idx + Re_idx(i, q))/ Res (i, q) &
1254+ Re_L(i) = qL_prim_rs${XYZ}$_vf(j, k, l, E_idx + Re_idx(i, q))/ Res_gs (i, q) &
12551255 + Re_L(i)
1256- Re_R(i) = qR_prim_rs${XYZ}$_vf(j + 1 , k, l, E_idx + Re_idx(i, q))/ Res (i, q) &
1256+ Re_R(i) = qR_prim_rs${XYZ}$_vf(j + 1 , k, l, E_idx + Re_idx(i, q))/ Res_gs (i, q) &
12571257 + Re_R(i)
12581258 end do
12591259 Re_L(i) = 1._wp / max (Re_L(i), sgm_eps)
@@ -1923,9 +1923,9 @@ contains
19231923
19241924 $:GPU_LOOP(parallelism= ' [seq]' )
19251925 do q = 1 , Re_size(i)
1926- Re_L(i) = (1._wp - qL_prim_rs${XYZ}$_vf(j, k, l, E_idx + Re_idx(i, q)))/ Res (i, q) &
1926+ Re_L(i) = (1._wp - qL_prim_rs${XYZ}$_vf(j, k, l, E_idx + Re_idx(i, q)))/ Res_gs (i, q) &
19271927 + Re_L(i)
1928- Re_R(i) = (1._wp - qR_prim_rs${XYZ}$_vf(j + 1 , k, l, E_idx + Re_idx(i, q)))/ Res (i, q) &
1928+ Re_R(i) = (1._wp - qR_prim_rs${XYZ}$_vf(j + 1 , k, l, E_idx + Re_idx(i, q)))/ Res_gs (i, q) &
19291929 + Re_R(i)
19301930 end do
19311931
@@ -2363,9 +2363,9 @@ contains
23632363
23642364 $:GPU_LOOP(parallelism= ' [seq]' )
23652365 do q = 1 , Re_size(i)
2366- Re_L(i) = qL_prim_rs${XYZ}$_vf(j, k, l, E_idx + Re_idx(i, q))/ Res (i, q) &
2366+ Re_L(i) = qL_prim_rs${XYZ}$_vf(j, k, l, E_idx + Re_idx(i, q))/ Res_gs (i, q) &
23672367 + Re_L(i)
2368- Re_R(i) = qR_prim_rs${XYZ}$_vf(j + 1 , k, l, E_idx + Re_idx(i, q))/ Res (i, q) &
2368+ Re_R(i) = qR_prim_rs${XYZ}$_vf(j + 1 , k, l, E_idx + Re_idx(i, q))/ Res_gs (i, q) &
23692369 + Re_R(i)
23702370 end do
23712371 Re_L(i) = 1._wp / max (Re_L(i), sgm_eps)
0 commit comments