From a351b17b852f927b67097d117e0ed6cecb5fb67f Mon Sep 17 00:00:00 2001
From: BS <brendan.shanahan@ipp.mpg.de>
Date: Wed, 12 May 2021 10:21:36 +0200
Subject: [PATCH] Use Bxy.ydown() for div_a_laplace_perp

---
 src/mesh/fv_ops.cxx | 5 +++++
 1 file changed, 5 insertions(+)

diff --git a/src/mesh/fv_ops.cxx b/src/mesh/fv_ops.cxx
index 118b8139db..1640c1bf0b 100644
--- a/src/mesh/fv_ops.cxx
+++ b/src/mesh/fv_ops.cxx
@@ -166,6 +166,11 @@ namespace FV {
           yzresult(i, j, k) = fout / (dyc(i, j, k) * Jc(i, j, k));
 
           // Calculate flux between j and j-1
+	  coef =
+              0.5
+              * (g_23c(i, j, k) / SQ(Jc(i, j, k) * Bxyc(i, j, k))
+                 + g_23down(i, j - 1, k) / SQ(Jdown(i, j - 1, k) * Bxydown(i, j - 1, k)));
+
           dfdz =
               0.5
               * (fc(i, j, kp) - fc(i, j, km) + fdown(i, j - 1, kp) - fdown(i, j - 1, km))