Contents

RZIP vs Free-boundary simulation with controller

This example shows how to perform a closed-loop simulation for magnetic position and current (R,Z,Ip) control. As well as control of coil currents orthogonal to those used for R,Z,Ip control.

We perform the simulation both with rzp (rigid body model), with fge (free boundary model), and with fgel (linearized model of fge) and plot the differencess

[+MEQ MatlabEQuilibrium Toolbox+] Swiss Plasma Center EPFL Lausanne 2022. All rights reserved.

% Set some basic parameters
dt = 1e-4; % step time time for simulation
tok = 'ana'; % anamak tokamak
shot = 2;  % anamak shot 2 is elongated
t0 = 0;    % start time
debuglevel = 0; % debugging level - set to 0 for publication

% concatenate parameters, add convergence tolerance and CDE
params = {tok,shot,t0,'izgrid',true,...
  'tol',1e-4,'tolF',1e-4,... % low tolerance both for fbt and fge
  'cde','OhmTor_rigid',...
  'ilim',3,'icsint',true}; % spline interpolation for limiter and axis finding

% Compute FGE and RZP models
disp('Computing FGE linearization')
[L,LX] = fge(params{:});

disp('Computing RZP linearization')
Lz  = rzp(params{:});
Lz  = rzpl(Lz,LX); % linearization for this equilibrium
Computing FGE linearization
Computing RZP linearization

This equilibrium is unstable, look at the eigenmode

figure(1);
subplot(121);fgeploteig(L) ; title(['FGEL',getfield(get(gca,'Title'),'String')])
subplot(122);fgeploteig(Lz); title(['RZP ',getfield(get(gca,'Title'),'String')])
% FGE and RZIP show similar growth rates and eigenmode structures

Compute the time-varying Va and Iadot vector to sustain plasma current

Since we model the resistive plasma current decay, a loop voltage is needed to maintain the plasma current. This loop voltage is generated by varying the poloidal field currents in time along a particular vector. Also a time-varying Va (coil voltage) is needed to sustain these time-varying voltages.

t  = t0+(0:dt:0.1); % time grid
LX = fgex([],t,L,LX);

Compute controller parameters

Use meqctrl function that implements a simple R,Z,Ip controller

help meqctrl

debugdesign = debuglevel; % optional debugging of controller design, set to 1 to see bode diagrams

% Inspect meqctrl for details of what this function does.
[~,~,ctrlpar] = meqctrl('setup',[],[],L,LX,'debug',debugdesign,'dt',dt);
  [cstate,Va,ctrlpar,ctdebug] = meqctrl(mode,cstate,ctrlpar,L,LY,varargin)
 
  Generic magnetic controller for meq - implements a simple R,Z,Ip
  controller for anamak and provides extension options for other tokamaks.
  See tutorials/anamak_rzp_control.m for example of its use for anamak.
  Also serves as template for possible design of other controllers.
 
  Inputs:
    mode: 'setup','init','step' or 'reset' (see later)
    cstate: vector with the controller state at the previous time step.
    ctrlpar: controller parameter structure as input. For details, inspect this function
    L,LY: MEQ L and LY (equilibrium) structures
    varargin: optional parameter-value pairs to pass into ctrlpar
 
  Outputs:
    cstate: vector with the controller state at the next time step
    Va is the controller output voltage (including feedforward part
    ctrlpar are the controller parameters, returned in the 'setup' mode
    ctdebug are optional debugging information.
 
  
  MODES:
  'setup' mode:
    Performs controller design and returns ctrlpar. Calls a tokamak-specific function 
    meqctrlp<tok> to load the controller parameters.
    * First, a mutual decoupling and resistive compensation term of the form 
      MR=Maa+Ra/s is placed in front of the open-loop plant.
    * A vertical control loop (KzQ) is closed on the voltages directly. By default
      the control direction is chosen to maximize the coupling to the
      most unstable eigenmode, but custom directions can also be chosen.
    * A coil current control loop (KIo) orthogonal to the vertical control loop 
      is closed as well. This gives a stable plant with inputs being the
      ortogonal coil currents and the vertical position reference.
    * A radial and current control direction is chosen using a pseudoinverse
      of the plant gain at low frequencies, to control rIp,Ip while not
      perturbing zIp
    * The rIp, Ip loops are closed with their own controllers (KrQ, KIQ)
      
    A note on nomenclature, the suffix 'Q' refers to [Piras, PhD thesis EPFL 2010] 
    and denotes the 'common' mode (mnemonic: 'Q'ommon), to distinguish if 
    from the symmetric mode 'S' for doublets.
 
  'init' mode: Initializes the state-space matrices for the controller as 
    persistent variables and returns the initial state and initial command voltage 
    (called by fget in the initialization phase)
 
  'step' mode: Executes controller step, receives previous state, parameters, 
    and LY from current time step and runs one step of the controller, 
    returning the control voltages and next controller state.
    (Called by fget in the time step)
  
  'reset' mode: Clears the persistent variables
 
  Setting debug=true plots some transfer functions in the design process

Check controller closed-loop behaviour

ctrlpar.Kmain = 1;
[cstate,V,ctrlpar,ctdebug] = meqctrl('init',[],ctrlpar);
KKd = ctdebug.KKd;

z=tf('z',dt);
KKd=KKd/z; % add unit delay

% Linear system
Ts=dt; % continuous time
sysd = fgess(Lz,Ts,{'Ia','zIp','rIp','Ip'});
sysd = sysd(:,'Va'); % Va only

syscl=feedback(sysd*KKd,eye(size(sysd,1)));

% Plot transfer function for z,r,Ip closed-loop step-responses
izrp = [sysd.OutputGroup.zIp;sysd.OutputGroup.rIp;sysd.OutputGroup.Ip];
step(syscl(izrp,izrp),0.2)
Warning: The "hold" command is not supported for multi-axes plots. Clearing the
previous plot. 

Define targets for closed-loop simulation

ctrlpar.debug = debuglevel;
% perturb Z/R ref
ctrlpar.zQref = timeseries((LX.zIp/LX.Ip(:,1) + interp1([0 0.005 0.01 0.02],+[0 0 0.02 0.02],LX.t,'linear','extrap'))',...
    LX.t);
ctrlpar.rQref = timeseries((LX.rIp/LX.Ip(:,1) + interp1([0 0.02 0.05 0.1],+[0 0.00 0.02 0.02],LX.t,'linear','extrap'))',...
    LX.t);
ctrlpar.IQref = timeseries((LX.Ip + interp1([0 0.06 0.07 0.1],-[0 0.00 20e3 20e3],LX.t,'linear','extrap'))',...
    LX.t);

% perturb bp as a disturbance
LX.bp = LX.bp(1) + interp1([0 0.07 0.09 0.1],[0 0 1 1],LX.t,'linear','extrap');

% Set control function handle and parameter
L.P.ctrlfct  = @meqctrl;
L.P.ctrlpar = ctrlpar;

% show debugging info during sim
L.P.debug=debuglevel;

Run FGE simulation

LY = fget(L,LX,'debug',1);
 *** FGE ***
  Tokamak: anamak, shot#2, t=0.000
  Solver: all-nl + jfnkGMRES: mgsaya
  prec:1, mkryl:200, kmax:50   tolF:1.000e-04
  Vessel:v, nu:200
  Doublet: 0
  Constraints: 
    Fun    iD    LXfield 
   --------------------- 
    {@Ip}    {[1]}    {'Ip'}    {[1]}
    {@bp}    {[1]}    {'bp'}    {[1]}
    {@qA}    {[1]}    {'qA'}    {[1]}




        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.00010 |   2 |  200.00 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      5 |    1 | 5.81e-07
  0.00020 |   3 |  200.00 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 5.92e-07
  0.00030 |   4 |  200.00 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 6.05e-07
  0.00040 |   5 |  200.00 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 6.49e-07
  0.00050 |   6 |  200.00 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 6.30e-07
  0.00060 |   7 |  200.00 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 6.43e-07
  0.00070 |   8 |  200.00 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 6.52e-07
  0.00080 |   9 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 6.54e-07
  0.00090 |  10 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 6.43e-07
  0.00100 |  11 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 6.16e-07

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.00110 |  12 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 6.19e-07
  0.00120 |  13 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 5.93e-07
  0.00130 |  14 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 5.87e-07
  0.00140 |  15 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 5.36e-07
  0.00150 |  16 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 5.07e-07
  0.00160 |  17 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 5.01e-07
  0.00170 |  18 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 4.64e-07
  0.00180 |  19 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 4.10e-07
  0.00190 |  20 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 3.88e-07
  0.00200 |  21 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 3.66e-07

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.00210 |  22 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 3.28e-07
  0.00220 |  23 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 2.80e-07
  0.00230 |  24 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 2.58e-07
  0.00240 |  25 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 2.22e-07
  0.00250 |  26 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 1.98e-07
  0.00260 |  27 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 1.52e-07
  0.00270 |  28 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 1.24e-07
  0.00280 |  29 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 8.58e-08
  0.00290 |  30 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 4.31e-08
  0.00300 |  31 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 2.88e-08

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.00310 |  32 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 2.92e-08
  0.00320 |  33 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 5.19e-08
  0.00330 |  34 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 3.61e-08
  0.00340 |  35 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 7.25e-08
  0.00350 |  36 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 1.05e-07
  0.00360 |  37 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 9.15e-08
  0.00370 |  38 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 1.24e-07
  0.00380 |  39 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 1.43e-07
  0.00390 |  40 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 1.60e-07
  0.00400 |  41 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 1.36e-07

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.00410 |  42 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 1.81e-07
  0.00420 |  43 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 1.86e-07
  0.00430 |  44 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 1.72e-07
  0.00440 |  45 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 1.71e-07
  0.00450 |  46 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 2.08e-07
  0.00460 |  47 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 1.81e-07
  0.00470 |  48 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      7 |    1 | 1.71e-07
  0.00480 |  49 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      7 |    1 | 2.06e-07
  0.00490 |  50 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      7 |    1 | 1.88e-07
  0.00500 |  51 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      7 |    1 | 1.92e-07

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.00510 |  52 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      7 |    1 | 1.75e-07
  0.00520 |  53 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      6 |    1 | 1.88e-07
  0.00530 |  54 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      7 |    1 | 2.51e-07
  0.00540 |  55 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      7 |    1 | 4.29e-07
  0.00550 |  56 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      7 |    1 | 7.09e-06
  0.00560 |  57 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      7 |    1 | 2.36e-05
  0.00570 |  58 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      7 |    1 | 5.15e-05
  0.00580 |  59 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     1 |      7 |    1 | 8.88e-05
  0.00590 |  60 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     2 |     14 |    1 | 2.39e-09
  0.00600 |  61 |  199.99 |  0.059 |  1.112 |  1.000 |  1.000 |     2 |     15 |    1 | 5.14e-09

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.00610 |  62 |  200.00 |  0.060 |  1.112 |  1.000 |  1.000 |     2 |     15 |    1 | 7.72e-09
  0.00620 |  63 |  200.00 |  0.060 |  1.112 |  1.000 |  1.000 |     2 |     16 |    1 | 2.32e-08
  0.00630 |  64 |  200.00 |  0.060 |  1.112 |  1.000 |  1.000 |     2 |     16 |    1 | 1.42e-08
  0.00640 |  65 |  200.00 |  0.061 |  1.112 |  1.000 |  1.000 |     2 |     16 |    1 | 1.67e-08
  0.00650 |  66 |  200.00 |  0.061 |  1.112 |  1.000 |  1.000 |     2 |     16 |    1 | 1.93e-08
  0.00660 |  67 |  200.00 |  0.061 |  1.113 |  1.000 |  1.000 |     2 |     18 |    1 | 2.21e-08
  0.00670 |  68 |  200.00 |  0.061 |  1.113 |  1.000 |  1.000 |     2 |     18 |    1 | 3.06e-08
  0.00680 |  69 |  200.00 |  0.062 |  1.113 |  1.000 |  1.000 |     2 |     18 |    1 | 3.42e-08
  0.00690 |  70 |  200.00 |  0.062 |  1.113 |  1.000 |  1.000 |     2 |     18 |    1 | 5.92e-08
  0.00700 |  71 |  200.01 |  0.063 |  1.113 |  1.000 |  1.000 |     2 |     18 |    1 | 4.82e-08

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.00710 |  72 |  200.01 |  0.063 |  1.113 |  1.000 |  1.000 |     2 |     18 |    1 | 4.70e-08
  0.00720 |  73 |  200.01 |  0.063 |  1.113 |  1.000 |  1.000 |     2 |     18 |    1 | 4.56e-08
  0.00730 |  74 |  200.02 |  0.064 |  1.113 |  1.000 |  1.000 |     2 |     20 |    1 | 3.24e-09
  0.00740 |  75 |  200.03 |  0.064 |  1.113 |  1.000 |  1.000 |     2 |     20 |    1 | 6.86e-08
  0.00750 |  76 |  200.04 |  0.065 |  1.113 |  1.000 |  1.000 |     2 |     20 |    1 | 6.30e-08
  0.00760 |  77 |  200.05 |  0.065 |  1.113 |  1.000 |  1.000 |     2 |     20 |    1 | 5.86e-08
  0.00770 |  78 |  200.05 |  0.066 |  1.113 |  1.000 |  1.000 |     2 |     20 |    1 | 5.61e-08
  0.00780 |  79 |  200.06 |  0.066 |  1.113 |  1.000 |  1.000 |     2 |     20 |    1 | 5.13e-08
  0.00790 |  80 |  200.06 |  0.066 |  1.114 |  1.000 |  1.000 |     2 |     20 |    1 | 4.70e-08
  0.00800 |  81 |  200.06 |  0.067 |  1.114 |  1.000 |  1.000 |     2 |     20 |    1 | 5.04e-08

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.00810 |  82 |  200.06 |  0.067 |  1.114 |  1.000 |  1.000 |     2 |     20 |    1 | 4.31e-08
  0.00820 |  83 |  200.06 |  0.068 |  1.114 |  1.000 |  1.000 |     2 |     20 |    1 | 4.50e-08
  0.00830 |  84 |  200.06 |  0.068 |  1.114 |  1.000 |  1.000 |     2 |     20 |    1 | 4.03e-08
  0.00840 |  85 |  200.06 |  0.068 |  1.114 |  1.000 |  1.000 |     2 |     20 |    1 | 3.58e-08
  0.00850 |  86 |  200.06 |  0.069 |  1.114 |  1.000 |  1.000 |     2 |     20 |    1 | 3.52e-08
  0.00860 |  87 |  200.05 |  0.069 |  1.114 |  1.000 |  1.000 |     2 |     20 |    1 | 3.28e-08
  0.00870 |  88 |  200.05 |  0.070 |  1.114 |  1.000 |  1.000 |     2 |     20 |    1 | 2.54e-08
  0.00880 |  89 |  200.05 |  0.070 |  1.114 |  1.000 |  1.000 |     2 |     21 |    1 | 4.78e-08
  0.00890 |  90 |  200.04 |  0.071 |  1.114 |  1.000 |  1.000 |     2 |     21 |    1 | 4.43e-08
  0.00900 |  91 |  200.04 |  0.071 |  1.114 |  1.000 |  1.000 |     2 |     21 |    1 | 5.07e-08

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.00910 |  92 |  200.03 |  0.072 |  1.114 |  1.000 |  1.000 |     2 |     21 |    1 | 4.88e-08
  0.00920 |  93 |  200.03 |  0.072 |  1.114 |  1.000 |  1.000 |     2 |     22 |    1 | 1.04e-07
  0.00930 |  94 |  200.02 |  0.072 |  1.114 |  1.000 |  1.000 |     2 |     21 |    1 | 4.67e-08
  0.00940 |  95 |  200.02 |  0.073 |  1.114 |  1.000 |  1.000 |     2 |     21 |    1 | 4.42e-08
  0.00950 |  96 |  200.01 |  0.073 |  1.114 |  1.000 |  1.000 |     2 |     21 |    1 | 3.81e-08
  0.00960 |  97 |  200.00 |  0.074 |  1.114 |  1.000 |  1.000 |     2 |     22 |    1 | 2.29e-09
  0.00970 |  98 |  200.00 |  0.074 |  1.114 |  1.000 |  1.000 |     2 |     22 |    1 | 1.55e-07
  0.00980 |  99 |  199.99 |  0.075 |  1.114 |  1.000 |  1.000 |     2 |     22 |    1 | 4.16e-08
  0.00990 | 100 |  199.99 |  0.075 |  1.114 |  1.000 |  1.000 |     2 |     22 |    1 | 4.11e-08
  0.01000 | 101 |  199.98 |  0.076 |  1.114 |  1.000 |  1.000 |     2 |     24 |    1 | 8.24e-07

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.01010 | 102 |  199.97 |  0.076 |  1.114 |  1.000 |  1.000 |     2 |     23 |    1 | 7.39e-08
  0.01020 | 103 |  199.96 |  0.077 |  1.114 |  1.000 |  1.000 |     2 |     23 |    1 | 3.18e-08
  0.01030 | 104 |  199.96 |  0.077 |  1.114 |  1.000 |  1.000 |     2 |     23 |    1 | 2.71e-08
  0.01040 | 105 |  199.95 |  0.077 |  1.114 |  1.000 |  1.000 |     2 |     23 |    1 | 4.04e-08
  0.01050 | 106 |  199.94 |  0.078 |  1.114 |  1.000 |  1.000 |     2 |     23 |    1 | 2.62e-10
  0.01060 | 107 |  199.94 |  0.078 |  1.114 |  1.000 |  1.000 |     2 |     23 |    1 | 1.12e-08
  0.01070 | 108 |  199.93 |  0.078 |  1.114 |  1.000 |  1.000 |     2 |     23 |    1 | 6.91e-09
  0.01080 | 109 |  199.93 |  0.079 |  1.114 |  1.000 |  1.000 |     2 |     23 |    1 | 3.92e-09
  0.01090 | 110 |  199.93 |  0.079 |  1.114 |  1.000 |  1.000 |     2 |     23 |    1 | 1.82e-09
  0.01100 | 111 |  199.93 |  0.079 |  1.114 |  1.000 |  1.000 |     2 |     23 |    1 | 8.29e-10

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.01110 | 112 |  199.92 |  0.079 |  1.114 |  1.000 |  1.000 |     2 |     23 |    1 | 7.49e-10
  0.01120 | 113 |  199.92 |  0.080 |  1.114 |  1.000 |  1.000 |     2 |     23 |    1 | 3.53e-10
  0.01130 | 114 |  199.92 |  0.080 |  1.114 |  1.000 |  1.000 |     2 |     23 |    1 | 6.36e-11
  0.01140 | 115 |  199.92 |  0.080 |  1.114 |  1.000 |  1.000 |     2 |     23 |    1 | 1.53e-10
  0.01150 | 116 |  199.92 |  0.080 |  1.113 |  1.000 |  1.000 |     2 |     23 |    1 | 6.19e-11
  0.01160 | 117 |  199.92 |  0.080 |  1.113 |  1.000 |  1.000 |     2 |     23 |    1 | 5.55e-11
  0.01170 | 118 |  199.93 |  0.080 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 9.81e-05
  0.01180 | 119 |  199.93 |  0.080 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 7.37e-05
  0.01190 | 120 |  199.93 |  0.080 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 5.39e-05
  0.01200 | 121 |  199.93 |  0.080 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 3.77e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.01210 | 122 |  199.93 |  0.080 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 2.48e-05
  0.01220 | 123 |  199.93 |  0.081 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 1.45e-05
  0.01230 | 124 |  199.94 |  0.081 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 6.50e-06
  0.01240 | 125 |  199.94 |  0.081 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 1.07e-06
  0.01250 | 126 |  199.94 |  0.081 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 3.73e-06
  0.01260 | 127 |  199.94 |  0.081 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 6.76e-06
  0.01270 | 128 |  199.94 |  0.081 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 8.62e-06
  0.01280 | 129 |  199.95 |  0.080 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 9.48e-06
  0.01290 | 130 |  199.95 |  0.080 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 9.32e-06
  0.01300 | 131 |  199.95 |  0.080 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 8.81e-06

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.01310 | 132 |  199.96 |  0.080 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 7.89e-06
  0.01320 | 133 |  199.96 |  0.080 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 6.44e-06
  0.01330 | 134 |  199.97 |  0.080 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 4.40e-06
  0.01340 | 135 |  199.97 |  0.080 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 2.45e-06
  0.01350 | 136 |  199.97 |  0.080 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 2.13e-07
  0.01360 | 137 |  199.98 |  0.080 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 2.03e-06
  0.01370 | 138 |  199.98 |  0.080 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 4.26e-06
  0.01380 | 139 |  199.99 |  0.080 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 6.49e-06
  0.01390 | 140 |  199.99 |  0.080 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 8.55e-06
  0.01400 | 141 |  200.00 |  0.080 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 1.07e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.01410 | 142 |  200.00 |  0.080 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 1.26e-05
  0.01420 | 143 |  200.01 |  0.080 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 1.46e-05
  0.01430 | 144 |  200.01 |  0.080 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 1.61e-05
  0.01440 | 145 |  200.02 |  0.080 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 1.76e-05
  0.01450 | 146 |  200.02 |  0.080 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 1.92e-05
  0.01460 | 147 |  200.03 |  0.080 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 2.04e-05
  0.01470 | 148 |  200.03 |  0.080 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 2.14e-05
  0.01480 | 149 |  200.04 |  0.079 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 2.23e-05
  0.01490 | 150 |  200.04 |  0.079 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 2.31e-05
  0.01500 | 151 |  200.05 |  0.079 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 2.36e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.01510 | 152 |  200.05 |  0.079 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 2.41e-05
  0.01520 | 153 |  200.06 |  0.079 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 2.42e-05
  0.01530 | 154 |  200.06 |  0.079 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 2.45e-05
  0.01540 | 155 |  200.07 |  0.079 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 2.45e-05
  0.01550 | 156 |  200.07 |  0.079 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 2.44e-05
  0.01560 | 157 |  200.08 |  0.079 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 2.42e-05
  0.01570 | 158 |  200.08 |  0.079 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 2.40e-05
  0.01580 | 159 |  200.09 |  0.079 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 2.37e-05
  0.01590 | 160 |  200.09 |  0.079 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 2.33e-05
  0.01600 | 161 |  200.10 |  0.079 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 2.29e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.01610 | 162 |  200.10 |  0.079 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 2.25e-05
  0.01620 | 163 |  200.11 |  0.079 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 2.20e-05
  0.01630 | 164 |  200.11 |  0.079 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 2.14e-05
  0.01640 | 165 |  200.12 |  0.079 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 2.09e-05
  0.01650 | 166 |  200.12 |  0.078 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 2.03e-05
  0.01660 | 167 |  200.12 |  0.078 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.97e-05
  0.01670 | 168 |  200.13 |  0.078 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.92e-05
  0.01680 | 169 |  200.13 |  0.078 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.86e-05
  0.01690 | 170 |  200.13 |  0.078 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.81e-05
  0.01700 | 171 |  200.14 |  0.078 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.76e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.01710 | 172 |  200.14 |  0.078 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.70e-05
  0.01720 | 173 |  200.14 |  0.078 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.64e-05
  0.01730 | 174 |  200.15 |  0.078 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.59e-05
  0.01740 | 175 |  200.15 |  0.078 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.54e-05
  0.01750 | 176 |  200.15 |  0.078 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.49e-05
  0.01760 | 177 |  200.15 |  0.078 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.45e-05
  0.01770 | 178 |  200.15 |  0.078 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.40e-05
  0.01780 | 179 |  200.16 |  0.078 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.35e-05
  0.01790 | 180 |  200.16 |  0.078 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.31e-05
  0.01800 | 181 |  200.16 |  0.078 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.27e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.01810 | 182 |  200.16 |  0.078 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.23e-05
  0.01820 | 183 |  200.16 |  0.078 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.19e-05
  0.01830 | 184 |  200.16 |  0.078 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.15e-05
  0.01840 | 185 |  200.16 |  0.078 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.12e-05
  0.01850 | 186 |  200.17 |  0.078 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.09e-05
  0.01860 | 187 |  200.17 |  0.078 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.06e-05
  0.01870 | 188 |  200.17 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.02e-05
  0.01880 | 189 |  200.17 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 9.92e-06
  0.01890 | 190 |  200.17 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 9.69e-06
  0.01900 | 191 |  200.17 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 9.40e-06

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.01910 | 192 |  200.17 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 8.55e-06
  0.01920 | 193 |  200.17 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 8.28e-06
  0.01930 | 194 |  200.17 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 8.01e-06
  0.01940 | 195 |  200.17 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 7.83e-06
  0.01950 | 196 |  200.17 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 7.63e-06
  0.01960 | 197 |  200.17 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 7.42e-06
  0.01970 | 198 |  200.17 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 7.28e-06
  0.01980 | 199 |  200.16 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 7.12e-06
  0.01990 | 200 |  200.16 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 1.06e-05
  0.02000 | 201 |  200.16 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 6.94e-06

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.02010 | 202 |  200.16 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 6.75e-06
  0.02020 | 203 |  200.16 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 6.66e-06
  0.02030 | 204 |  200.16 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 6.48e-06
  0.02040 | 205 |  200.16 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 6.33e-06
  0.02050 | 206 |  200.16 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 6.21e-06
  0.02060 | 207 |  200.16 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 6.16e-06
  0.02070 | 208 |  200.16 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 6.12e-06
  0.02080 | 209 |  200.15 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 6.11e-06
  0.02090 | 210 |  200.15 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 6.11e-06
  0.02100 | 211 |  200.15 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 6.19e-06

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.02110 | 212 |  200.15 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 6.23e-06
  0.02120 | 213 |  200.15 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 6.29e-06
  0.02130 | 214 |  200.15 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     12 |    1 | 6.27e-06
  0.02140 | 215 |  200.15 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 6.22e-06
  0.02150 | 216 |  200.15 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 6.05e-06
  0.02160 | 217 |  200.15 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 5.78e-06
  0.02170 | 218 |  200.15 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 5.41e-06
  0.02180 | 219 |  200.15 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 4.87e-06
  0.02190 | 220 |  200.15 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 4.14e-06
  0.02200 | 221 |  200.15 |  0.077 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 2.71e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.02210 | 222 |  200.15 |  0.076 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 3.06e-06
  0.02220 | 223 |  200.15 |  0.076 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 1.89e-06
  0.02230 | 224 |  200.15 |  0.076 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 5.86e-07
  0.02240 | 225 |  200.15 |  0.076 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 8.12e-07
  0.02250 | 226 |  200.16 |  0.076 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 2.41e-06
  0.02260 | 227 |  200.16 |  0.076 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 4.16e-06
  0.02270 | 228 |  200.16 |  0.076 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 5.93e-06
  0.02280 | 229 |  200.16 |  0.076 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 7.85e-06
  0.02290 | 230 |  200.16 |  0.076 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 1.00e-05
  0.02300 | 231 |  200.17 |  0.076 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 1.21e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.02310 | 232 |  200.17 |  0.076 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 1.40e-05
  0.02320 | 233 |  200.17 |  0.076 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 1.63e-05
  0.02330 | 234 |  200.17 |  0.076 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 1.82e-05
  0.02340 | 235 |  200.17 |  0.076 |  1.111 |  1.000 |  1.000 |     1 |     13 |    1 | 2.04e-05
  0.02350 | 236 |  200.18 |  0.076 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 2.22e-05
  0.02360 | 237 |  200.18 |  0.076 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 2.42e-05
  0.02370 | 238 |  200.18 |  0.076 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 2.59e-05
  0.02380 | 239 |  200.18 |  0.076 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 2.77e-05
  0.02390 | 240 |  200.18 |  0.076 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 2.93e-05
  0.02400 | 241 |  200.19 |  0.076 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 3.09e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.02410 | 242 |  200.19 |  0.076 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 3.24e-05
  0.02420 | 243 |  200.19 |  0.076 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 3.36e-05
  0.02430 | 244 |  200.19 |  0.076 |  1.112 |  1.000 |  1.000 |     2 |     23 |    1 | 1.34e-09
  0.02440 | 245 |  200.19 |  0.076 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 3.64e-05
  0.02450 | 246 |  200.19 |  0.076 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 3.75e-05
  0.02460 | 247 |  200.19 |  0.076 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 3.84e-05
  0.02470 | 248 |  200.19 |  0.076 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 3.90e-05
  0.02480 | 249 |  200.20 |  0.076 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 3.97e-05
  0.02490 | 250 |  200.20 |  0.076 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 4.04e-05
  0.02500 | 251 |  200.20 |  0.076 |  1.112 |  1.000 |  1.000 |     1 |     13 |    1 | 4.08e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.02510 | 252 |  200.20 |  0.076 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 4.13e-05
  0.02520 | 253 |  200.20 |  0.076 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 4.17e-05
  0.02530 | 254 |  200.20 |  0.076 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 4.20e-05
  0.02540 | 255 |  200.20 |  0.076 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 9.49e-05
  0.02550 | 256 |  200.20 |  0.076 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 4.26e-05
  0.02560 | 257 |  200.20 |  0.076 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 4.27e-05
  0.02570 | 258 |  200.19 |  0.076 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 4.28e-05
  0.02580 | 259 |  200.19 |  0.076 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 4.31e-05
  0.02590 | 260 |  200.19 |  0.076 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 4.32e-05
  0.02600 | 261 |  200.19 |  0.076 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 4.32e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.02610 | 262 |  200.19 |  0.076 |  1.113 |  1.000 |  1.000 |     2 |     23 |    1 | 7.46e-09
  0.02620 | 263 |  200.19 |  0.076 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 4.90e-05
  0.02630 | 264 |  200.19 |  0.076 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 4.92e-05
  0.02640 | 265 |  200.19 |  0.076 |  1.113 |  1.000 |  1.000 |     1 |     13 |    1 | 4.95e-05
  0.02650 | 266 |  200.19 |  0.076 |  1.114 |  1.000 |  1.000 |     1 |     13 |    1 | 4.99e-05
  0.02660 | 267 |  200.19 |  0.076 |  1.114 |  1.000 |  1.000 |     1 |     13 |    1 | 5.00e-05
  0.02670 | 268 |  200.19 |  0.076 |  1.114 |  1.000 |  1.000 |     1 |     13 |    1 | 5.05e-05
  0.02680 | 269 |  200.18 |  0.076 |  1.114 |  1.000 |  1.000 |     1 |     13 |    1 | 5.07e-05
  0.02690 | 270 |  200.18 |  0.076 |  1.114 |  1.000 |  1.000 |     1 |     13 |    1 | 5.05e-05
  0.02700 | 271 |  200.18 |  0.076 |  1.114 |  1.000 |  1.000 |     1 |     13 |    1 | 5.09e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.02710 | 272 |  200.18 |  0.076 |  1.114 |  1.000 |  1.000 |     1 |     13 |    1 | 5.06e-05
  0.02720 | 273 |  200.18 |  0.076 |  1.114 |  1.000 |  1.000 |     1 |     13 |    1 | 5.07e-05
  0.02730 | 274 |  200.18 |  0.076 |  1.114 |  1.000 |  1.000 |     1 |     13 |    1 | 5.06e-05
  0.02740 | 275 |  200.17 |  0.076 |  1.114 |  1.000 |  1.000 |     1 |     13 |    1 | 5.05e-05
  0.02750 | 276 |  200.17 |  0.076 |  1.114 |  1.000 |  1.000 |     1 |     13 |    1 | 5.37e-05
  0.02760 | 277 |  200.17 |  0.076 |  1.114 |  1.000 |  1.000 |     1 |     13 |    1 | 5.47e-05
  0.02770 | 278 |  200.17 |  0.076 |  1.114 |  1.000 |  1.000 |     1 |     13 |    1 | 5.44e-05
  0.02780 | 279 |  200.17 |  0.076 |  1.115 |  1.000 |  1.000 |     1 |     13 |    1 | 5.44e-05
  0.02790 | 280 |  200.17 |  0.076 |  1.115 |  1.000 |  1.000 |     1 |     13 |    1 | 5.43e-05
  0.02800 | 281 |  200.16 |  0.076 |  1.115 |  1.000 |  1.000 |     1 |     13 |    1 | 5.43e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.02810 | 282 |  200.16 |  0.076 |  1.115 |  1.000 |  1.000 |     1 |     13 |    1 | 5.42e-05
  0.02820 | 283 |  200.16 |  0.076 |  1.115 |  1.000 |  1.000 |     1 |     13 |    1 | 5.42e-05
  0.02830 | 284 |  200.16 |  0.076 |  1.115 |  1.000 |  1.000 |     1 |     13 |    1 | 5.40e-05
  0.02840 | 285 |  200.16 |  0.076 |  1.115 |  1.000 |  1.000 |     1 |     13 |    1 | 5.40e-05
  0.02850 | 286 |  200.15 |  0.076 |  1.115 |  1.000 |  1.000 |     2 |     23 |    1 | 5.92e-09
  0.02860 | 287 |  200.15 |  0.076 |  1.115 |  1.000 |  1.000 |     1 |     13 |    1 | 6.19e-05
  0.02870 | 288 |  200.15 |  0.076 |  1.115 |  1.000 |  1.000 |     1 |     13 |    1 | 6.16e-05
  0.02880 | 289 |  200.15 |  0.076 |  1.115 |  1.000 |  1.000 |     1 |     13 |    1 | 6.19e-05
  0.02890 | 290 |  200.15 |  0.076 |  1.115 |  1.000 |  1.000 |     1 |     13 |    1 | 6.14e-05
  0.02900 | 291 |  200.14 |  0.076 |  1.116 |  1.000 |  1.000 |     1 |     13 |    1 | 6.11e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.02910 | 292 |  200.14 |  0.076 |  1.116 |  1.000 |  1.000 |     1 |     13 |    1 | 6.06e-05
  0.02920 | 293 |  200.14 |  0.076 |  1.116 |  1.000 |  1.000 |     1 |     13 |    1 | 6.05e-05
  0.02930 | 294 |  200.14 |  0.076 |  1.116 |  1.000 |  1.000 |     1 |     13 |    1 | 5.98e-05
  0.02940 | 295 |  200.14 |  0.076 |  1.116 |  1.000 |  1.000 |     1 |     13 |    1 | 5.95e-05
  0.02950 | 296 |  200.13 |  0.076 |  1.116 |  1.000 |  1.000 |     1 |     13 |    1 | 6.01e-05
  0.02960 | 297 |  200.13 |  0.076 |  1.116 |  1.000 |  1.000 |     1 |     13 |    1 | 5.92e-05
  0.02970 | 298 |  200.13 |  0.076 |  1.116 |  1.000 |  1.000 |     1 |     13 |    1 | 5.87e-05
  0.02980 | 299 |  200.13 |  0.076 |  1.116 |  1.000 |  1.000 |     1 |     13 |    1 | 5.84e-05
  0.02990 | 300 |  200.12 |  0.076 |  1.116 |  1.000 |  1.000 |     1 |     13 |    1 | 5.74e-05
  0.03000 | 301 |  200.12 |  0.076 |  1.116 |  1.000 |  1.000 |     1 |     13 |    1 | 5.68e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.03010 | 302 |  200.12 |  0.076 |  1.116 |  1.000 |  1.000 |     1 |     13 |    1 | 5.62e-05
  0.03020 | 303 |  200.11 |  0.076 |  1.116 |  1.000 |  1.000 |     1 |     13 |    1 | 5.56e-05
  0.03030 | 304 |  200.11 |  0.076 |  1.117 |  1.000 |  1.000 |     1 |     13 |    1 | 5.50e-05
  0.03040 | 305 |  200.11 |  0.076 |  1.117 |  1.000 |  1.000 |     1 |     13 |    1 | 5.46e-05
  0.03050 | 306 |  200.11 |  0.076 |  1.117 |  1.000 |  1.000 |     1 |     13 |    1 | 5.41e-05
  0.03060 | 307 |  200.10 |  0.076 |  1.117 |  1.000 |  1.000 |     1 |     13 |    1 | 5.34e-05
  0.03070 | 308 |  200.10 |  0.076 |  1.117 |  1.000 |  1.000 |     1 |     13 |    1 | 5.63e-05
  0.03080 | 309 |  200.10 |  0.076 |  1.117 |  1.000 |  1.000 |     1 |     13 |    1 | 5.88e-05
  0.03090 | 310 |  200.09 |  0.076 |  1.117 |  1.000 |  1.000 |     1 |     13 |    1 | 5.80e-05
  0.03100 | 311 |  200.09 |  0.076 |  1.117 |  1.000 |  1.000 |     1 |     13 |    1 | 5.80e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.03110 | 312 |  200.09 |  0.076 |  1.117 |  1.000 |  1.000 |     1 |     13 |    1 | 5.77e-05
  0.03120 | 313 |  200.09 |  0.076 |  1.117 |  1.000 |  1.000 |     1 |     13 |    1 | 5.74e-05
  0.03130 | 314 |  200.08 |  0.076 |  1.117 |  1.000 |  1.000 |     1 |     13 |    1 | 5.71e-05
  0.03140 | 315 |  200.08 |  0.076 |  1.117 |  1.000 |  1.000 |     1 |     13 |    1 | 5.66e-05
  0.03150 | 316 |  200.08 |  0.076 |  1.118 |  1.000 |  1.000 |     1 |     13 |    1 | 5.64e-05
  0.03160 | 317 |  200.07 |  0.076 |  1.118 |  1.000 |  1.000 |     1 |     13 |    1 | 5.59e-05
  0.03170 | 318 |  200.07 |  0.076 |  1.118 |  1.000 |  1.000 |     1 |     13 |    1 | 5.57e-05
  0.03180 | 319 |  200.07 |  0.076 |  1.118 |  1.000 |  1.000 |     1 |     13 |    1 | 5.53e-05
  0.03190 | 320 |  200.07 |  0.076 |  1.118 |  1.000 |  1.000 |     1 |     13 |    1 | 5.53e-05
  0.03200 | 321 |  200.06 |  0.076 |  1.118 |  1.000 |  1.000 |     1 |     13 |    1 | 5.49e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.03210 | 322 |  200.06 |  0.076 |  1.118 |  1.000 |  1.000 |     1 |     13 |    1 | 5.47e-05
  0.03220 | 323 |  200.06 |  0.076 |  1.118 |  1.000 |  1.000 |     1 |     13 |    1 | 5.45e-05
  0.03230 | 324 |  200.05 |  0.076 |  1.118 |  1.000 |  1.000 |     1 |     13 |    1 | 5.44e-05
  0.03240 | 325 |  200.05 |  0.076 |  1.118 |  1.000 |  1.000 |     1 |     13 |    1 | 5.42e-05
  0.03250 | 326 |  200.05 |  0.076 |  1.118 |  1.000 |  1.000 |     1 |     13 |    1 | 5.40e-05
  0.03260 | 327 |  200.04 |  0.076 |  1.118 |  1.000 |  1.000 |     1 |     13 |    1 | 5.79e-05
  0.03270 | 328 |  200.04 |  0.076 |  1.118 |  1.000 |  1.000 |     1 |     13 |    1 | 5.87e-05
  0.03280 | 329 |  200.04 |  0.076 |  1.119 |  1.000 |  1.000 |     1 |     13 |    1 | 5.91e-05
  0.03290 | 330 |  200.04 |  0.076 |  1.119 |  1.000 |  1.000 |     1 |     13 |    1 | 5.89e-05
  0.03300 | 331 |  200.04 |  0.076 |  1.119 |  1.000 |  1.000 |     1 |     13 |    1 | 5.87e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.03310 | 332 |  200.03 |  0.076 |  1.119 |  1.000 |  1.000 |     1 |     13 |    1 | 5.86e-05
  0.03320 | 333 |  200.03 |  0.076 |  1.119 |  1.000 |  1.000 |     1 |     13 |    1 | 5.86e-05
  0.03330 | 334 |  200.03 |  0.076 |  1.119 |  1.000 |  1.000 |     1 |     13 |    1 | 5.85e-05
  0.03340 | 335 |  200.03 |  0.076 |  1.119 |  1.000 |  1.000 |     1 |     13 |    1 | 6.52e-05
  0.03350 | 336 |  200.03 |  0.076 |  1.119 |  1.000 |  1.000 |     1 |     13 |    1 | 6.56e-05
  0.03360 | 337 |  200.02 |  0.076 |  1.119 |  1.000 |  1.000 |     1 |     13 |    1 | 6.57e-05
  0.03370 | 338 |  200.02 |  0.076 |  1.119 |  1.000 |  1.000 |     1 |     13 |    1 | 6.56e-05
  0.03380 | 339 |  200.02 |  0.076 |  1.119 |  1.000 |  1.000 |     1 |     13 |    1 | 6.54e-05
  0.03390 | 340 |  200.02 |  0.076 |  1.119 |  1.000 |  1.000 |     1 |     13 |    1 | 6.60e-05
  0.03400 | 341 |  200.02 |  0.076 |  1.119 |  1.000 |  1.000 |     1 |     13 |    1 | 6.71e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.03410 | 342 |  200.02 |  0.076 |  1.120 |  1.000 |  1.000 |     1 |     13 |    1 | 6.67e-05
  0.03420 | 343 |  200.02 |  0.076 |  1.120 |  1.000 |  1.000 |     1 |     13 |    1 | 6.64e-05
  0.03430 | 344 |  200.02 |  0.076 |  1.120 |  1.000 |  1.000 |     1 |     13 |    1 | 6.61e-05
  0.03440 | 345 |  200.02 |  0.076 |  1.120 |  1.000 |  1.000 |     1 |     13 |    1 | 6.56e-05
  0.03450 | 346 |  200.02 |  0.076 |  1.120 |  1.000 |  1.000 |     1 |     13 |    1 | 6.50e-05
  0.03460 | 347 |  200.01 |  0.076 |  1.120 |  1.000 |  1.000 |     1 |     13 |    1 | 6.43e-05
  0.03470 | 348 |  200.01 |  0.076 |  1.120 |  1.000 |  1.000 |     1 |     13 |    1 | 6.36e-05
  0.03480 | 349 |  200.01 |  0.076 |  1.120 |  1.000 |  1.000 |     1 |     13 |    1 | 6.29e-05
  0.03490 | 350 |  200.01 |  0.076 |  1.120 |  1.000 |  1.000 |     1 |     13 |    1 | 6.24e-05
  0.03500 | 351 |  200.01 |  0.076 |  1.120 |  1.000 |  1.000 |     1 |     13 |    1 | 6.19e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.03510 | 352 |  200.01 |  0.076 |  1.120 |  1.000 |  1.000 |     1 |     13 |    1 | 6.11e-05
  0.03520 | 353 |  200.01 |  0.077 |  1.120 |  1.000 |  1.000 |     1 |     13 |    1 | 6.06e-05
  0.03530 | 354 |  200.01 |  0.077 |  1.121 |  1.000 |  1.000 |     1 |     13 |    1 | 6.00e-05
  0.03540 | 355 |  200.00 |  0.077 |  1.121 |  1.000 |  1.000 |     1 |     13 |    1 | 5.94e-05
  0.03550 | 356 |  200.00 |  0.077 |  1.121 |  1.000 |  1.000 |     1 |     13 |    1 | 5.91e-05
  0.03560 | 357 |  200.00 |  0.077 |  1.121 |  1.000 |  1.000 |     1 |     13 |    1 | 5.87e-05
  0.03570 | 358 |  200.00 |  0.077 |  1.121 |  1.000 |  1.000 |     1 |     13 |    1 | 6.26e-05
  0.03580 | 359 |  200.00 |  0.077 |  1.121 |  1.000 |  1.000 |     1 |     13 |    1 | 6.61e-05
  0.03590 | 360 |  200.00 |  0.077 |  1.121 |  1.000 |  1.000 |     1 |     13 |    1 | 6.57e-05
  0.03600 | 361 |  200.00 |  0.077 |  1.121 |  1.000 |  1.000 |     1 |     13 |    1 | 9.37e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.03610 | 362 |  200.00 |  0.077 |  1.121 |  1.000 |  1.000 |     1 |     13 |    1 | 6.51e-05
  0.03620 | 363 |  200.00 |  0.077 |  1.121 |  1.000 |  1.000 |     1 |     13 |    1 | 6.53e-05
  0.03630 | 364 |  200.00 |  0.077 |  1.121 |  1.000 |  1.000 |     1 |     13 |    1 | 6.47e-05
  0.03640 | 365 |  200.00 |  0.077 |  1.121 |  1.000 |  1.000 |     1 |     13 |    1 | 6.45e-05
  0.03650 | 366 |  200.00 |  0.077 |  1.121 |  1.000 |  1.000 |     1 |     13 |    1 | 6.43e-05
  0.03660 | 367 |  199.99 |  0.077 |  1.122 |  1.000 |  1.000 |     1 |     13 |    1 | 6.39e-05
  0.03670 | 368 |  199.99 |  0.077 |  1.122 |  1.000 |  1.000 |     2 |     23 |    1 | 9.89e-09
  0.03680 | 369 |  199.99 |  0.077 |  1.122 |  1.000 |  1.000 |     1 |     13 |    1 | 6.91e-05
  0.03690 | 370 |  199.99 |  0.077 |  1.122 |  1.000 |  1.000 |     1 |     13 |    1 | 6.86e-05
  0.03700 | 371 |  200.00 |  0.077 |  1.122 |  1.000 |  1.000 |     1 |     13 |    1 | 7.54e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.03710 | 372 |  200.00 |  0.077 |  1.122 |  1.000 |  1.000 |     1 |     13 |    1 | 7.63e-05
  0.03720 | 373 |  200.00 |  0.077 |  1.122 |  1.000 |  1.000 |     1 |     13 |    1 | 7.60e-05
  0.03730 | 374 |  200.00 |  0.077 |  1.122 |  1.000 |  1.000 |     1 |     13 |    1 | 7.56e-05
  0.03740 | 375 |  200.00 |  0.077 |  1.122 |  1.000 |  1.000 |     1 |     13 |    1 | 8.33e-05
  0.03750 | 376 |  200.01 |  0.077 |  1.122 |  1.000 |  1.000 |     1 |     13 |    1 | 8.36e-05
  0.03760 | 377 |  200.01 |  0.077 |  1.122 |  1.000 |  1.000 |     1 |     13 |    1 | 8.32e-05
  0.03770 | 378 |  200.01 |  0.077 |  1.122 |  1.000 |  1.000 |     1 |     13 |    1 | 8.26e-05
  0.03780 | 379 |  200.01 |  0.077 |  1.123 |  1.000 |  1.000 |     1 |     13 |    1 | 8.18e-05
  0.03790 | 380 |  200.02 |  0.077 |  1.123 |  1.000 |  1.000 |     1 |     13 |    1 | 8.10e-05
  0.03800 | 381 |  200.02 |  0.077 |  1.123 |  1.000 |  1.000 |     1 |     13 |    1 | 8.01e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.03810 | 382 |  200.02 |  0.077 |  1.123 |  1.000 |  1.000 |     1 |     13 |    1 | 7.88e-05
  0.03820 | 383 |  200.02 |  0.077 |  1.123 |  1.000 |  1.000 |     1 |     13 |    1 | 7.72e-05
  0.03830 | 384 |  200.02 |  0.077 |  1.123 |  1.000 |  1.000 |     1 |     13 |    1 | 7.56e-05
  0.03840 | 385 |  200.02 |  0.077 |  1.123 |  1.000 |  1.000 |     1 |     13 |    1 | 7.43e-05
  0.03850 | 386 |  200.02 |  0.077 |  1.123 |  1.000 |  1.000 |     1 |     13 |    1 | 7.27e-05
  0.03860 | 387 |  200.02 |  0.077 |  1.123 |  1.000 |  1.000 |     1 |     13 |    1 | 7.09e-05
  0.03870 | 388 |  200.03 |  0.077 |  1.123 |  1.000 |  1.000 |     1 |     13 |    1 | 6.89e-05
  0.03880 | 389 |  200.03 |  0.077 |  1.123 |  1.000 |  1.000 |     1 |     13 |    1 | 6.79e-05
  0.03890 | 390 |  200.03 |  0.077 |  1.123 |  1.000 |  1.000 |     1 |     13 |    1 | 6.62e-05
  0.03900 | 391 |  200.02 |  0.077 |  1.124 |  1.000 |  1.000 |     1 |     13 |    1 | 6.44e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.03910 | 392 |  200.02 |  0.077 |  1.124 |  1.000 |  1.000 |     1 |     13 |    1 | 6.27e-05
  0.03920 | 393 |  200.02 |  0.077 |  1.124 |  1.000 |  1.000 |     1 |     13 |    1 | 6.14e-05
  0.03930 | 394 |  200.02 |  0.077 |  1.124 |  1.000 |  1.000 |     1 |     13 |    1 | 5.98e-05
  0.03940 | 395 |  200.02 |  0.077 |  1.124 |  1.000 |  1.000 |     1 |     13 |    1 | 5.84e-05
  0.03950 | 396 |  200.02 |  0.077 |  1.124 |  1.000 |  1.000 |     1 |     13 |    1 | 5.75e-05
  0.03960 | 397 |  200.02 |  0.077 |  1.124 |  1.000 |  1.000 |     1 |     13 |    1 | 5.61e-05
  0.03970 | 398 |  200.02 |  0.077 |  1.124 |  1.000 |  1.000 |     1 |     13 |    1 | 5.48e-05
  0.03980 | 399 |  200.02 |  0.077 |  1.124 |  1.000 |  1.000 |     1 |     13 |    1 | 5.40e-05
  0.03990 | 400 |  200.01 |  0.077 |  1.124 |  1.000 |  1.000 |     1 |     14 |    1 | 5.61e-05
  0.04000 | 401 |  200.01 |  0.077 |  1.124 |  1.000 |  1.000 |     1 |     14 |    1 | 5.84e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.04010 | 402 |  200.01 |  0.077 |  1.124 |  1.000 |  1.000 |     1 |     14 |    1 | 5.77e-05
  0.04020 | 403 |  200.01 |  0.077 |  1.124 |  1.000 |  1.000 |     1 |     14 |    1 | 5.71e-05
  0.04030 | 404 |  200.01 |  0.077 |  1.125 |  1.000 |  1.000 |     1 |     14 |    1 | 5.61e-05
  0.04040 | 405 |  200.01 |  0.077 |  1.125 |  1.000 |  1.000 |     1 |     14 |    1 | 5.55e-05
  0.04050 | 406 |  200.01 |  0.077 |  1.125 |  1.000 |  1.000 |     1 |     14 |    1 | 5.62e-05
  0.04060 | 407 |  200.01 |  0.077 |  1.125 |  1.000 |  1.000 |     1 |     14 |    1 | 5.68e-05
  0.04070 | 408 |  200.01 |  0.077 |  1.125 |  1.000 |  1.000 |     1 |     14 |    1 | 5.63e-05
  0.04080 | 409 |  200.01 |  0.077 |  1.125 |  1.000 |  1.000 |     1 |     14 |    1 | 5.55e-05
  0.04090 | 410 |  200.01 |  0.077 |  1.125 |  1.000 |  1.000 |     1 |     14 |    1 | 5.51e-05
  0.04100 | 411 |  200.01 |  0.077 |  1.125 |  1.000 |  1.000 |     1 |     14 |    1 | 5.43e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.04110 | 412 |  200.01 |  0.077 |  1.125 |  1.000 |  1.000 |     1 |     14 |    1 | 5.34e-05
  0.04120 | 413 |  200.00 |  0.077 |  1.125 |  1.000 |  1.000 |     1 |     14 |    1 | 5.51e-05
  0.04130 | 414 |  200.01 |  0.077 |  1.125 |  1.000 |  1.000 |     1 |     14 |    1 | 5.80e-05
  0.04140 | 415 |  200.01 |  0.077 |  1.125 |  1.000 |  1.000 |     1 |     14 |    1 | 5.75e-05
  0.04150 | 416 |  200.01 |  0.077 |  1.125 |  1.000 |  1.000 |     1 |     14 |    1 | 5.68e-05
  0.04160 | 417 |  200.01 |  0.077 |  1.126 |  1.000 |  1.000 |     1 |     14 |    1 | 5.61e-05
  0.04170 | 418 |  200.01 |  0.077 |  1.126 |  1.000 |  1.000 |     1 |     14 |    1 | 5.56e-05
  0.04180 | 419 |  200.01 |  0.077 |  1.126 |  1.000 |  1.000 |     1 |     14 |    1 | 5.49e-05
  0.04190 | 420 |  200.01 |  0.077 |  1.126 |  1.000 |  1.000 |     1 |     14 |    1 | 5.42e-05
  0.04200 | 421 |  200.01 |  0.077 |  1.126 |  1.000 |  1.000 |     1 |     14 |    1 | 5.38e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.04210 | 422 |  200.01 |  0.077 |  1.126 |  1.000 |  1.000 |     2 |     25 |    1 | 2.73e-09
  0.04220 | 423 |  200.01 |  0.077 |  1.126 |  1.000 |  1.000 |     1 |     14 |    1 | 5.89e-05
  0.04230 | 424 |  200.01 |  0.077 |  1.126 |  1.000 |  1.000 |     1 |     14 |    1 | 5.85e-05
  0.04240 | 425 |  200.01 |  0.077 |  1.126 |  1.000 |  1.000 |     1 |     14 |    1 | 5.75e-05
  0.04250 | 426 |  200.01 |  0.077 |  1.126 |  1.000 |  1.000 |     1 |     14 |    1 | 5.71e-05
  0.04260 | 427 |  200.01 |  0.077 |  1.126 |  1.000 |  1.000 |     1 |     14 |    1 | 5.63e-05
  0.04270 | 428 |  200.01 |  0.077 |  1.126 |  1.000 |  1.000 |     1 |     14 |    1 | 5.58e-05
  0.04280 | 429 |  200.01 |  0.077 |  1.126 |  1.000 |  1.000 |     1 |     14 |    1 | 5.49e-05
  0.04290 | 430 |  200.01 |  0.077 |  1.127 |  1.000 |  1.000 |     1 |     14 |    1 | 5.41e-05
  0.04300 | 431 |  200.01 |  0.077 |  1.127 |  1.000 |  1.000 |     1 |     14 |    1 | 5.34e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.04310 | 432 |  200.01 |  0.077 |  1.127 |  1.000 |  1.000 |     1 |     14 |    1 | 5.25e-05
  0.04320 | 433 |  200.01 |  0.077 |  1.127 |  1.000 |  1.000 |     1 |     14 |    1 | 5.16e-05
  0.04330 | 434 |  200.01 |  0.077 |  1.127 |  1.000 |  1.000 |     1 |     14 |    1 | 6.15e-05
  0.04340 | 435 |  200.01 |  0.077 |  1.127 |  1.000 |  1.000 |     1 |     14 |    1 | 4.99e-05
  0.04350 | 436 |  200.01 |  0.077 |  1.127 |  1.000 |  1.000 |     1 |     14 |    1 | 4.89e-05
  0.04360 | 437 |  200.01 |  0.077 |  1.127 |  1.000 |  1.000 |     1 |     14 |    1 | 4.82e-05
  0.04370 | 438 |  200.01 |  0.077 |  1.127 |  1.000 |  1.000 |     1 |     14 |    1 | 4.72e-05
  0.04380 | 439 |  200.01 |  0.077 |  1.127 |  1.000 |  1.000 |     1 |     14 |    1 | 4.66e-05
  0.04390 | 440 |  200.01 |  0.077 |  1.127 |  1.000 |  1.000 |     1 |     14 |    1 | 4.57e-05
  0.04400 | 441 |  200.01 |  0.077 |  1.127 |  1.000 |  1.000 |     1 |     14 |    1 | 4.50e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.04410 | 442 |  200.00 |  0.077 |  1.127 |  1.000 |  1.000 |     1 |     14 |    1 | 4.44e-05
  0.04420 | 443 |  200.00 |  0.077 |  1.128 |  1.000 |  1.000 |     1 |     14 |    1 | 4.37e-05
  0.04430 | 444 |  200.00 |  0.077 |  1.128 |  1.000 |  1.000 |     1 |     14 |    1 | 4.29e-05
  0.04440 | 445 |  200.00 |  0.077 |  1.128 |  1.000 |  1.000 |     1 |     14 |    1 | 4.24e-05
  0.04450 | 446 |  200.00 |  0.077 |  1.128 |  1.000 |  1.000 |     1 |     14 |    1 | 4.20e-05
  0.04460 | 447 |  200.00 |  0.077 |  1.128 |  1.000 |  1.000 |     1 |     14 |    1 | 4.17e-05
  0.04470 | 448 |  200.00 |  0.077 |  1.128 |  1.000 |  1.000 |     1 |     14 |    1 | 4.12e-05
  0.04480 | 449 |  200.00 |  0.077 |  1.128 |  1.000 |  1.000 |     1 |     14 |    1 | 4.09e-05
  0.04490 | 450 |  199.99 |  0.077 |  1.128 |  1.000 |  1.000 |     1 |     14 |    1 | 4.06e-05
  0.04500 | 451 |  199.99 |  0.077 |  1.128 |  1.000 |  1.000 |     1 |     14 |    1 | 4.07e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.04510 | 452 |  199.99 |  0.077 |  1.128 |  1.000 |  1.000 |     1 |     14 |    1 | 4.02e-05
  0.04520 | 453 |  199.99 |  0.077 |  1.128 |  1.000 |  1.000 |     1 |     14 |    1 | 4.04e-05
  0.04530 | 454 |  199.99 |  0.077 |  1.128 |  1.000 |  1.000 |     1 |     14 |    1 | 4.02e-05
  0.04540 | 455 |  199.99 |  0.077 |  1.128 |  1.000 |  1.000 |     1 |     14 |    1 | 3.99e-05
  0.04550 | 456 |  199.98 |  0.077 |  1.128 |  1.000 |  1.000 |     1 |     14 |    1 | 4.01e-05
  0.04560 | 457 |  199.98 |  0.077 |  1.129 |  1.000 |  1.000 |     1 |     14 |    1 | 8.66e-05
  0.04570 | 458 |  199.98 |  0.077 |  1.129 |  1.000 |  1.000 |     1 |     14 |    1 | 4.19e-05
  0.04580 | 459 |  199.98 |  0.077 |  1.129 |  1.000 |  1.000 |     1 |     14 |    1 | 4.18e-05
  0.04590 | 460 |  199.98 |  0.077 |  1.129 |  1.000 |  1.000 |     1 |     14 |    1 | 4.19e-05
  0.04600 | 461 |  199.98 |  0.077 |  1.129 |  1.000 |  1.000 |     1 |     14 |    1 | 4.20e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.04610 | 462 |  199.98 |  0.077 |  1.129 |  1.000 |  1.000 |     1 |     14 |    1 | 4.22e-05
  0.04620 | 463 |  199.98 |  0.077 |  1.129 |  1.000 |  1.000 |     1 |     14 |    1 | 4.25e-05
  0.04630 | 464 |  199.98 |  0.077 |  1.129 |  1.000 |  1.000 |     1 |     14 |    1 | 4.28e-05
  0.04640 | 465 |  199.97 |  0.077 |  1.129 |  1.000 |  1.000 |     1 |     14 |    1 | 4.32e-05
  0.04650 | 466 |  199.97 |  0.077 |  1.129 |  1.000 |  1.000 |     1 |     14 |    1 | 4.33e-05
  0.04660 | 467 |  199.97 |  0.077 |  1.129 |  1.000 |  1.000 |     1 |     14 |    1 | 4.37e-05
  0.04670 | 468 |  199.97 |  0.077 |  1.129 |  1.000 |  1.000 |     1 |     14 |    1 | 4.39e-05
  0.04680 | 469 |  199.97 |  0.077 |  1.129 |  1.000 |  1.000 |     1 |     14 |    1 | 4.44e-05
  0.04690 | 470 |  199.97 |  0.077 |  1.129 |  1.000 |  1.000 |     1 |     14 |    1 | 4.45e-05
  0.04700 | 471 |  199.97 |  0.077 |  1.130 |  1.000 |  1.000 |     1 |     14 |    1 | 4.46e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.04710 | 472 |  199.97 |  0.077 |  1.130 |  1.000 |  1.000 |     2 |     25 |    1 | 6.52e-09
  0.04720 | 473 |  199.97 |  0.077 |  1.130 |  1.000 |  1.000 |     1 |     14 |    1 | 4.78e-05
  0.04730 | 474 |  199.97 |  0.077 |  1.130 |  1.000 |  1.000 |     1 |     14 |    1 | 4.80e-05
  0.04740 | 475 |  199.97 |  0.077 |  1.130 |  1.000 |  1.000 |     1 |     14 |    1 | 4.79e-05
  0.04750 | 476 |  199.97 |  0.077 |  1.130 |  1.000 |  1.000 |     1 |     14 |    1 | 4.83e-05
  0.04760 | 477 |  199.97 |  0.077 |  1.130 |  1.000 |  1.000 |     1 |     14 |    1 | 4.88e-05
  0.04770 | 478 |  199.97 |  0.077 |  1.130 |  1.000 |  1.000 |     1 |     14 |    1 | 4.90e-05
  0.04780 | 479 |  199.97 |  0.077 |  1.130 |  1.000 |  1.000 |     1 |     14 |    1 | 4.92e-05
  0.04790 | 480 |  199.97 |  0.077 |  1.130 |  1.000 |  1.000 |     1 |     14 |    1 | 4.93e-05
  0.04800 | 481 |  199.97 |  0.077 |  1.130 |  1.000 |  1.000 |     1 |     14 |    1 | 4.96e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.04810 | 482 |  199.97 |  0.077 |  1.130 |  1.000 |  1.000 |     1 |     14 |    1 | 5.01e-05
  0.04820 | 483 |  199.96 |  0.077 |  1.130 |  1.000 |  1.000 |     1 |     14 |    1 | 4.99e-05
  0.04830 | 484 |  199.96 |  0.077 |  1.131 |  1.000 |  1.000 |     1 |     14 |    1 | 5.03e-05
  0.04840 | 485 |  199.96 |  0.077 |  1.131 |  1.000 |  1.000 |     1 |     14 |    1 | 5.02e-05
  0.04850 | 486 |  199.96 |  0.077 |  1.131 |  1.000 |  1.000 |     1 |     14 |    1 | 5.06e-05
  0.04860 | 487 |  199.96 |  0.077 |  1.131 |  1.000 |  1.000 |     1 |     14 |    1 | 5.06e-05
  0.04870 | 488 |  199.96 |  0.077 |  1.131 |  1.000 |  1.000 |     1 |     14 |    1 | 5.06e-05
  0.04880 | 489 |  199.96 |  0.077 |  1.131 |  1.000 |  1.000 |     1 |     14 |    1 | 5.04e-05
  0.04890 | 490 |  199.96 |  0.077 |  1.131 |  1.000 |  1.000 |     1 |     14 |    1 | 5.06e-05
  0.04900 | 491 |  199.96 |  0.077 |  1.131 |  1.000 |  1.000 |     1 |     14 |    1 | 5.04e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.04910 | 492 |  199.96 |  0.077 |  1.131 |  1.000 |  1.000 |     1 |     14 |    1 | 5.08e-05
  0.04920 | 493 |  199.96 |  0.077 |  1.131 |  1.000 |  1.000 |     1 |     14 |    1 | 5.17e-05
  0.04930 | 494 |  199.96 |  0.077 |  1.131 |  1.000 |  1.000 |     1 |     14 |    1 | 5.14e-05
  0.04940 | 495 |  199.96 |  0.077 |  1.131 |  1.000 |  1.000 |     1 |     14 |    1 | 5.15e-05
  0.04950 | 496 |  199.96 |  0.077 |  1.131 |  1.000 |  1.000 |     1 |     14 |    1 | 5.57e-05
  0.04960 | 497 |  199.96 |  0.077 |  1.131 |  1.000 |  1.000 |     1 |     14 |    1 | 5.71e-05
  0.04970 | 498 |  199.96 |  0.077 |  1.132 |  1.000 |  1.000 |     1 |     14 |    1 | 5.71e-05
  0.04980 | 499 |  199.96 |  0.077 |  1.132 |  1.000 |  1.000 |     1 |     14 |    1 | 5.76e-05
  0.04990 | 500 |  199.96 |  0.077 |  1.132 |  1.000 |  1.000 |     1 |     14 |    1 | 5.74e-05
  0.05000 | 501 |  199.96 |  0.077 |  1.132 |  1.000 |  1.000 |     1 |     14 |    1 | 5.79e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.05010 | 502 |  199.96 |  0.077 |  1.132 |  1.000 |  1.000 |     1 |     14 |    1 | 5.81e-05
  0.05020 | 503 |  199.96 |  0.077 |  1.132 |  1.000 |  1.000 |     1 |     14 |    1 | 5.83e-05
  0.05030 | 504 |  199.96 |  0.077 |  1.132 |  1.000 |  1.000 |     1 |     14 |    1 | 5.81e-05
  0.05040 | 505 |  199.96 |  0.077 |  1.132 |  1.000 |  1.000 |     1 |     14 |    1 | 5.79e-05
  0.05050 | 506 |  199.96 |  0.077 |  1.132 |  1.000 |  1.000 |     1 |     14 |    1 | 5.76e-05
  0.05060 | 507 |  199.96 |  0.077 |  1.132 |  1.000 |  1.000 |     1 |     14 |    1 | 5.71e-05
  0.05070 | 508 |  199.96 |  0.077 |  1.132 |  1.000 |  1.000 |     1 |     14 |    1 | 5.63e-05
  0.05080 | 509 |  199.96 |  0.077 |  1.132 |  1.000 |  1.000 |     1 |     14 |    1 | 5.50e-05
  0.05090 | 510 |  199.96 |  0.077 |  1.132 |  1.000 |  1.000 |     1 |     14 |    1 | 5.34e-05
  0.05100 | 511 |  199.96 |  0.077 |  1.133 |  1.000 |  1.000 |     1 |     14 |    1 | 5.12e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.05110 | 512 |  199.96 |  0.077 |  1.133 |  1.000 |  1.000 |     1 |     14 |    1 | 4.88e-05
  0.05120 | 513 |  199.96 |  0.077 |  1.133 |  1.000 |  1.000 |     1 |     14 |    1 | 5.16e-05
  0.05130 | 514 |  199.96 |  0.077 |  1.133 |  1.000 |  1.000 |     1 |     14 |    1 | 4.93e-05
  0.05140 | 515 |  199.96 |  0.077 |  1.133 |  1.000 |  1.000 |     1 |     14 |    1 | 4.61e-05
  0.05150 | 516 |  199.96 |  0.077 |  1.133 |  1.000 |  1.000 |     1 |     14 |    1 | 4.28e-05
  0.05160 | 517 |  199.96 |  0.077 |  1.133 |  1.000 |  1.000 |     1 |     14 |    1 | 3.94e-05
  0.05170 | 518 |  199.95 |  0.077 |  1.133 |  1.000 |  1.000 |     1 |     14 |    1 | 3.60e-05
  0.05180 | 519 |  199.95 |  0.077 |  1.133 |  1.000 |  1.000 |     1 |     14 |    1 | 3.33e-05
  0.05190 | 520 |  199.95 |  0.077 |  1.133 |  1.000 |  1.000 |     1 |     14 |    1 | 2.99e-05
  0.05200 | 521 |  199.95 |  0.077 |  1.133 |  1.000 |  1.000 |     1 |     14 |    1 | 2.67e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.05210 | 522 |  199.94 |  0.077 |  1.133 |  1.000 |  1.000 |     1 |     14 |    1 | 2.37e-05
  0.05220 | 523 |  199.94 |  0.077 |  1.133 |  1.000 |  1.000 |     1 |     14 |    1 | 2.09e-05
  0.05230 | 524 |  199.94 |  0.077 |  1.133 |  1.000 |  1.000 |     1 |     14 |    1 | 1.79e-05
  0.05240 | 525 |  199.93 |  0.077 |  1.133 |  1.000 |  1.000 |     1 |     14 |    1 | 1.55e-05
  0.05250 | 526 |  199.93 |  0.077 |  1.133 |  1.000 |  1.000 |     1 |     14 |    1 | 1.30e-05
  0.05260 | 527 |  199.93 |  0.077 |  1.133 |  1.000 |  1.000 |     1 |     14 |    1 | 1.10e-05
  0.05270 | 528 |  199.92 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 9.15e-06
  0.05280 | 529 |  199.92 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 7.34e-06
  0.05290 | 530 |  199.91 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 5.99e-06
  0.05300 | 531 |  199.91 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 4.55e-06

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.05310 | 532 |  199.90 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 3.57e-06
  0.05320 | 533 |  199.90 |  0.077 |  1.134 |  1.000 |  1.000 |     2 |     25 |    1 | 2.45e-09
  0.05330 | 534 |  199.89 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 2.52e-06
  0.05340 | 535 |  199.89 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 1.79e-06
  0.05350 | 536 |  199.88 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 1.23e-06
  0.05360 | 537 |  199.88 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 6.20e-07
  0.05370 | 538 |  199.87 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 3.43e-07
  0.05380 | 539 |  199.87 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 2.20e-07
  0.05390 | 540 |  199.87 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 1.94e-07
  0.05400 | 541 |  199.86 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 3.75e-07

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.05410 | 542 |  199.86 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 5.86e-07
  0.05420 | 543 |  199.85 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 6.18e-07
  0.05430 | 544 |  199.85 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 7.57e-07
  0.05440 | 545 |  199.84 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 6.54e-07
  0.05450 | 546 |  199.84 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 7.21e-07
  0.05460 | 547 |  199.84 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 6.60e-07
  0.05470 | 548 |  199.83 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 6.76e-07
  0.05480 | 549 |  199.83 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 6.07e-07
  0.05490 | 550 |  199.83 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 6.08e-07
  0.05500 | 551 |  199.82 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 5.06e-07

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.05510 | 552 |  199.82 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 4.93e-07
  0.05520 | 553 |  199.82 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 4.17e-07
  0.05530 | 554 |  199.82 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 3.54e-07
  0.05540 | 555 |  199.81 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 2.59e-07
  0.05550 | 556 |  199.81 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 2.51e-07
  0.05560 | 557 |  199.81 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 2.08e-07
  0.05570 | 558 |  199.81 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 1.59e-07
  0.05580 | 559 |  199.81 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 8.35e-08
  0.05590 | 560 |  199.80 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 4.96e-08
  0.05600 | 561 |  199.80 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 2.91e-08

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.05610 | 562 |  199.80 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 2.34e-08
  0.05620 | 563 |  199.80 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 6.45e-08
  0.05630 | 564 |  199.80 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 9.30e-08
  0.05640 | 565 |  199.80 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 1.33e-07
  0.05650 | 566 |  199.80 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 1.33e-07
  0.05660 | 567 |  199.80 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 1.84e-07
  0.05670 | 568 |  199.80 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 1.85e-07
  0.05680 | 569 |  199.80 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 2.10e-07
  0.05690 | 570 |  199.80 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 2.25e-07
  0.05700 | 571 |  199.80 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 2.31e-07

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.05710 | 572 |  199.80 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 2.65e-07
  0.05720 | 573 |  199.80 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 2.84e-07
  0.05730 | 574 |  199.80 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 2.92e-07
  0.05740 | 575 |  199.80 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 2.98e-07
  0.05750 | 576 |  199.80 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 3.05e-07
  0.05760 | 577 |  199.80 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 3.05e-07
  0.05770 | 578 |  199.80 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 3.19e-07
  0.05780 | 579 |  199.81 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 3.15e-07
  0.05790 | 580 |  199.81 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 3.17e-07
  0.05800 | 581 |  199.81 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 3.17e-07

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.05810 | 582 |  199.81 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 3.15e-07
  0.05820 | 583 |  199.81 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 3.23e-07
  0.05830 | 584 |  199.81 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 3.19e-07
  0.05840 | 585 |  199.81 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 3.05e-07
  0.05850 | 586 |  199.81 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 6.30e-07
  0.05860 | 587 |  199.82 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 2.86e-07
  0.05870 | 588 |  199.82 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 2.87e-07
  0.05880 | 589 |  199.82 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 2.92e-07
  0.05890 | 590 |  199.82 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 2.69e-07
  0.05900 | 591 |  199.82 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 2.55e-07

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.05910 | 592 |  199.82 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 2.45e-07
  0.05920 | 593 |  199.82 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 2.33e-07
  0.05930 | 594 |  199.83 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 2.22e-07
  0.05940 | 595 |  199.83 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     13 |    1 | 2.11e-07
  0.05950 | 596 |  199.83 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 2.01e-07
  0.05960 | 597 |  199.83 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 1.87e-07
  0.05970 | 598 |  199.83 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 1.66e-07
  0.05980 | 599 |  199.84 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 1.65e-07
  0.05990 | 600 |  199.84 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 1.44e-07
  0.06000 | 601 |  199.84 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 1.27e-07

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.06010 | 602 |  199.84 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 1.16e-07
  0.06020 | 603 |  199.84 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 1.04e-07
  0.06030 | 604 |  199.84 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 9.27e-08
  0.06040 | 605 |  199.85 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 3.90e-08
  0.06050 | 606 |  199.85 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 8.42e-08
  0.06060 | 607 |  199.85 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 3.01e-07
  0.06070 | 608 |  199.84 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 4.00e-07
  0.06080 | 609 |  199.84 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 1.13e-07
  0.06090 | 610 |  199.83 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 2.12e-06
  0.06100 | 611 |  199.82 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 6.84e-06

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.06110 | 612 |  199.80 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 1.58e-05
  0.06120 | 613 |  199.77 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 3.08e-05
  0.06130 | 614 |  199.74 |  0.077 |  1.134 |  1.000 |  1.000 |     1 |     14 |    1 | 5.36e-05
  0.06140 | 615 |  199.69 |  0.077 |  1.134 |  1.000 |  1.000 |     2 |     26 |    1 | 1.93e-09
  0.06150 | 616 |  199.64 |  0.077 |  1.134 |  1.000 |  1.000 |     2 |     24 |    1 | 3.66e-11
  0.06160 | 617 |  199.58 |  0.077 |  1.134 |  1.000 |  1.000 |     2 |     24 |    1 | 1.31e-11
  0.06170 | 618 |  199.52 |  0.077 |  1.134 |  1.000 |  1.000 |     2 |     24 |    1 | 1.52e-10
  0.06180 | 619 |  199.44 |  0.077 |  1.134 |  1.000 |  1.000 |     2 |     25 |    1 | 1.96e-10
  0.06190 | 620 |  199.36 |  0.078 |  1.134 |  1.000 |  1.000 |     2 |     25 |    1 | 1.48e-10
  0.06200 | 621 |  199.27 |  0.078 |  1.134 |  1.000 |  1.000 |     2 |     25 |    1 | 3.09e-10

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.06210 | 622 |  199.17 |  0.078 |  1.133 |  1.000 |  1.000 |     2 |     25 |    1 | 2.75e-10
  0.06220 | 623 |  199.06 |  0.078 |  1.133 |  1.000 |  1.000 |     2 |     25 |    1 | 4.51e-10
  0.06230 | 624 |  198.95 |  0.078 |  1.133 |  1.000 |  1.000 |     2 |     26 |    1 | 8.29e-09
  0.06240 | 625 |  198.83 |  0.078 |  1.133 |  1.000 |  1.000 |     2 |     25 |    1 | 7.78e-10
  0.06250 | 626 |  198.70 |  0.078 |  1.133 |  1.000 |  1.000 |     2 |     25 |    1 | 5.50e-10
  0.06260 | 627 |  198.56 |  0.078 |  1.133 |  1.000 |  1.000 |     2 |     26 |    1 | 2.70e-10
  0.06270 | 628 |  198.42 |  0.078 |  1.133 |  1.000 |  1.000 |     2 |     27 |    1 | 4.01e-09
  0.06280 | 629 |  198.28 |  0.078 |  1.132 |  1.000 |  1.000 |     2 |     25 |    1 | 1.45e-09
  0.06290 | 630 |  198.13 |  0.078 |  1.132 |  1.000 |  1.000 |     2 |     25 |    1 | 1.19e-09
  0.06300 | 631 |  197.97 |  0.078 |  1.132 |  1.000 |  1.000 |     2 |     25 |    1 | 1.34e-09

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.06310 | 632 |  197.82 |  0.078 |  1.132 |  1.000 |  1.000 |     2 |     26 |    1 | 1.86e-08
  0.06320 | 633 |  197.66 |  0.078 |  1.132 |  1.000 |  1.000 |     2 |     26 |    1 | 1.39e-08
  0.06330 | 634 |  197.49 |  0.079 |  1.132 |  1.000 |  1.000 |     2 |     25 |    1 | 1.29e-10
  0.06340 | 635 |  197.33 |  0.079 |  1.131 |  1.000 |  1.000 |     2 |     25 |    1 | 1.30e-09
  0.06350 | 636 |  197.15 |  0.079 |  1.131 |  1.000 |  1.000 |     2 |     25 |    1 | 9.52e-10
  0.06360 | 637 |  196.98 |  0.079 |  1.131 |  1.000 |  1.000 |     2 |     25 |    1 | 2.45e-09
  0.06370 | 638 |  196.80 |  0.079 |  1.131 |  1.000 |  1.000 |     2 |     25 |    1 | 3.80e-09
  0.06380 | 639 |  196.62 |  0.079 |  1.130 |  1.000 |  1.000 |     2 |     25 |    1 | 4.48e-09
  0.06390 | 640 |  196.44 |  0.079 |  1.130 |  1.000 |  1.000 |     2 |     25 |    1 | 3.95e-09
  0.06400 | 641 |  196.25 |  0.079 |  1.130 |  1.000 |  1.000 |     2 |     25 |    1 | 6.55e-09

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.06410 | 642 |  196.06 |  0.079 |  1.130 |  1.000 |  1.000 |     2 |     25 |    1 | 7.30e-09
  0.06420 | 643 |  195.87 |  0.079 |  1.129 |  1.000 |  1.000 |     2 |     25 |    1 | 7.41e-09
  0.06430 | 644 |  195.68 |  0.079 |  1.129 |  1.000 |  1.000 |     2 |     25 |    1 | 9.22e-09
  0.06440 | 645 |  195.48 |  0.079 |  1.129 |  1.000 |  1.000 |     2 |     25 |    1 | 1.16e-08
  0.06450 | 646 |  195.28 |  0.079 |  1.129 |  1.000 |  1.000 |     2 |     25 |    1 | 1.00e-08
  0.06460 | 647 |  195.08 |  0.079 |  1.128 |  1.000 |  1.000 |     2 |     25 |    1 | 3.26e-09
  0.06470 | 648 |  194.88 |  0.079 |  1.128 |  1.000 |  1.000 |     2 |     25 |    1 | 1.41e-08
  0.06480 | 649 |  194.68 |  0.079 |  1.128 |  1.000 |  1.000 |     2 |     26 |    1 | 7.69e-08
  0.06490 | 650 |  194.48 |  0.079 |  1.128 |  1.000 |  1.000 |     2 |     25 |    1 | 1.23e-08
  0.06500 | 651 |  194.29 |  0.079 |  1.127 |  1.000 |  1.000 |     2 |     25 |    1 | 1.36e-08

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.06510 | 652 |  194.09 |  0.079 |  1.127 |  1.000 |  1.000 |     2 |     25 |    1 | 1.32e-08
  0.06520 | 653 |  193.89 |  0.079 |  1.127 |  1.000 |  1.000 |     2 |     25 |    1 | 1.37e-08
  0.06530 | 654 |  193.69 |  0.079 |  1.127 |  1.000 |  1.000 |     2 |     25 |    1 | 1.55e-08
  0.06540 | 655 |  193.49 |  0.079 |  1.126 |  1.000 |  1.000 |     2 |     25 |    1 | 1.25e-08
  0.06550 | 656 |  193.29 |  0.078 |  1.126 |  1.000 |  1.000 |     2 |     25 |    1 | 1.67e-08
  0.06560 | 657 |  193.09 |  0.078 |  1.126 |  1.000 |  1.000 |     2 |     25 |    1 | 1.54e-08
  0.06570 | 658 |  192.89 |  0.078 |  1.126 |  1.000 |  1.000 |     2 |     25 |    1 | 1.63e-08
  0.06580 | 659 |  192.69 |  0.078 |  1.126 |  1.000 |  1.000 |     2 |     25 |    1 | 1.68e-08
  0.06590 | 660 |  192.49 |  0.078 |  1.125 |  1.000 |  1.000 |     2 |     25 |    1 | 6.93e-08
  0.06600 | 661 |  192.31 |  0.078 |  1.125 |  1.000 |  1.000 |     2 |     25 |    1 | 1.20e-08

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.06610 | 662 |  192.12 |  0.078 |  1.125 |  1.000 |  1.000 |     2 |     25 |    1 | 1.26e-08
  0.06620 | 663 |  191.94 |  0.078 |  1.125 |  1.000 |  1.000 |     2 |     25 |    1 | 1.31e-08
  0.06630 | 664 |  191.75 |  0.078 |  1.125 |  1.000 |  1.000 |     2 |     25 |    1 | 3.86e-09
  0.06640 | 665 |  191.57 |  0.078 |  1.124 |  1.000 |  1.000 |     2 |     25 |    1 | 6.49e-08
  0.06650 | 666 |  191.39 |  0.078 |  1.124 |  1.000 |  1.000 |     2 |     25 |    1 | 1.22e-08
  0.06660 | 667 |  191.20 |  0.078 |  1.124 |  1.000 |  1.000 |     2 |     25 |    1 | 1.28e-08
  0.06670 | 668 |  191.02 |  0.078 |  1.124 |  1.000 |  1.000 |     2 |     25 |    1 | 1.29e-08
  0.06680 | 669 |  190.83 |  0.078 |  1.124 |  1.000 |  1.000 |     2 |     25 |    1 | 1.24e-08
  0.06690 | 670 |  190.64 |  0.078 |  1.124 |  1.000 |  1.000 |     2 |     24 |    1 | 1.23e-08
  0.06700 | 671 |  190.45 |  0.078 |  1.124 |  1.000 |  1.000 |     2 |     24 |    1 | 1.12e-08

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.06710 | 672 |  190.26 |  0.077 |  1.123 |  1.000 |  1.000 |     2 |     24 |    1 | 1.29e-08
  0.06720 | 673 |  190.07 |  0.077 |  1.123 |  1.000 |  1.000 |     2 |     24 |    1 | 1.33e-08
  0.06730 | 674 |  189.88 |  0.077 |  1.123 |  1.000 |  1.000 |     2 |     25 |    1 | 1.29e-08
  0.06740 | 675 |  189.69 |  0.077 |  1.123 |  1.000 |  1.000 |     2 |     25 |    1 | 1.21e-08
  0.06750 | 676 |  189.49 |  0.077 |  1.123 |  1.000 |  1.000 |     2 |     23 |    1 | 1.10e-08
  0.06760 | 677 |  189.30 |  0.077 |  1.123 |  1.000 |  1.000 |     2 |     23 |    1 | 1.19e-08
  0.06770 | 678 |  189.11 |  0.077 |  1.123 |  1.000 |  1.000 |     2 |     25 |    1 | 1.16e-08
  0.06780 | 679 |  188.91 |  0.077 |  1.122 |  1.000 |  1.000 |     2 |     25 |    1 | 1.21e-08
  0.06790 | 680 |  188.71 |  0.077 |  1.122 |  1.000 |  1.000 |     2 |     25 |    1 | 1.32e-08
  0.06800 | 681 |  188.52 |  0.077 |  1.122 |  1.000 |  1.000 |     2 |     23 |    1 | 1.21e-08

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.06810 | 682 |  188.32 |  0.077 |  1.122 |  1.000 |  1.000 |     2 |     23 |    1 | 1.02e-08
  0.06820 | 683 |  188.12 |  0.077 |  1.122 |  1.000 |  1.000 |     2 |     24 |    1 | 1.08e-08
  0.06830 | 684 |  187.92 |  0.077 |  1.122 |  1.000 |  1.000 |     2 |     23 |    1 | 1.00e-08
  0.06840 | 685 |  187.72 |  0.077 |  1.122 |  1.000 |  1.000 |     2 |     23 |    1 | 1.03e-08
  0.06850 | 686 |  187.52 |  0.077 |  1.122 |  1.000 |  1.000 |     2 |     23 |    1 | 9.95e-09
  0.06860 | 687 |  187.32 |  0.077 |  1.121 |  1.000 |  1.000 |     2 |     23 |    1 | 9.68e-09
  0.06870 | 688 |  187.12 |  0.077 |  1.121 |  1.000 |  1.000 |     2 |     23 |    1 | 9.74e-09
  0.06880 | 689 |  186.92 |  0.077 |  1.121 |  1.000 |  1.000 |     2 |     23 |    1 | 9.31e-09
  0.06890 | 690 |  186.72 |  0.077 |  1.121 |  1.000 |  1.000 |     2 |     23 |    1 | 9.07e-09
  0.06900 | 691 |  186.52 |  0.076 |  1.121 |  1.000 |  1.000 |     2 |     23 |    1 | 8.97e-09

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.06910 | 692 |  186.32 |  0.076 |  1.121 |  1.000 |  1.000 |     2 |     23 |    1 | 9.23e-09
  0.06920 | 693 |  186.12 |  0.076 |  1.121 |  1.000 |  1.000 |     2 |     23 |    1 | 9.67e-09
  0.06930 | 694 |  185.92 |  0.076 |  1.121 |  1.000 |  1.000 |     2 |     24 |    1 | 7.82e-09
  0.06940 | 695 |  185.72 |  0.076 |  1.120 |  1.000 |  1.000 |     2 |     24 |    1 | 8.49e-09
  0.06950 | 696 |  185.51 |  0.076 |  1.120 |  1.000 |  1.000 |     2 |     23 |    1 | 5.92e-09
  0.06960 | 697 |  185.31 |  0.076 |  1.120 |  1.000 |  1.000 |     2 |     25 |    1 | 1.63e-07
  0.06970 | 698 |  185.11 |  0.076 |  1.120 |  1.000 |  1.000 |     2 |     23 |    1 | 1.00e-08
  0.06980 | 699 |  184.91 |  0.076 |  1.120 |  1.000 |  1.000 |     2 |     23 |    1 | 1.15e-08
  0.06990 | 700 |  184.71 |  0.076 |  1.120 |  1.000 |  1.000 |     2 |     25 |    1 | 9.25e-09
  0.07000 | 701 |  184.51 |  0.076 |  1.120 |  1.000 |  1.000 |     2 |     25 |    1 | 1.11e-08

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.07010 | 702 |  184.38 |  0.076 |  1.120 |  1.005 |  1.000 |     2 |     25 |    1 | 3.01e-08
  0.07020 | 703 |  184.25 |  0.076 |  1.120 |  1.010 |  1.000 |     2 |     25 |    1 | 1.76e-09
  0.07030 | 704 |  184.13 |  0.076 |  1.120 |  1.015 |  1.000 |     2 |     25 |    1 | 4.60e-10
  0.07040 | 705 |  184.00 |  0.076 |  1.121 |  1.020 |  1.000 |     2 |     25 |    1 | 1.32e-09
  0.07050 | 706 |  183.87 |  0.076 |  1.121 |  1.025 |  1.000 |     2 |     25 |    1 | 9.01e-10
  0.07060 | 707 |  183.74 |  0.076 |  1.121 |  1.030 |  1.000 |     2 |     25 |    1 | 3.81e-10
  0.07070 | 708 |  183.61 |  0.076 |  1.121 |  1.035 |  1.000 |     2 |     25 |    1 | 1.18e-09
  0.07080 | 709 |  183.48 |  0.076 |  1.122 |  1.040 |  1.000 |     2 |     25 |    1 | 1.57e-09
  0.07090 | 710 |  183.36 |  0.076 |  1.122 |  1.045 |  1.000 |     2 |     25 |    1 | 7.60e-10
  0.07100 | 711 |  183.23 |  0.076 |  1.122 |  1.050 |  1.000 |     2 |     25 |    1 | 1.25e-09

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.07110 | 712 |  183.11 |  0.076 |  1.122 |  1.055 |  1.000 |     2 |     25 |    1 | 2.33e-10
  0.07120 | 713 |  182.99 |  0.076 |  1.123 |  1.060 |  1.000 |     2 |     25 |    1 | 1.14e-09
  0.07130 | 714 |  182.88 |  0.076 |  1.123 |  1.065 |  1.000 |     2 |     25 |    1 | 2.21e-09
  0.07140 | 715 |  182.77 |  0.076 |  1.123 |  1.070 |  1.000 |     2 |     25 |    1 | 2.22e-09
  0.07150 | 716 |  182.66 |  0.075 |  1.123 |  1.075 |  1.000 |     2 |     25 |    1 | 2.37e-10
  0.07160 | 717 |  182.55 |  0.075 |  1.123 |  1.080 |  1.000 |     2 |     25 |    1 | 2.22e-09
  0.07170 | 718 |  182.45 |  0.075 |  1.124 |  1.085 |  1.000 |     2 |     25 |    1 | 1.31e-09
  0.07180 | 719 |  182.36 |  0.075 |  1.124 |  1.090 |  1.000 |     2 |     25 |    1 | 1.59e-09
  0.07190 | 720 |  182.26 |  0.075 |  1.124 |  1.095 |  1.000 |     1 |     14 |    1 | 8.65e-05
  0.07200 | 721 |  182.17 |  0.075 |  1.125 |  1.100 |  1.000 |     2 |     25 |    1 | 1.74e-08

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.07210 | 722 |  182.09 |  0.075 |  1.125 |  1.105 |  1.000 |     1 |     14 |    1 | 6.32e-05
  0.07220 | 723 |  182.01 |  0.075 |  1.125 |  1.110 |  1.000 |     2 |     25 |    1 | 5.49e-10
  0.07230 | 724 |  181.93 |  0.075 |  1.125 |  1.115 |  1.000 |     2 |     25 |    1 | 2.63e-09
  0.07240 | 725 |  181.85 |  0.075 |  1.126 |  1.120 |  1.000 |     2 |     25 |    1 | 2.54e-09
  0.07250 | 726 |  181.78 |  0.075 |  1.126 |  1.125 |  1.000 |     2 |     25 |    1 | 1.16e-09
  0.07260 | 727 |  181.71 |  0.075 |  1.126 |  1.130 |  1.000 |     2 |     25 |    1 | 1.59e-09
  0.07270 | 728 |  181.64 |  0.075 |  1.126 |  1.135 |  1.000 |     2 |     25 |    1 | 8.64e-10
  0.07280 | 729 |  181.58 |  0.075 |  1.127 |  1.140 |  1.000 |     2 |     25 |    1 | 3.26e-09
  0.07290 | 730 |  181.52 |  0.075 |  1.127 |  1.145 |  1.000 |     2 |     25 |    1 | 7.02e-10
  0.07300 | 731 |  181.46 |  0.075 |  1.127 |  1.150 |  1.000 |     2 |     25 |    1 | 1.10e-09

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.07310 | 732 |  181.41 |  0.075 |  1.127 |  1.155 |  1.000 |     2 |     25 |    1 | 1.98e-09
  0.07320 | 733 |  181.36 |  0.075 |  1.128 |  1.160 |  1.000 |     2 |     25 |    1 | 3.38e-09
  0.07330 | 734 |  181.31 |  0.075 |  1.128 |  1.165 |  1.000 |     2 |     25 |    1 | 1.20e-09
  0.07340 | 735 |  181.26 |  0.075 |  1.128 |  1.170 |  1.000 |     2 |     25 |    1 | 3.50e-08
  0.07350 | 736 |  181.21 |  0.075 |  1.128 |  1.175 |  1.000 |     2 |     25 |    1 | 1.16e-10
  0.07360 | 737 |  181.17 |  0.075 |  1.129 |  1.180 |  1.000 |     2 |     25 |    1 | 3.05e-09
  0.07370 | 738 |  181.12 |  0.075 |  1.129 |  1.185 |  1.000 |     2 |     25 |    1 | 1.31e-09
  0.07380 | 739 |  181.08 |  0.075 |  1.129 |  1.190 |  1.000 |     2 |     25 |    1 | 2.72e-10
  0.07390 | 740 |  181.04 |  0.075 |  1.130 |  1.195 |  1.000 |     2 |     25 |    1 | 1.05e-09
  0.07400 | 741 |  181.01 |  0.075 |  1.130 |  1.200 |  1.000 |     2 |     25 |    1 | 1.51e-09

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.07410 | 742 |  180.97 |  0.075 |  1.130 |  1.205 |  1.000 |     2 |     25 |    1 | 2.08e-10
  0.07420 | 743 |  180.93 |  0.075 |  1.130 |  1.210 |  1.000 |     2 |     25 |    1 | 3.23e-09
  0.07430 | 744 |  180.90 |  0.075 |  1.131 |  1.215 |  1.000 |     2 |     25 |    1 | 1.24e-09
  0.07440 | 745 |  180.87 |  0.075 |  1.131 |  1.220 |  1.000 |     2 |     25 |    1 | 3.14e-09
  0.07450 | 746 |  180.84 |  0.075 |  1.131 |  1.225 |  1.000 |     2 |     25 |    1 | 1.09e-09
  0.07460 | 747 |  180.81 |  0.075 |  1.131 |  1.230 |  1.000 |     2 |     26 |    1 | 3.31e-09
  0.07470 | 748 |  180.78 |  0.075 |  1.132 |  1.235 |  1.000 |     2 |     25 |    1 | 1.26e-10
  0.07480 | 749 |  180.75 |  0.075 |  1.132 |  1.240 |  1.000 |     2 |     26 |    1 | 1.94e-09
  0.07490 | 750 |  180.72 |  0.075 |  1.132 |  1.245 |  1.000 |     2 |     26 |    1 | 8.60e-10
  0.07500 | 751 |  180.70 |  0.075 |  1.132 |  1.250 |  1.000 |     2 |     27 |    1 | 1.22e-09

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.07510 | 752 |  180.67 |  0.075 |  1.133 |  1.255 |  1.000 |     2 |     27 |    1 | 1.02e-09
  0.07520 | 753 |  180.65 |  0.075 |  1.133 |  1.260 |  1.000 |     2 |     27 |    1 | 3.12e-09
  0.07530 | 754 |  180.63 |  0.075 |  1.133 |  1.265 |  1.000 |     2 |     27 |    1 | 6.97e-09
  0.07540 | 755 |  180.60 |  0.075 |  1.133 |  1.270 |  1.000 |     2 |     27 |    1 | 4.03e-08
  0.07550 | 756 |  180.58 |  0.075 |  1.134 |  1.275 |  1.000 |     2 |     27 |    1 | 8.57e-08
  0.07560 | 757 |  180.56 |  0.075 |  1.134 |  1.280 |  1.000 |     2 |     28 |    1 | 1.59e-07
  0.07570 | 758 |  180.54 |  0.075 |  1.134 |  1.285 |  1.000 |     2 |     27 |    1 | 2.12e-09
  0.07580 | 759 |  180.52 |  0.075 |  1.134 |  1.290 |  1.000 |     2 |     27 |    1 | 9.89e-10
  0.07590 | 760 |  180.50 |  0.075 |  1.134 |  1.295 |  1.000 |     2 |     27 |    1 | 1.35e-09
  0.07600 | 761 |  180.48 |  0.075 |  1.135 |  1.300 |  1.000 |     2 |     27 |    1 | 1.88e-09

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.07610 | 762 |  180.46 |  0.076 |  1.135 |  1.305 |  1.000 |     2 |     27 |    1 | 7.86e-10
  0.07620 | 763 |  180.45 |  0.076 |  1.135 |  1.310 |  1.000 |     2 |     27 |    1 | 1.25e-09
  0.07630 | 764 |  180.43 |  0.076 |  1.135 |  1.315 |  1.000 |     2 |     27 |    1 | 2.09e-09
  0.07640 | 765 |  180.41 |  0.076 |  1.135 |  1.320 |  1.000 |     2 |     27 |    1 | 2.27e-09
  0.07650 | 766 |  180.40 |  0.076 |  1.136 |  1.325 |  1.000 |     2 |     28 |    1 | 9.75e-08
  0.07660 | 767 |  180.38 |  0.076 |  1.136 |  1.330 |  1.000 |     2 |     28 |    1 | 4.80e-08
  0.07670 | 768 |  180.37 |  0.076 |  1.136 |  1.335 |  1.000 |     2 |     27 |    1 | 5.16e-09
  0.07680 | 769 |  180.36 |  0.076 |  1.136 |  1.340 |  1.000 |     2 |     27 |    1 | 3.47e-09
  0.07690 | 770 |  180.34 |  0.076 |  1.136 |  1.345 |  1.000 |     2 |     27 |    1 | 2.27e-09
  0.07700 | 771 |  180.33 |  0.076 |  1.137 |  1.350 |  1.000 |     2 |     27 |    1 | 2.00e-09

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.07710 | 772 |  180.32 |  0.076 |  1.137 |  1.355 |  1.000 |     2 |     27 |    1 | 6.98e-10
  0.07720 | 773 |  180.31 |  0.076 |  1.137 |  1.360 |  1.000 |     2 |     27 |    1 | 9.31e-10
  0.07730 | 774 |  180.29 |  0.076 |  1.137 |  1.365 |  1.000 |     2 |     27 |    1 | 3.85e-09
  0.07740 | 775 |  180.28 |  0.076 |  1.137 |  1.370 |  1.000 |     2 |     27 |    1 | 1.36e-09
  0.07750 | 776 |  180.27 |  0.076 |  1.137 |  1.375 |  1.000 |     2 |     27 |    1 | 3.18e-09
  0.07760 | 777 |  180.26 |  0.076 |  1.138 |  1.380 |  1.000 |     2 |     27 |    1 | 3.46e-09
  0.07770 | 778 |  180.25 |  0.076 |  1.138 |  1.385 |  1.000 |     2 |     27 |    1 | 1.54e-09
  0.07780 | 779 |  180.24 |  0.076 |  1.138 |  1.390 |  1.000 |     2 |     27 |    1 | 1.48e-09
  0.07790 | 780 |  180.23 |  0.076 |  1.138 |  1.395 |  1.000 |     2 |     27 |    1 | 2.49e-09
  0.07800 | 781 |  180.22 |  0.076 |  1.138 |  1.400 |  1.000 |     2 |     27 |    1 | 3.42e-09

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.07810 | 782 |  180.21 |  0.076 |  1.138 |  1.405 |  1.000 |     2 |     27 |    1 | 2.18e-09
  0.07820 | 783 |  180.20 |  0.076 |  1.138 |  1.410 |  1.000 |     2 |     27 |    1 | 1.43e-09
  0.07830 | 784 |  180.19 |  0.076 |  1.139 |  1.415 |  1.000 |     2 |     27 |    1 | 3.52e-09
  0.07840 | 785 |  180.19 |  0.076 |  1.139 |  1.420 |  1.000 |     2 |     28 |    1 | 9.17e-10
  0.07850 | 786 |  180.18 |  0.076 |  1.139 |  1.425 |  1.000 |     2 |     28 |    1 | 4.20e-10
  0.07860 | 787 |  180.17 |  0.076 |  1.139 |  1.430 |  1.000 |     2 |     28 |    1 | 1.78e-09
  0.07870 | 788 |  180.16 |  0.076 |  1.139 |  1.435 |  1.000 |     2 |     28 |    1 | 2.26e-09
  0.07880 | 789 |  180.16 |  0.076 |  1.139 |  1.440 |  1.000 |     2 |     28 |    1 | 5.88e-10
  0.07890 | 790 |  180.15 |  0.076 |  1.139 |  1.445 |  1.000 |     2 |     28 |    1 | 2.58e-09
  0.07900 | 791 |  180.14 |  0.076 |  1.139 |  1.450 |  1.000 |     2 |     28 |    1 | 8.36e-10

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.07910 | 792 |  180.14 |  0.076 |  1.140 |  1.455 |  1.000 |     2 |     28 |    1 | 3.03e-09
  0.07920 | 793 |  180.13 |  0.077 |  1.140 |  1.460 |  1.000 |     2 |     29 |    1 | 9.63e-09
  0.07930 | 794 |  180.13 |  0.077 |  1.140 |  1.465 |  1.000 |     2 |     28 |    1 | 1.71e-09
  0.07940 | 795 |  180.12 |  0.077 |  1.140 |  1.470 |  1.000 |     2 |     28 |    1 | 3.30e-09
  0.07950 | 796 |  180.11 |  0.077 |  1.140 |  1.475 |  1.000 |     2 |     28 |    1 | 2.19e-09
  0.07960 | 797 |  180.11 |  0.077 |  1.140 |  1.480 |  1.000 |     2 |     28 |    1 | 1.55e-09
  0.07970 | 798 |  180.10 |  0.077 |  1.140 |  1.485 |  1.000 |     2 |     28 |    1 | 2.27e-09
  0.07980 | 799 |  180.10 |  0.077 |  1.140 |  1.490 |  1.000 |     2 |     28 |    1 | 2.59e-09
  0.07990 | 800 |  180.09 |  0.077 |  1.140 |  1.495 |  1.000 |     2 |     28 |    1 | 2.04e-09
  0.08000 | 801 |  180.09 |  0.077 |  1.140 |  1.500 |  1.000 |     2 |     28 |    1 | 1.36e-09

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.08010 | 802 |  180.08 |  0.077 |  1.141 |  1.505 |  1.000 |     2 |     28 |    1 | 1.39e-09
  0.08020 | 803 |  180.08 |  0.077 |  1.141 |  1.510 |  1.000 |     2 |     29 |    1 | 1.94e-09
  0.08030 | 804 |  180.07 |  0.077 |  1.141 |  1.515 |  1.000 |     2 |     29 |    1 | 1.56e-09
  0.08040 | 805 |  180.06 |  0.077 |  1.141 |  1.520 |  1.000 |     2 |     29 |    1 | 2.35e-09
  0.08050 | 806 |  180.06 |  0.077 |  1.141 |  1.525 |  1.000 |     2 |     28 |    1 | 2.44e-09
  0.08060 | 807 |  180.06 |  0.077 |  1.141 |  1.530 |  1.000 |     2 |     28 |    1 | 3.16e-09
  0.08070 | 808 |  180.06 |  0.077 |  1.141 |  1.535 |  1.000 |     2 |     29 |    1 | 6.57e-10
  0.08080 | 809 |  180.06 |  0.077 |  1.141 |  1.540 |  1.000 |     2 |     29 |    1 | 1.22e-10
  0.08090 | 810 |  180.06 |  0.077 |  1.141 |  1.545 |  1.000 |     2 |     29 |    1 | 1.84e-09
  0.08100 | 811 |  180.06 |  0.077 |  1.141 |  1.550 |  1.000 |     2 |     28 |    1 | 1.85e-09

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.08110 | 812 |  180.05 |  0.077 |  1.141 |  1.555 |  1.000 |     2 |     28 |    1 | 1.12e-09
  0.08120 | 813 |  180.05 |  0.077 |  1.141 |  1.560 |  1.000 |     2 |     28 |    1 | 4.88e-10
  0.08130 | 814 |  180.05 |  0.077 |  1.141 |  1.565 |  1.000 |     2 |     28 |    1 | 2.73e-09
  0.08140 | 815 |  180.05 |  0.077 |  1.142 |  1.570 |  1.000 |     2 |     28 |    1 | 6.01e-10
  0.08150 | 816 |  180.05 |  0.077 |  1.142 |  1.575 |  1.000 |     2 |     28 |    1 | 1.56e-09
  0.08160 | 817 |  180.05 |  0.077 |  1.142 |  1.580 |  1.000 |     2 |     28 |    1 | 1.69e-09
  0.08170 | 818 |  180.05 |  0.077 |  1.142 |  1.585 |  1.000 |     2 |     28 |    1 | 9.56e-10
  0.08180 | 819 |  180.05 |  0.077 |  1.142 |  1.590 |  1.000 |     2 |     28 |    1 | 1.57e-09
  0.08190 | 820 |  180.05 |  0.077 |  1.142 |  1.595 |  1.000 |     2 |     29 |    1 | 1.97e-09
  0.08200 | 821 |  180.05 |  0.077 |  1.142 |  1.600 |  1.000 |     2 |     29 |    1 | 3.80e-10

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.08210 | 822 |  180.05 |  0.077 |  1.142 |  1.605 |  1.000 |     2 |     29 |    1 | 7.71e-09
  0.08220 | 823 |  180.05 |  0.077 |  1.142 |  1.610 |  1.000 |     2 |     29 |    1 | 2.99e-10
  0.08230 | 824 |  180.05 |  0.077 |  1.142 |  1.615 |  1.000 |     2 |     30 |    1 | 2.46e-09
  0.08240 | 825 |  180.05 |  0.077 |  1.142 |  1.620 |  1.000 |     2 |     30 |    1 | 2.04e-09
  0.08250 | 826 |  180.05 |  0.077 |  1.142 |  1.625 |  1.000 |     2 |     30 |    1 | 5.41e-10
  0.08260 | 827 |  180.06 |  0.077 |  1.142 |  1.630 |  1.000 |     2 |     30 |    1 | 2.00e-09
  0.08270 | 828 |  180.06 |  0.077 |  1.142 |  1.635 |  1.000 |     2 |     30 |    1 | 6.98e-10
  0.08280 | 829 |  180.06 |  0.077 |  1.142 |  1.640 |  1.000 |     2 |     29 |    1 | 1.57e-09
  0.08290 | 830 |  180.06 |  0.077 |  1.142 |  1.645 |  1.000 |     2 |     29 |    1 | 2.40e-10
  0.08300 | 831 |  180.07 |  0.077 |  1.143 |  1.650 |  1.000 |     2 |     29 |    1 | 2.11e-09

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.08310 | 832 |  180.07 |  0.077 |  1.143 |  1.655 |  1.000 |     2 |     29 |    1 | 1.97e-09
  0.08320 | 833 |  180.08 |  0.077 |  1.143 |  1.660 |  1.000 |     2 |     30 |    1 | 1.18e-09
  0.08330 | 834 |  180.08 |  0.077 |  1.143 |  1.665 |  1.000 |     2 |     30 |    1 | 8.02e-10
  0.08340 | 835 |  180.08 |  0.077 |  1.143 |  1.670 |  1.000 |     2 |     30 |    1 | 3.22e-10
  0.08350 | 836 |  180.09 |  0.077 |  1.143 |  1.675 |  1.000 |     2 |     30 |    1 | 2.36e-10
  0.08360 | 837 |  180.09 |  0.077 |  1.143 |  1.680 |  1.000 |     2 |     30 |    1 | 1.13e-09
  0.08370 | 838 |  180.10 |  0.077 |  1.143 |  1.685 |  1.000 |     2 |     30 |    1 | 2.45e-09
  0.08380 | 839 |  180.10 |  0.077 |  1.143 |  1.690 |  1.000 |     2 |     30 |    1 | 1.64e-09
  0.08390 | 840 |  180.11 |  0.077 |  1.143 |  1.695 |  1.000 |     2 |     30 |    1 | 2.76e-09
  0.08400 | 841 |  180.12 |  0.077 |  1.143 |  1.700 |  1.000 |     2 |     30 |    1 | 1.32e-09

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.08410 | 842 |  180.12 |  0.077 |  1.143 |  1.705 |  1.000 |     2 |     30 |    1 | 5.62e-10
  0.08420 | 843 |  180.13 |  0.077 |  1.143 |  1.710 |  1.000 |     2 |     30 |    1 | 2.04e-09
  0.08430 | 844 |  180.14 |  0.077 |  1.143 |  1.715 |  1.000 |     2 |     30 |    1 | 3.47e-09
  0.08440 | 845 |  180.16 |  0.077 |  1.143 |  1.720 |  1.000 |     2 |     30 |    1 | 3.16e-09
  0.08450 | 846 |  180.17 |  0.077 |  1.143 |  1.725 |  1.000 |     2 |     30 |    1 | 1.46e-09
  0.08460 | 847 |  180.18 |  0.077 |  1.143 |  1.730 |  1.000 |     2 |     30 |    1 | 2.29e-09
  0.08470 | 848 |  180.19 |  0.077 |  1.143 |  1.735 |  1.000 |     2 |     30 |    1 | 7.79e-10
  0.08480 | 849 |  180.20 |  0.077 |  1.144 |  1.740 |  1.000 |     2 |     30 |    1 | 2.11e-09
  0.08490 | 850 |  180.22 |  0.077 |  1.144 |  1.745 |  1.000 |     2 |     31 |    1 | 6.40e-08
  0.08500 | 851 |  180.24 |  0.077 |  1.144 |  1.750 |  1.000 |     2 |     30 |    1 | 4.08e-09

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.08510 | 852 |  180.25 |  0.077 |  1.144 |  1.755 |  1.000 |     2 |     30 |    1 | 8.15e-10
  0.08520 | 853 |  180.26 |  0.077 |  1.144 |  1.760 |  1.000 |     2 |     30 |    1 | 2.78e-09
  0.08530 | 854 |  180.28 |  0.077 |  1.144 |  1.765 |  1.000 |     2 |     30 |    1 | 1.45e-09
  0.08540 | 855 |  180.29 |  0.077 |  1.144 |  1.770 |  1.000 |     2 |     30 |    1 | 3.95e-09
  0.08550 | 856 |  180.31 |  0.077 |  1.144 |  1.775 |  1.000 |     2 |     30 |    1 | 3.19e-09
  0.08560 | 857 |  180.32 |  0.077 |  1.144 |  1.780 |  1.000 |     2 |     30 |    1 | 2.89e-09
  0.08570 | 858 |  180.33 |  0.077 |  1.144 |  1.785 |  1.000 |     2 |     30 |    1 | 7.01e-10
  0.08580 | 859 |  180.34 |  0.077 |  1.144 |  1.790 |  1.000 |     2 |     30 |    1 | 8.91e-10
  0.08590 | 860 |  180.35 |  0.077 |  1.144 |  1.795 |  1.000 |     2 |     30 |    1 | 5.09e-09
  0.08600 | 861 |  180.36 |  0.077 |  1.144 |  1.800 |  1.000 |     2 |     30 |    1 | 4.40e-09

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.08610 | 862 |  180.37 |  0.077 |  1.144 |  1.805 |  1.000 |     2 |     30 |    1 | 6.43e-10
  0.08620 | 863 |  180.37 |  0.077 |  1.144 |  1.810 |  1.000 |     2 |     30 |    1 | 3.27e-09
  0.08630 | 864 |  180.38 |  0.077 |  1.144 |  1.815 |  1.000 |     2 |     30 |    1 | 1.07e-09
  0.08640 | 865 |  180.39 |  0.077 |  1.144 |  1.820 |  1.000 |     2 |     30 |    1 | 1.94e-09
  0.08650 | 866 |  180.39 |  0.077 |  1.145 |  1.825 |  1.000 |     2 |     30 |    1 | 4.73e-09
  0.08660 | 867 |  180.40 |  0.077 |  1.145 |  1.830 |  1.000 |     2 |     30 |    1 | 6.79e-10
  0.08670 | 868 |  180.40 |  0.077 |  1.145 |  1.835 |  1.000 |     2 |     31 |    1 | 2.26e-09
  0.08680 | 869 |  180.40 |  0.077 |  1.145 |  1.840 |  1.000 |     2 |     31 |    1 | 9.47e-10
  0.08690 | 870 |  180.40 |  0.077 |  1.145 |  1.845 |  1.000 |     2 |     32 |    1 | 6.42e-10
  0.08700 | 871 |  180.40 |  0.077 |  1.145 |  1.850 |  1.000 |     2 |     32 |    1 | 3.24e-09

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.08710 | 872 |  180.40 |  0.077 |  1.145 |  1.855 |  1.000 |     1 |     18 |    1 | 9.29e-05
  0.08720 | 873 |  180.40 |  0.077 |  1.145 |  1.860 |  1.000 |     2 |     32 |    1 | 4.79e-09
  0.08730 | 874 |  180.41 |  0.077 |  1.145 |  1.865 |  1.000 |     2 |     32 |    1 | 9.16e-10
  0.08740 | 875 |  180.41 |  0.077 |  1.145 |  1.870 |  1.000 |     2 |     32 |    1 | 1.26e-10
  0.08750 | 876 |  180.41 |  0.077 |  1.145 |  1.875 |  1.000 |     2 |     31 |    1 | 3.02e-09
  0.08760 | 877 |  180.42 |  0.077 |  1.145 |  1.880 |  1.000 |     2 |     32 |    1 | 4.53e-09
  0.08770 | 878 |  180.42 |  0.077 |  1.145 |  1.885 |  1.000 |     2 |     32 |    1 | 1.77e-09
  0.08780 | 879 |  180.43 |  0.077 |  1.145 |  1.890 |  1.000 |     2 |     32 |    1 | 9.47e-10
  0.08790 | 880 |  180.43 |  0.077 |  1.145 |  1.895 |  1.000 |     2 |     32 |    1 | 3.64e-10
  0.08800 | 881 |  180.44 |  0.077 |  1.145 |  1.900 |  1.000 |     2 |     32 |    1 | 1.56e-09

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.08810 | 882 |  180.44 |  0.077 |  1.145 |  1.905 |  1.000 |     2 |     31 |    1 | 1.14e-09
  0.08820 | 883 |  180.45 |  0.077 |  1.145 |  1.910 |  1.000 |     2 |     31 |    1 | 1.12e-09
  0.08830 | 884 |  180.45 |  0.077 |  1.145 |  1.915 |  1.000 |     2 |     31 |    1 | 9.03e-10
  0.08840 | 885 |  180.45 |  0.077 |  1.145 |  1.920 |  1.000 |     2 |     32 |    1 | 6.40e-10
  0.08850 | 886 |  180.45 |  0.077 |  1.145 |  1.925 |  1.000 |     2 |     32 |    1 | 2.76e-09
  0.08860 | 887 |  180.45 |  0.077 |  1.145 |  1.930 |  1.000 |     2 |     32 |    1 | 3.73e-09
  0.08870 | 888 |  180.45 |  0.077 |  1.145 |  1.935 |  1.000 |     2 |     32 |    1 | 6.06e-10
  0.08880 | 889 |  180.45 |  0.077 |  1.145 |  1.940 |  1.000 |     2 |     31 |    1 | 2.68e-09
  0.08890 | 890 |  180.46 |  0.077 |  1.145 |  1.945 |  1.000 |     2 |     32 |    1 | 3.20e-09
  0.08900 | 891 |  180.46 |  0.077 |  1.145 |  1.950 |  1.000 |     2 |     32 |    1 | 5.62e-09

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.08910 | 892 |  180.46 |  0.077 |  1.145 |  1.955 |  1.000 |     2 |     32 |    1 | 2.43e-09
  0.08920 | 893 |  180.46 |  0.077 |  1.145 |  1.960 |  1.000 |     2 |     32 |    1 | 1.30e-09
  0.08930 | 894 |  180.46 |  0.077 |  1.145 |  1.965 |  1.000 |     2 |     32 |    1 | 5.21e-10
  0.08940 | 895 |  180.46 |  0.077 |  1.145 |  1.970 |  1.000 |     2 |     32 |    1 | 2.20e-09
  0.08950 | 896 |  180.47 |  0.077 |  1.145 |  1.975 |  1.000 |     2 |     32 |    1 | 2.68e-09
  0.08960 | 897 |  180.48 |  0.077 |  1.145 |  1.980 |  1.000 |     2 |     32 |    1 | 1.98e-09
  0.08970 | 898 |  180.49 |  0.077 |  1.145 |  1.985 |  1.000 |     2 |     32 |    1 | 2.86e-09
  0.08980 | 899 |  180.50 |  0.077 |  1.146 |  1.990 |  1.000 |     2 |     32 |    1 | 3.44e-09
  0.08990 | 900 |  180.51 |  0.077 |  1.146 |  1.995 |  1.000 |     2 |     32 |    1 | 5.28e-09
  0.09000 | 901 |  180.52 |  0.077 |  1.146 |  2.000 |  1.000 |     2 |     32 |    1 | 7.34e-09

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.09010 | 902 |  180.45 |  0.077 |  1.145 |  2.000 |  1.000 |     2 |     33 |    1 | 1.50e-08
  0.09020 | 903 |  180.39 |  0.077 |  1.145 |  2.000 |  1.000 |     2 |     33 |    1 | 1.47e-08
  0.09030 | 904 |  180.32 |  0.077 |  1.145 |  2.000 |  1.000 |     2 |     33 |    1 | 2.10e-08
  0.09040 | 905 |  180.25 |  0.077 |  1.145 |  2.000 |  1.000 |     2 |     33 |    1 | 2.66e-08
  0.09050 | 906 |  180.18 |  0.077 |  1.144 |  2.000 |  1.000 |     2 |     33 |    1 | 2.75e-08
  0.09060 | 907 |  180.11 |  0.077 |  1.144 |  2.000 |  1.000 |     2 |     33 |    1 | 3.64e-08
  0.09070 | 908 |  180.04 |  0.077 |  1.144 |  2.000 |  1.000 |     2 |     33 |    1 | 3.66e-08
  0.09080 | 909 |  179.96 |  0.077 |  1.143 |  2.000 |  1.000 |     2 |     33 |    1 | 4.52e-08
  0.09090 | 910 |  179.89 |  0.077 |  1.143 |  2.000 |  1.000 |     2 |     33 |    1 | 5.28e-08
  0.09100 | 911 |  179.81 |  0.077 |  1.143 |  2.000 |  1.000 |     2 |     33 |    1 | 6.71e-08

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.09110 | 912 |  179.74 |  0.077 |  1.142 |  2.000 |  1.000 |     2 |     34 |    1 | 8.11e-08
  0.09120 | 913 |  179.68 |  0.077 |  1.142 |  2.000 |  1.000 |     2 |     33 |    1 | 3.33e-08
  0.09130 | 914 |  179.61 |  0.077 |  1.142 |  2.000 |  1.000 |     2 |     33 |    1 | 4.40e-08
  0.09140 | 915 |  179.55 |  0.077 |  1.142 |  2.000 |  1.000 |     2 |     33 |    1 | 4.73e-08
  0.09150 | 916 |  179.49 |  0.077 |  1.141 |  2.000 |  1.000 |     2 |     33 |    1 | 6.33e-08
  0.09160 | 917 |  179.44 |  0.077 |  1.141 |  2.000 |  1.000 |     2 |     34 |    1 | 7.94e-08
  0.09170 | 918 |  179.38 |  0.077 |  1.141 |  2.000 |  1.000 |     2 |     35 |    1 | 1.07e-07
  0.09180 | 919 |  179.32 |  0.076 |  1.140 |  2.000 |  1.000 |     2 |     35 |    1 | 1.32e-07
  0.09190 | 920 |  179.27 |  0.076 |  1.140 |  2.000 |  1.000 |     2 |     35 |    1 | 5.19e-09
  0.09200 | 921 |  179.23 |  0.076 |  1.140 |  2.000 |  1.000 |     2 |     34 |    1 | 6.64e-08

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.09210 | 922 |  179.20 |  0.076 |  1.140 |  2.000 |  1.000 |     2 |     35 |    1 | 2.48e-06
  0.09220 | 923 |  179.17 |  0.076 |  1.139 |  2.000 |  1.000 |     2 |     33 |    1 | 8.61e-08
  0.09230 | 924 |  179.15 |  0.076 |  1.139 |  2.000 |  1.000 |     2 |     34 |    1 | 2.85e-07
  0.09240 | 925 |  179.13 |  0.075 |  1.139 |  2.000 |  1.000 |     2 |     33 |    1 | 8.21e-08
  0.09250 | 926 |  179.11 |  0.075 |  1.139 |  2.000 |  1.000 |     2 |     33 |    1 | 8.10e-08
  0.09260 | 927 |  179.09 |  0.075 |  1.139 |  2.000 |  1.000 |     2 |     33 |    1 | 8.68e-08
  0.09270 | 928 |  179.08 |  0.075 |  1.138 |  2.000 |  1.000 |     2 |     33 |    1 | 8.16e-08
  0.09280 | 929 |  179.06 |  0.075 |  1.138 |  2.000 |  1.000 |     2 |     33 |    1 | 7.49e-08
  0.09290 | 930 |  179.05 |  0.074 |  1.138 |  2.000 |  1.000 |     2 |     33 |    1 | 6.44e-08
  0.09300 | 931 |  179.04 |  0.074 |  1.138 |  2.000 |  1.000 |     2 |     33 |    1 | 6.71e-08

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.09310 | 932 |  179.04 |  0.074 |  1.138 |  2.000 |  1.000 |     2 |     33 |    1 | 3.92e-05
  0.09320 | 933 |  179.05 |  0.074 |  1.138 |  2.000 |  1.000 |     2 |     33 |    1 | 1.74e-08
  0.09330 | 934 |  179.06 |  0.074 |  1.138 |  2.000 |  1.000 |     1 |     18 |    1 | 6.05e-05
  0.09340 | 935 |  179.07 |  0.074 |  1.137 |  2.000 |  1.000 |     2 |     34 |    1 | 7.33e-08
  0.09350 | 936 |  179.08 |  0.074 |  1.137 |  2.000 |  1.000 |     1 |     18 |    1 | 5.21e-05
  0.09360 | 937 |  179.10 |  0.074 |  1.137 |  2.000 |  1.000 |     1 |     18 |    1 | 6.53e-05
  0.09370 | 938 |  179.12 |  0.074 |  1.137 |  2.000 |  1.000 |     1 |     18 |    1 | 7.29e-05
  0.09380 | 939 |  179.15 |  0.073 |  1.137 |  2.000 |  1.000 |     1 |     18 |    1 | 7.66e-05
  0.09390 | 940 |  179.17 |  0.073 |  1.137 |  2.000 |  1.000 |     1 |     18 |    1 | 7.59e-05
  0.09400 | 941 |  179.20 |  0.073 |  1.137 |  2.000 |  1.000 |     1 |     18 |    1 | 7.26e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.09410 | 942 |  179.22 |  0.073 |  1.137 |  2.000 |  1.000 |     1 |     18 |    1 | 6.76e-05
  0.09420 | 943 |  179.25 |  0.073 |  1.137 |  2.000 |  1.000 |     1 |     18 |    1 | 6.16e-05
  0.09430 | 944 |  179.28 |  0.073 |  1.137 |  2.000 |  1.000 |     1 |     18 |    1 | 5.51e-05
  0.09440 | 945 |  179.30 |  0.073 |  1.137 |  2.000 |  1.000 |     1 |     18 |    1 | 4.91e-05
  0.09450 | 946 |  179.33 |  0.073 |  1.137 |  2.000 |  1.000 |     1 |     18 |    1 | 4.38e-05
  0.09460 | 947 |  179.36 |  0.073 |  1.137 |  2.000 |  1.000 |     1 |     18 |    1 | 3.96e-05
  0.09470 | 948 |  179.39 |  0.073 |  1.137 |  2.000 |  1.000 |     1 |     18 |    1 | 3.62e-05
  0.09480 | 949 |  179.42 |  0.073 |  1.137 |  2.000 |  1.000 |     1 |     18 |    1 | 3.41e-05
  0.09490 | 950 |  179.45 |  0.073 |  1.137 |  2.000 |  1.000 |     1 |     18 |    1 | 3.30e-05
  0.09500 | 951 |  179.47 |  0.073 |  1.137 |  2.000 |  1.000 |     1 |     18 |    1 | 3.29e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.09510 | 952 |  179.50 |  0.073 |  1.137 |  2.000 |  1.000 |     1 |     18 |    1 | 3.35e-05
  0.09520 | 953 |  179.53 |  0.073 |  1.137 |  2.000 |  1.000 |     1 |     18 |    1 | 3.49e-05
  0.09530 | 954 |  179.55 |  0.073 |  1.137 |  2.000 |  1.000 |     1 |     18 |    1 | 3.64e-05
  0.09540 | 955 |  179.58 |  0.073 |  1.136 |  2.000 |  1.000 |     1 |     18 |    1 | 3.82e-05
  0.09550 | 956 |  179.60 |  0.073 |  1.136 |  2.000 |  1.000 |     1 |     18 |    1 | 4.01e-05
  0.09560 | 957 |  179.63 |  0.073 |  1.136 |  2.000 |  1.000 |     1 |     18 |    1 | 4.19e-05
  0.09570 | 958 |  179.65 |  0.073 |  1.136 |  2.000 |  1.000 |     1 |     18 |    1 | 4.36e-05
  0.09580 | 959 |  179.67 |  0.073 |  1.136 |  2.000 |  1.000 |     1 |     18 |    1 | 4.48e-05
  0.09590 | 960 |  179.70 |  0.074 |  1.136 |  2.000 |  1.000 |     1 |     18 |    1 | 4.57e-05
  0.09600 | 961 |  179.72 |  0.074 |  1.136 |  2.000 |  1.000 |     1 |     18 |    1 | 4.63e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.09610 | 962 |  179.73 |  0.074 |  1.136 |  2.000 |  1.000 |     1 |     18 |    1 | 4.62e-05
  0.09620 | 963 |  179.75 |  0.074 |  1.136 |  2.000 |  1.000 |     1 |     18 |    1 | 4.56e-05
  0.09630 | 964 |  179.77 |  0.074 |  1.136 |  2.000 |  1.000 |     1 |     18 |    1 | 4.43e-05
  0.09640 | 965 |  179.79 |  0.074 |  1.136 |  2.000 |  1.000 |     1 |     18 |    1 | 5.08e-05
  0.09650 | 966 |  179.81 |  0.074 |  1.136 |  2.000 |  1.000 |     1 |     18 |    1 | 5.03e-05
  0.09660 | 967 |  179.82 |  0.074 |  1.136 |  2.000 |  1.000 |     1 |     18 |    1 | 4.51e-05
  0.09670 | 968 |  179.84 |  0.074 |  1.136 |  2.000 |  1.000 |     1 |     18 |    1 | 4.19e-05
  0.09680 | 969 |  179.85 |  0.074 |  1.136 |  2.000 |  1.000 |     1 |     18 |    1 | 3.96e-05
  0.09690 | 970 |  179.87 |  0.074 |  1.136 |  2.000 |  1.000 |     1 |     18 |    1 | 3.69e-05
  0.09700 | 971 |  179.88 |  0.074 |  1.136 |  2.000 |  1.000 |     1 |     18 |    1 | 3.40e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.09710 | 972 |  179.90 |  0.074 |  1.136 |  2.000 |  1.000 |     1 |     18 |    1 | 3.07e-05
  0.09720 | 973 |  179.91 |  0.074 |  1.136 |  2.000 |  1.000 |     1 |     18 |    1 | 2.78e-05
  0.09730 | 974 |  179.92 |  0.074 |  1.136 |  2.000 |  1.000 |     1 |     18 |    1 | 2.46e-05
  0.09740 | 975 |  179.93 |  0.074 |  1.135 |  2.000 |  1.000 |     1 |     18 |    1 | 2.10e-05
  0.09750 | 976 |  179.94 |  0.074 |  1.135 |  2.000 |  1.000 |     1 |     18 |    1 | 1.73e-05
  0.09760 | 977 |  179.95 |  0.074 |  1.135 |  2.000 |  1.000 |     1 |     18 |    1 | 1.36e-05
  0.09770 | 978 |  179.96 |  0.074 |  1.135 |  2.000 |  1.000 |     1 |     18 |    1 | 9.89e-06
  0.09780 | 979 |  179.97 |  0.074 |  1.135 |  2.000 |  1.000 |     1 |     18 |    1 | 6.51e-06
  0.09790 | 980 |  179.98 |  0.074 |  1.135 |  2.000 |  1.000 |     1 |     18 |    1 | 2.46e-06
  0.09800 | 981 |  179.99 |  0.074 |  1.135 |  2.000 |  1.000 |     1 |     18 |    1 | 1.27e-06

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.09810 | 982 |  179.99 |  0.074 |  1.135 |  2.000 |  1.000 |     1 |     18 |    1 | 4.44e-06
  0.09820 | 983 |  180.00 |  0.074 |  1.135 |  2.000 |  1.000 |     1 |     18 |    1 | 8.17e-06
  0.09830 | 984 |  180.01 |  0.074 |  1.135 |  2.000 |  1.000 |     1 |     18 |    1 | 1.14e-05
  0.09840 | 985 |  180.01 |  0.074 |  1.135 |  2.000 |  1.000 |     1 |     18 |    1 | 1.48e-05
  0.09850 | 986 |  180.02 |  0.074 |  1.135 |  2.000 |  1.000 |     1 |     18 |    1 | 1.82e-05
  0.09860 | 987 |  180.03 |  0.074 |  1.135 |  2.000 |  1.000 |     1 |     18 |    1 | 2.11e-05
  0.09870 | 988 |  180.03 |  0.074 |  1.135 |  2.000 |  1.000 |     1 |     18 |    1 | 2.38e-05
  0.09880 | 989 |  180.04 |  0.074 |  1.135 |  2.000 |  1.000 |     1 |     18 |    1 | 2.65e-05
  0.09890 | 990 |  180.04 |  0.074 |  1.134 |  2.000 |  1.000 |     1 |     18 |    1 | 2.88e-05
  0.09900 | 991 |  180.04 |  0.074 |  1.134 |  2.000 |  1.000 |     1 |     18 |    1 | 3.16e-05

        t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ==============================================================================================
  0.09910 | 992 |  180.05 |  0.075 |  1.134 |  2.000 |  1.000 |     1 |     18 |    1 | 3.38e-05
  0.09920 | 993 |  180.05 |  0.075 |  1.134 |  2.000 |  1.000 |     1 |     18 |    1 | 3.54e-05
  0.09930 | 994 |  180.06 |  0.075 |  1.134 |  2.000 |  1.000 |     2 |     34 |    1 | 5.41e-08
  0.09940 | 995 |  180.06 |  0.075 |  1.134 |  2.000 |  1.000 |     1 |     18 |    1 | 2.00e-05
  0.09950 | 996 |  180.07 |  0.075 |  1.134 |  2.000 |  1.000 |     1 |     18 |    1 | 2.14e-05
  0.09960 | 997 |  180.07 |  0.075 |  1.134 |  2.000 |  1.000 |     1 |     18 |    1 | 2.25e-05
  0.09970 | 998 |  180.08 |  0.075 |  1.134 |  2.000 |  1.000 |     1 |     18 |    1 | 2.31e-05
  0.09980 | 999 |  180.08 |  0.075 |  1.134 |  2.000 |  1.000 |     1 |     18 |    1 | 2.35e-05
  0.09990 | 1000 |  180.09 |  0.075 |  1.134 |  2.000 |  1.000 |     1 |     18 |    1 | 2.42e-05
  0.10000 | 1001 |  180.10 |  0.075 |  1.134 |  2.000 |  1.000 |     1 |     18 |    1 | 1.07e-05
Time elapsed: 1303.81[s]

Run RZP simulation

Lz.P.ctrlfct = @meqctrl;
Lz.P.ctrlpar = ctrlpar;
Lz.P.debug=0;
LYz = rzpt(Lz,LX);
[..................................................] 1001/1001

Run FGEL simulation

just switch L.P.lin to true

Ll=L;
Ll.P.lin = true;
LYl = fget(Ll,LX);
[..................................................] 1001/1001

Throw an error if simulation diverges

assert(max(abs(LYl.zIp./LYl.Ip - ctrlpar.zQref.Data'))<1e-2,'too large z error');

Compare time evolutions

clf;
meqplotevo(L,LY,LYl,LYz,'argplot',{'Ia','Va','Ip','zIp','rIp','Bm','bp','li','qA'});
legend('FGE','FGEL','RZP');
shg

Post-simulation movie from FGE simulation

if debuglevel
  clf;
  meqmovie(L,LY,'decimate',5,'dosave',false)
end