Skip to content

Replace existing .geqdsk and .gpeckf files with TokaMaker DIII-D-like synthetic equilibrium#245

Open
d-burg wants to merge 7 commits intodevelopfrom
D3D-like_examples
Open

Replace existing .geqdsk and .gpeckf files with TokaMaker DIII-D-like synthetic equilibrium#245
d-burg wants to merge 7 commits intodevelopfrom
D3D-like_examples

Conversation

@d-burg
Copy link
Collaborator

@d-burg d-burg commented Feb 27, 2026

All .geqdsk and .gpeckf files have been replaced, and equil.in etol has been set as 1e-07 to allow successful DCON runs.

Copilot AI review requested due to automatic review settings February 27, 2026 20:47
Copy link
Contributor

Copilot AI left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Pull request overview

Updates the DIII-D example inputs to use a new TokaMaker DIII-D-like synthetic equilibrium (GEQDSK) and corresponding kinetic profiles (GPECKF), removing the legacy KEFIT .kin profile tables and tightening equilibrium reconstruction tolerances for DCON in several examples.

Changes:

  • Switch example equil.in files to reference TkMkr_D3Dlike_inexp=1.0.geqdsk (and add etol = 1e-7 in several cases).
  • Switch example pentrc.in files to reference TkMkr_D3Dlike_inexp=1.0.gpeckf instead of g147131.02300_DIIID_KEFIT.kin.
  • Remove legacy .kin profile files and add at least one new .gpeckf profile file.

Reviewed changes

Copilot reviewed 18 out of 26 changed files in this pull request and generated 1 comment.

Show a summary per file
File Description
docs/examples/run_kinetic_example/pentrc.in Points PENTRC to the new .gpeckf kinetic profiles.
docs/examples/run_kinetic_example/g147131.02300_DIIID_KEFIT.kin Removes legacy KEFIT kinetic profile table.
docs/examples/run_kinetic_example/equil.in Points equilibrium input to the new synthetic .geqdsk.
docs/examples/run_ideal_example/pentrc.in Points PENTRC to the new .gpeckf kinetic profiles.
docs/examples/run_ideal_example/g147131.02300_DIIID_KEFIT.kin Removes legacy KEFIT kinetic profile table.
docs/examples/run_ideal_example/equil.in Points equilibrium input to the new synthetic .geqdsk.
docs/examples/run_ideal_example/TkMkr_D3Dlike_inexp=1.0.gpeckf Adds the new synthetic kinetic profile table.
docs/examples/DIIID_resistive_example/pentrc.in Updates resistive example to reference the new kinetic profiles (via ideal example path).
docs/examples/DIIID_resistive_example/equil.in Updates resistive example equilibrium reference and adds etol = 1e-7.
docs/examples/DIIID_kinetic_example/pentrc.in Points PENTRC to the new .gpeckf kinetic profiles.
docs/examples/DIIID_kinetic_example/g147131.02300_DIIID_KEFIT.kin Removes legacy KEFIT kinetic profile table.
docs/examples/DIIID_kinetic_example/equil.in Updates equilibrium reference and adds etol = 1e-7.
docs/examples/DIIID_ideal_example/pentrc.in Points PENTRC to the new .gpeckf kinetic profiles.
docs/examples/DIIID_ideal_example/g147131.02300_DIIID_KEFIT.kin Removes legacy KEFIT kinetic profile table.
docs/examples/DIIID_ideal_example/equil.in Updates equilibrium reference and adds etol = 1e-7.

💡 Add Copilot custom instructions for smarter, more guided reviews. Learn how to get started.

Comment on lines 1 to 5
&EQUIL_CONTROL
eq_type="efit" ! Type of the input 2D equilibrium file. Accepts efit, chease, fluxgrid, transp, jsolver, lar, sol, etc.
eq_filename="g147131.02300_DIIID_KEFIT" ! Path to input file (geqdsk, etc.)
eq_filename="TkMkr_D3Dlike_inexp=1.0.geqdsk" ! Path to input file (geqdsk, etc.)

jac_type="hamada" ! Working coordinate system for all DCON and GPEC calculations. Overrides individual powers. Accepts hamada, pest, boozer, equal_arc
Copy link

Copilot AI Feb 27, 2026

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

PR description says equil.in etol was set to 1e-07 to enable successful DCON runs, but in this example equil.in still has etol=1e-6 (later in the file). If this example is also expected to run with the new synthetic equilibrium, consider updating etol here as well (or clarify in the PR description which examples require the tighter tolerance).

Copilot uses AI. Check for mistakes.
…ontained geqdsk and kinetic files

This is how the others work. It is duplication of files, but it enables easy copying of the examples for new users to copy into their own sandbox to play around with.
@logan-nc
Copy link
Contributor

logan-nc commented Mar 1, 2026

@d-burg a plot of the old and new equilibria would be nice to have on the record of this PR! Any table of typical values to compare that is easy to make would be appreciated too. Help everyone understand what is changing!

@d-burg
Copy link
Collaborator Author

d-burg commented Mar 1, 2026

@logan-nc original geqdsk:

CURRENT = 1.187 MA
BCENTR = -1.723 T
RCENTR = 1.696 m
RMAXIS = 1.759 m

Original geqdsk summary

New geqdsk:

CURRENT = 1.200 MA
BCENTR = 2.024 T
RCENTR = 1.680 m
RMAXIS = 1.744 m

TkMkr_D3Dlike_Hmode geqdsk summary

@logan-nc
Copy link
Contributor

logan-nc commented Mar 2, 2026

Bt should be negative to be DIII-D like. Standard DIII-D operation has Bt clockwise and Ip counterclockwise when looking from above.

Is this easy enough to flip in tokamaker? is there a reason you made what would be called a "reverse Bt" equilibrium?

When you update it, it would be nice to see the equilibria over plotted. It is also important to document the results thoroughly. Perhaps just pasting the terminal output and the DCON and GPEC summary plots from OMFIT would be good. Note that you can use markdown to put he terminal output in a collapsible code block so it doesn't dominate the pr comments.

Thanks again!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants