Skip to content

Commit

Permalink
use MXG from grid file
Browse files Browse the repository at this point in the history
  • Loading branch information
dschwoerer committed Jan 7, 2025
1 parent c29d2bb commit 8a81d6b
Showing 1 changed file with 8 additions and 4 deletions.
12 changes: 8 additions & 4 deletions tests/integrated/test-fci-boundary/runtest
Original file line number Diff line number Diff line change
Expand Up @@ -29,15 +29,19 @@ mthread = 2
directory = "data"

with DataFile("grid.fci.nc") as grid:
xfwd = grid.read("forward_xt_prime")[2:-2]
xbwd = grid.read("backward_xt_prime")[2:-2]
try:
MXG = grid.read("MXG")
except:
MXG = 1
xfwd = grid.read("forward_xt_prime")[MXG:-MXG]
xbwd = grid.read("backward_xt_prime")[MXG:-MXG]

nx = xfwd.shape[0]

regions = {
"xin_fwd": xfwd < 1,
"xin_fwd": xfwd < MXG,
"xout_fwd": xfwd > nx,
"xin_bwd": xbwd < 1,
"xin_bwd": xbwd < MXG,
"xout_bwd": xbwd > nx,
}
regions = {k: v.astype(int) for k, v in regions.items()}
Expand Down

0 comments on commit 8a81d6b

Please sign in to comment.