From 9d68a3df854ca61510e0a0612d5f76ba749c9fe0 Mon Sep 17 00:00:00 2001 From: Jeremy Lore Date: Mon, 10 Aug 2026 13:30:26 -0400 Subject: [PATCH] Use positive face-normal projection for fcQgam --- src/preprocessing/b2agmt.F | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/preprocessing/b2agmt.F b/src/preprocessing/b2agmt.F index 48ba7fca9..1fd26e626 100644 --- a/src/preprocessing/b2agmt.F +++ b/src/preprocessing/b2agmt.F @@ -529,10 +529,10 @@ subroutine b2agmt (g, m, isymm, artificial_slab) ! With this method, numerical error on the 5-point contribution on ! gradient calculation is significantly reduced for general grids - fcQgam1(iFc,0) = sinang(ec1fc(0),ec1fc(1),fn(0),fn(1)) + fcQgam1(iFc,0) = abs(sinang(ec1fc(0),ec1fc(1),fn(0),fn(1))) fcQgam1(iFc,1) = cosang(ec1fc(0),ec1fc(1),fn(0),fn(1)) - fcQgam2(iFc,0) = sinang(efcc2(0),efcc2(1),fn(0),fn(1)) + fcQgam2(iFc,0) = abs(sinang(efcc2(0),efcc2(1),fn(0),fn(1))) fcQgam2(iFc,1) = cosang(efcc2(0),efcc2(1),fn(0),fn(1)) g%fcQgam(iFc,0) = ( fcQgam1(iFc,0)*g%fcHc(iFc,1) +