using EFIT, Equilibrium, GuidingCenterOrbits
#Read in equilibrium
M, wall = read_geqdsk(EFIT.test_gfile, clockwise_phi=false)
#Define initial conditions
gcp = GCDeuteron(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)