Skip to content

Commit a36cc6c

Browse files
committed
Fix in 2D set_eb_data that gave incorrect vcent and vfrac
1 parent 772868a commit a36cc6c

1 file changed

Lines changed: 2 additions & 2 deletions

File tree

Src/EB/AMReX_EB2_2D_C.cpp

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -101,14 +101,14 @@ void set_eb_data (const int i, const int j,
101101
}
102102
vcent(i,j,0,1) = 0.0_rt;
103103
} else {
104-
Real aa = nxabs/ny*dx[0]/dx[1];
104+
Real aa = nxabs/ny;
105105
const Real dxx = x_ym - x_yp;
106106
const Real dx2 = dxx * (x_ym + x_yp);
107107
const Real dx3 = dxx * (x_ym*x_ym + x_ym*x_yp + x_yp*x_yp);
108108
const Real af1 = 0.5_rt*(axm+axp)*dx[0] + aa*0.5_rt*dx2;
109109
vcent(i,j,0,0) = -0.125_rt*daxp*dx[0]*dx[0] + aa*(1._rt/6._rt)*dx3;
110110

111-
aa = nyabs/nx*dx[1]/dx[0];
111+
aa = nyabs/nx;
112112
const Real dy = y_xm - y_xp;
113113
const Real dy2 = dy * (y_xm + y_xp);
114114
const Real dy3 = dy * (y_xm*y_xm + y_xm*y_xp + y_xp*y_xp);

0 commit comments

Comments
 (0)