using EFIT, Equilibrium, GuidingCenterOrbits
#Read in equilibrium
M, wall = read_geqdsk(EFIT.test_gfile)
#Define initial conditions
gcp = GCParticle(80.0,0.2,1.9,0.0)
#Calculate trajectory
path, stat = integrate(M, gcp, tmax=1000.0)
#Calculate orbit trajectory for single poloidal transit
o = get_orbit(M,gcp)
forked from JuliaFusion/GuidingCenterOrbits.jl
-
Notifications
You must be signed in to change notification settings - Fork 0
Julia package for calculating guiding-center orbits in magnetic equilibrium
License
StuartBenjamin/GuidingCenterOrbits.jl
Folders and files
| Name | Name | Last commit message | Last commit date | |
|---|---|---|---|---|
Repository files navigation
About
Julia package for calculating guiding-center orbits in magnetic equilibrium
Resources
License
Stars
Watchers
Forks
Releases
No releases published
Packages 0
No packages published
Languages
- Julia 100.0%