Replace existing .geqdsk and .gpeckf files with TokaMaker DIII-D-like synthetic equilibrium#245
Replace existing .geqdsk and .gpeckf files with TokaMaker DIII-D-like synthetic equilibrium#245
Conversation
There was a problem hiding this comment.
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.infiles to referenceTkMkr_D3Dlike_inexp=1.0.geqdsk(and addetol = 1e-7in several cases). - Switch example
pentrc.infiles to referenceTkMkr_D3Dlike_inexp=1.0.gpeckfinstead ofg147131.02300_DIIID_KEFIT.kin. - Remove legacy
.kinprofile files and add at least one new.gpeckfprofile 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.
| &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 |
There was a problem hiding this comment.
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).
…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.
|
@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! |
|
@logan-nc original geqdsk: CURRENT = 1.187 MA
New geqdsk: CURRENT = 1.200 MA
|
|
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! |


All .geqdsk and .gpeckf files have been replaced, and
equil.inetolhas been set as1e-07to allow successful DCON runs.