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 080b16c commit e23c832
Showing 1 changed file with 9 additions and 0 deletions.
9 changes: 9 additions & 0 deletions hypnotoad/core/mesh.py
Original file line number Diff line number Diff line change
Expand Up @@ -3861,6 +3861,15 @@ 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 e23c832

Please sign in to comment.