pysimple

Python interface to SIMPLE.

pysimple is a thin layer over the Fortran module simple_api. A run is configured by the &config namelist of simple.in; every keyword argument of init() is a &config name (see examples/simple_full.in):

import pysimple

pysimple.init(config="simple.in")             # exactly as simple.x
pysimple.init("wout.nc", trace_time=1e-3)     # equilibrium plus settings
pysimple.init("wout.nc", config="simple.in")  # file with overrides

start = pysimple.sample()                     # as configured by startmode
result = pysimple.trace_parallel(start)

States are (s, theta, phi, v/v0, v_par/v) in reference coordinates and have shape (5,) for one particle or (5, n) for n particles. Angles are radians. pysimple.params exposes the Fortran params module for reading derived values such as params.v0.