Populate grid subsets for faces and key points (like X-points, strike points)

[eldond]

• x_points: 6 (nodes defining X-points)
• core_cut: 7 (y-aligned faces inside the separatrix connecting to the active x-point)
  • This tells us where the cut defined by b2fgmtry was
  • If extrapolated, this curve would pass through the magnetic axis
• PFR_cut: 8 (y-aligned faces inside the private flux region connecting to the active x-point)
  • Another b2fgmtry cut
  • This is the boundary between inner and outer divertors
• outer_throat: 9 (y-aligned faces in the outer SOL connecting to the active X-point)
  • This is the boundary between SOL and outer divertor
• inner_throat: 10
• outer_midplane: 11 (y-aligned faces connecting to the node closest to the outer midplane on the separatrix)
• inner_midplane: 12
• outer_target: 13 (y-aligned faces defining the outer target)
• inner_target: 14
• core_boundary: 15 (innermost x-aligned faces)
  • the inner edge of the mesh; boundary between computation region and blank spot in the core
• separatrix: 16 (x-aligned faces defining the active separatrix)
• outer_mid_plane_separatrix: 101 (point on active separatrix at outer midplane)
  • Intersection of 11 with 16
• inner_mid_plane_separatrix: 102
  • Intersection of 12 with 16
• outer_target_separatrix: 103 (point on active separatrix at outer [divertor] target)
  • Intersection of 13 with 16
• inner_target_separatrix: 104
  • Intersection of 14 with 16

[anchal-physics]

a994b4d Added x_points

[anchal-physics]

8e399d0 Added faces and core_cut

[anchal-physics]

760ef7b Added pfr_cut, throats, targets, fixed speed issue

[anchal-physics]

3d5972b Finished all but midplane subsets

[anchal-physics]

62371f2 Added inner and outer midplane subsets

[anchal-physics]

With new visualitation through GGDUtils, it has come to light that some subsets are not created properly:

SOLPSgridVisualization.pdf

Issues to fix:

• Inner throat edge is wrongly identified.
• Outer midplane are y-aligned edges but is that what we really wanted?
• Inner midplane did not get populated.

@eldond Please add more issues to fix if you can spot them.

[anchal-physics]

Actually, the inner midplane did not get populated because jxi index is not provided in our sample SOLPS files. So I'm marking this resolved. It will probably show the same feature as y-aligned (poliodal direction aligned) edges which is what is written in the OMAS data structure schema webpage. Is that a typo there?

[anchal-physics]

3062359 Inner throat edge issue fixed.

[anchal-physics]

3523eda

That was my bad. I got confused about y-aligned edge definition. I have corrected the outer midplane and inner midplane generation. The inner midplane would be populated if jxi information is provided in b2mn.dat. We can close this issue again now.

[eldond]

Beautiful.