Skip to content

Commit

Permalink
Fix ixseps* indices for divertor-leg-only grid
Browse files Browse the repository at this point in the history
It makes most sense to set ixseps1 and ixseps2 so that the divertor-leg
grid is treated as a SOL-only grid, where neither y-boundary connects to
anything.
  • Loading branch information
johnomotani committed Nov 13, 2024
1 parent 6319293 commit 4124957
Showing 1 changed file with 6 additions and 0 deletions.
6 changes: 6 additions & 0 deletions hypnotoad/core/mesh.py
Original file line number Diff line number Diff line change
Expand Up @@ -3860,6 +3860,12 @@ def writeGridfile(self, filename):
# separatrices in the same radial location
ixseps2 = ixseps1

if "single_region" in self.equilibrium.user_options and "divertor" in self.equilibrium.user_options.single_region:
# Need to hack topology indices to get something sensible for
# single_region case.
ixseps1 = -1
ixseps2 = -1

f.write("ixseps1", ixseps1)
f.write("ixseps2", ixseps2)
f.write("jyseps1_1", jyseps1_1)
Expand Down

0 comments on commit 4124957

Please sign in to comment.