Skip to content

GC RK: VMEC vs coils comparison example #195

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Closed
krystophny wants to merge 18 commits into from
Closed

GC RK: VMEC vs coils comparison example #195

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
18 commits
Select commit Hold shift + click to select a range
4d4e8d0
Add NCSX Meiss vs VMEC Cash–Karp comparison script
krystophny Dec 1, 2025
9299885
Add VMEC vs coils RK guiding-center comparison example
krystophny Dec 1, 2025
2dc8fa1
Wrap GC comparison angles mod 2pi
krystophny Dec 1, 2025
53b634b
Add magfie_coils backend and use coils standalone in GC RK example
krystophny Dec 1, 2025
93934a8
Add fixed-step RK4 GC integrator (integmode=-2) and NCSX VMEC+coils d…
krystophny Dec 1, 2025
e3fac85
Make GC RK comparison reconstruct time axis explicitly
krystophny Dec 1, 2025
d816020
Fix uninitialized final timestep in orbit NetCDF output
krystophny Dec 1, 2025
e03a92f
Lower GC RK comparison energy via facE_al=1000
krystophny Dec 1, 2025
b101654
Use STELLOPT DIAGNO_TEST free-boundary wout_ncsx in golden tests
krystophny Dec 1, 2025
2ac206e
Compare GC NCSX VMEC vs coils for RK4 and Cash-Karp, parallelized in …
krystophny Dec 1, 2025
4851fa2
Improve line styles/colors in NCSX GC RK comparison plot
krystophny Dec 1, 2025
18efda0
Add Meiss-based GC coils orbits to VMEC vs coils RK comparison
krystophny Dec 1, 2025
5ffd6b5
Fix direct coils GC to use Biot-Savart and clean init_field
krystophny Dec 1, 2025
abf21e6
Add Cartesian GC mode for coils field with VMEC output
krystophny Dec 1, 2025
9322aa0
Fix VMEC output mapping for Cartesian coils GC orbits
krystophny Dec 1, 2025
17eede4
Wire Cash–Karp GC integrator for integmode=-1 (including coils Cartes…
krystophny Dec 1, 2025
f5f51dd
Add NCSX vacuum equilibrium and GC comparison examples
krystophny Dec 2, 2025
004ceb7
Tidy GC coils RK comparison plot and test runtime
krystophny Dec 2, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
1 change: 1 addition & 0 deletions .gitignore
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -7,6 +7,7 @@ fort.*
*~
*.nc
*.nc.*
!golden_record/test_data/wout_ncsx_vacuum.nc
*.x
*.png
*.o
Expand Down
1 change: 1 addition & 0 deletions examples/coils_vacuum/coils.simple
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
/home/ert/code/SIMPLE/golden_record/test_data/coils.c09r00.simple
15 changes: 15 additions & 0 deletions examples/coils_vacuum/simple.in
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,15 @@
&config
trace_time = 0.0001
sbeg = 0.6d0
ntestpart = 1
ntimstep = 1000
netcdffile = 'wout.nc'
isw_field_type = 5 ! Coils magfie backend (Cartesian GC)
integmode = -1 ! Cash-Karp RK5(4)
npoiper2 = 128
deterministic = .True.
facE_al = 1000.0
field_input = 'coils.simple'
contr_pp = -1000d0
output_orbits_macrostep = .True.
/
Loading
Loading