From b6b3175a5a49e4477096e44ebbe7ec6a22795c89 Mon Sep 17 00:00:00 2001
From: Andreas Holm <60451789+holm10@users.noreply.github.com>
Date: Tue, 14 Jan 2025 23:16:16 -0800
Subject: [PATCH] Revert "Uses precomputed logpg"

This reverts commit fd7fef7477143279b1997994c6d1aba689ca5535.
---
 bbb/oderhs.m | 24 ++++++++++++------------
 1 file changed, 12 insertions(+), 12 deletions(-)

diff --git a/bbb/oderhs.m b/bbb/oderhs.m
index 4703b7b8..392323fa 100755
--- a/bbb/oderhs.m
+++ b/bbb/oderhs.m
@@ -6371,16 +6371,16 @@ call fd2tra (nx,ny,floxg,floyg,conxg,conyg,
      .               (ix==ixrb(jx).and.ixmxbcl==1) ) isxyfl = .false.
                enddo
                if (methgx .eq. 6) then  # log interpolation
-               grdnv =(   ( fym (ix,iy,1)*logpg(ix2,iy1 ,igsp) + 
-     .                      fy0 (ix,iy,1)*logpg(ix2,iy  ,igsp) +
-     .                      fyp (ix,iy,1)*logpg(ix2,iy+1,igsp) + 
-     .                      fymx(ix,iy,1)*logpg(ix ,iy1 ,igsp) +
-     .                      fypx(ix,iy,1)*logpg(ix, iy+1,igsp) )
-     .                   -( fym (ix,iy,0)*logpg(ix ,iy1 ,igsp) +
-     .                      fy0 (ix,iy,0)*logpg(ix ,iy  ,igsp) +
-     .                      fyp (ix,iy,0)*logpg(ix ,iy+1,igsp) +
-     .                      fymx(ix,iy,0)*logpg(ix4,iy1 ,igsp) + 
-     .                      fypx(ix,iy,0)*logpg(ix6,iy+1,igsp) ) )/ 
+               grdnv =(   ( fym (ix,iy,1)*log(pg(ix2,iy1 ,igsp)) + 
+     .                      fy0 (ix,iy,1)*log(pg(ix2,iy  ,igsp)) +
+     .                      fyp (ix,iy,1)*log(pg(ix2,iy+1,igsp)) + 
+     .                      fymx(ix,iy,1)*log(pg(ix ,iy1 ,igsp)) +
+     .                      fypx(ix,iy,1)*log(pg(ix, iy+1,igsp)) )
+     .                   -( fym (ix,iy,0)*log(pg(ix ,iy1 ,igsp)) +
+     .                      fy0 (ix,iy,0)*log(pg(ix ,iy  ,igsp)) +
+     .                      fyp (ix,iy,0)*log(pg(ix ,iy+1,igsp)) +
+     .                      fymx(ix,iy,0)*log(pg(ix4,iy1 ,igsp)) + 
+     .                      fypx(ix,iy,0)*log(pg(ix6,iy+1,igsp)) ) )/ 
      .                                                  dxnog(ix,iy)
                elseif (methgx .eq. 7) then # inverse interpolation
                grdnv =( 1/(fym (ix,iy,1)/pg(ix2,iy1 ,igsp) + 
@@ -6413,9 +6413,9 @@ call fd2tra (nx,ny,floxg,floyg,conxg,conyg,
      .                           0.5*(nuiz(ix,iy,igsp)+nuiz(ix2,iy,igsp))
               if (methgx .eq. 6) then
                fngxy(ix,iy,igsp) =  exp( 0.5*
-     .                     (logpg(ix2,iy,igsp)+logpg(ix,iy,igsp)) )*
+     .                     (log(pg(ix2,iy,igsp))+log(pg(ix,iy,igsp))) )*
      .                               difgx2*(grdnv/cosangfx(ix,iy) -
-     .                     (logpg(ix2,iy,igsp) - logpg(ix,iy,igsp))*
+     .                     (log(pg(ix2,iy,igsp)) - log(pg(ix,iy,igsp)))*
      .                                 gxf(ix,iy) ) * sx(ix,iy)
               else
                fngxy(ix,iy,igsp) = difgx2*(grdnv/cosangfx(ix,iy) -