Fix scalar diffusion rhs.

Author: skang67

Source/Diffusion/ERF_DiffusionSrcForMom_EB.cpp

@@ -318,8 +318,6 @@ DiffusionSrcForMom_EB (const MFIter& mfi,
 
                     rho_w_rhs(i,j,k) -= mu_eff * barea * dwdn / (vol * w_volfrac(i,j,k));
 
-                } else if (l_surface_layer) {
-                    // SK TODO: implement surface layer BC contribution to diffusion
                 }
             }
         }

Source/Diffusion/ERF_DiffusionSrcForState_EB.cpp

@@ -243,9 +243,8 @@ DiffusionSrcForState_EB (const Box& bx, const Box& domain,
             ParallelFor(bx,[=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept
             {
                 if (cfg_arr(i,j,k).isSingleValued()) {
-
-                    cell_rhs(i,j,k,qty_index) += barea_arr(i,j,k) / hfx_EB(i,j,k) / (vol * detJ(i,j,k));
-
+                                       
+                    cell_rhs(i,j,k,qty_index) += barea_arr(i,j,k) * hfx_EB(i,j,k) / (vol * detJ(i,j,k));
                 }
             });
         }