Skip to content

Commit

Permalink
Revert "Use precomputed logte"
Browse files Browse the repository at this point in the history 
This reverts commit 8f779d1.
  • Loading branch information
@holm10
holm10 committed Jan 16, 2025
1 parent 4dad69a commit 882443a
Showing 1 changed file with 12 additions and 12 deletions.
24 changes: 12 additions & 12 deletions bbb/oderhs.m
View file
Original file line number Diff line number Diff line change
Expand Up @@ -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.
Expand Down

0 comments on commit 882443a

Please sign in to comment.