Example: Pulse Design Tool
This notebook demonstrates how to use the FreeGSNKE Pulse Design Tool (FPDT).
The FPDT combines the evolutive equilibrium solver introduced in previous notebooks with a virtual plasma control system (PCS) that governs the time evolution of the plasma. It enables closed-loop plasma scenario and control design studies.
The virtual PCS generates voltage requests for the active poloidal field coils, which are then passed to the evolutive solver. At present, FPDT is primarily intended for use during the flat-top phase of a discharge. While it may function during ramp-up or ramp-down phases, its performance in those regimes is not guaranteed.
The controllers implemented in the virtual PCS are inspired by the MAST-U PCS and allow control of the following quantities:
Plasma current.
Geometric shape parameters (derived from the equilibrium).
Vertical plasma position.
Coil activation times.
Coil current and voltage limits.
For more detailed information on these controllers and features, please refer to the FPDT paper on arXiv.
The FPDT can be used to:
Plan challenging new plasma scenarios.
Test new control schemes and parameter settings.
Serve as a training tool for plasma experimentalists.
In the following sections, we demonstrate how to set up the FPDT for the controlled evolution of a MAST-U–like plasma.
Generate a starting equilibrium
As in prior notebooks, we need a starting equilibrium.
This equilibrium will be the initial condition from which the FPDT will evolve the plasma.
# package
import numpy as np
import matplotlib.pyplot as plt
import time
# build machine
from freegsnke import build_machine
tokamak = build_machine.tokamak(
active_coils_path=f"../machine_configs/MAST-U/MAST-U_like_active_coils.pickle",
passive_coils_path=f"../machine_configs/MAST-U/MAST-U_like_passive_coils.pickle",
limiter_path=f"../machine_configs/MAST-U/MAST-U_like_limiter.pickle",
wall_path=f"../machine_configs/MAST-U/MAST-U_like_wall.pickle",
)
# initialise equilibrium
from freegsnke import equilibrium_update
eq = equilibrium_update.Equilibrium(
tokamak=tokamak,
Rmin=0.1, Rmax=2.0, # Radial range
Zmin=-2.2, Zmax=2.2, # Vertical range
nx=65, # Number of grid points in the radial direction (needs to be of the form (2**n + 1) with n being an integer)
ny=65, # Number of grid points in the vertical direction (needs to be of the form (2**n + 1) with n being an integer)
)
# initialise profiles
from freegsnke.jtor_update import Lao85
profiles = Lao85(
eq=eq, # equilibrium object
Ip=7.59e5, # plasma current
fvac=-0.52, # fvac = rB_{tor}
alpha=[362685, 17696], # Lao profiles parameters
beta=[0.103, 0.753], # Lao profiles parameters
alpha_logic=True,
beta_logic=True,
)
# set initial coil currents (passive structure currents start at zero here)
coil_currents = np.array([2766.10782056, 193.12768191, 4019.86873871, 4857.42357466,
-726.99282376, -1696.92521166, -77.7252396 , 266.38207493,
-73.01483716, -4169.86874271, -4518.28147417, 5.2 ])
for i, key in enumerate(tokamak.coil_names[0:12]):
eq.tokamak.set_coil_current(coil_label=key, current_value=coil_currents[i])
# initialise the static solver
from freegsnke import GSstaticsolver
GSStaticSolver = GSstaticsolver.NKGSsolver(eq)
# solve for the equilbirium
GSStaticSolver.solve(
eq=eq,
profiles=profiles,
constrain=None,
target_relative_tolerance=1e-8,
)
# plot the equilibrium
fig1, ax1 = plt.subplots(1, 1, figsize=(4, 8), dpi=70)
ax1.grid(True, which='both')
eq.plot(axis=ax1, show=False)
eq.tokamak.plot(axis=ax1, show=False)
ax1.set_xlim(0.1, 2.15)
ax1.set_ylim(-2.25, 2.25)
plt.tight_layout()
Active coils --> built from pickle file.
Passive structures --> built from pickle file.
Limiter --> built from pickle file.
Wall --> built from pickle file.
Magnetic probes --> none provided.
Resistance (R) and inductance (M) matrices --> built using actives (and passives if present).
Tokamak built.
Forward static solve SUCCESS. Tolerance 5.72e-09 (vs. requested 1.00e-08) reached in 29/100 iterations.
Prepare inputs for the PCS class
Here, we will set up the inputs (lists and dictionaries) for each of the (internal) controller classes within the virtual PCS (e.g. plasma, shape, virtual circuits).
We start by defining which coils have which purpose.
# COIL NAMES AND SHAPE TARGET NAMES
# active coils: must match the order they appear in FreeGSNKE machine description
active_coils = ['Solenoid', 'PX', 'D1', 'D2', 'D3', 'Dp', 'D5', 'D6', 'D7', 'P4', 'P5', 'P6']
# control coils: a subset of the active coils to be used for plasma and shape parameter control (i.e. we exclude vertical control coil P6 here)
ctrl_coils = ['Solenoid', 'PX', 'D1', 'D2', 'D3', 'Dp', 'D5', 'D6', 'D7', 'P4', 'P5']
# vertical control coil
vertical_coils=["P6"]
# names of the plasma targets (this will be clear later, equal to the number of solenoid coils)
plasma_targets = ['plasma']
All of the following inputs are going to be specified as waveforms, i.e. they will be a dictionary of times and vals.
For each array/list of times (in seconds), we must specify a corresponding list of list/arrays in vals.
Most waveform quantities (e.g. *_ref, *_ff, blend) will use linear interpolation inside the PCS class internally, to enable querying at a given time \(t\) within the simulation.
Some other quantities (e.g. PID gains, virtual circuits, and other arrays), will use previous-value (or step) interpolation where the value defined at time tmin is used over the interval \( [t_\mathrm{min},\, t_\mathrm{next})\), after which the value defined at tnext is applied.
We will denote which quantity uses which type of interpolation with linear or step.
The following two cells provide an example of how to set these two different waveforms types and what they look like when plotted.
# waveform dictionary structure (for linearly interpolated quantity)
# e.g. this could be the desired radial position of the X-point over time
waveform = {
"times": np.array([0.0, 0.05, 0.15, 0.35, 0.45, 0.5]),
"vals": np.array([0.57, 0.57, 0.5, 0.5, 0.57, 0.57]),
}
# plot what it looks like
fig, ax = plt.subplots(
nrows=1,
ncols=1,
figsize=(8, 4),
dpi=80
)
ax.plot(waveform['times'], waveform['vals'], color='k', linewidth=1, linestyle="--", marker="x", markersize=7, label="linearly interpolated waveform")
ax.set_xlabel(r"Shot time [$s$]")
ax.set_ylabel("Scalar quantity [units]")
ax.grid()
ax.legend()
plt.tight_layout()
plt.show()
# waveform dictionary structure (for step interpolated quantity)
# e.g. this could be a proportional PID gain over time
waveform = {
"times": np.array([0.0, 0.15, 0.3, 0.5]),
"vals": np.array([100, 150, 75, 100]),
}
# plot what it looks like
fig, ax = plt.subplots(
nrows=1,
ncols=1,
figsize=(8, 4),
dpi=80
)
ax.step(waveform['times'], waveform['vals'], where='post',
color='k', linewidth=1, linestyle="--", marker="x", markersize=7,
label="step interpolated waveform")
ax.set_xlabel(r"Shot time [$s$]")
ax.set_ylabel("Scalar quantity [units]")
ax.grid()
ax.legend()
plt.tight_layout()
plt.show()
Before we define each of the controller settings, let us define some global times, i.e. the start and end of the simulation say.
# choose some min amd max simulation times
tmin = 0.0
tmax = 0.25
First, we define the desired control settings for the plasma category.
These waveforms govern how the plasma current, \(I_p\), is controlled—either via feedback (FB), feedforward (FF), or a combination of both.
The following quantities must be specified:
ip_ref: Plasma current feedback reference waveform [A] (linear).vloop_ff: Loop voltage feedforward reference waveform [V·s⁻¹] (linear).blend: Control blending parameter [dimensionless] (linear):1→ purely FB control0→ purely FF controlvalues in
(0, 1)→ mixed FB/FF control.
k_prop: Proportional gain waveform for the FB PID controller [s⁻¹] (step).k_int: Integral gain waveform for the FB PID controller [s⁻²] (step).k_deriv: Derivative gain waveform for the FB PID controller (no units) (step).M_solenoid: Mutual inductance between the plasma and solenoid (required for FF control) [V·s·A⁻¹] (step).
In this particular example, we will tell the PCS to hold the plasma current constant at it’s starting value.
# PLASMA CATEGORY DATA
# build the dictionary which will be passed to the virtual PCS (all keys below are required)
plasma_data = {}
# plasma current feedback (FB) reference waveform
# here we hold the value for the initial equilibrium constant
plasma_data["ip_ref"] = {
'times': np.array([tmin, tmax]),
'vals': np.array([eq._profiles.Ip, eq._profiles.Ip])
}
# loop voltage feedforward (FF) reference waveform (can specify the desired loop voltage on the plasma)
# (here we don't use a loop voltage as we will use the FB reference waveform above)
plasma_data["vloop_ff"] = {
'times': np.array([tmin, tmax]),
'vals': np.array([0, 0])
}
# blending waveform: tells controller to use purely FB control (1), purely FF control (0), or a mixture (between 0 and 1)
# we are doing pure FB control here
plasma_data["ip_blend"] = {
'times': np.array([tmin, tmax]),
'vals': np.array([1, 1])
}
# proportional gain for the FB PID controller
plasma_data["k_prop"] = {
'times': np.array([tmin]),
'vals': np.array([-5*4])
}
# integral gain for the FB PID controller
plasma_data["k_int"] = {
'times': np.array([tmin]),
'vals': np.array([-50*2])
}
# derivative gain for the FB PID controller
plasma_data["k_deriv"] = {
'times': np.array([tmin]),
'vals': np.array([0])
}
# mutual inductance between plasma and solenoid: required for FF control
# not used here as we do not set a loop voltage
plasma_data["M_solenoid"] = {
"times": np.array([0]),
"vals": np.array([1.0])
}
plasma_data.keys()
dict_keys(['ip_ref', 'vloop_ff', 'ip_blend', 'k_prop', 'k_int', 'k_deriv', 'M_solenoid'])
Next, we define the desired control settings for the shape category.
First, we define a function that extracts the relevant shape parameters from the equilibrium object (here referred to as plasma_descriptors). In this example, we extract the following quantities:
Rin: Inboard midplane radius.Rout: Outboard midplane radius.Rx: Radial position of the lower X-point.Zx: Vertical position of the lower X-point.
These are taken directly from the equilibrium object (without measurement noise). You could instead build synthetic diagnostics (e.g. magnetic probes) to take readings around the equilibrium (with measurement noise) and then reconstruct the plasma shape based on these.
# define the descriptors function (it should return an array of values and take in an eq object)
# outputs ("measurements") from this function will be passed to the PCS later on during simulation
def plasma_descriptors(eq):
# inboard/outboard midplane radii
RinRout = eq.innerOuterSeparatrix()
# find lower X-point
# define a "box" in which to search for the lower X-point
XPT_BOX = [[0.33, -0.88], [0.95, -1.38]]
# mask those points
xpt_mask = (
(eq.xpt[:, 0] >= XPT_BOX[0][0])
& (eq.xpt[:, 0] <= XPT_BOX[1][0])
& (eq.xpt[:, 1] <= XPT_BOX[0][1])
& (eq.xpt[:, 1] >= XPT_BOX[1][1])
)
xpts = eq.xpt[xpt_mask, 0:2].squeeze()
if xpts.ndim > 1 and xpts.shape[0] > 1:
opt = eq.opt[0, 0:2]
dists = np.linalg.norm(xpts - opt, axis=1)
idx = np.argmin(dists) # index of closest point
Rx, Zx = xpts[idx, :]
else:
Rx, Zx = xpts
return np.array([RinRout[0], RinRout[1], Rx, Zx])
# give these descriptors some names for use in the PCS
ctrl_targets = ['Rin','Rout', 'Rx', 'Zx']
Next, we specify the waveforms used to control these shape parameters throughout the simulation.
Not all shape parameters need to be defined in units of metres—these are simply illustrative examples.
For each shape parameter, the following waveforms must be specified:
ff: Feedforward (FF) waveform (linear).
Not used in this example, but can be employed to impose a fixed offset on the parameter [m].ref: Feedback (FB) reference waveform (linear).
These are the desired target values of the parameter over the course of the simulation [m].blend: Control blending waveform specifying the relative contribution of FF and FB control dimensionless.0→ purely FF control1→ purely FB controlvalues in
(0, 1)→ mixed FF/FB control.
k_prop: Proportional gain waveform for the FB PID controller [s⁻¹] (step).k_int: Integral gain waveform for the FB PID controller (not used here) [s⁻²] (step).k_deriv: Derivative gain waveform for the FB PID controller (no units) (step).damping: Damping waveform (not used here) [dimensionless] (step).
In the following, we will tell the controllers to:
hold both
RinandRoutconstant.hold
Rxconstant for a short period, before ramping down (to increase plasma triangularity), then ramp back up to the original value.allow
Zxto be uncontrolled.
# SHAPE/DIVERTOR CATEGORY
# build the dictionary which will be passed to the virtual PCS (all keys below are required)
shape_data = {}
# waveforms for Rin
shape_data["Rin"] = {}
shape_data["Rin"]["ff"] = {"times": np.array([tmin, tmax]), "vals": np.array([0.0, 0.0])}
shape_data["Rin"]["ref"] = {"times": np.array([tmin, tmax]), "vals": np.array([0.28, 0.28])}
shape_data["Rin"]["blend"] = {"times": np.array([tmin, tmax]), "vals": np.array([1.0, 1.0])}
shape_data["Rin"]["k_prop"] = {"times": np.array([tmin]), "vals": np.array([120])}
shape_data["Rin"]["k_int"] = {"times": np.array([tmin]), "vals": np.array([0.0])}
shape_data["Rin"]["k_deriv"] = {"times": np.array([tmin]), "vals": np.array([0.0])}
# waveforms for Rout
shape_data["Rout"] = {}
shape_data["Rout"]["ff"] = {"times": np.array([tmin, tmax]), "vals": np.array([0.0, 0.0])}
shape_data["Rout"]["ref"] = {"times": np.array([tmin, tmax]), "vals": np.array([1.35, 1.35])}
shape_data["Rout"]["blend"] = {"times": np.array([tmin, tmax]), "vals": np.array([1.0, 1.0])}
shape_data["Rout"]["k_prop"] = {"times": np.array([tmin]), "vals": np.array([170])}
shape_data["Rout"]["k_int"] = {"times": np.array([tmin]), "vals": np.array([0.0])}
shape_data["Rout"]["k_deriv"] = {"times": np.array([tmin]), "vals": np.array([0.0])}
# waveforms for Rx
shape_data["Rx"] = {}
shape_data["Rx"]["ff"] = {"times": np.array([tmin, tmax]), "vals": np.array([0.0, 0.0])}
shape_data["Rx"]["ref"] = {"times": np.array([tmin, 0.025, 0.075, 0.175, 0.225, tmax]), "vals": np.array([0.57, 0.57, 0.5, 0.5, 0.57, 0.57])} # we vary Rx
shape_data["Rx"]["blend"] = {"times": np.array([tmin, tmax]), "vals": np.array([1.0, 1.0])}
shape_data["Rx"]["k_prop"] = {"times": np.array([tmin]), "vals": np.array([120])}
shape_data["Rx"]["k_int"] = {"times": np.array([tmin]), "vals": np.array([0.0])}
shape_data["Rx"]["k_deriv"] = {"times": np.array([tmin]), "vals": np.array([0.0])}
# waveforms for Zx
shape_data["Zx"] = {}
shape_data["Zx"]["ff"] = {"times": np.array([tmin, tmax]), "vals": np.array([0.0, 0.0])}
shape_data["Zx"]["ref"] = {"times": np.array([tmin, tmax]), "vals": np.array([-1.24, -1.24])}
shape_data["Zx"]["blend"] = {"times": np.array([tmin, tmax]), "vals": np.array([0.0, 0.0])}
shape_data["Zx"]["k_prop"] = {"times": np.array([tmin]), "vals": np.array([120])}
shape_data["Zx"]["k_int"] = {"times": np.array([tmin]), "vals": np.array([0.0])}
shape_data["Zx"]["k_deriv"] = {"times": np.array([tmin]), "vals": np.array([0.0])}
shape_data.keys()
dict_keys(['Rin', 'Rout', 'Rx', 'Zx'])
Next, we define the schedule of virtual circuits (VCs) in the virtual circuits category.
For each shape parameter in the shape category, a VC is required to map the discrepancy between the target value (from the shape controller) and the measured value (from the equilibrium) to requests for the poloidal field (PF) coil currents. The VCs used here were obtained using the procedure described in the example notebook on VC construction. Recall that VCs are defined with units [m/A].
Each VC is an array with length equal to the number of ctrl_coils, with entries ordered consistently with that list. When specified in the waveform dictionary, VCs are previous-value interpolated by the PCS class. For example, for the Rin controller, the VC defined at time tmin is used over the interval \( [t_\mathrm{min},\, t_\mathrm{max})\), after which the VC defined at tmax is applied. Note here that we use only one VC for each shape parameter as the shape changes are not that significant. If more complex shape changes are required, then more VCs (linearised around the expected plasma shape), will be required.
Note that the first entry of each shape-parameter VC is zero. This ensures that the solenoid does not contribute to shape control and is reserved exclusively for plasma current, \(I_p\), control.
In addition to the shape VCs, a final VC is defined for the plasma category. This VC maps the requested change in \(I_p\) from the plasma controller to PF coil current requests (well in this case, only a request to the solenoid).
Finally, optional pre-programmed feedforward (FF) requests for the PF coil currents may be specified. These allow the coil currents to be set directly if desired. In this example, no such FF requests are used, and the coil currents are driven purely via feedback control.
# VIRTUAL CIRCUITS CATEGORY (provide VCs calculated elsewhere)
# build the dictionary which will be passed to the virtual PCS (all keys below are required)
circuits_data = {}
circuits_data["Rin"] = {
"times": np.array([tmin]),
"vals": [
np.array([0, 3.1394e+04, 8.5990e+03, -8.8300e+02, -1.2460e+03, -6.6690e+03, 1.2660e+03, 2.9770e+03, 3.5350e+03, 8.3110e+03, -7.8520e+03]),
]
}
circuits_data["Rout"] = {
"times": np.array([tmin]),
"vals": [
np.array([0, -9.9700e+02, 1.1790e+03, 3.9700e+02, -2.0000e+02, -1.8970e+03, -5.5500e+02, -2.2640e+03, -1.5860e+03, -1.9750e+03, 5.5800e+03]),
]
}
circuits_data["Rx"] = {
"times": np.array([tmin]),
"vals": [
np.array([0, -3.0021e+04, 3.3980e+03, 8.6900e+03, 5.8310e+03, 1.5858e+04, 2.4000e+01, -1.2100e+02, -2.2160e+03, -9.5290e+03, 1.3570e+03]),
]
}
circuits_data["Zx"] = {
"times": np.array([tmin]),
"vals": [
np.zeros(11),
]
}
circuits_data["plasma"] = {
"times": np.array([tmin]),
"vals": [
np.array([1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]),
]
}
# store the coil order for future reference
circuits_data["coil_order"] = ctrl_coils
# define any feedforward coil current drives on each control coil (no drives used here)
zeros_dict = {"times": np.array([tmin, tmax]), "vals": np.array([0.0, 0.0])}
for coil in ctrl_coils: # linearly interpolated
circuits_data[coil+"_ref"] = zeros_dict
circuits_data.keys()
dict_keys(['Rin', 'Rout', 'Rx', 'Zx', 'plasma', 'coil_order', 'Solenoid_ref', 'PX_ref', 'D1_ref', 'D2_ref', 'D3_ref', 'Dp_ref', 'D5_ref', 'D6_ref', 'D7_ref', 'P4_ref', 'P5_ref'])
The cell above shows how to explicitly provide your own VCs for each shape parameter (and the plasma).
The following cell shows how you can generate a VC schedule using a set of FreeGSNKE equilbiria that you might have to hand.
Here, we build the VCs around the initial equilibrium generated above (it essentially builds almost the same VCs as we hard coded above).
Here, the VCGenerator class is initialised with:
a
solverobject (i.e. the static GS solver).a
target_calculatorto compute the plasma shape parameters (plasma_descriptors())the list of corresponding shape parameter names (
target_names).the full list of coils available for shape control (
coils), only needed ifgenerate_fixed_scheduleis going to be used (as below)the list of coils and targets to use in a VC computation / matrix inversion (
targets_calc,coils_calc)the list of shape targets that will use the new computed VCs, overwriting any pre-scheduled ones (
targets_ctrl)
The method generate_fixed_schedule then takes a list of eq and profile objects (at chosen times), computing VCs in the required waveform format above. It uses the targets_calc, targets_ctrl, coils_calc, and coils choices already given to VCGenerator above.
The schedule dictionary containts the VC waveforms only for the shape parameters (and the coil ordering), these need to be combined with the plasma VC and any FF drive waveforms.
Refer back to the virtual_circuits notebook for more details on VC calculations.
# from freegsnke.control_loop.vc_provider import VCGenerator
# # VC calculator object
# VCG = VCGenerator(
# solver=GSStaticSolver, # solver object
# target_calculator=plasma_descriptors, # for calculating shape parameters
# target_names=ctrl_targets, # names of shape parameters
# coils=ctrl_coils, # coils available for shape control (only needed for generate_fixed_schedule)
# coils_calc=ctrl_coils[1:], # subset of coils used for VC computation (solenoid excluded here)
# targets_calc=ctrl_targets, # shapes included in the VC calculation
# targets_ctrl=["Rin", "Rout", "Rx"], # shapes to control (others have VCs set to zero)
# vc_update_rate=0.05, # how often (in seconds) to recompute the VCs during a simulation (if initialised in the PCS class, see comment input to PCS class below)
# n_vc_workers=4, # number of works used to build the VCs
# )
# # list of start times, equilibria, and profiles used to generate the schedule (here we use the initial condition only, i.e one set of VCs)
# vc_times = [0.0]
# vc_eq_list = [eq]
# vc_prof_list = [profiles]
# # generate schedule (shape-target VCs only -- the plasma VC and feedforward
# # coil drives are not produced here and are already defined separately above)
# circuits_data = VCG.generate_fixed_schedule(
# times=vc_times,
# eq_list=vc_eq_list,
# profile_list=vc_prof_list,
# tikhonov_lambda=1e-12*np.ones(len(ctrl_coils[1:])), # optional regularisation of the sensitivity matrix inversion (for calculating VCs)
# verbose=True, # print output
# )
# # build the rest of the dictionary which will be passed to the virtual PCS (all keys below are required)
# circuits_data["plasma"] = {
# "times": np.array([tmin]),
# "vals": [
# np.array([1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]),
# ]
# }
# # define any feedforward coil current drives on each control coil (no drives used here)
# zeros_dict = {"times": np.array([tmin, tmax]), "vals": np.array([0.0, 0.0])}
# for coil in ctrl_coils: # linearly interpolated
# circuits_data[coil+"_ref"] = zeros_dict
# circuits_data.keys()
Next, we define the inputs for the systems category.
This controller specifies operational limits for the coils listed in ctrl_coils, including:
min_coil_curr_lims: Minimum allowable coil currents [A] (step).max_coil_curr_lims: Maximum allowable coil currents [A] (step).max_coil_curr_ramp_lims: Maximum allowable coil current ramp rates [A·s⁻¹] (step).
There is also the option to define feedforward (FF) perturbations (linear) to the control coil currents—similar to the FF drives available in the virtual circuits category. These are not used here.
# SYSTEMS CATEGORY
# build the dictionary which will be passed to the virtual PCS (all keys below are required)
systems_data = {}
# define the min coil current limits and ramp rate limits
systems_data["min_coil_curr_lims"] = {
'times': [0.0],
'vals': [np.array([-10000, -7000, -9000, -9000, -9000, -9000, -9000, -9000, -9000, -10000, -10000])],
}
# define the max coil current limits
systems_data["max_coil_curr_lims"] = {
'times': [0.0],
'vals': [np.array([10000, 7000, 9000, 9000, 9000, 9000, 9000, 9000, 9000, 0, 0])],
}
# define the max ramp rate limits for the coils
systems_data["max_coil_curr_ramp_lims"] = {
'times': [0.0],
'vals': [1e12*np.ones(len(ctrl_coils))],
}
# define the control coil current perturbations
for name in ctrl_coils: # linearly interpolated
systems_data[name+"_pert"] = zeros_dict
systems_data.keys()
dict_keys(['min_coil_curr_lims', 'max_coil_curr_lims', 'max_coil_curr_ramp_lims', 'Solenoid_pert', 'PX_pert', 'D1_pert', 'D2_pert', 'D3_pert', 'Dp_pert', 'D5_pert', 'D6_pert', 'D7_pert', 'P4_pert', 'P5_pert'])
Next, we define the inputs for the PF category.
This category specifies the electrical properties, gains, and operational limits of the poloidal field (PF) coils. These inputs are used by the virtual PCS to convert coil current requests into physically consistent voltage commands.
The following quantities must be provided:
R_matrix: Coil resistance array for the PF coils [Ω] (step).
Defined as a time-dependent waveform and restricted to the coils listed inctrl_coils.M_FF_matrix: Coil mutual inductance matrix used in feedforward (FF) terms [Vs/A] (step).
This accounts for inductive coupling between PF coils when computing FF voltage requests.M_FB_matrix: Coil mutual inductance matrix used in feedback (FB) terms [Vs/A] (step). This accounts for inductive coupling between PF coils when computing FF voltage requests. Often it is identical toM_FF_matrix, but is provided separately for flexibility.coil_gains: Feedback gain applied to each PF coil [s] (step).
These gains scale the feedback voltage contributions on a per-coil basis.coil_voltage_lims: Maximum allowable voltages for each PF coil [V] (step).
These limits are enforced by the PCS when generating coil voltage commands.coil_voltage_slew_lims: Maximum allowable voltage ramp rates for each PF coil [V·s⁻¹] (step).
In this example, very large values are used, effectively disabling slew-rate limiting.
All quantities are defined as time-dependent waveforms, even when constant, to allow future extension to time-varying PF system parameters.
# PF CATEGORY
# build the dictionary which will be passed to the virtual PCS (all keys below are required)
pf_data = {}
# coil resistances
pf_data["R_matrix"] = {
'times': [0.0],
'vals': [tokamak.coil_resist[0:11]],
}
# coil mutual inductances (on feedforward terms)
pf_data["M_FF_matrix"] = {
'times': [0.0],
'vals': [tokamak.coil_self_ind[0:11,0:11]],
}
# coil mutual inductances (on feedback terms)
pf_data["M_FB_matrix"] = {
'times': [0.0],
'vals': [tokamak.coil_self_ind[0:11,0:11]],
}
# gains on the coils (for feedback term)
pf_data["coil_gains"] = {
'times': [0.0],
'vals': [0.015*np.ones(len(ctrl_coils))],
}
# limits on voltages in the coils
pf_data["coil_voltage_lims"] = {
"times": np.array([0.0]),
'vals': [[2000, 750, 750, 750, 750, 750, 750, 750, 750, 750, 750]],
}
# limits on voltage ramp rates in the coils (we set no limits here)
pf_data["coil_voltage_slew_lims"] = {
'times': [0.0],
'vals': [1e12*np.ones(len(ctrl_coils))],
}
pf_data.keys()
dict_keys(['R_matrix', 'M_FF_matrix', 'M_FB_matrix', 'coil_gains', 'coil_voltage_lims', 'coil_voltage_slew_lims'])
Next, we define the inputs for the vertical category.
This category configures the feedback controller responsible for stabilising and controlling the vertical position of the plasma column using the vertical control coil (here the P6 coil).
The following quantities must be specified:
z_ref: Plasma vertical position reference waveform [m] (linear).
This defines the desired vertical position of the plasma as a function of time.k_prop: Proportional gain waveform for the vertical position feedback controller [s⁻¹] (step).
This gain determines the strength of the restoring force in response to vertical displacement.k_deriv: Derivative gain waveform for the vertical position feedback controller [s⁻²] (step). This term provides damping by responding to the rate of change of the vertical position.
All inputs are defined as time-dependent waveforms, allowing for time-varying vertical control behaviour if required.
In this example, we steadily increase the vertical position to before holding it constant thereafter.
# VERTICAL CATEGORY DATA
# build the dictionary which will be passed to the virtual PCS (all keys below are required)
vertical_data = {}
# plasma vertical position reference
vertical_data["z_ref"] = {"times": np.array([tmin, tmin+0.125, tmax]), "vals": np.array([0.0, 0.01, 0.01])}
# proportional gain
vertical_data["k_prop"] = {"times": np.array([tmin]), "vals": np.array([0.025])}
# derivative gain
vertical_data["k_deriv"] = {"times": np.array([tmin]), "vals": np.array([5e-7])}
vertical_data.keys()
dict_keys(['z_ref', 'k_prop', 'k_deriv'])
Finally, we define the inputs for the coil activations category.
This category specifies when each active poloidal field (PF) coil is enabled or disabled during the simulation. Coil activations are represented as time-dependent waveforms and are applied on a per-coil basis.
For each coil listed in active_coils, the following entry is defined:
<coil_name>_activation: Coil activation waveform (step).
A value of:1.0indicates that the coil is active and available to the PCS.0.0indicates that the coil is inactive and excluded from PCS.
In this example, all coils are activated for the entire simulation interval \([t_\mathrm{min},\, t_\mathrm{max}]\), using constant activation waveforms. This structure allows coils to be selectively enabled or disabled at different times in more advanced scenarios.
# COIL ACTIVATIONS CATEGORY
# build the dictionary which will be passed to the virtual PCS (all keys below are required)
coil_activation_data = {}
# coil activation times
default_coil_dict = {"times": np.array([tmin]), "vals": np.array([1.0])}
for name in active_coils:
coil_activation_data[name+"_activation"] = default_coil_dict
coil_activation_data.keys()
dict_keys(['Solenoid_activation', 'PX_activation', 'D1_activation', 'D2_activation', 'D3_activation', 'Dp_activation', 'D5_activation', 'D6_activation', 'D7_activation', 'P4_activation', 'P5_activation', 'P6_activation'])
Initialise the virtual PCS class
Now that the inputs have been prepared we can initialise the class and use in-built functions to view or plot the waveform data in each control category. This will build internal classes for each of the above listed controllers.
# FIRE UP THE PLASMA CONTROL SYSTEM
from freegsnke.control_loop.pcs import PlasmaControlSystem
PCS = PlasmaControlSystem(
plasma_data=plasma_data,
shape_data=shape_data,
shape_control_mode="PID",
circuits_data=circuits_data,
systems_data=systems_data,
pf_data=pf_data,
vertical_data=vertical_data,
coil_activation_data=coil_activation_data,
active_coils=active_coils,
ctrl_coils=ctrl_coils,
vertical_coils=vertical_coils,
ctrl_targets=ctrl_targets,
plasma_target=plasma_targets,
# vc_generator=VCG, # optional: pass a VCGenerator (see commented-out cell above) to recompute VCs during the simulation below, else use pre-set schedule
)
PCS.CoilActivationController.data
{'Solenoid_activation': {'times': array([0.]), 'vals': array([1.])},
'PX_activation': {'times': array([0.]), 'vals': array([1.])},
'D1_activation': {'times': array([0.]), 'vals': array([1.])},
'D2_activation': {'times': array([0.]), 'vals': array([1.])},
'D3_activation': {'times': array([0.]), 'vals': array([1.])},
'Dp_activation': {'times': array([0.]), 'vals': array([1.])},
'D5_activation': {'times': array([0.]), 'vals': array([1.])},
'D6_activation': {'times': array([0.]), 'vals': array([1.])},
'D7_activation': {'times': array([0.]), 'vals': array([1.])},
'P4_activation': {'times': array([0.]), 'vals': array([1.])},
'P5_activation': {'times': array([0.]), 'vals': array([1.])},
'P6_activation': {'times': array([0.]), 'vals': array([1.])}}
# each controller should contain a copy of the waveform dictionaries internally
PCS.PlasmaController.data
{'ip_ref': {'times': array([0. , 0.25]), 'vals': array([759000., 759000.])},
'vloop_ff': {'times': array([0. , 0.25]), 'vals': array([0, 0])},
'ip_blend': {'times': array([0. , 0.25]), 'vals': array([1, 1])},
'k_prop': {'times': array([0.]), 'vals': array([-20])},
'k_int': {'times': array([0.]), 'vals': array([-100])},
'k_deriv': {'times': array([0.]), 'vals': array([0])},
'M_solenoid': {'times': array([0]), 'vals': array([1.])}}
# it will also contain a dictionary of functions that has linearly (or previous value) interpolated the various waveforms
PCS.PlasmaController.interpolants
{'ip_ref': <freegsnke.control_loop.useful_functions.ConstantInterpolant at 0x7f7ca98d36a0>,
'vloop_ff': <freegsnke.control_loop.useful_functions.ConstantInterpolant at 0x7f7ca98d06d0>,
'ip_blend': <freegsnke.control_loop.useful_functions.ConstantInterpolant at 0x7f7ca98d0640>,
'k_prop': <freegsnke.control_loop.useful_functions.ConstantInterpolant at 0x7f7ca98d2cb0>,
'k_int': <freegsnke.control_loop.useful_functions.ConstantInterpolant at 0x7f7ca98d09d0>,
'k_deriv': <freegsnke.control_loop.useful_functions.ConstantInterpolant at 0x7f7ca98d2470>,
'M_solenoid': <freegsnke.control_loop.useful_functions.ConstantInterpolant at 0x7f7ca98d1c90>}
# # it should also be possible to plot the data (uncomment following cell)
# # green shading indiciates times that FB control is ON, yellow that FF is ON, a mix that both are ON, and white that there is no control
# PCS.PlasmaController.plot_data(tmin=tmin, tmax=tmax)
# the "run_control" method will be used internally later on but can also be called explicitly if required
PCS.PlasmaController.run_control
<bound method PlasmaController.run_control of <freegsnke.control_loop.plasma_category.PlasmaController object at 0x7f7ca98d3250>>
Finally, suppose you wish to modify the shot setup data after initialising the PCS class.
To do this, simply edit the waveforms required directly in the .data attribute of the controller of interest and then call .update_interpolants(). This has to be done to ensure that the controller refreshes any stale interpolants from a prior initialisation.
Uncomment the code below to see how it works for an example in the ShapeController.
# # update the data entry with your new waveform
# new_waveform = {"times": np.array([tmin, 0.05, 0.15, 0.35, 0.45, tmax]), "vals": np.array([0.28, 0.28, 0.30, 0.30, 0.28, 0.28])} # we vary Rin
# PCS.ShapeController.data["Rin"]["ref"] = new_waveform
# # call the update function to refresh inteprolants (essential)
# PCS.ShapeController.update_interpolants()
# # plot to see new waveform
# PCS.ShapeController.plot_data(targ="Rin", tmin=tmin, tmax=tmax)
Initialise the nonlinear solver object
Now choose the desired settings for the nonlinear solver object - recall the prior example notebook on this.
The simulation timestep was chosen based on the vertical instability timescale in the output below (you need to always choose a value that is 5-10x smaller than this timescale) - this keeps the simulation numerically stable.
This example explicitly uses coupling-based mode selection. The coupling ranking is calculated about the initial equilibrium and can evolve as the plasma moves or changes shape. The 30-mode reduction used here is therefore specific to this trajectory; for larger shape changes, validate the retained modes across representative equilibria or use mode_selection="timescale".
The simulation timestep was chosen based on the vertical instability timescale in the output below (you need to always choose a value that is 5-10x smaller than this timescale) - this keeps the simulation numerically stable.
The new n_linearization_workers argument controls the number of independent worker processes used to build both the current-response (dIydI) and profile-parameter-response (dIydtheta) Jacobian columns. 1 retains serial construction; set it higher to parallelise both the initial linearisation and later relinearisations. This setting is scoped to Jacobian construction and does not limit the inner threading of the GS solves used during ordinary evolution steps.
from freegsnke import nonlinear_solve
stepping = nonlinear_solve.nl_solver(
eq=eq,
profiles=profiles,
GSStaticSolver=GSStaticSolver,
full_timestep=5e-4,
plasma_resistivity=1e-7,
mode_selection="coupling",
fix_n_vessel_modes=30,
n_linearization_workers=4, # set >1 to parallelise Jacobian construction
plasma_descriptor_function=plasma_descriptors,
)
-----
Checking that the provided 'eq' and 'profiles' are a GS solution...
Forward static solve SUCCESS. Tolerance 5.72e-09 (vs. requested 1.00e-08) reached in 0/100 iterations.
-----
Instantiating nonlinear solver objects...
done.
-----
Identifying mode selection criteria...
'fix_n_vessel_modes' option selected --> passive structure modes that couple most to the strongest passive structure mode are being selected.
-----
Initial mode selection:
Active coils
total selected = 12 (out of 12)
Passive structures
30 selected using 'fix_n_vessel_modes'
Total number of modes = 42 (12 active coils + 30 passive structures)
(Note: some additional modes may be removed after Jacobian calculation if 'mode_removal=True')
-----
Building the 1513 x 43 Jacobian (dIy/dI) of plasma current density (inside the LCFS) with respect to all metal currents and the total plasma current.
Building the 1513 x 4 Jacobian (dIy/dtheta) of plasma current density (inside the LCFS) with respect to all plasma current density profile parameters within Jtor.
Built the 4 x 43 Jacobian (ds/dI) of plasma descriptors with respect to all metal currents and the total plasma current.
Built the 4 x 4 Jacobian (ds/dtheta) of plasma descriptors with respect to all plasma current density profile parameters within Jtor.
-----
Stability paramters:
Deformable plasma metrics:
Growth rate = [24.57695603] [1/s]
Instability timescale = [0.04068852] [s]
Inductive stability margin = [1.07876721+0.j]
Rigid plasma metrics:
Leuer parameter (ratio of stabilsing to de-stabilising force gradients):
between all metals and all metals = 2.0951019313803987
between all metals and active metals = 2.0951019313803987
between passive metals and active metals = 2.0950840014147594
-----
Evolutive solver timestep:
Solver timestep 'dt_step' has been set to 0.0005 as requested.
Ensure it is smaller than the growth rate else you may find numerical instability in any subsequent evoltuive simulations!
-----
Define FPDT simulation parameters
Next set the key simulation parameters for the FPDT.
# FPDT SETUP
# number of simulation time steps
n = 500
# simulation time step size (must be smaller than the vertical instability timescale)
dt = stepping.dt_step
# PCS time step (this must be whole fraction of the simulation time step)
# it governs the frequency at which the PCS is called
dt_PCS = dt/2
# # starting time (leave as it is)
# tmin = 0.0
# automatically calculates time array and sets time step in solver object
t_end = tmin + n*dt
times = np.arange(tmin, t_end, dt)
# (re-)initialise the dynamic solver with the initial eq and profiles
stepping.initialize_from_ICs(eq, profiles)
Before moving forward, we note that there are two time-dependent quantities that have not been set explicitly here. In this example, we hold the:
plasma resistivity constant.
plasma current density profiles constant.
For modelling real plasma discharges, both of these quantities should passed as time-dependent inputs to the time-stepping loop below. For further details, see the earlier example notebooks where time-dependent resistivity and profile evolution are demonstrated.
FPDT Simulation
In the following cell, we initialise lists/arrays to store any equilbirium related data we wish to view after the simulation. See the FreeGSNKE example notebook (example03 - extracting_equilibrium_quantites) for a list of these and how to extract them.
Note that some quantities can be computationally costly to extract and have not been optimised for speed yet!
# equilibrium-related data storage
dynamic_psi = np.zeros((stepping.eq1.psi().shape[0], stepping.eq1.psi().shape[1], len(times))) # total poloidal flux
dynamic_limiter_flag = np.zeros(len(times)) # flag if plasma is limited
dynamic_psi_boundary = np.zeros(len(times)) # poloidal flux on plasma boundary
dynamic_xpts = [] # list of X-point locations and associated poloidal flux
dynamic_opts = [] # list of O-point locations and associated poloidal flux
dynamic_currents = np.zeros((len(stepping.vessel_currents_vec),len(times))) # PF coil and vessel eigenmode currents
dynamic_ip = np.zeros(len(times)) # total plasma current
dynamic_shape_targets = np.zeros((len(ctrl_targets),len(times))) # shape parameters
dynamic_triangularity = np.zeros(len(times)) # plasma triangularity
linear_growth_rate = np.zeros(len(times)) # linearised (deformable) growth rate estimate
# PCS-related data storage
V_approved = np.zeros((len(ctrl_coils)+1,len(times))) # final PF coil voltages from the PCS class (passed into solver)
ip_hist = np.zeros(len(times)) # integral term feeding into Plasma Category PID FB controller (at prior time step)
ip_err = np.zeros(len(times)) # (error) for derivative term feeding into Plasma Category PID FB controller (at prior time step)
T_err = np.zeros((len(ctrl_targets),len(times))) # proportional term feeding into Shape Category PID FB controller (at prior time step)
T_hist = np.zeros((len(ctrl_targets),len(times))) # integral term feeding into Shape Category PID FB controller (at prior time step)
I_approved = np.zeros((len(ctrl_coils),len(times))) # approved PF coil currents from prior time step feeding into System Category
coil_resists = np.zeros((len(active_coils),len(times))) # PF coil resistances (to tell solver if a coil is switched on or off)
z_current = np.zeros(len(times)) # vertical position of the plasma (average jtor position)
dynamic_timings = np.zeros(len(times)) # solver runtime at each time step
dynamic_jtor_norm = np.zeros(len(times)) # norm change in jtor between time steps (used to trigger relinearisation)
threshold = 0.05 # relative jtor threshold (since last linearisation) above which relinearisation is triggered (here not used)
# extract any initial values from the initial equilibrium
dynamic_psi[:,:,0] = stepping.eq1.psi()
dynamic_limiter_flag[0] = stepping.eq1._profiles.flag_limiter
dynamic_psi_boundary[0] = stepping.eq1._profiles.psi_bndry
dynamic_xpts.append(stepping.eq1.xpt)
dynamic_opts.append(stepping.eq1.opt)
dynamic_currents[:,0] = stepping.vessel_currents_vec.copy()
dynamic_ip[0] = stepping.profiles1.Ip
dynamic_shape_targets[:,0] = plasma_descriptors(eq=stepping.eq1)
dynamic_triangularity[0] = eq.triangularity()
linear_growth_rate[0] = np.max(np.real(stepping.linearised_sol.all_timescales))
z_current[0] = stepping.eq1.Zcurrent()
Before running the following cell (which takes a few minutes), be sure to familiarise yourself with all of the steps.
During each time step:
The primary method in the
PCSclass,calculate_ctrl_voltages, is called. Using “measurements” from the current equilibrium (plasma current, coil currents, and shape parameters), it computes the voltages to be applied to the evolutive solver, as well as the coil resistances (to determine whether any coils are switched off). These voltages enact the plasma control.The plasma current density profile parameters are placed into a dictionary for use by the evolutive solver.
Note: they are constant here, but they can be made time-dependent at this stage.The evolutive solver is then called using these parameters (along with additional ones). At this stage, several choices can be made:
Specify a time-dependent plasma resistivity (usually required for resimulation of prior discharges).
Select, via
linear_only, either a linear or fully nonlinear solution of the circuit, plasma, and Grad–Shafranov equations.Set a relinearisation threshold when using linear mode (see example notebook 05c).
Choose not to solve the Grad–Shafranov equation fully (when using linear mode), but instead solve only for the shape parameters specified in
plasma_descriptors(see example notebooks 05b and 05c).
The required data are stored following successful completion of the time step.
# RUN FPDT SIMULATION
for i, t in enumerate(times[0:-1]):
print("-----")
print(f"t = {np.round(t,5)}s (step {i+1}/{n-1})")
# start timer
start_time = time.time()
# initialise any historical quantities for PCS PID controllers
if i == 0:
ip_hist_prev = 0.0
ip_err_prev = 0.0
T_err_prev = np.zeros(len(ctrl_targets))
T_hist_prev = np.zeros(len(ctrl_targets))
I_approved_prev = eq.tokamak.getCurrentsVec()[0:11] # must use starting currents here
V_approved_prev = np.zeros(len(ctrl_coils))
zipv_meas = 0.0
else:
ip_hist_prev=ip_hist[i-1].copy() # integral term from FB PID Plasma controller at prior step
ip_err_prev=ip_err[i-1].copy() # error term from FB PID Plasma controller at prior step
T_err_prev=T_err[:,i-1].copy() # proportional term from FB PID Shape controller at prior step
T_hist_prev=T_hist[:,i-1].copy() # integral term from FB PID Shape controller at prior step
I_approved_prev=I_approved[:,i-1].copy() # approved PF coil currents from prior step (ctrl_coils only)
V_approved_prev=V_approved[0:-1,i-1].copy() # approved PF coil voltages from prior step (ctrl_coils only)
zipv_meas = ((z_current[i]-z_current[i-1])/dt)*dynamic_ip[i] # rate of change of vertical plasma position
# # if you wish to calculate new VCs on the fly uncomment below
# vcg_inputs = VCG.get_inputs_from_eq(eq=stepping.eq1, profiles=stepping.profiles1)
# call the PCS class to attain PF coil voltages (on ctrl_coils and vertical coil)
V_approved[:,i], ip_hist[i], ip_err[i], T_err[:,i], T_hist[:,i], I_approved[:,i], coil_resists[:,i] = PCS.calculate_ctrl_voltages(
t=t, # current simulation time
dt=dt_PCS, # PCS time step
dt_simulator=dt, # simulator (solver) time step
ip_meas=dynamic_ip[i], # measured plasma current
ip_hist_prev=ip_hist_prev, # plasma controller integral term history
ip_err_prev=ip_err_prev, # plasma controller error term history (for derivative term)
T_meas=dynamic_shape_targets[:,i].copy(), # measured shape parameters
T_err_prev=T_err_prev, # shape controller proportional term history
T_hist_prev=T_hist_prev, # shape controller integral term history
I_approved_prev=I_approved_prev, # approved PF currents from prior timestep
I_meas=dynamic_currents[0:11,i].copy(), # measured PF coil currents
V_approved_prev=V_approved_prev, # approved PF voltages from prior timestep
zip_meas=z_current[i]*dynamic_ip[i], # measured vertical position x measured plasma current
zipv_meas=zipv_meas, # derivative of above
active_coil_resists=tokamak.coil_resist[0:12].copy(), # PF coil resistances (these are constant)
verbose=False, # print some output?
# vcg_inputs=vcg_inputs, # pass in eq and profiles at each time to re-calculate VCs after each time period (vc_update_rate)
)
# extract plasma current density profile parameters (these are constant but can be time-dependent if desired)
profile_params = {
"alpha": profiles.alpha[0:2],
"beta": profiles.beta[0:2],
}
# run freegsnke over the time step with the calculated voltages
stepping.nlstepper(
plasma_resistivity=1e-7, # assign plasma resistivity (chosen to be constant here)
active_voltage_vec=V_approved[:,i], # assign approved PF coil voltages from the PCS
profiles_parameters=profile_params, # assign profile parameters
custom_active_coil_resistances=coil_resists[:,i], # assign active coil resistances from PCS (tells solver if any coils are switched off)
linear_only=True, # linear or nonlinear solve?
target_relative_tol_currents=1e-2, # relative tolerance in the currents required for convergence
target_relative_tol_GS=1e-2, # relative tolerance in the plasma flux required for convergence
working_relative_tol_GS=(1e-2)/2, # tolerance used when solving GS equation, expressed in terms of the change in the plasma flux due to one timestep of evolution (must be smaller tolerance above)
max_solving_iterations=20, # stop after this many iterations
verbose=False, # print some output?
relinearise_threshold=threshold, # if the relative jtor norm change from last linearisation is above threshold, relinearise around current equilibrium
no_GS=False, # do not solve GS at each time?
)
# stop timer
end_time = time.time()
# extract and store relevant data
dynamic_psi[:,:,i+1] = stepping.eq1.psi()
dynamic_limiter_flag[i+1] = stepping.eq1._profiles.flag_limiter
dynamic_psi_boundary[i+1] = stepping.eq1._profiles.psi_bndry
dynamic_xpts.append(stepping.eq1.xpt)
dynamic_opts.append(stepping.eq1.opt)
dynamic_currents[:,i+1] = stepping.vessel_currents_vec.copy()
dynamic_ip[i+1] = stepping.currents_vec[-1] * stepping.plasma_norm_factor
dynamic_shape_targets[:,i+1] = plasma_descriptors(stepping.eq1)
dynamic_timings[i] = end_time - start_time
dynamic_triangularity[i] = stepping.eq1.triangularity()
z_current[i+1] = stepping.eq1.Zcurrent()
dynamic_jtor_norm[i+1] = stepping.relinearise_criteria
stepping.linearised_sol.calculate_linear_growth_rate()
linear_growth_rate[i+1] = np.max(np.real(stepping.linearised_sol.all_timescales))
# print some stuff to track solve
print(f" PLASMA QUANTITIES")
print(f" Ip = {np.round(dynamic_ip[i+1]/1000,1)} [kA]")
print(f" Shape parameters {ctrl_targets} = {np.round(dynamic_shape_targets[:,i+1],3)}")
print(f" Z position = {np.round(z_current[i+1],5)} [m]")
print(f" Linearised (deformable) growth rate = {np.round(1/linear_growth_rate[i+1],3)} [1/s]")
print(f" SOLVER STATUS")
print(f" Relinearisation status = {np.round(dynamic_jtor_norm[i+1],3)} (threshold = {threshold}) ")
# break if growth rate timescale exceeds timestep multiple
if dt > linear_growth_rate[i+1]:
print(" BREAKING LOOP - INSTABILITY TIMESCALE EXCEEDS TIMESTEP SIZE")
break
-----
t = 0.0s (step 1/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.277 1.347 0.567 -1.245]
Z position = -0.00021 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.0 (threshold = 0.05)
-----
t = 0.0005s (step 2/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.277 1.347 0.567 -1.245]
Z position = 0.00013 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.001 (threshold = 0.05)
-----
t = 0.001s (step 3/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.277 1.347 0.567 -1.245]
Z position = 0.00031 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.002 (threshold = 0.05)
-----
t = 0.0015s (step 4/499)
PLASMA QUANTITIES
Ip = 758.2 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.277 1.347 0.568 -1.245]
Z position = 0.00036 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.002 (threshold = 0.05)
-----
t = 0.002s (step 5/499)
PLASMA QUANTITIES
Ip = 758.1 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.277 1.347 0.568 -1.245]
Z position = 0.00037 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.002 (threshold = 0.05)
-----
t = 0.0025s (step 6/499)
PLASMA QUANTITIES
Ip = 757.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.278 1.347 0.568 -1.245]
Z position = 0.00036 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.003 (threshold = 0.05)
-----
t = 0.003s (step 7/499)
PLASMA QUANTITIES
Ip = 757.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.278 1.347 0.568 -1.245]
Z position = 0.00036 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.003 (threshold = 0.05)
-----
t = 0.0035s (step 8/499)
PLASMA QUANTITIES
Ip = 757.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.278 1.348 0.568 -1.245]
Z position = 0.00037 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.004 (threshold = 0.05)
-----
t = 0.004s (step 9/499)
PLASMA QUANTITIES
Ip = 757.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.278 1.348 0.569 -1.245]
Z position = 0.00039 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.004 (threshold = 0.05)
-----
t = 0.0045s (step 10/499)
PLASMA QUANTITIES
Ip = 757.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.279 1.348 0.569 -1.245]
Z position = 0.00042 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.005 (threshold = 0.05)
-----
t = 0.005s (step 11/499)
PLASMA QUANTITIES
Ip = 757.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.279 1.348 0.569 -1.245]
Z position = 0.00046 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.006 (threshold = 0.05)
-----
t = 0.0055s (step 12/499)
PLASMA QUANTITIES
Ip = 757.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.279 1.348 0.569 -1.245]
Z position = 0.0005 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.007 (threshold = 0.05)
-----
t = 0.006s (step 13/499)
PLASMA QUANTITIES
Ip = 757.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.279 1.349 0.569 -1.245]
Z position = 0.00052 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.008 (threshold = 0.05)
-----
t = 0.0065s (step 14/499)
PLASMA QUANTITIES
Ip = 757.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.279 1.349 0.57 -1.245]
Z position = 0.00055 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.009 (threshold = 0.05)
-----
t = 0.007s (step 15/499)
PLASMA QUANTITIES
Ip = 757.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.28 1.349 0.57 -1.245]
Z position = 0.00059 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.01 (threshold = 0.05)
-----
t = 0.0075s (step 16/499)
PLASMA QUANTITIES
Ip = 757.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.28 1.349 0.57 -1.245]
Z position = 0.00063 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.011 (threshold = 0.05)
-----
t = 0.008s (step 17/499)
PLASMA QUANTITIES
Ip = 757.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.28 1.35 0.57 -1.245]
Z position = 0.00068 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.012 (threshold = 0.05)
-----
t = 0.0085s (step 18/499)
PLASMA QUANTITIES
Ip = 757.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.28 1.35 0.57 -1.245]
Z position = 0.00073 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.009s (step 19/499)
PLASMA QUANTITIES
Ip = 757.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.28 1.351 0.571 -1.245]
Z position = 0.00077 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.0095s (step 20/499)
PLASMA QUANTITIES
Ip = 757.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.571 -1.245]
Z position = 0.00078 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.01s (step 21/499)
PLASMA QUANTITIES
Ip = 757.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.571 -1.245]
Z position = 0.00082 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.016 (threshold = 0.05)
-----
t = 0.0105s (step 22/499)
PLASMA QUANTITIES
Ip = 757.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.571 -1.245]
Z position = 0.00087 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.017 (threshold = 0.05)
-----
t = 0.011s (step 23/499)
PLASMA QUANTITIES
Ip = 757.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.571 -1.245]
Z position = 0.00091 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.019 (threshold = 0.05)
-----
t = 0.0115s (step 24/499)
PLASMA QUANTITIES
Ip = 757.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.571 -1.245]
Z position = 0.00095 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.02 (threshold = 0.05)
-----
t = 0.012s (step 25/499)
PLASMA QUANTITIES
Ip = 757.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.353 0.572 -1.245]
Z position = 0.00098 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.021 (threshold = 0.05)
-----
t = 0.0125s (step 26/499)
PLASMA QUANTITIES
Ip = 757.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.353 0.572 -1.245]
Z position = 0.00101 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.022 (threshold = 0.05)
-----
t = 0.013s (step 27/499)
PLASMA QUANTITIES
Ip = 757.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.353 0.572 -1.245]
Z position = 0.00107 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.023 (threshold = 0.05)
-----
t = 0.0135s (step 28/499)
PLASMA QUANTITIES
Ip = 757.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.572 -1.245]
Z position = 0.00112 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.025 (threshold = 0.05)
-----
t = 0.014s (step 29/499)
PLASMA QUANTITIES
Ip = 757.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.572 -1.245]
Z position = 0.00114 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.026 (threshold = 0.05)
-----
t = 0.0145s (step 30/499)
PLASMA QUANTITIES
Ip = 758.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.572 -1.245]
Z position = 0.00118 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.027 (threshold = 0.05)
-----
t = 0.015s (step 31/499)
PLASMA QUANTITIES
Ip = 758.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.572 -1.245]
Z position = 0.00123 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.027 (threshold = 0.05)
-----
t = 0.0155s (step 32/499)
PLASMA QUANTITIES
Ip = 758.1 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.355 0.572 -1.245]
Z position = 0.00128 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.028 (threshold = 0.05)
-----
t = 0.016s (step 33/499)
PLASMA QUANTITIES
Ip = 758.1 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.355 0.572 -1.245]
Z position = 0.00131 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.029 (threshold = 0.05)
-----
t = 0.0165s (step 34/499)
PLASMA QUANTITIES
Ip = 758.1 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.355 0.573 -1.245]
Z position = 0.00135 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.03 (threshold = 0.05)
-----
t = 0.017s (step 35/499)
PLASMA QUANTITIES
Ip = 758.2 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.355 0.573 -1.245]
Z position = 0.00139 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.031 (threshold = 0.05)
-----
t = 0.0175s (step 36/499)
PLASMA QUANTITIES
Ip = 758.2 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.355 0.573 -1.245]
Z position = 0.00144 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.032 (threshold = 0.05)
-----
t = 0.018s (step 37/499)
PLASMA QUANTITIES
Ip = 758.3 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.356 0.573 -1.246]
Z position = 0.00149 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.033 (threshold = 0.05)
-----
t = 0.0185s (step 38/499)
PLASMA QUANTITIES
Ip = 758.3 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.356 0.573 -1.246]
Z position = 0.00154 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.033 (threshold = 0.05)
-----
t = 0.019s (step 39/499)
PLASMA QUANTITIES
Ip = 758.3 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.356 0.573 -1.246]
Z position = 0.00158 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.034 (threshold = 0.05)
-----
t = 0.0195s (step 40/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.356 0.573 -1.246]
Z position = 0.00161 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.035 (threshold = 0.05)
-----
t = 0.02s (step 41/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.356 0.573 -1.246]
Z position = 0.00165 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.035 (threshold = 0.05)
-----
t = 0.0205s (step 42/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.356 0.573 -1.246]
Z position = 0.00168 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.036 (threshold = 0.05)
-----
t = 0.021s (step 43/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.356 0.573 -1.246]
Z position = 0.00172 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.036 (threshold = 0.05)
-----
t = 0.0215s (step 44/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.356 0.573 -1.246]
Z position = 0.00176 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.036 (threshold = 0.05)
-----
t = 0.022s (step 45/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.284 1.356 0.573 -1.246]
Z position = 0.0018 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.037 (threshold = 0.05)
-----
t = 0.0225s (step 46/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.284 1.356 0.573 -1.246]
Z position = 0.00184 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.037 (threshold = 0.05)
-----
t = 0.023s (step 47/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.284 1.356 0.573 -1.246]
Z position = 0.00188 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.037 (threshold = 0.05)
-----
t = 0.0235s (step 48/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.284 1.356 0.573 -1.247]
Z position = 0.00192 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.024s (step 49/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.284 1.356 0.573 -1.247]
Z position = 0.00197 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.0245s (step 50/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.284 1.356 0.573 -1.247]
Z position = 0.00201 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.025s (step 51/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.284 1.356 0.573 -1.247]
Z position = 0.00205 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.0255s (step 52/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.284 1.356 0.573 -1.247]
Z position = 0.00208 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.026s (step 53/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.284 1.356 0.573 -1.247]
Z position = 0.00211 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.0265s (step 54/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.284 1.356 0.573 -1.247]
Z position = 0.00215 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.027s (step 55/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.284 1.356 0.573 -1.247]
Z position = 0.00218 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.0275s (step 56/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.284 1.356 0.572 -1.247]
Z position = 0.00222 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.028s (step 57/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.356 0.572 -1.247]
Z position = 0.00226 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.0285s (step 58/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.356 0.572 -1.247]
Z position = 0.0023 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.029s (step 59/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.356 0.572 -1.248]
Z position = 0.00235 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.0295s (step 60/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.356 0.571 -1.248]
Z position = 0.00239 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.03s (step 61/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.356 0.571 -1.248]
Z position = 0.00244 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.0305s (step 62/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.356 0.571 -1.248]
Z position = 0.00248 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.039 (threshold = 0.05)
-----
t = 0.031s (step 63/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.356 0.57 -1.248]
Z position = 0.00253 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.039 (threshold = 0.05)
-----
t = 0.0315s (step 64/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.356 0.57 -1.248]
Z position = 0.00257 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.039 (threshold = 0.05)
-----
t = 0.032s (step 65/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.356 0.57 -1.248]
Z position = 0.00261 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.039 (threshold = 0.05)
-----
t = 0.0325s (step 66/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.356 0.569 -1.248]
Z position = 0.00265 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.039 (threshold = 0.05)
-----
t = 0.033s (step 67/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.356 0.569 -1.248]
Z position = 0.00269 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.039 (threshold = 0.05)
-----
t = 0.0335s (step 68/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.356 0.568 -1.248]
Z position = 0.00273 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.039 (threshold = 0.05)
-----
t = 0.034s (step 69/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.355 0.568 -1.248]
Z position = 0.00277 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.039 (threshold = 0.05)
-----
t = 0.0345s (step 70/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.355 0.567 -1.248]
Z position = 0.00281 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.039 (threshold = 0.05)
-----
t = 0.035s (step 71/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.355 0.567 -1.248]
Z position = 0.00285 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.0355s (step 72/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.355 0.566 -1.248]
Z position = 0.00289 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.036s (step 73/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.355 0.565 -1.248]
Z position = 0.00293 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.0365s (step 74/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.355 0.565 -1.248]
Z position = 0.00297 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.037s (step 75/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.355 0.564 -1.249]
Z position = 0.00301 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.0375s (step 76/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.355 0.564 -1.249]
Z position = 0.00305 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.038s (step 77/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.355 0.563 -1.249]
Z position = 0.0031 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.0385s (step 78/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.355 0.563 -1.249]
Z position = 0.00314 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.039s (step 79/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.355 0.562 -1.249]
Z position = 0.00318 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.0395s (step 80/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.355 0.561 -1.249]
Z position = 0.00322 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.04s (step 81/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.561 -1.249]
Z position = 0.00326 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.0405s (step 82/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.561 -1.249]
Z position = 0.0033 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.041s (step 83/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.56 -1.249]
Z position = 0.00343 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.0415s (step 84/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.559 -1.249]
Z position = 0.00342 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.042s (step 85/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.559 -1.249]
Z position = 0.00328 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.0425s (step 86/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.558 -1.249]
Z position = 0.00337 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.043s (step 87/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.557 -1.249]
Z position = 0.00347 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.0435s (step 88/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.556 -1.249]
Z position = 0.00355 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.038 (threshold = 0.05)
-----
t = 0.044s (step 89/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.556 -1.249]
Z position = 0.00361 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.039 (threshold = 0.05)
-----
t = 0.0445s (step 90/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.555 -1.249]
Z position = 0.00364 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.039 (threshold = 0.05)
-----
t = 0.045s (step 91/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.554 -1.249]
Z position = 0.00367 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.039 (threshold = 0.05)
-----
t = 0.0455s (step 92/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.553 -1.25 ]
Z position = 0.0037 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.039 (threshold = 0.05)
-----
t = 0.046s (step 93/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.553 -1.25 ]
Z position = 0.00374 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.039 (threshold = 0.05)
-----
t = 0.0465s (step 94/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.552 -1.25 ]
Z position = 0.00377 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.039 (threshold = 0.05)
-----
t = 0.047s (step 95/499)
/opt/hostedtoolcache/Python/3.10.22/x64/lib/python3.10/site-packages/freegsnke/critical.py:1094: UserWarning: Theta grid too close to X-point, shifting by half-step
warnings.warn("Theta grid too close to X-point, shifting by half-step")
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.551 -1.25 ]
Z position = 0.00381 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.04 (threshold = 0.05)
-----
t = 0.0475s (step 96/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.551 -1.25 ]
Z position = 0.00385 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.04 (threshold = 0.05)
-----
t = 0.048s (step 97/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.55 -1.25 ]
Z position = 0.00389 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.04 (threshold = 0.05)
-----
t = 0.0485s (step 98/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.549 -1.25 ]
Z position = 0.00394 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.04 (threshold = 0.05)
-----
t = 0.049s (step 99/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.548 -1.25 ]
Z position = 0.00398 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.041 (threshold = 0.05)
-----
t = 0.0495s (step 100/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.548 -1.25 ]
Z position = 0.00402 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.041 (threshold = 0.05)
-----
t = 0.05s (step 101/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.547 -1.25 ]
Z position = 0.00406 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.041 (threshold = 0.05)
-----
t = 0.0505s (step 102/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.546 -1.25 ]
Z position = 0.0041 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.041 (threshold = 0.05)
-----
t = 0.051s (step 103/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.546 -1.25 ]
Z position = 0.00414 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.041 (threshold = 0.05)
-----
t = 0.0515s (step 104/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.545 -1.25 ]
Z position = 0.00418 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.042 (threshold = 0.05)
-----
t = 0.052s (step 105/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.544 -1.249]
Z position = 0.00422 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.042 (threshold = 0.05)
-----
t = 0.0525s (step 106/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.543 -1.249]
Z position = 0.00426 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.042 (threshold = 0.05)
-----
t = 0.053s (step 107/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.543 -1.249]
Z position = 0.0043 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.043 (threshold = 0.05)
-----
t = 0.0535s (step 108/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.542 -1.249]
Z position = 0.00434 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.043 (threshold = 0.05)
-----
t = 0.054s (step 109/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.541 -1.249]
Z position = 0.00438 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.043 (threshold = 0.05)
-----
t = 0.0545s (step 110/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.54 -1.249]
Z position = 0.00442 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.043 (threshold = 0.05)
-----
t = 0.055s (step 111/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.54 -1.249]
Z position = 0.00447 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.044 (threshold = 0.05)
-----
t = 0.0555s (step 112/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.539 -1.249]
Z position = 0.00452 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.044 (threshold = 0.05)
-----
t = 0.056s (step 113/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.538 -1.249]
Z position = 0.00455 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.044 (threshold = 0.05)
-----
t = 0.0565s (step 114/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.538 -1.249]
Z position = 0.00459 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.045 (threshold = 0.05)
-----
t = 0.057s (step 115/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.537 -1.249]
Z position = 0.00463 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.045 (threshold = 0.05)
-----
t = 0.0575s (step 116/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.536 -1.249]
Z position = 0.00466 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.045 (threshold = 0.05)
-----
t = 0.058s (step 117/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.535 -1.249]
Z position = 0.0047 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.045 (threshold = 0.05)
-----
t = 0.0585s (step 118/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.535 -1.249]
Z position = 0.00474 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.046 (threshold = 0.05)
-----
t = 0.059s (step 119/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.534 -1.249]
Z position = 0.00478 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.046 (threshold = 0.05)
-----
t = 0.0595s (step 120/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.533 -1.249]
Z position = 0.00482 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.046 (threshold = 0.05)
-----
t = 0.06s (step 121/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.533 -1.249]
Z position = 0.00486 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.047 (threshold = 0.05)
-----
t = 0.0605s (step 122/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.532 -1.249]
Z position = 0.0049 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.047 (threshold = 0.05)
-----
t = 0.061s (step 123/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.531 -1.249]
Z position = 0.00494 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.047 (threshold = 0.05)
-----
t = 0.0615s (step 124/499)
PLASMA QUANTITIES
Ip = 759.1 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.531 -1.25 ]
Z position = 0.00498 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.048 (threshold = 0.05)
-----
t = 0.062s (step 125/499)
PLASMA QUANTITIES
Ip = 759.1 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.53 -1.25 ]
Z position = 0.00508 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.048 (threshold = 0.05)
-----
t = 0.0625s (step 126/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.53 -1.25 ]
Z position = 0.00508 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.048 (threshold = 0.05)
-----
t = 0.063s (step 127/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.529 -1.25 ]
Z position = 0.00508 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.048 (threshold = 0.05)
-----
t = 0.0635s (step 128/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.528 -1.25 ]
Z position = 0.00497 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.048 (threshold = 0.05)
-----
t = 0.064s (step 129/499)
PLASMA QUANTITIES
Ip = 759.1 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.527 -1.25 ]
Z position = 0.0051 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.048 (threshold = 0.05)
-----
t = 0.0645s (step 130/499)
PLASMA QUANTITIES
Ip = 759.1 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.526 -1.25 ]
Z position = 0.00522 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.049 (threshold = 0.05)
-----
t = 0.065s (step 131/499)
PLASMA QUANTITIES
Ip = 759.1 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.526 -1.25 ]
Z position = 0.0053 [m]
Linearised (deformable) growth rate = 24.577 [1/s]
SOLVER STATUS
Relinearisation status = 0.05 (threshold = 0.05)
-----
t = 0.0655s (step 132/499)
Re-linearising around current equilibrium!
Relative relinearisation criteria change = 5.038% (threshold = 5.0%)
Building the 1513 x 43 Jacobian (dIy/dI) of plasma current density (inside the LCFS) with respect to all metal currents and the total plasma current.
Building the 1513 x 4 Jacobian (dIy/dtheta) of plasma current density (inside the LCFS) with respect to all plasma current density profile parameters within Jtor.
Built the 4 x 43 Jacobian (ds/dI) of plasma descriptors with respect to all metal currents and the total plasma current.
Built the 4 x 4 Jacobian (ds/dtheta) of plasma descriptors with respect to all plasma current density profile parameters within Jtor.
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.354 0.525 -1.25 ]
Z position = 0.00533 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.05 (threshold = 0.05)
-----
t = 0.066s (step 133/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.524 -1.249]
Z position = 0.00537 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.001 (threshold = 0.05)
-----
t = 0.0665s (step 134/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.523 -1.249]
Z position = 0.0054 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.001 (threshold = 0.05)
-----
t = 0.067s (step 135/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.522 -1.249]
Z position = 0.00544 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.002 (threshold = 0.05)
-----
t = 0.0675s (step 136/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.353 0.521 -1.249]
Z position = 0.00548 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.003 (threshold = 0.05)
-----
t = 0.068s (step 137/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.353 0.521 -1.249]
Z position = 0.00552 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.003 (threshold = 0.05)
-----
t = 0.0685s (step 138/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.52 -1.249]
Z position = 0.00556 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.004 (threshold = 0.05)
-----
t = 0.069s (step 139/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.519 -1.249]
Z position = 0.0056 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.004 (threshold = 0.05)
-----
t = 0.0695s (step 140/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.518 -1.249]
Z position = 0.00563 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.004 (threshold = 0.05)
-----
t = 0.07s (step 141/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.518 -1.249]
Z position = 0.00567 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.005 (threshold = 0.05)
-----
t = 0.0705s (step 142/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.517 -1.249]
Z position = 0.00571 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.005 (threshold = 0.05)
-----
t = 0.071s (step 143/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.516 -1.248]
Z position = 0.00575 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.005 (threshold = 0.05)
-----
t = 0.0715s (step 144/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.515 -1.248]
Z position = 0.00579 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.006 (threshold = 0.05)
-----
t = 0.072s (step 145/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.515 -1.248]
Z position = 0.00583 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.006 (threshold = 0.05)
-----
t = 0.0725s (step 146/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.514 -1.248]
Z position = 0.00587 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.006 (threshold = 0.05)
-----
t = 0.073s (step 147/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.513 -1.248]
Z position = 0.00591 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.006 (threshold = 0.05)
-----
t = 0.0735s (step 148/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.513 -1.248]
Z position = 0.00595 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.006 (threshold = 0.05)
-----
t = 0.074s (step 149/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.512 -1.248]
Z position = 0.00598 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.006 (threshold = 0.05)
-----
t = 0.0745s (step 150/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.511 -1.248]
Z position = 0.00602 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.006 (threshold = 0.05)
-----
t = 0.075s (step 151/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.511 -1.248]
Z position = 0.00605 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.006 (threshold = 0.05)
-----
t = 0.0755s (step 152/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.51 -1.247]
Z position = 0.0061 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.007 (threshold = 0.05)
-----
t = 0.076s (step 153/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.509 -1.247]
Z position = 0.00614 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.007 (threshold = 0.05)
-----
t = 0.0765s (step 154/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.509 -1.247]
Z position = 0.00618 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.007 (threshold = 0.05)
-----
t = 0.077s (step 155/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.508 -1.247]
Z position = 0.00622 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.007 (threshold = 0.05)
-----
t = 0.0775s (step 156/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.508 -1.247]
Z position = 0.00626 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.007 (threshold = 0.05)
-----
t = 0.078s (step 157/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.507 -1.247]
Z position = 0.0063 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.007 (threshold = 0.05)
-----
t = 0.0785s (step 158/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.507 -1.247]
Z position = 0.00634 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.007 (threshold = 0.05)
-----
t = 0.079s (step 159/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.506 -1.247]
Z position = 0.00638 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.008 (threshold = 0.05)
-----
t = 0.0795s (step 160/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.506 -1.247]
Z position = 0.00641 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.008 (threshold = 0.05)
-----
t = 0.08s (step 161/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.506 -1.247]
Z position = 0.00645 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.008 (threshold = 0.05)
-----
t = 0.0805s (step 162/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.505 -1.247]
Z position = 0.00649 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.008 (threshold = 0.05)
-----
t = 0.081s (step 163/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.505 -1.247]
Z position = 0.00653 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.008 (threshold = 0.05)
-----
t = 0.0815s (step 164/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.505 -1.247]
Z position = 0.00657 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.008 (threshold = 0.05)
-----
t = 0.082s (step 165/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.505 -1.247]
Z position = 0.00661 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.008 (threshold = 0.05)
-----
t = 0.0825s (step 166/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.504 -1.247]
Z position = 0.00665 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.009 (threshold = 0.05)
-----
t = 0.083s (step 167/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.504 -1.247]
Z position = 0.00669 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.009 (threshold = 0.05)
-----
t = 0.0835s (step 168/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.504 -1.247]
Z position = 0.00674 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.009 (threshold = 0.05)
-----
t = 0.084s (step 169/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.504 -1.246]
Z position = 0.00678 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.009 (threshold = 0.05)
-----
t = 0.0845s (step 170/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.503 -1.246]
Z position = 0.00682 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.009 (threshold = 0.05)
-----
t = 0.085s (step 171/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.503 -1.246]
Z position = 0.00686 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.009 (threshold = 0.05)
-----
t = 0.0855s (step 172/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.283 1.352 0.503 -1.246]
Z position = 0.0069 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.009 (threshold = 0.05)
-----
t = 0.086s (step 173/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.503 -1.246]
Z position = 0.00694 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.009 (threshold = 0.05)
-----
t = 0.0865s (step 174/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.503 -1.246]
Z position = 0.00698 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.01 (threshold = 0.05)
-----
t = 0.087s (step 175/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.503 -1.246]
Z position = 0.00702 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.01 (threshold = 0.05)
-----
t = 0.0875s (step 176/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.00706 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.01 (threshold = 0.05)
-----
t = 0.088s (step 177/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.0071 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.01 (threshold = 0.05)
-----
t = 0.0885s (step 178/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.00714 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.01 (threshold = 0.05)
-----
t = 0.089s (step 179/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.00718 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.01 (threshold = 0.05)
-----
t = 0.0895s (step 180/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.00722 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.01 (threshold = 0.05)
-----
t = 0.09s (step 181/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.00726 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.011 (threshold = 0.05)
-----
t = 0.0905s (step 182/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.0073 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.011 (threshold = 0.05)
-----
t = 0.091s (step 183/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.00734 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.011 (threshold = 0.05)
-----
t = 0.0915s (step 184/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.00738 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.011 (threshold = 0.05)
-----
t = 0.092s (step 185/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.00742 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.011 (threshold = 0.05)
-----
t = 0.0925s (step 186/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.00746 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.011 (threshold = 0.05)
-----
t = 0.093s (step 187/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.0075 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.011 (threshold = 0.05)
-----
t = 0.0935s (step 188/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.00754 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.011 (threshold = 0.05)
-----
t = 0.094s (step 189/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.00758 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.011 (threshold = 0.05)
-----
t = 0.0945s (step 190/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.502 -1.247]
Z position = 0.00762 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.012 (threshold = 0.05)
-----
t = 0.095s (step 191/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.502 -1.247]
Z position = 0.00766 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.012 (threshold = 0.05)
-----
t = 0.0955s (step 192/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.502 -1.247]
Z position = 0.0077 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.012 (threshold = 0.05)
-----
t = 0.096s (step 193/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.502 -1.247]
Z position = 0.00774 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.012 (threshold = 0.05)
-----
t = 0.0965s (step 194/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.502 -1.247]
Z position = 0.00778 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.012 (threshold = 0.05)
-----
t = 0.097s (step 195/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.502 -1.247]
Z position = 0.00782 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.012 (threshold = 0.05)
-----
t = 0.0975s (step 196/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.247]
Z position = 0.00786 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.012 (threshold = 0.05)
-----
t = 0.098s (step 197/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.247]
Z position = 0.0079 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.012 (threshold = 0.05)
-----
t = 0.0985s (step 198/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.247]
Z position = 0.00794 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.012 (threshold = 0.05)
-----
t = 0.099s (step 199/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.247]
Z position = 0.00798 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.012 (threshold = 0.05)
-----
t = 0.0995s (step 200/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.247]
Z position = 0.00802 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.012 (threshold = 0.05)
-----
t = 0.1s (step 201/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.247]
Z position = 0.00806 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.012 (threshold = 0.05)
-----
t = 0.1005s (step 202/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.247]
Z position = 0.0081 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.101s (step 203/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00814 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1015s (step 204/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00818 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.102s (step 205/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00822 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1025s (step 206/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00826 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.103s (step 207/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.0083 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1035s (step 208/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00834 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.104s (step 209/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00838 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1045s (step 210/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00842 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.105s (step 211/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00846 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1055s (step 212/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.0085 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.106s (step 213/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00854 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1065s (step 214/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00858 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.107s (step 215/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00862 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1075s (step 216/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00866 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.108s (step 217/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.0087 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1085s (step 218/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00874 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.109s (step 219/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00878 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1095s (step 220/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00882 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.11s (step 221/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00886 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1105s (step 222/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.0089 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.111s (step 223/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00894 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1115s (step 224/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00898 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.112s (step 225/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00902 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1125s (step 226/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00906 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.113s (step 227/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.0091 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1135s (step 228/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00914 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.114s (step 229/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00918 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1145s (step 230/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00922 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.115s (step 231/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00926 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1155s (step 232/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.0093 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.116s (step 233/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00934 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1165s (step 234/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00938 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.117s (step 235/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00942 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1175s (step 236/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00946 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.118s (step 237/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.0095 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1185s (step 238/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.501 -1.246]
Z position = 0.00954 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.119s (step 239/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.246]
Z position = 0.00958 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1195s (step 240/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.246]
Z position = 0.00962 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.12s (step 241/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.246]
Z position = 0.00966 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1205s (step 242/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.246]
Z position = 0.0097 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.121s (step 243/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.246]
Z position = 0.00974 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1215s (step 244/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.246]
Z position = 0.00978 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.122s (step 245/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.246]
Z position = 0.00982 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1225s (step 246/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.245]
Z position = 0.00986 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.123s (step 247/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.245]
Z position = 0.0099 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1235s (step 248/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.245]
Z position = 0.00994 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.124s (step 249/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.245]
Z position = 0.00998 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1245s (step 250/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.245]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.125s (step 251/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.245]
Z position = 0.01005 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1255s (step 252/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.245]
Z position = 0.01006 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.126s (step 253/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.245]
Z position = 0.01006 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1265s (step 254/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.245]
Z position = 0.01004 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.127s (step 255/499)
PLASMA QUANTITIES
Ip = 758.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.245]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1275s (step 256/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.245]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.128s (step 257/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1285s (step 258/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.129s (step 259/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1295s (step 260/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.13s (step 261/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1305s (step 262/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.131s (step 263/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.1315s (step 264/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.132s (step 265/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1325s (step 266/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.133s (step 267/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1335s (step 268/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.134s (step 269/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1345s (step 270/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.135s (step 271/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.501 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1355s (step 272/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.136s (step 273/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1365s (step 274/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.137s (step 275/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1375s (step 276/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.138s (step 277/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1385s (step 278/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.139s (step 279/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1395s (step 280/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.14s (step 281/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1405s (step 282/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.141s (step 283/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1415s (step 284/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.142s (step 285/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1425s (step 286/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.143s (step 287/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1435s (step 288/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.144s (step 289/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1445s (step 290/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.145s (step 291/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1455s (step 292/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.146s (step 293/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1465s (step 294/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.147s (step 295/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1475s (step 296/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.148s (step 297/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1485s (step 298/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.149s (step 299/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1495s (step 300/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.15s (step 301/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1505s (step 302/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.151s (step 303/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1515s (step 304/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.152s (step 305/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1525s (step 306/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.153s (step 307/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1535s (step 308/499)
PLASMA QUANTITIES
Ip = 758.5 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.154s (step 309/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1545s (step 310/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.155s (step 311/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1555s (step 312/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.156s (step 313/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.246]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1565s (step 314/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.157s (step 315/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1575s (step 316/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.158s (step 317/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1585s (step 318/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.159s (step 319/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1595s (step 320/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.16s (step 321/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1605s (step 322/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.161s (step 323/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1615s (step 324/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.162s (step 325/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1625s (step 326/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.163s (step 327/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1635s (step 328/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.164s (step 329/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1645s (step 330/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.165s (step 331/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1655s (step 332/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.166s (step 333/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1665s (step 334/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.167s (step 335/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1675s (step 336/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.168s (step 337/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1685s (step 338/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.169s (step 339/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1695s (step 340/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.17s (step 341/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1705s (step 342/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.171s (step 343/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1715s (step 344/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.172s (step 345/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1725s (step 346/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.173s (step 347/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1735s (step 348/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.174s (step 349/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1745s (step 350/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.175s (step 351/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1755s (step 352/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.176s (step 353/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1765s (step 354/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.177s (step 355/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1775s (step 356/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.178s (step 357/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1785s (step 358/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.179s (step 359/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.502 -1.247]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1795s (step 360/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.503 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.18s (step 361/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.352 0.503 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1805s (step 362/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.503 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.181s (step 363/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.504 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1815s (step 364/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.504 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.182s (step 365/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.505 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1825s (step 366/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.351 0.505 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.183s (step 367/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.506 -1.247]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1835s (step 368/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.506 -1.248]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.184s (step 369/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.507 -1.248]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1845s (step 370/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.507 -1.248]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.185s (step 371/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.508 -1.248]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1855s (step 372/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.508 -1.248]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.186s (step 373/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.509 -1.248]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1865s (step 374/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.509 -1.248]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.187s (step 375/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.51 -1.248]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1875s (step 376/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.511 -1.248]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.188s (step 377/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.511 -1.248]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1885s (step 378/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.512 -1.248]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.189s (step 379/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.512 -1.249]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1895s (step 380/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.513 -1.249]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.19s (step 381/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.514 -1.249]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1905s (step 382/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.514 -1.249]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.191s (step 383/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.515 -1.249]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1915s (step 384/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.516 -1.249]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.192s (step 385/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.516 -1.249]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1925s (step 386/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.517 -1.249]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.193s (step 387/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.518 -1.249]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1935s (step 388/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.518 -1.249]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.194s (step 389/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.519 -1.249]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.1945s (step 390/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.52 -1.249]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.195s (step 391/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.52 -1.25 ]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1955s (step 392/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.521 -1.25 ]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.196s (step 393/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.522 -1.25 ]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1965s (step 394/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.522 -1.25 ]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.197s (step 395/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.523 -1.25 ]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1975s (step 396/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.524 -1.25 ]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.198s (step 397/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.524 -1.25 ]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1985s (step 398/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.525 -1.25 ]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.199s (step 399/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.526 -1.25 ]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.1995s (step 400/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.527 -1.25 ]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.2s (step 401/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.527 -1.25 ]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.2005s (step 402/499)
PLASMA QUANTITIES
Ip = 758.6 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.528 -1.25 ]
Z position = 0.01018 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.201s (step 403/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.529 -1.25 ]
Z position = 0.01009 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.2015s (step 404/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.529 -1.25 ]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.202s (step 405/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.53 -1.251]
Z position = 0.00998 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.2025s (step 406/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.531 -1.251]
Z position = 0.00993 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.203s (step 407/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.532 -1.251]
Z position = 0.00995 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.2035s (step 408/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.531 -1.25 ]
Z position = 0.00998 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.204s (step 409/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.532 -1.25 ]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.2045s (step 410/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.533 -1.25 ]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.205s (step 411/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.533 -1.25 ]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.2055s (step 412/499)
PLASMA QUANTITIES
Ip = 758.7 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.534 -1.25 ]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.206s (step 413/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.535 -1.25 ]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.2065s (step 414/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.535 -1.251]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.207s (step 415/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.536 -1.251]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.2075s (step 416/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.537 -1.251]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.208s (step 417/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.538 -1.251]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.2085s (step 418/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.538 -1.251]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.209s (step 419/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.539 -1.251]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.2095s (step 420/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.54 -1.251]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.21s (step 421/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.54 -1.251]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.2105s (step 422/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.541 -1.251]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.211s (step 423/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.542 -1.251]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.2115s (step 424/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.543 -1.252]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.212s (step 425/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.543 -1.252]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.2125s (step 426/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.353 0.544 -1.252]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.213s (step 427/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.353 0.545 -1.252]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.2135s (step 428/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.353 0.545 -1.252]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.214s (step 429/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.353 0.546 -1.252]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.2145s (step 430/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.353 0.547 -1.252]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.215s (step 431/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.353 0.548 -1.252]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.013 (threshold = 0.05)
-----
t = 0.2155s (step 432/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.353 0.548 -1.252]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.216s (step 433/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.353 0.549 -1.252]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.2165s (step 434/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.353 0.55 -1.253]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.217s (step 435/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.353 0.55 -1.253]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.2175s (step 436/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.353 0.551 -1.253]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.218s (step 437/499)
PLASMA QUANTITIES
Ip = 759.1 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.353 0.552 -1.253]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.2185s (step 438/499)
PLASMA QUANTITIES
Ip = 759.1 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.353 0.553 -1.253]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.014 (threshold = 0.05)
-----
t = 0.219s (step 439/499)
PLASMA QUANTITIES
Ip = 759.1 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.353 0.553 -1.253]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.2195s (step 440/499)
PLASMA QUANTITIES
Ip = 759.1 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.353 0.554 -1.253]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.22s (step 441/499)
PLASMA QUANTITIES
Ip = 759.1 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.555 -1.253]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.2205s (step 442/499)
PLASMA QUANTITIES
Ip = 759.1 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.556 -1.253]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.221s (step 443/499)
PLASMA QUANTITIES
Ip = 759.2 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.556 -1.253]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.015 (threshold = 0.05)
-----
t = 0.2215s (step 444/499)
PLASMA QUANTITIES
Ip = 759.2 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.557 -1.254]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.016 (threshold = 0.05)
-----
t = 0.222s (step 445/499)
PLASMA QUANTITIES
Ip = 759.2 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.558 -1.254]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.016 (threshold = 0.05)
-----
t = 0.2225s (step 446/499)
PLASMA QUANTITIES
Ip = 759.2 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.558 -1.254]
Z position = 0.01022 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.016 (threshold = 0.05)
-----
t = 0.223s (step 447/499)
PLASMA QUANTITIES
Ip = 759.3 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.559 -1.254]
Z position = 0.01011 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.017 (threshold = 0.05)
-----
t = 0.2235s (step 448/499)
PLASMA QUANTITIES
Ip = 759.3 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.56 -1.254]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.017 (threshold = 0.05)
-----
t = 0.224s (step 449/499)
PLASMA QUANTITIES
Ip = 759.3 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.56 -1.254]
Z position = 0.00992 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.017 (threshold = 0.05)
-----
t = 0.2245s (step 450/499)
PLASMA QUANTITIES
Ip = 759.3 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.561 -1.254]
Z position = 0.00994 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.017 (threshold = 0.05)
-----
t = 0.225s (step 451/499)
PLASMA QUANTITIES
Ip = 759.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.561 -1.254]
Z position = 0.00997 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.018 (threshold = 0.05)
-----
t = 0.2255s (step 452/499)
PLASMA QUANTITIES
Ip = 759.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.562 -1.254]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.018 (threshold = 0.05)
-----
t = 0.226s (step 453/499)
PLASMA QUANTITIES
Ip = 759.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.562 -1.254]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.018 (threshold = 0.05)
-----
t = 0.2265s (step 454/499)
PLASMA QUANTITIES
Ip = 759.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.563 -1.254]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.018 (threshold = 0.05)
-----
t = 0.227s (step 455/499)
PLASMA QUANTITIES
Ip = 759.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.563 -1.254]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.018 (threshold = 0.05)
-----
t = 0.2275s (step 456/499)
PLASMA QUANTITIES
Ip = 759.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.564 -1.254]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.019 (threshold = 0.05)
-----
t = 0.228s (step 457/499)
PLASMA QUANTITIES
Ip = 759.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.564 -1.254]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.019 (threshold = 0.05)
-----
t = 0.2285s (step 458/499)
PLASMA QUANTITIES
Ip = 759.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.565 -1.254]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.019 (threshold = 0.05)
-----
t = 0.229s (step 459/499)
PLASMA QUANTITIES
Ip = 759.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.565 -1.254]
Z position = 0.01001 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.019 (threshold = 0.05)
-----
t = 0.2295s (step 460/499)
PLASMA QUANTITIES
Ip = 759.4 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.566 -1.255]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.019 (threshold = 0.05)
-----
t = 0.23s (step 461/499)
PLASMA QUANTITIES
Ip = 759.3 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.566 -1.255]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.019 (threshold = 0.05)
-----
t = 0.2305s (step 462/499)
PLASMA QUANTITIES
Ip = 759.3 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.566 -1.255]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.02 (threshold = 0.05)
-----
t = 0.231s (step 463/499)
PLASMA QUANTITIES
Ip = 759.3 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.567 -1.255]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.02 (threshold = 0.05)
-----
t = 0.2315s (step 464/499)
PLASMA QUANTITIES
Ip = 759.3 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.567 -1.255]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.02 (threshold = 0.05)
-----
t = 0.232s (step 465/499)
PLASMA QUANTITIES
Ip = 759.3 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.354 0.567 -1.255]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.02 (threshold = 0.05)
-----
t = 0.2325s (step 466/499)
PLASMA QUANTITIES
Ip = 759.3 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.353 0.568 -1.255]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.02 (threshold = 0.05)
-----
t = 0.233s (step 467/499)
PLASMA QUANTITIES
Ip = 759.3 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.353 0.568 -1.255]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.02 (threshold = 0.05)
-----
t = 0.2335s (step 468/499)
PLASMA QUANTITIES
Ip = 759.2 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.282 1.353 0.568 -1.255]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.02 (threshold = 0.05)
-----
t = 0.234s (step 469/499)
PLASMA QUANTITIES
Ip = 759.2 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.353 0.568 -1.255]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.02 (threshold = 0.05)
-----
t = 0.2345s (step 470/499)
PLASMA QUANTITIES
Ip = 759.2 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.353 0.568 -1.255]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.02 (threshold = 0.05)
-----
t = 0.235s (step 471/499)
PLASMA QUANTITIES
Ip = 759.2 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.353 0.569 -1.255]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.02 (threshold = 0.05)
-----
t = 0.2355s (step 472/499)
PLASMA QUANTITIES
Ip = 759.2 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.353 0.569 -1.255]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.021 (threshold = 0.05)
-----
t = 0.236s (step 473/499)
PLASMA QUANTITIES
Ip = 759.2 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.569 -1.255]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.021 (threshold = 0.05)
-----
t = 0.2365s (step 474/499)
PLASMA QUANTITIES
Ip = 759.1 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.569 -1.255]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.021 (threshold = 0.05)
-----
t = 0.237s (step 475/499)
PLASMA QUANTITIES
Ip = 759.1 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.569 -1.255]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.021 (threshold = 0.05)
-----
t = 0.2375s (step 476/499)
PLASMA QUANTITIES
Ip = 759.1 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.569 -1.255]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.021 (threshold = 0.05)
-----
t = 0.238s (step 477/499)
PLASMA QUANTITIES
Ip = 759.1 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.569 -1.255]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.021 (threshold = 0.05)
-----
t = 0.2385s (step 478/499)
PLASMA QUANTITIES
Ip = 759.1 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.352 0.57 -1.255]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.021 (threshold = 0.05)
-----
t = 0.239s (step 479/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.57 -1.255]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.021 (threshold = 0.05)
-----
t = 0.2395s (step 480/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.57 -1.255]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.022 (threshold = 0.05)
-----
t = 0.24s (step 481/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.57 -1.255]
Z position = 0.01004 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.022 (threshold = 0.05)
-----
t = 0.2405s (step 482/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.57 -1.255]
Z position = 0.01005 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.022 (threshold = 0.05)
-----
t = 0.241s (step 483/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.57 -1.255]
Z position = 0.01004 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.022 (threshold = 0.05)
-----
t = 0.2415s (step 484/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.57 -1.255]
Z position = 0.01004 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.022 (threshold = 0.05)
-----
t = 0.242s (step 485/499)
PLASMA QUANTITIES
Ip = 759.0 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.57 -1.255]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.022 (threshold = 0.05)
-----
t = 0.2425s (step 486/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.57 -1.255]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.022 (threshold = 0.05)
-----
t = 0.243s (step 487/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.57 -1.255]
Z position = 0.01003 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.022 (threshold = 0.05)
-----
t = 0.2435s (step 488/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.57 -1.255]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.022 (threshold = 0.05)
-----
t = 0.244s (step 489/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.351 0.57 -1.255]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.023 (threshold = 0.05)
-----
t = 0.2445s (step 490/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.571 -1.255]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.023 (threshold = 0.05)
-----
t = 0.245s (step 491/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.571 -1.255]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.023 (threshold = 0.05)
-----
t = 0.2455s (step 492/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.571 -1.255]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.023 (threshold = 0.05)
-----
t = 0.246s (step 493/499)
PLASMA QUANTITIES
Ip = 758.9 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.571 -1.255]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.023 (threshold = 0.05)
-----
t = 0.2465s (step 494/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.571 -1.255]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.023 (threshold = 0.05)
-----
t = 0.247s (step 495/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.571 -1.255]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.023 (threshold = 0.05)
-----
t = 0.2475s (step 496/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.571 -1.255]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.023 (threshold = 0.05)
-----
t = 0.248s (step 497/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.571 -1.255]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.023 (threshold = 0.05)
-----
t = 0.2485s (step 498/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.571 -1.255]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.023 (threshold = 0.05)
-----
t = 0.249s (step 499/499)
PLASMA QUANTITIES
Ip = 758.8 [kA]
Shape parameters ['Rin', 'Rout', 'Rx', 'Zx'] = [ 0.281 1.35 0.571 -1.255]
Z position = 0.01002 [m]
Linearised (deformable) growth rate = 23.139 [1/s]
SOLVER STATUS
Relinearisation status = 0.023 (threshold = 0.05)
Plot some results
Let’s now plot some of the output. We can compare the evolution of the controlled parameters (plasma current, shape parameters, and vertical position) against their FB reference waveforms. In these plots:
the black dashed indicated the FB reference waveform.
the solid blue the simulated quantity.
green background shading indicates when FB control is ON.
yellow background shading indicates when FF control is ON (not used).
white background shading indicates no control is used.
# PLASMA CURRENT EVOLUTION
fig, ax = plt.subplots(
nrows=1,
ncols=1,
figsize=(8, 4),
dpi=80
)
# --- references and masks ---
FB_reference = PCS.PlasmaController.interpolants['ip_ref'](times)
FF_reference = PCS.PlasmaController.interpolants['vloop_ff'](times)
blend = PCS.PlasmaController.interpolants['ip_blend'](times)
FB_mask = (blend > 0) & (np.abs(FB_reference) > 0)
FF_mask = (blend < 1) & (np.abs(FF_reference) > 0)
# --- shade FB regions (green) ---
on_regions = np.where(np.diff(FB_mask.astype(int)) != 0)[0] + 1
for seg_t, seg_state in zip(np.split(times, on_regions),
np.split(FB_mask, on_regions)):
if np.all(seg_state):
ax.axvspan(seg_t[0], seg_t[-1],
color='green', alpha=0.25,
label="FB active")
# --- shade FF regions (yellow) ---
on_regions = np.where(np.diff(FF_mask.astype(int)) != 0)[0] + 1
for seg_t, seg_state in zip(np.split(times, on_regions),
np.split(FF_mask, on_regions)):
if np.all(seg_state):
ax.axvspan(seg_t[0], seg_t[-1],
color='gold', alpha=0.25,
label="FF active")
# --- FreeGSNKE ---
ax.plot(times, dynamic_ip,
color='navy', linewidth=1,
marker='x', markersize=0,
label="FreeGSNKE")
# --- FB reference ---
ax.plot(times[FB_mask], FB_reference[FB_mask],
color='k', linestyle='--',
linewidth=1.5,
label="FB reference")
ax.set_xlabel(r"Shot time [$s$]")
ax.set_ylabel("Plasma current [A]")
ax.grid()
# deduplicate legend entries
handles, labels = ax.get_legend_handles_labels()
ax.legend(dict(zip(labels, handles)).values(),
dict(zip(labels, handles)).keys())
ax.set_ylim([7.54e5, 7.64e5])
fig.suptitle("Plasma Current", y=0.98)
plt.tight_layout()
plt.show()
# PLASMA VERTICAL POSITION EVOLUTION
fig, ax = plt.subplots(
nrows=1,
ncols=1,
figsize=(8, 4),
dpi=80
)
# --- FB reference ---
FB_reference = PCS.VerticalController.interpolants['z_ref'](times)
# --- shade FB-active region (always on here) ---
ax.axvspan(times[0], times[-1],
color='green', alpha=0.2,
label="FB active")
# --- references ---
ax.plot(times, FB_reference,
color='k', linestyle='--',
linewidth=1.5,
label="FB reference")
# --- FreeGSNKE ---
ax.plot(times, z_current,
color='navy', linewidth=1,
marker='x', markersize=0,
label="FreeGSNKE")
ax.set_xlabel(r"Shot time [$s$]")
ax.set_ylabel("Vertical position [m]")
ax.grid()
ax.legend()
fig.suptitle("Vertical Position", y=0.98)
plt.tight_layout()
plt.show()
# SHAPE PARAMETER EVOLUTION
# the ramp down in Rx here is designed to show how one might increase triangularity of the plasma (see next plot for this effect)
ntarg = len(ctrl_targets)
fig, axes = plt.subplots(
nrows=ntarg,
ncols=1,
figsize=(8, 12),
dpi=80,
sharex=True
)
for j, targ in enumerate(ctrl_targets):
ax = axes[j]
# --- references and masks ---
FF_reference = PCS.ShapeController.interpolants[targ]['ff'](times)
FB_reference = PCS.ShapeController.interpolants[targ]['ref'](times)
FF_mask = (
(PCS.ShapeController.interpolants[targ]['blend'](times) < 1)
& (np.abs(PCS.ShapeController.interpolants[targ]['ff'].derivative()(times)) > 0)
)
FB_mask = (
(PCS.ShapeController.interpolants[targ]['blend'](times) > 0)
& (np.abs(FB_reference) > 0)
)
# --- shade FB regions (green) ---
on_regions = np.where(np.diff(FB_mask.astype(int)) != 0)[0] + 1
for seg_t, seg_state in zip(np.split(times, on_regions),
np.split(FB_mask, on_regions)):
if np.all(seg_state):
ax.axvspan(seg_t[0], seg_t[-1], color='green', alpha=0.25,
label="FB active" if j == 0 else None)
# --- shade FF regions (yellow) ---
on_regions = np.where(np.diff(FF_mask.astype(int)) != 0)[0] + 1
for seg_t, seg_state in zip(np.split(times, on_regions),
np.split(FF_mask, on_regions)):
if np.all(seg_state):
ax.axvspan(seg_t[0], seg_t[-1], color='gold', alpha=0.25,
label="FF active" if j == 0 else None)
# --- references ---
ax.plot(times[FB_mask], FB_reference[FB_mask],
color='k', linestyle='--', linewidth=1.5,
label="FB reference" if j == 0 else None)
if np.any(FF_mask):
ax.plot(times[FF_mask], FF_reference[FF_mask],
color='r', linestyle='--', linewidth=1.5,
label="FF reference" if j == 0 else None)
# --- FreeGSNKE ---
ax.plot(times, dynamic_shape_targets[j, :],
color='navy', linewidth=1,
marker='x', markersize=0,
label="FreeGSNKE" if j == 0 else None)
ax.set_ylabel(f"{ctrl_targets[j]} [m]")
ax.grid()
ax.set_ylim([np.min(FB_reference)-0.02, np.max(FB_reference)+0.02])
# shared x label
axes[-1].set_xlabel(r"Shot time [$s$]")
# legend once
axes[0].legend(ncol=4, loc="upper right")
fig.suptitle("Shape parameters", y=0.995, fontsize=14)
plt.tight_layout()
plt.show()
# PLASMA TRIANGULARITY EVOLUTION
fig, ax = plt.subplots(
nrows=1,
ncols=1,
figsize=(8, 4),
dpi=80
)
# --- FreeGSNKE ---
ax.plot(times[0:-1], dynamic_triangularity[0:-1],
color='navy', linewidth=1,
marker='x', markersize=0,
label="FreeGSNKE")
ax.set_xlabel(r"Shot time [$s$]")
ax.set_ylabel("Triangularity")
ax.grid()
ax.legend()
plt.tight_layout()
plt.show()
We can also show the voltages produced by the PCS and the subsequent evolution of the currents for each of the PF coils. To display the voltage/current limits, uncomment the relevant sections in the next cell.
# PF COIL VOLTAGES AND CURRENTS
ncoils = len(active_coils)
fig, axes = plt.subplots(
nrows=ncoils,
ncols=2,
figsize=(16, 45),
dpi=80,
sharex=True
)
for j, coil in enumerate(active_coils):
axV = axes[j, 0] # voltage column
axI = axes[j, 1] # current column
# VOLTAGES (left)
# # voltage limits
# if coil not in ["P6"]:
# vlim = PCS.PFController.data['coil_voltage_lims']['vals'][0][j]
# axV.hlines([-vlim, vlim], times[0], times[-1],
# colors='k', linestyles='--', linewidth=1.2,
# label="Coil limits" if j == 0 else None)
axV.plot(times[0:-1], V_approved[j, 0:-1],
color='navy', linewidth=1,
marker='x', markersize=0,
label="FreeGSNKE")
axV.set_ylabel(f'{coil} voltage [V]')
axV.grid()
if j == 0:
axV.legend()
# CURRENTS (right)
# # current limits
# if coil not in ["P6"]:
# imin = PCS.SystemsController.data['min_coil_curr_lims']['vals'][0][j]
# imax = PCS.SystemsController.data['max_coil_curr_lims']['vals'][0][j]
# axI.hlines([imin, imax], times[0], times[-1],
# colors='k', linestyles='--', linewidth=1.2,
# label="Coil limits" if j == 0 else None)
axI.plot(times, dynamic_currents[j, :],
color='navy', linewidth=1,
marker='x', markersize=0,
label="FreeGSNKE")
axI.set_ylabel(f'{coil} current [A]')
axI.grid()
if j == 0:
axI.legend()
# x-labels only on bottom row
axes[-1, 0].set_xlabel(r'Shot time [$s$]')
axes[-1, 1].set_xlabel(r'Shot time [$s$]')
# column titles
axes[0, 0].set_title("PF Coil Voltages")
axes[0, 1].set_title("PF Coil Currents")
plt.tight_layout()
plt.show()
Given that we have passive conducting structures in the machine, we can also plot the evolution of their currents. Note how they all start from zero given that we did not initialise them in this example.
# PASSIVE CURRENTS EVOLUTION
# enter the names of some passives here
passives_to_plot = ["vessel_1", "centrecolumn_1", "gas_baffle_upper_1", "p5_case_upper_1"]
# indices of passives_to_plot in tokamak.coil_names
indxs = [tokamak.coil_names.index(name) for name in passives_to_plot]
fig, ax = plt.subplots(
nrows=1,
ncols=1,
figsize=(8, 4),
dpi=80
)
# --- FreeGSNKE ---
for name, idx in zip(passives_to_plot, indxs):
ax.plot(
times, dynamic_currents[idx, :],
linewidth=1,
marker='x', markersize=0,
label=name,
)
ax.set_xlabel(r"Shot time [$s$]")
ax.set_ylabel("Currents [A]")
ax.grid()
ax.legend()
fig.suptitle("Passive currents evolution", y=0.98)
plt.tight_layout()
plt.show()
Here we see the evolution of the vertical instability timescale (i.e. reciprocal of the growth rate). We can see that the plasma remains stable, with the timescale only changing after the internal Jacobians are re-linearised (around 0.06 seconds).
# LINEAR INSTABILITY TIMESCALE EVOLUTION (1/GROWTH_RATE)
fig, ax = plt.subplots(
nrows=1,
ncols=1,
figsize=(8, 4),
dpi=80
)
# plot results
ax.plot(times[0:len(linear_growth_rate)], linear_growth_rate, color='navy', linewidth=1, marker=None, markersize=3, label="FPDT")
ax.hlines(dt, xmin=times[0], xmax=times[-1], color='black', linestyle='--', linewidth=1.5, label=f"dt")
ax.set_ylim([1e-4, 1e-1])
ax.set_yscale('log')
# ax.legend()
ax.grid()
ax.set_xlabel(r'Shot time [$s$]')
ax.set_ylabel('Linearised instability timescale [s]')
fig.suptitle("Vertical instability timescale (linearised deformable)", y=0.98)
plt.tight_layout()
plt.show()
The following cell can be uncommented to make an animation of the equilibrium in the machine over time. Might take a few seconds to run and will save the output .mp4 in the data directory.
# # SAVE AN ANIMATION OF THE SIMULATION
# %matplotlib inline
# import matplotlib.animation as animation
# plt.rcParams['figure.dpi'] = 90
# fig1, ax1 = plt.subplots(1, 1, figsize=(8, 8))
# # Plot static wall / tokamak background
# eq.tokamak.plot(axis=ax1, show=False)
# ax1.plot(tokamak.wall.R, tokamak.wall.Z, color='k', linewidth=1.2)
# ax1.grid(True, which='both', alpha=0.5)
# ax1.set_aspect('equal')
# ax1.set_xlabel(r'Major radius, $R$ $[m]$')
# ax1.set_ylabel(r'Height, $Z$ $[m]$')
# ax1.set_xlim(0.05, 2.15)
# ax1.set_ylim(-2.25, 2.25)
# # Determine contour levels from entire dynamic_psi
# min_psi = np.min(dynamic_psi)
# max_psi = np.max(dynamic_psi)
# levels = np.linspace(min_psi, max_psi, 40)
# # --- Storage for dynamic artists ---
# contour_artists = []
# scatter_artists = []
# # --- Update function ---
# def update(i):
# global contour_artists, scatter_artists
# # Remove previous dynamic artists
# for c in contour_artists + scatter_artists:
# if isinstance(c, list):
# for coll in c:
# coll.remove()
# else:
# c.remove()
# contour_artists = []
# scatter_artists = []
# ax1.set_title(rf"$t$ = {np.round(times[i],3)}")
# # Main psi contours
# c1 = ax1.contour(eq.R, eq.Z, dynamic_psi[:,:,i], levels=levels, alpha=0.8, cmap='viridis')
# contour_artists.append(c1)
# # plasma boundary
# c2 = ax1.contour(eq.R, eq.Z, dynamic_psi[:,:,i],
# levels=[dynamic_psi_boundary[i]],
# linestyles="-", colors='r', linewidths=1.4)
# contour_artists.append(c2)
# # adds separatrix of primary X-point if plasma limited
# if dynamic_limiter_flag[i]:
# c3 = ax1.contour(eq.R, eq.Z, dynamic_psi[:,:,i],
# levels=[dynamic_xpts[i][0,2]],
# linestyles="--", colors='k', linewidths=1.4)
# contour_artists.append(c3)
# # X-points and O-points
# sc1 = ax1.scatter(dynamic_xpts[i][:,0], dynamic_xpts[i][:,1], color='r', marker='x', s=30)
# sc2 = ax1.scatter(dynamic_opts[i][:,0], dynamic_opts[i][:,1], color='g', marker='2', s=30)
# # indicate which is the primary X-point more clearly
# sc3 = ax1.scatter(dynamic_xpts[i][0,0], dynamic_xpts[i][0,1], color='r', marker='x', s=60)
# scatter_artists.extend([sc1, sc2, sc3])
# return contour_artists + scatter_artists
# # --- Animation ---
# num_frames = len(times)
# frames_to_plot = range(0, num_frames, 5)
# fps = int(len(frames_to_plot)/10)
# ani = animation.FuncAnimation(fig1, update, frames=frames_to_plot, interval=1, blit=True)
# # save to video (much faster & smaller than GIF)
# ani.save(f"data/animated_equilibrium.mp4", writer="ffmpeg", fps=fps)