From 882443ae4fae388f85ab941db53965d1657c43a5 Mon Sep 17 00:00:00 2001 From: Andreas Holm <60451789+holm10@users.noreply.github.com> Date: Thu, 16 Jan 2025 13:11:25 -0800 Subject: [PATCH] Revert "Use precomputed logte" This reverts commit 8f779d1c9c4d77b5070b0820aa658a76de8cfbbb. --- bbb/oderhs.m | 24 ++++++++++++------------ 1 file changed, 12 insertions(+), 12 deletions(-) diff --git a/bbb/oderhs.m b/bbb/oderhs.m index d506cd59..9cb3b834 100755 --- a/bbb/oderhs.m +++ b/bbb/oderhs.m @@ -4330,23 +4330,23 @@ c if (ifld .ne. iigsp) then ix4 = ixp1(ix,iy1) ix5 = ixm1(ix,iy+1) ix6 = ixp1(ix,iy+1) - grdnv=( ( fym (ix,iy,1)*logte(ix2,iy1 ) + - . fy0 (ix,iy,1)*logte(ix2,iy ) + - . fyp (ix,iy,1)*logte(ix2,iy+1) + - . fymx(ix,iy,1)*logte(ix ,iy1 ) + - . fypx(ix,iy,1)*logte(ix ,iy+1) ) - . -( fym (ix,iy,0)*logte(ix ,iy1 ) + - . fy0 (ix,iy,0)*logte(ix ,iy ) + - . fyp (ix,iy,0)*logte(ix ,iy+1) + - . fymx(ix,iy,0)*logte(ix4,iy1 ) + - . fypx(ix,iy,0)*logte(ix6,iy+1) ) ) / + grdnv=( ( fym (ix,iy,1)*log(te(ix2,iy1 )) + + . fy0 (ix,iy,1)*log(te(ix2,iy )) + + . fyp (ix,iy,1)*log(te(ix2,iy+1)) + + . fymx(ix,iy,1)*log(te(ix ,iy1 )) + + . fypx(ix,iy,1)*log(te(ix ,iy+1)) ) + . -( fym (ix,iy,0)*log(te(ix ,iy1 )) + + . fy0 (ix,iy,0)*log(te(ix ,iy )) + + . fyp (ix,iy,0)*log(te(ix ,iy+1)) + + . fymx(ix,iy,0)*log(te(ix4,iy1 )) + + . fypx(ix,iy,0)*log(te(ix6,iy+1)) ) ) / . dxnog(ix,iy) feexy(ix,iy) = exp( 0.5* - . (logte(ix2,iy) + logte(ix,iy)) )* + . (log(te(ix2,iy)) + log(te(ix,iy))) )* . (fcdif*kye+kye_use(ix,iy))*0.5* . (ne(ix2,iy)+ne(ix,iy))* . (grdnv/cosangfx(ix,iy) - - . (logte(ix2,iy) - logte(ix,iy))* + . (log(te(ix2,iy)) - log(te(ix,iy)))* . gxf(ix,iy))*sx(ix,iy) c... Now do the Ti equation.