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)