From 705857e157bf720036c705420e1bc7088266f926 Mon Sep 17 00:00:00 2001
From: jsbamboo <zhang73@llnl.gov>
Date: Mon, 16 Sep 2024 11:18:44 -0700
Subject: [PATCH] remove the exact integration for condition (abs(yseg(1)
 -yseg(2))<fuzzy_width)

---
 cube_to_target/remap.F90 | 6 ++++--
 1 file changed, 4 insertions(+), 2 deletions(-)

diff --git a/cube_to_target/remap.F90 b/cube_to_target/remap.F90
index 65b001c..9b6eed3 100644
--- a/cube_to_target/remap.F90
+++ b/cube_to_target/remap.F90
@@ -1093,7 +1093,8 @@ subroutine side_integral(&
 
             if (xseg(1).EQ.xseg(2))then
               slope = bignum
-            else if (abs(yseg(1) -yseg(2))<fuzzy_width) then
+           !else if (abs(yseg(1) -yseg(2))<fuzzy_width) then !remove the exact integration
+            else if (.false.) then
               slope = 0.0
             else
               slope    = (yseg(2)-yseg(1))/(xseg(2)-xseg(1))
@@ -1196,7 +1197,8 @@ subroutine get_weights_gauss(weights,xseg,yseg,nreconstruction,ngauss,gauss_weig
 !    if (fuzzy(abs(xseg(1) -xseg(2)),fuzzy_width)==0)then
     if (xseg(1).EQ.xseg(2))then
       weights = 0.0D0
-    else if (abs(yseg(1) -yseg(2))<fuzzy_width) then
+    !else if (abs(yseg(1) -yseg(2))<fuzzy_width) then !remove the exact integration
+    else if (.false.) then
       !
       ! line segment parallel to latitude - compute weights exactly
       !