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-5,'tolF',1e-4,... % low tolerance both for fbt and fge
  'cde','OhmTor_rigid_0D',...
  'ilim',3,'icsint',true}; % spline interpolation for limiter and axis finding

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

disp('Computing RZP linearization')
[Lz,LXz] = rzp(params{:},'agcon',{'bp', 'li'});
Lz  = rzpl(Lz,LXz); % 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);
LXz = rzpx([],t,L,LX); % Start from FGE initial guess but set RZP ag values

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 + jfnk
  GMRES: aya
  prec:1, mkryl:200, kmax:50
  tolF:1.000e-04
  Vessel:v, nu:200
  Doublet: 0
  CDE: OhmTor_rigid_0D
  Constraints: 
    Fun    iD       LXfield     ag index 
   ----------------------------------- 
    {@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.087 |  1.092 |  1.000 |  1.000 |     1 |      5 |    1 | 7.51e-08
   0.00020 |   3 |  200.00 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      5 |    1 | 9.16e-08
   0.00030 |   4 |  200.00 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 1.37e-07
   0.00040 |   5 |  200.00 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 1.59e-07
   0.00050 |   6 |  200.00 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 1.23e-07
   0.00060 |   7 |  200.00 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 1.40e-07
   0.00070 |   8 |  200.00 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 5.96e-08
   0.00080 |   9 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 2.65e-07
   0.00090 |  10 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 5.60e-08
   0.00100 |  11 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 8.46e-08

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.00110 |  12 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 3.22e-08
   0.00120 |  13 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 2.12e-08
   0.00130 |  14 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 1.11e-07
   0.00140 |  15 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 1.04e-07
   0.00150 |  16 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 1.35e-08
   0.00160 |  17 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 7.33e-09
   0.00170 |  18 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 3.68e-08
   0.00180 |  19 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 3.38e-08
   0.00190 |  20 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 4.64e-08
   0.00200 |  21 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 7.65e-08

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.00210 |  22 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 1.21e-08
   0.00220 |  23 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 1.40e-08
   0.00230 |  24 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 2.16e-07
   0.00240 |  25 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 3.29e-08
   0.00250 |  26 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 1.12e-07
   0.00260 |  27 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 1.87e-07
   0.00270 |  28 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 2.38e-08
   0.00280 |  29 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 1.81e-08
   0.00290 |  30 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 5.87e-08
   0.00300 |  31 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 5.17e-08

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.00310 |  32 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 1.82e-07
   0.00320 |  33 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 1.74e-07
   0.00330 |  34 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 3.04e-08
   0.00340 |  35 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 8.71e-08
   0.00350 |  36 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 3.75e-07
   0.00360 |  37 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      6 |    1 | 1.00e-07
   0.00370 |  38 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      8 |    1 | 9.81e-08
   0.00380 |  39 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      7 |    1 | 2.72e-08
   0.00390 |  40 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      7 |    1 | 3.12e-08
   0.00400 |  41 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      8 |    1 | 1.86e-07

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.00410 |  42 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      7 |    1 | 1.85e-08
   0.00420 |  43 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      8 |    1 | 3.08e-07
   0.00430 |  44 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      7 |    1 | 1.15e-07
   0.00440 |  45 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      8 |    1 | 3.22e-07
   0.00450 |  46 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      8 |    1 | 3.60e-08
   0.00460 |  47 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      8 |    1 | 6.83e-08
   0.00470 |  48 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      8 |    1 | 4.70e-08
   0.00480 |  49 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      8 |    1 | 2.20e-08
   0.00490 |  50 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      8 |    1 | 2.12e-07
   0.00500 |  51 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      8 |    1 | 1.12e-07

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.00510 |  52 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      8 |    1 | 2.83e-08
   0.00520 |  53 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      8 |    1 | 7.41e-08
   0.00530 |  54 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      8 |    1 | 5.31e-08
   0.00540 |  55 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      8 |    1 | 3.51e-06
   0.00550 |  56 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      8 |    1 | 1.17e-05
   0.00560 |  57 |  199.99 |  0.087 |  1.092 |  1.000 |  1.000 |     1 |      9 |    1 | 4.32e-05
   0.00570 |  58 |  199.99 |  0.088 |  1.092 |  1.000 |  1.000 |     1 |      9 |    1 | 8.77e-05
   0.00580 |  59 |  199.99 |  0.088 |  1.092 |  1.000 |  1.000 |     2 |     16 |    1 | 1.87e-08
   0.00590 |  60 |  199.99 |  0.088 |  1.092 |  1.000 |  1.000 |     2 |     16 |    1 | 1.10e-08
   0.00600 |  61 |  199.99 |  0.088 |  1.092 |  1.000 |  1.000 |     2 |     16 |    1 | 1.15e-08

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.00610 |  62 |  199.99 |  0.088 |  1.092 |  1.000 |  1.000 |     2 |     18 |    1 | 5.39e-07
   0.00620 |  63 |  199.99 |  0.088 |  1.092 |  1.000 |  1.000 |     2 |     18 |    1 | 6.61e-07
   0.00630 |  64 |  200.00 |  0.089 |  1.092 |  1.000 |  1.000 |     2 |     18 |    1 | 4.20e-08
   0.00640 |  65 |  200.00 |  0.089 |  1.092 |  1.000 |  1.000 |     2 |     18 |    1 | 5.65e-08
   0.00650 |  66 |  200.00 |  0.089 |  1.092 |  1.000 |  1.000 |     2 |     18 |    1 | 4.24e-07
   0.00660 |  67 |  200.00 |  0.090 |  1.092 |  1.000 |  1.000 |     2 |     18 |    1 | 2.97e-07
   0.00670 |  68 |  200.00 |  0.090 |  1.092 |  1.000 |  1.000 |     2 |     19 |    1 | 1.27e-06
   0.00680 |  69 |  200.01 |  0.090 |  1.092 |  1.000 |  1.000 |     2 |     19 |    1 | 4.58e-07
   0.00690 |  70 |  200.01 |  0.091 |  1.092 |  1.000 |  1.000 |     2 |     19 |    1 | 6.13e-08
   0.00700 |  71 |  200.02 |  0.091 |  1.092 |  1.000 |  1.000 |     2 |     19 |    1 | 1.27e-07

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.00710 |  72 |  200.02 |  0.091 |  1.092 |  1.000 |  1.000 |     2 |     20 |    1 | 1.74e-07
   0.00720 |  73 |  200.03 |  0.092 |  1.092 |  1.000 |  1.000 |     2 |     20 |    1 | 1.50e-07
   0.00730 |  74 |  200.03 |  0.092 |  1.092 |  1.000 |  1.000 |     2 |     20 |    1 | 8.12e-08
   0.00740 |  75 |  200.03 |  0.092 |  1.092 |  1.000 |  1.000 |     2 |     20 |    1 | 6.84e-08
   0.00750 |  76 |  200.03 |  0.093 |  1.092 |  1.000 |  1.000 |     2 |     20 |    1 | 2.62e-07
   0.00760 |  77 |  200.03 |  0.093 |  1.092 |  1.000 |  1.000 |     2 |     20 |    1 | 1.17e-06
   0.00770 |  78 |  200.04 |  0.094 |  1.092 |  1.000 |  1.000 |     2 |     20 |    1 | 2.87e-07
   0.00780 |  79 |  200.04 |  0.094 |  1.092 |  1.000 |  1.000 |     2 |     20 |    1 | 1.94e-07
   0.00790 |  80 |  200.04 |  0.094 |  1.092 |  1.000 |  1.000 |     2 |     20 |    1 | 9.36e-08
   0.00800 |  81 |  200.04 |  0.095 |  1.092 |  1.000 |  1.000 |     2 |     22 |    1 | 3.52e-07

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.00810 |  82 |  200.04 |  0.095 |  1.093 |  1.000 |  1.000 |     2 |     21 |    1 | 4.16e-07
   0.00820 |  83 |  200.03 |  0.096 |  1.093 |  1.000 |  1.000 |     2 |     21 |    1 | 2.50e-07
   0.00830 |  84 |  200.03 |  0.096 |  1.093 |  1.000 |  1.000 |     2 |     21 |    1 | 5.52e-08
   0.00840 |  85 |  200.03 |  0.096 |  1.093 |  1.000 |  1.000 |     2 |     21 |    1 | 2.24e-07
   0.00850 |  86 |  200.03 |  0.097 |  1.093 |  1.000 |  1.000 |     2 |     22 |    1 | 2.27e-07
   0.00860 |  87 |  200.02 |  0.097 |  1.093 |  1.000 |  1.000 |     2 |     21 |    1 | 7.30e-08
   0.00870 |  88 |  200.02 |  0.098 |  1.093 |  1.000 |  1.000 |     2 |     22 |    1 | 5.18e-07
   0.00880 |  89 |  200.01 |  0.098 |  1.093 |  1.000 |  1.000 |     2 |     21 |    1 | 2.89e-07
   0.00890 |  90 |  200.00 |  0.098 |  1.093 |  1.000 |  1.000 |     2 |     21 |    1 | 1.97e-07
   0.00900 |  91 |  200.00 |  0.099 |  1.093 |  1.000 |  1.000 |     2 |     22 |    1 | 1.78e-07

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.00910 |  92 |  199.99 |  0.099 |  1.093 |  1.000 |  1.000 |     2 |     22 |    1 | 4.32e-07
   0.00920 |  93 |  199.99 |  0.100 |  1.093 |  1.000 |  1.000 |     2 |     22 |    1 | 2.69e-07
   0.00930 |  94 |  199.98 |  0.100 |  1.093 |  1.000 |  1.000 |     2 |     22 |    1 | 3.16e-07
   0.00940 |  95 |  199.97 |  0.101 |  1.093 |  1.000 |  1.000 |     2 |     23 |    1 | 1.88e-07
   0.00950 |  96 |  199.97 |  0.101 |  1.093 |  1.000 |  1.000 |     2 |     23 |    1 | 1.05e-06
   0.00960 |  97 |  199.97 |  0.101 |  1.093 |  1.000 |  1.000 |     2 |     23 |    1 | 6.88e-07
   0.00970 |  98 |  199.96 |  0.102 |  1.093 |  1.000 |  1.000 |     2 |     23 |    1 | 8.37e-08
   0.00980 |  99 |  199.96 |  0.102 |  1.093 |  1.000 |  1.000 |     2 |     23 |    1 | 3.77e-07
   0.00990 | 100 |  199.95 |  0.103 |  1.093 |  1.000 |  1.000 |     2 |     23 |    1 | 2.14e-07
   0.01000 | 101 |  199.95 |  0.103 |  1.093 |  1.000 |  1.000 |     2 |     23 |    1 | 3.30e-07

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.01010 | 102 |  199.95 |  0.104 |  1.093 |  1.000 |  1.000 |     2 |     23 |    1 | 3.65e-07
   0.01020 | 103 |  199.95 |  0.104 |  1.093 |  1.000 |  1.000 |     2 |     24 |    1 | 1.13e-07
   0.01030 | 104 |  199.95 |  0.105 |  1.093 |  1.000 |  1.000 |     2 |     23 |    1 | 2.36e-07
   0.01040 | 105 |  199.95 |  0.105 |  1.093 |  1.000 |  1.000 |     2 |     23 |    1 | 1.34e-07
   0.01050 | 106 |  199.95 |  0.105 |  1.093 |  1.000 |  1.000 |     2 |     24 |    1 | 2.86e-07
   0.01060 | 107 |  199.95 |  0.106 |  1.093 |  1.000 |  1.000 |     2 |     23 |    1 | 4.86e-07
   0.01070 | 108 |  199.95 |  0.106 |  1.093 |  1.000 |  1.000 |     2 |     23 |    1 | 1.38e-07
   0.01080 | 109 |  199.95 |  0.106 |  1.093 |  1.000 |  1.000 |     2 |     23 |    1 | 2.01e-07
   0.01090 | 110 |  199.95 |  0.107 |  1.093 |  1.000 |  1.000 |     2 |     24 |    1 | 1.60e-06
   0.01100 | 111 |  199.95 |  0.107 |  1.093 |  1.000 |  1.000 |     2 |     23 |    1 | 2.71e-07

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.01110 | 112 |  199.95 |  0.107 |  1.093 |  1.000 |  1.000 |     2 |     23 |    1 | 3.42e-08
   0.01120 | 113 |  199.95 |  0.108 |  1.093 |  1.000 |  1.000 |     2 |     23 |    1 | 2.72e-08
   0.01130 | 114 |  199.95 |  0.108 |  1.093 |  1.000 |  1.000 |     2 |     23 |    1 | 8.49e-08
   0.01140 | 115 |  199.95 |  0.108 |  1.092 |  1.000 |  1.000 |     2 |     23 |    1 | 3.82e-08
   0.01150 | 116 |  199.96 |  0.108 |  1.092 |  1.000 |  1.000 |     2 |     24 |    1 | 3.85e-08
   0.01160 | 117 |  199.96 |  0.108 |  1.092 |  1.000 |  1.000 |     2 |     24 |    1 | 9.81e-09
   0.01170 | 118 |  199.96 |  0.108 |  1.092 |  1.000 |  1.000 |     2 |     24 |    1 | 2.01e-08
   0.01180 | 119 |  199.96 |  0.108 |  1.092 |  1.000 |  1.000 |     2 |     24 |    1 | 2.47e-08
   0.01190 | 120 |  199.95 |  0.109 |  1.092 |  1.000 |  1.000 |     2 |     24 |    1 | 2.52e-08
   0.01200 | 121 |  199.95 |  0.109 |  1.092 |  1.000 |  1.000 |     2 |     24 |    1 | 2.18e-07

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.01210 | 122 |  199.95 |  0.109 |  1.092 |  1.000 |  1.000 |     2 |     24 |    1 | 8.97e-09
   0.01220 | 123 |  199.95 |  0.109 |  1.092 |  1.000 |  1.000 |     1 |     14 |    1 | 9.42e-05
   0.01230 | 124 |  199.95 |  0.109 |  1.092 |  1.000 |  1.000 |     1 |     14 |    1 | 6.99e-05
   0.01240 | 125 |  199.95 |  0.109 |  1.092 |  1.000 |  1.000 |     1 |     14 |    1 | 7.08e-05
   0.01250 | 126 |  199.95 |  0.109 |  1.092 |  1.000 |  1.000 |     1 |     14 |    1 | 5.21e-05
   0.01260 | 127 |  199.95 |  0.109 |  1.092 |  1.000 |  1.000 |     1 |     14 |    1 | 4.55e-05
   0.01270 | 128 |  199.95 |  0.109 |  1.092 |  1.000 |  1.000 |     1 |     14 |    1 | 4.37e-05
   0.01280 | 129 |  199.95 |  0.109 |  1.092 |  1.000 |  1.000 |     1 |     14 |    1 | 4.14e-05
   0.01290 | 130 |  199.94 |  0.109 |  1.092 |  1.000 |  1.000 |     1 |     14 |    1 | 3.41e-05
   0.01300 | 131 |  199.94 |  0.109 |  1.092 |  1.000 |  1.000 |     1 |     14 |    1 | 1.55e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.01310 | 132 |  199.94 |  0.109 |  1.092 |  1.000 |  1.000 |     1 |     14 |    1 | 2.45e-05
   0.01320 | 133 |  199.94 |  0.109 |  1.092 |  1.000 |  1.000 |     1 |     14 |    1 | 2.56e-05
   0.01330 | 134 |  199.94 |  0.109 |  1.091 |  1.000 |  1.000 |     1 |     14 |    1 | 3.25e-05
   0.01340 | 135 |  199.94 |  0.109 |  1.091 |  1.000 |  1.000 |     1 |     14 |    1 | 9.38e-07
   0.01350 | 136 |  199.94 |  0.109 |  1.091 |  1.000 |  1.000 |     1 |     14 |    1 | 2.65e-05
   0.01360 | 137 |  199.94 |  0.109 |  1.091 |  1.000 |  1.000 |     1 |     14 |    1 | 2.76e-05
   0.01370 | 138 |  199.95 |  0.109 |  1.091 |  1.000 |  1.000 |     1 |     14 |    1 | 2.75e-05
   0.01380 | 139 |  199.95 |  0.109 |  1.091 |  1.000 |  1.000 |     1 |     14 |    1 | 1.73e-05
   0.01390 | 140 |  199.95 |  0.109 |  1.091 |  1.000 |  1.000 |     1 |     14 |    1 | 5.48e-06
   0.01400 | 141 |  199.95 |  0.109 |  1.091 |  1.000 |  1.000 |     1 |     14 |    1 | 4.69e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.01410 | 142 |  199.95 |  0.109 |  1.091 |  1.000 |  1.000 |     1 |     14 |    1 | 9.11e-06
   0.01420 | 143 |  199.95 |  0.109 |  1.091 |  1.000 |  1.000 |     1 |     14 |    1 | 1.73e-05
   0.01430 | 144 |  199.96 |  0.109 |  1.091 |  1.000 |  1.000 |     1 |     14 |    1 | 3.17e-05
   0.01440 | 145 |  199.96 |  0.109 |  1.091 |  1.000 |  1.000 |     1 |     14 |    1 | 6.69e-05
   0.01450 | 146 |  199.96 |  0.109 |  1.091 |  1.000 |  1.000 |     1 |     14 |    1 | 2.43e-05
   0.01460 | 147 |  199.96 |  0.109 |  1.091 |  1.000 |  1.000 |     1 |     14 |    1 | 4.35e-05
   0.01470 | 148 |  199.97 |  0.108 |  1.091 |  1.000 |  1.000 |     1 |     14 |    1 | 2.79e-05
   0.01480 | 149 |  199.97 |  0.108 |  1.091 |  1.000 |  1.000 |     1 |     14 |    1 | 4.97e-05
   0.01490 | 150 |  199.97 |  0.108 |  1.091 |  1.000 |  1.000 |     1 |     14 |    1 | 5.12e-05
   0.01500 | 151 |  199.98 |  0.108 |  1.091 |  1.000 |  1.000 |     1 |     14 |    1 | 5.29e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.01510 | 152 |  199.98 |  0.108 |  1.091 |  1.000 |  1.000 |     1 |     14 |    1 | 4.81e-05
   0.01520 | 153 |  199.99 |  0.108 |  1.091 |  1.000 |  1.000 |     1 |     14 |    1 | 3.44e-05
   0.01530 | 154 |  199.99 |  0.108 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.52e-05
   0.01540 | 155 |  199.99 |  0.108 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 3.38e-05
   0.01550 | 156 |  200.00 |  0.108 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 3.61e-05
   0.01560 | 157 |  200.00 |  0.108 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.87e-05
   0.01570 | 158 |  200.01 |  0.108 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 3.85e-05
   0.01580 | 159 |  200.01 |  0.108 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 3.16e-05
   0.01590 | 160 |  200.01 |  0.108 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 4.11e-05
   0.01600 | 161 |  200.02 |  0.108 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 4.01e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.01610 | 162 |  200.02 |  0.108 |  1.091 |  1.000 |  1.000 |     2 |     24 |    1 | 1.97e-07
   0.01620 | 163 |  200.03 |  0.108 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 3.14e-05
   0.01630 | 164 |  200.03 |  0.108 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 4.50e-05
   0.01640 | 165 |  200.04 |  0.108 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 3.57e-05
   0.01650 | 166 |  200.04 |  0.108 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 3.03e-05
   0.01660 | 167 |  200.04 |  0.108 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.59e-05
   0.01670 | 168 |  200.05 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.22e-05
   0.01680 | 169 |  200.05 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.16e-05
   0.01690 | 170 |  200.06 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.02e-05
   0.01700 | 171 |  200.06 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.65e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.01710 | 172 |  200.06 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 4.44e-05
   0.01720 | 173 |  200.07 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.56e-05
   0.01730 | 174 |  200.07 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.49e-05
   0.01740 | 175 |  200.07 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.00e-05
   0.01750 | 176 |  200.08 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.23e-05
   0.01760 | 177 |  200.08 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.46e-05
   0.01770 | 178 |  200.08 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.68e-05
   0.01780 | 179 |  200.08 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 3.27e-06
   0.01790 | 180 |  200.09 |  0.107 |  1.091 |  1.000 |  1.000 |     2 |     24 |    1 | 4.68e-07
   0.01800 | 181 |  200.09 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 6.62e-06

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.01810 | 182 |  200.09 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 6.43e-06
   0.01820 | 183 |  200.10 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.18e-05
   0.01830 | 184 |  200.10 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 3.06e-05
   0.01840 | 185 |  200.10 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 8.07e-06
   0.01850 | 186 |  200.10 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.84e-05
   0.01860 | 187 |  200.10 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.58e-05
   0.01870 | 188 |  200.10 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 7.98e-06
   0.01880 | 189 |  200.11 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.06e-05
   0.01890 | 190 |  200.11 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 9.62e-06
   0.01900 | 191 |  200.11 |  0.107 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.72e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.01910 | 192 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.80e-05
   0.01920 | 193 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 5.23e-06
   0.01930 | 194 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 7.07e-06
   0.01940 | 195 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.32e-05
   0.01950 | 196 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.64e-05
   0.01960 | 197 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.49e-05
   0.01970 | 198 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.83e-05
   0.01980 | 199 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.80e-05
   0.01990 | 200 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 5.97e-06
   0.02000 | 201 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 3.49e-06

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.02010 | 202 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.01e-05
   0.02020 | 203 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.30e-05
   0.02030 | 204 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.31e-05
   0.02040 | 205 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.22e-05
   0.02050 | 206 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.18e-05
   0.02060 | 207 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.57e-05
   0.02070 | 208 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 5.54e-06
   0.02080 | 209 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.34e-05
   0.02090 | 210 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 4.57e-06
   0.02100 | 211 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.28e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.02110 | 212 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.28e-05
   0.02120 | 213 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.44e-05
   0.02130 | 214 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 8.04e-06
   0.02140 | 215 |  200.11 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.07e-05
   0.02150 | 216 |  200.12 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.11e-05
   0.02160 | 217 |  200.12 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 8.18e-06
   0.02170 | 218 |  200.12 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.37e-05
   0.02180 | 219 |  200.12 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.52e-05
   0.02190 | 220 |  200.12 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.42e-05
   0.02200 | 221 |  200.12 |  0.106 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 1.96e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.02210 | 222 |  200.13 |  0.105 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.11e-05
   0.02220 | 223 |  200.13 |  0.105 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.34e-05
   0.02230 | 224 |  200.13 |  0.105 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.48e-05
   0.02240 | 225 |  200.13 |  0.105 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.03e-05
   0.02250 | 226 |  200.14 |  0.105 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.38e-05
   0.02260 | 227 |  200.14 |  0.105 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.52e-05
   0.02270 | 228 |  200.14 |  0.105 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.71e-05
   0.02280 | 229 |  200.15 |  0.105 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.73e-05
   0.02290 | 230 |  200.15 |  0.105 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.92e-05
   0.02300 | 231 |  200.16 |  0.105 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 3.00e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.02310 | 232 |  200.16 |  0.105 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 2.96e-05
   0.02320 | 233 |  200.16 |  0.105 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 3.39e-05
   0.02330 | 234 |  200.17 |  0.105 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 3.56e-05
   0.02340 | 235 |  200.17 |  0.105 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 3.09e-05
   0.02350 | 236 |  200.17 |  0.105 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 3.39e-05
   0.02360 | 237 |  200.18 |  0.105 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 3.76e-05
   0.02370 | 238 |  200.18 |  0.105 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 3.83e-05
   0.02380 | 239 |  200.19 |  0.105 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 3.21e-05
   0.02390 | 240 |  200.19 |  0.105 |  1.091 |  1.000 |  1.000 |     1 |     13 |    1 | 3.70e-05
   0.02400 | 241 |  200.19 |  0.105 |  1.092 |  1.000 |  1.000 |     1 |     13 |    1 | 4.06e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.02410 | 242 |  200.20 |  0.105 |  1.092 |  1.000 |  1.000 |     1 |     13 |    1 | 4.52e-05
   0.02420 | 243 |  200.20 |  0.105 |  1.092 |  1.000 |  1.000 |     1 |     13 |    1 | 4.84e-05
   0.02430 | 244 |  200.20 |  0.105 |  1.092 |  1.000 |  1.000 |     1 |     13 |    1 | 3.21e-05
   0.02440 | 245 |  200.21 |  0.105 |  1.092 |  1.000 |  1.000 |     1 |     13 |    1 | 3.84e-05
   0.02450 | 246 |  200.21 |  0.105 |  1.092 |  1.000 |  1.000 |     1 |     13 |    1 | 5.05e-05
   0.02460 | 247 |  200.21 |  0.105 |  1.092 |  1.000 |  1.000 |     1 |     13 |    1 | 4.83e-05
   0.02470 | 248 |  200.22 |  0.105 |  1.092 |  1.000 |  1.000 |     1 |     13 |    1 | 7.04e-05
   0.02480 | 249 |  200.22 |  0.105 |  1.092 |  1.000 |  1.000 |     1 |     13 |    1 | 4.13e-05
   0.02490 | 250 |  200.22 |  0.105 |  1.092 |  1.000 |  1.000 |     1 |     13 |    1 | 6.51e-05
   0.02500 | 251 |  200.22 |  0.105 |  1.092 |  1.000 |  1.000 |     1 |     13 |    1 | 4.70e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.02510 | 252 |  200.23 |  0.105 |  1.092 |  1.000 |  1.000 |     1 |     13 |    1 | 6.69e-05
   0.02520 | 253 |  200.23 |  0.105 |  1.092 |  1.000 |  1.000 |     1 |     13 |    1 | 5.37e-05
   0.02530 | 254 |  200.23 |  0.105 |  1.092 |  1.000 |  1.000 |     1 |     13 |    1 | 7.00e-05
   0.02540 | 255 |  200.23 |  0.105 |  1.092 |  1.000 |  1.000 |     1 |     13 |    1 | 5.21e-05
   0.02550 | 256 |  200.23 |  0.105 |  1.092 |  1.000 |  1.000 |     1 |     13 |    1 | 6.30e-05
   0.02560 | 257 |  200.23 |  0.105 |  1.093 |  1.000 |  1.000 |     1 |     13 |    1 | 3.41e-05
   0.02570 | 258 |  200.24 |  0.105 |  1.093 |  1.000 |  1.000 |     1 |     13 |    1 | 8.04e-05
   0.02580 | 259 |  200.24 |  0.105 |  1.093 |  1.000 |  1.000 |     1 |     14 |    1 | 6.25e-05
   0.02590 | 260 |  200.24 |  0.105 |  1.093 |  1.000 |  1.000 |     1 |     14 |    1 | 6.09e-05
   0.02600 | 261 |  200.24 |  0.105 |  1.093 |  1.000 |  1.000 |     1 |     14 |    1 | 5.96e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.02610 | 262 |  200.24 |  0.105 |  1.093 |  1.000 |  1.000 |     1 |     14 |    1 | 7.93e-05
   0.02620 | 263 |  200.24 |  0.105 |  1.093 |  1.000 |  1.000 |     1 |     14 |    1 | 5.88e-05
   0.02630 | 264 |  200.24 |  0.105 |  1.093 |  1.000 |  1.000 |     1 |     14 |    1 | 6.14e-05
   0.02640 | 265 |  200.24 |  0.105 |  1.093 |  1.000 |  1.000 |     1 |     14 |    1 | 8.45e-05
   0.02650 | 266 |  200.24 |  0.105 |  1.093 |  1.000 |  1.000 |     1 |     14 |    1 | 6.84e-05
   0.02660 | 267 |  200.25 |  0.105 |  1.093 |  1.000 |  1.000 |     1 |     14 |    1 | 6.42e-05
   0.02670 | 268 |  200.25 |  0.105 |  1.093 |  1.000 |  1.000 |     1 |     14 |    1 | 3.78e-05
   0.02680 | 269 |  200.25 |  0.105 |  1.093 |  1.000 |  1.000 |     1 |     14 |    1 | 3.90e-05
   0.02690 | 270 |  200.25 |  0.105 |  1.093 |  1.000 |  1.000 |     1 |     14 |    1 | 5.76e-05
   0.02700 | 271 |  200.25 |  0.105 |  1.093 |  1.000 |  1.000 |     1 |     14 |    1 | 5.74e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.02710 | 272 |  200.25 |  0.105 |  1.094 |  1.000 |  1.000 |     1 |     14 |    1 | 7.26e-05
   0.02720 | 273 |  200.25 |  0.105 |  1.094 |  1.000 |  1.000 |     1 |     14 |    1 | 9.79e-05
   0.02730 | 274 |  200.25 |  0.105 |  1.094 |  1.000 |  1.000 |     1 |     14 |    1 | 6.63e-05
   0.02740 | 275 |  200.25 |  0.105 |  1.094 |  1.000 |  1.000 |     1 |     14 |    1 | 3.35e-05
   0.02750 | 276 |  200.24 |  0.105 |  1.094 |  1.000 |  1.000 |     2 |     24 |    1 | 6.93e-09
   0.02760 | 277 |  200.24 |  0.105 |  1.094 |  1.000 |  1.000 |     1 |     14 |    1 | 9.06e-05
   0.02770 | 278 |  200.24 |  0.105 |  1.094 |  1.000 |  1.000 |     1 |     14 |    1 | 4.83e-06
   0.02780 | 279 |  200.24 |  0.105 |  1.094 |  1.000 |  1.000 |     2 |     24 |    1 | 6.16e-09
   0.02790 | 280 |  200.24 |  0.105 |  1.094 |  1.000 |  1.000 |     1 |     14 |    1 | 5.09e-05
   0.02800 | 281 |  200.24 |  0.105 |  1.094 |  1.000 |  1.000 |     1 |     14 |    1 | 8.97e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.02810 | 282 |  200.24 |  0.105 |  1.094 |  1.000 |  1.000 |     2 |     24 |    1 | 4.61e-07
   0.02820 | 283 |  200.24 |  0.105 |  1.094 |  1.000 |  1.000 |     1 |     14 |    1 | 5.20e-05
   0.02830 | 284 |  200.24 |  0.105 |  1.094 |  1.000 |  1.000 |     1 |     14 |    1 | 2.60e-05
   0.02840 | 285 |  200.24 |  0.105 |  1.095 |  1.000 |  1.000 |     1 |     14 |    1 | 8.19e-05
   0.02850 | 286 |  200.24 |  0.105 |  1.095 |  1.000 |  1.000 |     1 |     14 |    1 | 5.67e-05
   0.02860 | 287 |  200.24 |  0.105 |  1.095 |  1.000 |  1.000 |     1 |     14 |    1 | 6.18e-05
   0.02870 | 288 |  200.24 |  0.105 |  1.095 |  1.000 |  1.000 |     2 |     24 |    1 | 1.07e-08
   0.02880 | 289 |  200.23 |  0.105 |  1.095 |  1.000 |  1.000 |     2 |     24 |    1 | 1.45e-08
   0.02890 | 290 |  200.23 |  0.105 |  1.095 |  1.000 |  1.000 |     1 |     14 |    1 | 5.47e-05
   0.02900 | 291 |  200.23 |  0.105 |  1.095 |  1.000 |  1.000 |     1 |     14 |    1 | 7.30e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.02910 | 292 |  200.23 |  0.105 |  1.095 |  1.000 |  1.000 |     1 |     14 |    1 | 7.25e-05
   0.02920 | 293 |  200.23 |  0.105 |  1.095 |  1.000 |  1.000 |     1 |     14 |    1 | 4.53e-05
   0.02930 | 294 |  200.22 |  0.105 |  1.095 |  1.000 |  1.000 |     2 |     24 |    1 | 4.57e-07
   0.02940 | 295 |  200.22 |  0.105 |  1.095 |  1.000 |  1.000 |     1 |     14 |    1 | 6.41e-05
   0.02950 | 296 |  200.22 |  0.105 |  1.095 |  1.000 |  1.000 |     1 |     14 |    1 | 7.33e-05
   0.02960 | 297 |  200.22 |  0.105 |  1.095 |  1.000 |  1.000 |     1 |     14 |    1 | 8.78e-05
   0.02970 | 298 |  200.21 |  0.105 |  1.096 |  1.000 |  1.000 |     1 |     14 |    1 | 3.32e-05
   0.02980 | 299 |  200.21 |  0.105 |  1.096 |  1.000 |  1.000 |     1 |     14 |    1 | 6.00e-05
   0.02990 | 300 |  200.21 |  0.105 |  1.096 |  1.000 |  1.000 |     1 |     14 |    1 | 9.60e-05
   0.03000 | 301 |  200.21 |  0.105 |  1.096 |  1.000 |  1.000 |     1 |     14 |    1 | 6.86e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.03010 | 302 |  200.20 |  0.105 |  1.096 |  1.000 |  1.000 |     1 |     14 |    1 | 9.65e-06
   0.03020 | 303 |  200.20 |  0.105 |  1.096 |  1.000 |  1.000 |     2 |     24 |    1 | 8.99e-09
   0.03030 | 304 |  200.20 |  0.105 |  1.096 |  1.000 |  1.000 |     2 |     24 |    1 | 1.01e-08
   0.03040 | 305 |  200.20 |  0.105 |  1.096 |  1.000 |  1.000 |     1 |     14 |    1 | 7.51e-05
   0.03050 | 306 |  200.20 |  0.105 |  1.096 |  1.000 |  1.000 |     1 |     14 |    1 | 6.16e-05
   0.03060 | 307 |  200.19 |  0.105 |  1.096 |  1.000 |  1.000 |     1 |     14 |    1 | 4.07e-05
   0.03070 | 308 |  200.19 |  0.105 |  1.096 |  1.000 |  1.000 |     1 |     14 |    1 | 7.38e-05
   0.03080 | 309 |  200.19 |  0.105 |  1.096 |  1.000 |  1.000 |     1 |     14 |    1 | 5.62e-05
   0.03090 | 310 |  200.18 |  0.105 |  1.096 |  1.000 |  1.000 |     1 |     14 |    1 | 5.71e-05
   0.03100 | 311 |  200.18 |  0.105 |  1.097 |  1.000 |  1.000 |     1 |     14 |    1 | 1.37e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.03110 | 312 |  200.18 |  0.105 |  1.097 |  1.000 |  1.000 |     1 |     14 |    1 | 8.86e-05
   0.03120 | 313 |  200.18 |  0.105 |  1.097 |  1.000 |  1.000 |     2 |     24 |    1 | 6.55e-09
   0.03130 | 314 |  200.17 |  0.105 |  1.097 |  1.000 |  1.000 |     2 |     24 |    1 | 2.26e-07
   0.03140 | 315 |  200.17 |  0.105 |  1.097 |  1.000 |  1.000 |     1 |     14 |    1 | 8.18e-05
   0.03150 | 316 |  200.17 |  0.105 |  1.097 |  1.000 |  1.000 |     1 |     14 |    1 | 2.52e-05
   0.03160 | 317 |  200.16 |  0.105 |  1.097 |  1.000 |  1.000 |     1 |     14 |    1 | 3.38e-05
   0.03170 | 318 |  200.16 |  0.105 |  1.097 |  1.000 |  1.000 |     1 |     14 |    1 | 6.97e-05
   0.03180 | 319 |  200.16 |  0.105 |  1.097 |  1.000 |  1.000 |     1 |     14 |    1 | 3.94e-05
   0.03190 | 320 |  200.15 |  0.105 |  1.097 |  1.000 |  1.000 |     2 |     24 |    1 | 3.50e-09
   0.03200 | 321 |  200.15 |  0.105 |  1.097 |  1.000 |  1.000 |     1 |     14 |    1 | 5.81e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.03210 | 322 |  200.15 |  0.105 |  1.097 |  1.000 |  1.000 |     1 |     14 |    1 | 3.83e-05
   0.03220 | 323 |  200.14 |  0.105 |  1.097 |  1.000 |  1.000 |     1 |     14 |    1 | 9.19e-05
   0.03230 | 324 |  200.14 |  0.105 |  1.098 |  1.000 |  1.000 |     1 |     14 |    1 | 5.10e-05
   0.03240 | 325 |  200.14 |  0.105 |  1.098 |  1.000 |  1.000 |     1 |     14 |    1 | 7.18e-05
   0.03250 | 326 |  200.13 |  0.105 |  1.098 |  1.000 |  1.000 |     2 |     24 |    1 | 1.28e-08
   0.03260 | 327 |  200.13 |  0.105 |  1.098 |  1.000 |  1.000 |     1 |     14 |    1 | 6.32e-05
   0.03270 | 328 |  200.13 |  0.105 |  1.098 |  1.000 |  1.000 |     2 |     24 |    1 | 3.77e-07
   0.03280 | 329 |  200.12 |  0.105 |  1.098 |  1.000 |  1.000 |     1 |     14 |    1 | 7.69e-05
   0.03290 | 330 |  200.12 |  0.105 |  1.098 |  1.000 |  1.000 |     1 |     14 |    1 | 3.01e-05
   0.03300 | 331 |  200.12 |  0.105 |  1.098 |  1.000 |  1.000 |     1 |     14 |    1 | 3.62e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.03310 | 332 |  200.11 |  0.105 |  1.098 |  1.000 |  1.000 |     1 |     14 |    1 | 4.57e-05
   0.03320 | 333 |  200.11 |  0.105 |  1.098 |  1.000 |  1.000 |     1 |     14 |    1 | 5.54e-05
   0.03330 | 334 |  200.11 |  0.105 |  1.098 |  1.000 |  1.000 |     1 |     14 |    1 | 6.04e-05
   0.03340 | 335 |  200.10 |  0.105 |  1.098 |  1.000 |  1.000 |     1 |     14 |    1 | 1.82e-05
   0.03350 | 336 |  200.10 |  0.105 |  1.098 |  1.000 |  1.000 |     1 |     14 |    1 | 2.14e-05
   0.03360 | 337 |  200.10 |  0.105 |  1.098 |  1.000 |  1.000 |     1 |     14 |    1 | 6.74e-05
   0.03370 | 338 |  200.09 |  0.105 |  1.099 |  1.000 |  1.000 |     1 |     14 |    1 | 6.21e-05
   0.03380 | 339 |  200.09 |  0.105 |  1.099 |  1.000 |  1.000 |     1 |     14 |    1 | 5.38e-05
   0.03390 | 340 |  200.09 |  0.105 |  1.099 |  1.000 |  1.000 |     1 |     14 |    1 | 7.84e-05
   0.03400 | 341 |  200.09 |  0.105 |  1.099 |  1.000 |  1.000 |     2 |     24 |    1 | 3.96e-09

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.03410 | 342 |  200.08 |  0.105 |  1.099 |  1.000 |  1.000 |     1 |     14 |    1 | 8.28e-05
   0.03420 | 343 |  200.08 |  0.105 |  1.099 |  1.000 |  1.000 |     1 |     14 |    1 | 5.54e-05
   0.03430 | 344 |  200.08 |  0.105 |  1.099 |  1.000 |  1.000 |     2 |     24 |    1 | 7.86e-09
   0.03440 | 345 |  200.07 |  0.105 |  1.099 |  1.000 |  1.000 |     1 |     14 |    1 | 6.11e-05
   0.03450 | 346 |  200.07 |  0.105 |  1.099 |  1.000 |  1.000 |     1 |     14 |    1 | 8.05e-05
   0.03460 | 347 |  200.07 |  0.105 |  1.099 |  1.000 |  1.000 |     1 |     14 |    1 | 5.65e-05
   0.03470 | 348 |  200.06 |  0.105 |  1.099 |  1.000 |  1.000 |     2 |     24 |    1 | 3.82e-09
   0.03480 | 349 |  200.06 |  0.105 |  1.099 |  1.000 |  1.000 |     1 |     14 |    1 | 9.10e-05
   0.03490 | 350 |  200.06 |  0.105 |  1.099 |  1.000 |  1.000 |     1 |     14 |    1 | 3.60e-05
   0.03500 | 351 |  200.06 |  0.105 |  1.100 |  1.000 |  1.000 |     1 |     14 |    1 | 5.58e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.03510 | 352 |  200.05 |  0.105 |  1.100 |  1.000 |  1.000 |     1 |     14 |    1 | 5.82e-05
   0.03520 | 353 |  200.05 |  0.105 |  1.100 |  1.000 |  1.000 |     1 |     14 |    1 | 4.24e-05
   0.03530 | 354 |  200.05 |  0.105 |  1.100 |  1.000 |  1.000 |     1 |     14 |    1 | 8.76e-05
   0.03540 | 355 |  200.05 |  0.105 |  1.100 |  1.000 |  1.000 |     1 |     14 |    1 | 2.84e-05
   0.03550 | 356 |  200.04 |  0.105 |  1.100 |  1.000 |  1.000 |     1 |     14 |    1 | 5.62e-05
   0.03560 | 357 |  200.04 |  0.105 |  1.100 |  1.000 |  1.000 |     1 |     14 |    1 | 8.82e-05
   0.03570 | 358 |  200.04 |  0.105 |  1.100 |  1.000 |  1.000 |     2 |     24 |    1 | 2.24e-09
   0.03580 | 359 |  200.04 |  0.105 |  1.100 |  1.000 |  1.000 |     1 |     14 |    1 | 7.10e-05
   0.03590 | 360 |  200.03 |  0.105 |  1.100 |  1.000 |  1.000 |     1 |     14 |    1 | 2.41e-05
   0.03600 | 361 |  200.03 |  0.105 |  1.100 |  1.000 |  1.000 |     1 |     14 |    1 | 5.29e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.03610 | 362 |  200.03 |  0.105 |  1.100 |  1.000 |  1.000 |     1 |     14 |    1 | 3.34e-05
   0.03620 | 363 |  200.03 |  0.105 |  1.101 |  1.000 |  1.000 |     1 |     14 |    1 | 3.60e-05
   0.03630 | 364 |  200.03 |  0.105 |  1.101 |  1.000 |  1.000 |     1 |     14 |    1 | 7.82e-05
   0.03640 | 365 |  200.03 |  0.105 |  1.101 |  1.000 |  1.000 |     2 |     24 |    1 | 9.42e-09
   0.03650 | 366 |  200.02 |  0.105 |  1.101 |  1.000 |  1.000 |     1 |     14 |    1 | 4.04e-05
   0.03660 | 367 |  200.02 |  0.105 |  1.101 |  1.000 |  1.000 |     1 |     14 |    1 | 2.51e-06
   0.03670 | 368 |  200.02 |  0.105 |  1.101 |  1.000 |  1.000 |     1 |     14 |    1 | 2.89e-05
   0.03680 | 369 |  200.02 |  0.105 |  1.101 |  1.000 |  1.000 |     1 |     14 |    1 | 4.19e-05
   0.03690 | 370 |  200.02 |  0.105 |  1.101 |  1.000 |  1.000 |     1 |     14 |    1 | 4.89e-05
   0.03700 | 371 |  200.02 |  0.105 |  1.101 |  1.000 |  1.000 |     2 |     24 |    1 | 4.86e-09

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.03710 | 372 |  200.02 |  0.105 |  1.101 |  1.000 |  1.000 |     1 |     14 |    1 | 5.87e-05
   0.03720 | 373 |  200.02 |  0.105 |  1.101 |  1.000 |  1.000 |     1 |     14 |    1 | 7.70e-05
   0.03730 | 374 |  200.01 |  0.105 |  1.101 |  1.000 |  1.000 |     2 |     24 |    1 | 5.65e-09
   0.03740 | 375 |  200.01 |  0.105 |  1.101 |  1.000 |  1.000 |     1 |     14 |    1 | 5.76e-05
   0.03750 | 376 |  200.01 |  0.105 |  1.102 |  1.000 |  1.000 |     1 |     14 |    1 | 6.89e-05
   0.03760 | 377 |  200.01 |  0.105 |  1.102 |  1.000 |  1.000 |     1 |     14 |    1 | 7.49e-05
   0.03770 | 378 |  200.01 |  0.105 |  1.102 |  1.000 |  1.000 |     1 |     14 |    1 | 3.94e-05
   0.03780 | 379 |  200.01 |  0.105 |  1.102 |  1.000 |  1.000 |     1 |     14 |    1 | 9.67e-05
   0.03790 | 380 |  200.01 |  0.105 |  1.102 |  1.000 |  1.000 |     2 |     24 |    1 | 1.66e-08
   0.03800 | 381 |  200.00 |  0.105 |  1.102 |  1.000 |  1.000 |     1 |     14 |    1 | 6.30e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.03810 | 382 |  200.00 |  0.105 |  1.102 |  1.000 |  1.000 |     1 |     14 |    1 | 1.52e-05
   0.03820 | 383 |  200.00 |  0.105 |  1.102 |  1.000 |  1.000 |     1 |     14 |    1 | 8.06e-05
   0.03830 | 384 |  200.00 |  0.105 |  1.102 |  1.000 |  1.000 |     1 |     14 |    1 | 2.98e-05
   0.03840 | 385 |  200.00 |  0.105 |  1.102 |  1.000 |  1.000 |     1 |     14 |    1 | 4.07e-05
   0.03850 | 386 |  200.00 |  0.105 |  1.102 |  1.000 |  1.000 |     1 |     14 |    1 | 4.65e-05
   0.03860 | 387 |  200.00 |  0.105 |  1.102 |  1.000 |  1.000 |     1 |     14 |    1 | 5.73e-05
   0.03870 | 388 |  200.00 |  0.105 |  1.102 |  1.000 |  1.000 |     1 |     14 |    1 | 3.70e-05
   0.03880 | 389 |  200.00 |  0.105 |  1.103 |  1.000 |  1.000 |     2 |     24 |    1 | 7.85e-09
   0.03890 | 390 |  200.00 |  0.105 |  1.103 |  1.000 |  1.000 |     1 |     14 |    1 | 9.10e-05
   0.03900 | 391 |  200.00 |  0.105 |  1.103 |  1.000 |  1.000 |     2 |     25 |    1 | 3.90e-09

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.03910 | 392 |  200.00 |  0.105 |  1.103 |  1.000 |  1.000 |     2 |     24 |    1 | 6.85e-09
   0.03920 | 393 |  200.00 |  0.105 |  1.103 |  1.000 |  1.000 |     1 |     14 |    1 | 9.99e-05
   0.03930 | 394 |  200.00 |  0.105 |  1.103 |  1.000 |  1.000 |     1 |     14 |    1 | 3.26e-05
   0.03940 | 395 |  200.00 |  0.105 |  1.103 |  1.000 |  1.000 |     2 |     25 |    1 | 2.14e-08
   0.03950 | 396 |  200.01 |  0.105 |  1.103 |  1.000 |  1.000 |     2 |     25 |    1 | 1.79e-08
   0.03960 | 397 |  200.01 |  0.105 |  1.103 |  1.000 |  1.000 |     1 |     14 |    1 | 9.21e-05
   0.03970 | 398 |  200.01 |  0.105 |  1.103 |  1.000 |  1.000 |     2 |     25 |    1 | 3.94e-07
   0.03980 | 399 |  200.01 |  0.105 |  1.103 |  1.000 |  1.000 |     1 |     14 |    1 | 6.34e-05
   0.03990 | 400 |  200.01 |  0.105 |  1.104 |  1.000 |  1.000 |     2 |     25 |    1 | 4.58e-09
   0.04000 | 401 |  200.01 |  0.105 |  1.104 |  1.000 |  1.000 |     1 |     14 |    1 | 6.81e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.04010 | 402 |  200.01 |  0.105 |  1.104 |  1.000 |  1.000 |     2 |     25 |    1 | 2.43e-08
   0.04020 | 403 |  200.01 |  0.105 |  1.104 |  1.000 |  1.000 |     1 |     14 |    1 | 9.21e-05
   0.04030 | 404 |  200.01 |  0.105 |  1.104 |  1.000 |  1.000 |     2 |     24 |    1 | 1.44e-08
   0.04040 | 405 |  200.01 |  0.105 |  1.104 |  1.000 |  1.000 |     1 |     14 |    1 | 6.98e-05
   0.04050 | 406 |  200.01 |  0.105 |  1.104 |  1.000 |  1.000 |     1 |     14 |    1 | 6.98e-05
   0.04060 | 407 |  200.02 |  0.105 |  1.104 |  1.000 |  1.000 |     1 |     14 |    1 | 5.75e-05
   0.04070 | 408 |  200.01 |  0.105 |  1.104 |  1.000 |  1.000 |     1 |     14 |    1 | 8.13e-05
   0.04080 | 409 |  200.01 |  0.105 |  1.104 |  1.000 |  1.000 |     1 |     14 |    1 | 2.32e-05
   0.04090 | 410 |  200.01 |  0.105 |  1.104 |  1.000 |  1.000 |     1 |     14 |    1 | 7.48e-05
   0.04100 | 411 |  200.01 |  0.105 |  1.104 |  1.000 |  1.000 |     1 |     14 |    1 | 5.46e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.04110 | 412 |  200.01 |  0.105 |  1.104 |  1.000 |  1.000 |     1 |     14 |    1 | 8.99e-05
   0.04120 | 413 |  200.01 |  0.105 |  1.105 |  1.000 |  1.000 |     2 |     24 |    1 | 2.70e-08
   0.04130 | 414 |  200.01 |  0.105 |  1.105 |  1.000 |  1.000 |     2 |     24 |    1 | 1.25e-08
   0.04140 | 415 |  200.01 |  0.105 |  1.105 |  1.000 |  1.000 |     1 |     14 |    1 | 3.12e-05
   0.04150 | 416 |  200.01 |  0.105 |  1.105 |  1.000 |  1.000 |     2 |     24 |    1 | 2.38e-08
   0.04160 | 417 |  200.01 |  0.105 |  1.105 |  1.000 |  1.000 |     1 |     14 |    1 | 6.77e-05
   0.04170 | 418 |  200.01 |  0.105 |  1.105 |  1.000 |  1.000 |     1 |     14 |    1 | 6.78e-05
   0.04180 | 419 |  200.01 |  0.105 |  1.105 |  1.000 |  1.000 |     2 |     24 |    1 | 1.84e-08
   0.04190 | 420 |  200.01 |  0.105 |  1.105 |  1.000 |  1.000 |     1 |     14 |    1 | 7.86e-05
   0.04200 | 421 |  200.01 |  0.105 |  1.105 |  1.000 |  1.000 |     1 |     14 |    1 | 7.97e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.04210 | 422 |  200.01 |  0.105 |  1.105 |  1.000 |  1.000 |     1 |     14 |    1 | 6.24e-05
   0.04220 | 423 |  200.01 |  0.105 |  1.105 |  1.000 |  1.000 |     2 |     24 |    1 | 1.73e-08
   0.04230 | 424 |  200.01 |  0.105 |  1.105 |  1.000 |  1.000 |     2 |     25 |    1 | 3.09e-08
   0.04240 | 425 |  200.01 |  0.105 |  1.106 |  1.000 |  1.000 |     1 |     14 |    1 | 2.79e-05
   0.04250 | 426 |  200.01 |  0.105 |  1.106 |  1.000 |  1.000 |     2 |     24 |    1 | 1.53e-08
   0.04260 | 427 |  200.01 |  0.105 |  1.106 |  1.000 |  1.000 |     1 |     14 |    1 | 2.24e-05
   0.04270 | 428 |  200.01 |  0.105 |  1.106 |  1.000 |  1.000 |     1 |     14 |    1 | 9.04e-05
   0.04280 | 429 |  200.01 |  0.105 |  1.106 |  1.000 |  1.000 |     1 |     14 |    1 | 6.87e-05
   0.04290 | 430 |  200.01 |  0.105 |  1.106 |  1.000 |  1.000 |     1 |     14 |    1 | 7.98e-05
   0.04300 | 431 |  200.01 |  0.105 |  1.106 |  1.000 |  1.000 |     2 |     25 |    1 | 4.78e-09

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.04310 | 432 |  200.00 |  0.105 |  1.106 |  1.000 |  1.000 |     2 |     25 |    1 | 2.20e-08
   0.04320 | 433 |  200.00 |  0.105 |  1.106 |  1.000 |  1.000 |     1 |     14 |    1 | 9.98e-05
   0.04330 | 434 |  200.00 |  0.105 |  1.106 |  1.000 |  1.000 |     1 |     14 |    1 | 5.33e-05
   0.04340 | 435 |  200.00 |  0.105 |  1.106 |  1.000 |  1.000 |     2 |     24 |    1 | 9.73e-09
   0.04350 | 436 |  200.00 |  0.105 |  1.106 |  1.000 |  1.000 |     1 |     14 |    1 | 9.18e-05
   0.04360 | 437 |  200.00 |  0.105 |  1.106 |  1.000 |  1.000 |     2 |     24 |    1 | 1.94e-08
   0.04370 | 438 |  200.00 |  0.105 |  1.107 |  1.000 |  1.000 |     2 |     25 |    1 | 1.81e-08
   0.04380 | 439 |  200.00 |  0.105 |  1.107 |  1.000 |  1.000 |     1 |     14 |    1 | 8.05e-05
   0.04390 | 440 |  200.00 |  0.105 |  1.107 |  1.000 |  1.000 |     2 |     24 |    1 | 5.23e-09
   0.04400 | 441 |  200.00 |  0.105 |  1.107 |  1.000 |  1.000 |     1 |     14 |    1 | 1.59e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.04410 | 442 |  200.00 |  0.105 |  1.107 |  1.000 |  1.000 |     1 |     14 |    1 | 9.69e-05
   0.04420 | 443 |  200.00 |  0.105 |  1.107 |  1.000 |  1.000 |     1 |     14 |    1 | 4.73e-05
   0.04430 | 444 |  200.00 |  0.105 |  1.107 |  1.000 |  1.000 |     1 |     14 |    1 | 1.46e-05
   0.04440 | 445 |  200.00 |  0.105 |  1.107 |  1.000 |  1.000 |     1 |     14 |    1 | 5.15e-05
   0.04450 | 446 |  200.00 |  0.105 |  1.107 |  1.000 |  1.000 |     2 |     25 |    1 | 5.46e-09
   0.04460 | 447 |  199.99 |  0.105 |  1.107 |  1.000 |  1.000 |     1 |     14 |    1 | 4.91e-05
   0.04470 | 448 |  199.99 |  0.105 |  1.107 |  1.000 |  1.000 |     1 |     14 |    1 | 6.21e-05
   0.04480 | 449 |  199.99 |  0.105 |  1.107 |  1.000 |  1.000 |     2 |     24 |    1 | 6.32e-09
   0.04490 | 450 |  199.99 |  0.105 |  1.107 |  1.000 |  1.000 |     2 |     25 |    1 | 8.96e-09
   0.04500 | 451 |  199.99 |  0.105 |  1.108 |  1.000 |  1.000 |     1 |     14 |    1 | 6.93e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.04510 | 452 |  199.99 |  0.105 |  1.108 |  1.000 |  1.000 |     2 |     25 |    1 | 5.68e-09
   0.04520 | 453 |  199.99 |  0.105 |  1.108 |  1.000 |  1.000 |     2 |     24 |    1 | 4.45e-09
   0.04530 | 454 |  199.99 |  0.105 |  1.108 |  1.000 |  1.000 |     2 |     26 |    1 | 4.32e-07
   0.04540 | 455 |  199.99 |  0.105 |  1.108 |  1.000 |  1.000 |     2 |     25 |    1 | 1.46e-08
   0.04550 | 456 |  199.99 |  0.105 |  1.108 |  1.000 |  1.000 |     2 |     25 |    1 | 3.84e-09
   0.04560 | 457 |  199.99 |  0.105 |  1.108 |  1.000 |  1.000 |     1 |     14 |    1 | 1.95e-05
   0.04570 | 458 |  199.98 |  0.105 |  1.108 |  1.000 |  1.000 |     1 |     14 |    1 | 3.60e-05
   0.04580 | 459 |  199.98 |  0.105 |  1.108 |  1.000 |  1.000 |     1 |     14 |    1 | 8.36e-05
   0.04590 | 460 |  199.98 |  0.105 |  1.108 |  1.000 |  1.000 |     2 |     25 |    1 | 4.32e-09
   0.04600 | 461 |  199.98 |  0.105 |  1.108 |  1.000 |  1.000 |     1 |     14 |    1 | 7.01e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.04610 | 462 |  199.98 |  0.105 |  1.108 |  1.000 |  1.000 |     2 |     25 |    1 | 2.17e-08
   0.04620 | 463 |  199.98 |  0.105 |  1.108 |  1.000 |  1.000 |     2 |     25 |    1 | 1.80e-07
   0.04630 | 464 |  199.98 |  0.105 |  1.109 |  1.000 |  1.000 |     1 |     14 |    1 | 6.39e-05
   0.04640 | 465 |  199.98 |  0.105 |  1.109 |  1.000 |  1.000 |     1 |     14 |    1 | 2.41e-05
   0.04650 | 466 |  199.98 |  0.105 |  1.109 |  1.000 |  1.000 |     1 |     14 |    1 | 3.52e-05
   0.04660 | 467 |  199.98 |  0.105 |  1.109 |  1.000 |  1.000 |     2 |     24 |    1 | 1.25e-08
   0.04670 | 468 |  199.98 |  0.105 |  1.109 |  1.000 |  1.000 |     2 |     25 |    1 | 1.67e-08
   0.04680 | 469 |  199.97 |  0.105 |  1.109 |  1.000 |  1.000 |     2 |     25 |    1 | 1.12e-08
   0.04690 | 470 |  199.97 |  0.105 |  1.109 |  1.000 |  1.000 |     1 |     14 |    1 | 2.55e-05
   0.04700 | 471 |  199.97 |  0.105 |  1.109 |  1.000 |  1.000 |     1 |     14 |    1 | 5.25e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.04710 | 472 |  199.97 |  0.105 |  1.109 |  1.000 |  1.000 |     2 |     25 |    1 | 5.23e-09
   0.04720 | 473 |  199.97 |  0.105 |  1.109 |  1.000 |  1.000 |     2 |     25 |    1 | 2.75e-09
   0.04730 | 474 |  199.97 |  0.105 |  1.109 |  1.000 |  1.000 |     1 |     14 |    1 | 7.51e-05
   0.04740 | 475 |  199.97 |  0.105 |  1.109 |  1.000 |  1.000 |     2 |     25 |    1 | 2.98e-09
   0.04750 | 476 |  199.97 |  0.105 |  1.109 |  1.000 |  1.000 |     1 |     14 |    1 | 3.02e-05
   0.04760 | 477 |  199.97 |  0.105 |  1.110 |  1.000 |  1.000 |     1 |     14 |    1 | 3.93e-05
   0.04770 | 478 |  199.97 |  0.105 |  1.110 |  1.000 |  1.000 |     1 |     14 |    1 | 5.68e-05
   0.04780 | 479 |  199.97 |  0.105 |  1.110 |  1.000 |  1.000 |     1 |     14 |    1 | 1.59e-05
   0.04790 | 480 |  199.97 |  0.105 |  1.110 |  1.000 |  1.000 |     2 |     25 |    1 | 4.71e-09
   0.04800 | 481 |  199.97 |  0.105 |  1.110 |  1.000 |  1.000 |     2 |     25 |    1 | 1.52e-08

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.04810 | 482 |  199.97 |  0.105 |  1.110 |  1.000 |  1.000 |     1 |     14 |    1 | 7.70e-05
   0.04820 | 483 |  199.97 |  0.105 |  1.110 |  1.000 |  1.000 |     1 |     14 |    1 | 9.22e-05
   0.04830 | 484 |  199.97 |  0.105 |  1.110 |  1.000 |  1.000 |     1 |     14 |    1 | 7.25e-05
   0.04840 | 485 |  199.97 |  0.105 |  1.110 |  1.000 |  1.000 |     1 |     14 |    1 | 6.44e-05
   0.04850 | 486 |  199.97 |  0.105 |  1.110 |  1.000 |  1.000 |     1 |     14 |    1 | 9.44e-05
   0.04860 | 487 |  199.97 |  0.105 |  1.110 |  1.000 |  1.000 |     2 |     24 |    1 | 2.63e-09
   0.04870 | 488 |  199.97 |  0.105 |  1.110 |  1.000 |  1.000 |     1 |     14 |    1 | 3.57e-05
   0.04880 | 489 |  199.98 |  0.105 |  1.110 |  1.000 |  1.000 |     1 |     14 |    1 | 7.94e-05
   0.04890 | 490 |  199.98 |  0.105 |  1.111 |  1.000 |  1.000 |     2 |     24 |    1 | 1.21e-08
   0.04900 | 491 |  199.98 |  0.105 |  1.111 |  1.000 |  1.000 |     2 |     25 |    1 | 1.70e-08

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.04910 | 492 |  199.98 |  0.105 |  1.111 |  1.000 |  1.000 |     1 |     14 |    1 | 9.09e-05
   0.04920 | 493 |  199.98 |  0.105 |  1.111 |  1.000 |  1.000 |     2 |     25 |    1 | 2.25e-07
   0.04930 | 494 |  199.98 |  0.105 |  1.111 |  1.000 |  1.000 |     1 |     14 |    1 | 8.40e-05
   0.04940 | 495 |  199.98 |  0.105 |  1.111 |  1.000 |  1.000 |     1 |     14 |    1 | 7.74e-05
   0.04950 | 496 |  199.98 |  0.105 |  1.111 |  1.000 |  1.000 |     2 |     25 |    1 | 4.46e-09
   0.04960 | 497 |  199.98 |  0.106 |  1.111 |  1.000 |  1.000 |     1 |     14 |    1 | 2.73e-05
   0.04970 | 498 |  199.98 |  0.106 |  1.111 |  1.000 |  1.000 |     1 |     14 |    1 | 7.42e-05
   0.04980 | 499 |  199.98 |  0.106 |  1.111 |  1.000 |  1.000 |     1 |     14 |    1 | 8.52e-05
   0.04990 | 500 |  199.98 |  0.106 |  1.111 |  1.000 |  1.000 |     1 |     14 |    1 | 7.86e-06
   0.05000 | 501 |  199.99 |  0.106 |  1.111 |  1.000 |  1.000 |     1 |     14 |    1 | 5.68e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.05010 | 502 |  199.99 |  0.106 |  1.111 |  1.000 |  1.000 |     2 |     25 |    1 | 3.20e-08
   0.05020 | 503 |  199.99 |  0.106 |  1.112 |  1.000 |  1.000 |     2 |     24 |    1 | 3.73e-08
   0.05030 | 504 |  199.99 |  0.106 |  1.112 |  1.000 |  1.000 |     1 |     14 |    1 | 4.93e-05
   0.05040 | 505 |  199.99 |  0.106 |  1.112 |  1.000 |  1.000 |     2 |     25 |    1 | 2.86e-08
   0.05050 | 506 |  199.99 |  0.106 |  1.112 |  1.000 |  1.000 |     1 |     14 |    1 | 8.39e-05
   0.05060 | 507 |  199.99 |  0.106 |  1.112 |  1.000 |  1.000 |     1 |     14 |    1 | 3.48e-05
   0.05070 | 508 |  199.99 |  0.106 |  1.112 |  1.000 |  1.000 |     1 |     14 |    1 | 7.01e-05
   0.05080 | 509 |  200.00 |  0.106 |  1.112 |  1.000 |  1.000 |     1 |     14 |    1 | 4.19e-05
   0.05090 | 510 |  200.00 |  0.106 |  1.112 |  1.000 |  1.000 |     2 |     25 |    1 | 7.46e-09
   0.05100 | 511 |  200.00 |  0.106 |  1.112 |  1.000 |  1.000 |     2 |     25 |    1 | 2.17e-07

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.05110 | 512 |  200.00 |  0.106 |  1.112 |  1.000 |  1.000 |     1 |     14 |    1 | 2.54e-05
   0.05120 | 513 |  199.99 |  0.106 |  1.112 |  1.000 |  1.000 |     2 |     25 |    1 | 7.12e-09
   0.05130 | 514 |  199.99 |  0.106 |  1.112 |  1.000 |  1.000 |     1 |     14 |    1 | 9.41e-05
   0.05140 | 515 |  199.99 |  0.106 |  1.112 |  1.000 |  1.000 |     2 |     25 |    1 | 1.52e-08
   0.05150 | 516 |  199.99 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 2.38e-05
   0.05160 | 517 |  199.99 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 4.91e-05
   0.05170 | 518 |  199.98 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 5.51e-05
   0.05180 | 519 |  199.98 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 8.48e-05
   0.05190 | 520 |  199.98 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 2.10e-05
   0.05200 | 521 |  199.97 |  0.106 |  1.113 |  1.000 |  1.000 |     2 |     25 |    1 | 5.22e-09

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.05210 | 522 |  199.97 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 2.92e-05
   0.05220 | 523 |  199.96 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 2.74e-05
   0.05230 | 524 |  199.96 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 1.59e-06
   0.05240 | 525 |  199.95 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 3.98e-06
   0.05250 | 526 |  199.95 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 9.58e-05
   0.05260 | 527 |  199.94 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 3.11e-05
   0.05270 | 528 |  199.93 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 3.81e-05
   0.05280 | 529 |  199.93 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 1.36e-05
   0.05290 | 530 |  199.92 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 2.84e-05
   0.05300 | 531 |  199.92 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 4.63e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.05310 | 532 |  199.91 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 3.30e-05
   0.05320 | 533 |  199.90 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 1.39e-05
   0.05330 | 534 |  199.89 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 3.76e-05
   0.05340 | 535 |  199.89 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 3.82e-06
   0.05350 | 536 |  199.88 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 3.33e-05
   0.05360 | 537 |  199.87 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 8.92e-06
   0.05370 | 538 |  199.87 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 2.97e-05
   0.05380 | 539 |  199.86 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 5.79e-06
   0.05390 | 540 |  199.85 |  0.106 |  1.113 |  1.000 |  1.000 |     1 |     14 |    1 | 9.81e-06
   0.05400 | 541 |  199.85 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 2.10e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.05410 | 542 |  199.84 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 1.56e-06
   0.05420 | 543 |  199.84 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 3.15e-06
   0.05430 | 544 |  199.83 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 2.28e-06
   0.05440 | 545 |  199.83 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 4.59e-06
   0.05450 | 546 |  199.82 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 2.06e-05
   0.05460 | 547 |  199.82 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 1.16e-05
   0.05470 | 548 |  199.81 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 5.26e-06
   0.05480 | 549 |  199.81 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 1.41e-05
   0.05490 | 550 |  199.80 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 5.63e-06
   0.05500 | 551 |  199.80 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 4.40e-06

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.05510 | 552 |  199.79 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 9.66e-06
   0.05520 | 553 |  199.79 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 7.62e-06
   0.05530 | 554 |  199.79 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 3.09e-06
   0.05540 | 555 |  199.78 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 2.09e-06
   0.05550 | 556 |  199.78 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 4.14e-06
   0.05560 | 557 |  199.78 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 4.46e-06
   0.05570 | 558 |  199.77 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 2.47e-06
   0.05580 | 559 |  199.77 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 8.96e-06
   0.05590 | 560 |  199.77 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 1.50e-06
   0.05600 | 561 |  199.77 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 4.29e-06

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.05610 | 562 |  199.77 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 1.02e-05
   0.05620 | 563 |  199.76 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 1.17e-06
   0.05630 | 564 |  199.76 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 2.70e-06
   0.05640 | 565 |  199.76 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 1.25e-06
   0.05650 | 566 |  199.76 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 2.61e-06
   0.05660 | 567 |  199.76 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 3.27e-06
   0.05670 | 568 |  199.76 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 4.35e-06
   0.05680 | 569 |  199.76 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 2.04e-06
   0.05690 | 570 |  199.76 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 1.28e-06
   0.05700 | 571 |  199.76 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 4.12e-07

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.05710 | 572 |  199.76 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 2.09e-06
   0.05720 | 573 |  199.76 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 2.78e-06
   0.05730 | 574 |  199.76 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 4.07e-06
   0.05740 | 575 |  199.76 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 1.91e-06
   0.05750 | 576 |  199.76 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 1.96e-06
   0.05760 | 577 |  199.76 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 2.73e-06
   0.05770 | 578 |  199.76 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 1.47e-06
   0.05780 | 579 |  199.76 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 2.93e-06
   0.05790 | 580 |  199.76 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 1.40e-06
   0.05800 | 581 |  199.76 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 9.54e-07

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.05810 | 582 |  199.76 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 8.16e-07
   0.05820 | 583 |  199.76 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 9.84e-07
   0.05830 | 584 |  199.76 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 3.03e-07
   0.05840 | 585 |  199.77 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 8.96e-07
   0.05850 | 586 |  199.77 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 4.95e-07
   0.05860 | 587 |  199.77 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 8.19e-07
   0.05870 | 588 |  199.77 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 1.56e-07
   0.05880 | 589 |  199.77 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 2.30e-07
   0.05890 | 590 |  199.77 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 3.99e-07
   0.05900 | 591 |  199.77 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 2.01e-07

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.05910 | 592 |  199.78 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 2.99e-07
   0.05920 | 593 |  199.78 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 3.72e-07
   0.05930 | 594 |  199.78 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 5.33e-07
   0.05940 | 595 |  199.78 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 4.95e-07
   0.05950 | 596 |  199.78 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 3.85e-07
   0.05960 | 597 |  199.78 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 5.93e-07
   0.05970 | 598 |  199.78 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 1.75e-07
   0.05980 | 599 |  199.79 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 1.83e-06
   0.05990 | 600 |  199.79 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 3.97e-07
   0.06000 | 601 |  199.79 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 2.46e-06

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.06010 | 602 |  199.79 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 1.40e-06
   0.06020 | 603 |  199.79 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 5.69e-07
   0.06030 | 604 |  199.79 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 1.68e-06
   0.06040 | 605 |  199.80 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 2.20e-06
   0.06050 | 606 |  199.80 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     13 |    1 | 1.95e-06
   0.06060 | 607 |  199.80 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 1.52e-06
   0.06070 | 608 |  199.80 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 7.31e-07
   0.06080 | 609 |  199.79 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 1.23e-06
   0.06090 | 610 |  199.78 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 3.02e-06
   0.06100 | 611 |  199.77 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 1.14e-06

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.06110 | 612 |  199.75 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 5.96e-06
   0.06120 | 613 |  199.73 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 2.47e-05
   0.06130 | 614 |  199.70 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 2.33e-05
   0.06140 | 615 |  199.66 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     14 |    1 | 5.60e-05
   0.06150 | 616 |  199.61 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     15 |    1 | 6.64e-05
   0.06160 | 617 |  199.55 |  0.106 |  1.114 |  1.000 |  1.000 |     1 |     15 |    1 | 9.72e-05
   0.06170 | 618 |  199.49 |  0.106 |  1.114 |  1.000 |  1.000 |     2 |     26 |    1 | 3.12e-09
   0.06180 | 619 |  199.42 |  0.106 |  1.114 |  1.000 |  1.000 |     2 |     26 |    1 | 9.25e-09
   0.06190 | 620 |  199.34 |  0.106 |  1.113 |  1.000 |  1.000 |     2 |     26 |    1 | 2.01e-08
   0.06200 | 621 |  199.25 |  0.106 |  1.113 |  1.000 |  1.000 |     2 |     26 |    1 | 2.20e-08

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.06210 | 622 |  199.16 |  0.106 |  1.113 |  1.000 |  1.000 |     2 |     26 |    1 | 6.40e-09
   0.06220 | 623 |  199.05 |  0.106 |  1.113 |  1.000 |  1.000 |     2 |     26 |    1 | 1.74e-08
   0.06230 | 624 |  198.94 |  0.106 |  1.113 |  1.000 |  1.000 |     2 |     27 |    1 | 6.77e-09
   0.06240 | 625 |  198.83 |  0.106 |  1.113 |  1.000 |  1.000 |     2 |     27 |    1 | 2.06e-08
   0.06250 | 626 |  198.70 |  0.106 |  1.113 |  1.000 |  1.000 |     2 |     27 |    1 | 8.44e-08
   0.06260 | 627 |  198.57 |  0.106 |  1.112 |  1.000 |  1.000 |     2 |     27 |    1 | 9.60e-08
   0.06270 | 628 |  198.43 |  0.106 |  1.112 |  1.000 |  1.000 |     2 |     28 |    1 | 9.18e-07
   0.06280 | 629 |  198.29 |  0.107 |  1.112 |  1.000 |  1.000 |     2 |     27 |    1 | 7.60e-08
   0.06290 | 630 |  198.14 |  0.107 |  1.112 |  1.000 |  1.000 |     2 |     27 |    1 | 8.50e-08
   0.06300 | 631 |  197.99 |  0.107 |  1.112 |  1.000 |  1.000 |     2 |     27 |    1 | 9.76e-08

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.06310 | 632 |  197.84 |  0.107 |  1.111 |  1.000 |  1.000 |     2 |     27 |    1 | 2.76e-08
   0.06320 | 633 |  197.68 |  0.107 |  1.111 |  1.000 |  1.000 |     2 |     27 |    1 | 1.01e-07
   0.06330 | 634 |  197.51 |  0.107 |  1.111 |  1.000 |  1.000 |     2 |     27 |    1 | 2.73e-07
   0.06340 | 635 |  197.35 |  0.107 |  1.111 |  1.000 |  1.000 |     2 |     27 |    1 | 9.92e-08
   0.06350 | 636 |  197.18 |  0.107 |  1.111 |  1.000 |  1.000 |     2 |     27 |    1 | 3.90e-07
   0.06360 | 637 |  197.01 |  0.107 |  1.110 |  1.000 |  1.000 |     2 |     27 |    1 | 1.81e-07
   0.06370 | 638 |  196.85 |  0.107 |  1.110 |  1.000 |  1.000 |     2 |     27 |    1 | 9.29e-08
   0.06380 | 639 |  196.67 |  0.107 |  1.110 |  1.000 |  1.000 |     2 |     27 |    1 | 7.72e-08
   0.06390 | 640 |  196.50 |  0.107 |  1.110 |  1.000 |  1.000 |     2 |     27 |    1 | 3.13e-08
   0.06400 | 641 |  196.32 |  0.107 |  1.109 |  1.000 |  1.000 |     2 |     28 |    1 | 3.71e-06

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.06410 | 642 |  196.15 |  0.107 |  1.109 |  1.000 |  1.000 |     2 |     27 |    1 | 3.24e-07
   0.06420 | 643 |  195.97 |  0.107 |  1.109 |  1.000 |  1.000 |     2 |     27 |    1 | 4.92e-07
   0.06430 | 644 |  195.79 |  0.107 |  1.109 |  1.000 |  1.000 |     2 |     27 |    1 | 5.72e-07
   0.06440 | 645 |  195.61 |  0.107 |  1.108 |  1.000 |  1.000 |     2 |     27 |    1 | 4.34e-08
   0.06450 | 646 |  195.43 |  0.107 |  1.108 |  1.000 |  1.000 |     2 |     27 |    1 | 2.34e-07
   0.06460 | 647 |  195.24 |  0.107 |  1.108 |  1.000 |  1.000 |     2 |     27 |    1 | 3.64e-07
   0.06470 | 648 |  195.05 |  0.107 |  1.108 |  1.000 |  1.000 |     2 |     27 |    1 | 2.06e-07
   0.06480 | 649 |  194.87 |  0.107 |  1.107 |  1.000 |  1.000 |     2 |     28 |    1 | 1.17e-06
   0.06490 | 650 |  194.68 |  0.107 |  1.107 |  1.000 |  1.000 |     1 |     15 |    1 | 1.80e-05
   0.06500 | 651 |  194.50 |  0.107 |  1.107 |  1.000 |  1.000 |     2 |     27 |    1 | 1.22e-07

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.06510 | 652 |  194.31 |  0.107 |  1.107 |  1.000 |  1.000 |     2 |     27 |    1 | 1.36e-07
   0.06520 | 653 |  194.13 |  0.107 |  1.106 |  1.000 |  1.000 |     2 |     27 |    1 | 1.50e-07
   0.06530 | 654 |  193.94 |  0.107 |  1.106 |  1.000 |  1.000 |     2 |     27 |    1 | 1.10e-07
   0.06540 | 655 |  193.75 |  0.107 |  1.106 |  1.000 |  1.000 |     2 |     27 |    1 | 4.43e-07
   0.06550 | 656 |  193.56 |  0.107 |  1.106 |  1.000 |  1.000 |     2 |     27 |    1 | 4.82e-08
   0.06560 | 657 |  193.37 |  0.107 |  1.105 |  1.000 |  1.000 |     2 |     28 |    1 | 1.43e-06
   0.06570 | 658 |  193.19 |  0.107 |  1.105 |  1.000 |  1.000 |     2 |     26 |    1 | 3.21e-07
   0.06580 | 659 |  193.00 |  0.107 |  1.105 |  1.000 |  1.000 |     2 |     26 |    1 | 3.67e-08
   0.06590 | 660 |  192.82 |  0.107 |  1.105 |  1.000 |  1.000 |     2 |     25 |    1 | 1.48e-08
   0.06600 | 661 |  192.63 |  0.107 |  1.105 |  1.000 |  1.000 |     2 |     25 |    1 | 2.66e-08

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.06610 | 662 |  192.45 |  0.107 |  1.104 |  1.000 |  1.000 |     2 |     25 |    1 | 2.22e-07
   0.06620 | 663 |  192.26 |  0.107 |  1.104 |  1.000 |  1.000 |     2 |     25 |    1 | 1.44e-07
   0.06630 | 664 |  192.08 |  0.107 |  1.104 |  1.000 |  1.000 |     2 |     25 |    1 | 1.38e-07
   0.06640 | 665 |  191.89 |  0.107 |  1.104 |  1.000 |  1.000 |     2 |     25 |    1 | 1.32e-08
   0.06650 | 666 |  191.71 |  0.107 |  1.104 |  1.000 |  1.000 |     2 |     25 |    1 | 2.26e-07
   0.06660 | 667 |  191.52 |  0.106 |  1.103 |  1.000 |  1.000 |     2 |     25 |    1 | 7.06e-07
   0.06670 | 668 |  191.34 |  0.106 |  1.103 |  1.000 |  1.000 |     2 |     25 |    1 | 3.12e-07
   0.06680 | 669 |  191.15 |  0.106 |  1.103 |  1.000 |  1.000 |     2 |     25 |    1 | 3.76e-08
   0.06690 | 670 |  190.96 |  0.106 |  1.103 |  1.000 |  1.000 |     2 |     25 |    1 | 1.74e-08
   0.06700 | 671 |  190.77 |  0.106 |  1.103 |  1.000 |  1.000 |     2 |     25 |    1 | 1.39e-07

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.06710 | 672 |  190.58 |  0.106 |  1.103 |  1.000 |  1.000 |     2 |     25 |    1 | 1.28e-07
   0.06720 | 673 |  190.38 |  0.106 |  1.102 |  1.000 |  1.000 |     2 |     25 |    1 | 1.01e-06
   0.06730 | 674 |  190.19 |  0.106 |  1.102 |  1.000 |  1.000 |     2 |     25 |    1 | 1.60e-07
   0.06740 | 675 |  190.00 |  0.106 |  1.102 |  1.000 |  1.000 |     2 |     25 |    1 | 2.17e-07
   0.06750 | 676 |  189.80 |  0.106 |  1.102 |  1.000 |  1.000 |     2 |     25 |    1 | 5.98e-08
   0.06760 | 677 |  189.60 |  0.106 |  1.102 |  1.000 |  1.000 |     2 |     25 |    1 | 1.57e-07
   0.06770 | 678 |  189.40 |  0.106 |  1.101 |  1.000 |  1.000 |     2 |     25 |    1 | 6.78e-08
   0.06780 | 679 |  189.20 |  0.106 |  1.101 |  1.000 |  1.000 |     2 |     25 |    1 | 8.61e-08
   0.06790 | 680 |  189.00 |  0.106 |  1.101 |  1.000 |  1.000 |     2 |     25 |    1 | 1.98e-07
   0.06800 | 681 |  188.80 |  0.105 |  1.101 |  1.000 |  1.000 |     2 |     26 |    1 | 2.00e-06

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.06810 | 682 |  188.60 |  0.105 |  1.101 |  1.000 |  1.000 |     2 |     25 |    1 | 1.11e-07
   0.06820 | 683 |  188.40 |  0.105 |  1.101 |  1.000 |  1.000 |     2 |     25 |    1 | 1.71e-07
   0.06830 | 684 |  188.20 |  0.105 |  1.100 |  1.000 |  1.000 |     2 |     25 |    1 | 1.05e-07
   0.06840 | 685 |  188.00 |  0.105 |  1.100 |  1.000 |  1.000 |     2 |     25 |    1 | 3.35e-07
   0.06850 | 686 |  187.80 |  0.105 |  1.100 |  1.000 |  1.000 |     2 |     25 |    1 | 5.64e-07
   0.06860 | 687 |  187.60 |  0.105 |  1.100 |  1.000 |  1.000 |     2 |     25 |    1 | 5.61e-08
   0.06870 | 688 |  187.40 |  0.105 |  1.100 |  1.000 |  1.000 |     2 |     25 |    1 | 2.03e-08
   0.06880 | 689 |  187.20 |  0.105 |  1.100 |  1.000 |  1.000 |     2 |     25 |    1 | 3.16e-07
   0.06890 | 690 |  187.00 |  0.105 |  1.100 |  1.000 |  1.000 |     2 |     25 |    1 | 7.09e-08
   0.06900 | 691 |  186.80 |  0.105 |  1.099 |  1.000 |  1.000 |     2 |     25 |    1 | 6.38e-08

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.06910 | 692 |  186.60 |  0.105 |  1.099 |  1.000 |  1.000 |     2 |     25 |    1 | 7.20e-07
   0.06920 | 693 |  186.40 |  0.105 |  1.099 |  1.000 |  1.000 |     2 |     25 |    1 | 1.16e-07
   0.06930 | 694 |  186.21 |  0.105 |  1.099 |  1.000 |  1.000 |     2 |     25 |    1 | 3.40e-07
   0.06940 | 695 |  186.01 |  0.105 |  1.099 |  1.000 |  1.000 |     2 |     25 |    1 | 5.37e-08
   0.06950 | 696 |  185.81 |  0.105 |  1.099 |  1.000 |  1.000 |     2 |     25 |    1 | 3.14e-08
   0.06960 | 697 |  185.62 |  0.104 |  1.099 |  1.000 |  1.000 |     1 |     14 |    1 | 2.09e-05
   0.06970 | 698 |  185.42 |  0.104 |  1.098 |  1.000 |  1.000 |     2 |     25 |    1 | 2.79e-08
   0.06980 | 699 |  185.22 |  0.104 |  1.098 |  1.000 |  1.000 |     2 |     25 |    1 | 3.80e-08
   0.06990 | 700 |  185.03 |  0.104 |  1.098 |  1.000 |  1.000 |     2 |     25 |    1 | 1.02e-07
   0.07000 | 701 |  184.83 |  0.104 |  1.098 |  1.000 |  1.000 |     2 |     25 |    1 | 1.33e-07

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.07010 | 702 |  184.72 |  0.104 |  1.098 |  1.005 |  1.000 |     2 |     25 |    1 | 4.58e-08
   0.07020 | 703 |  184.60 |  0.104 |  1.099 |  1.010 |  1.000 |     1 |     14 |    1 | 6.70e-05
   0.07030 | 704 |  184.48 |  0.104 |  1.099 |  1.015 |  1.000 |     2 |     25 |    1 | 8.26e-07
   0.07040 | 705 |  184.37 |  0.104 |  1.099 |  1.020 |  1.000 |     2 |     25 |    1 | 2.89e-08
   0.07050 | 706 |  184.25 |  0.104 |  1.099 |  1.025 |  1.000 |     2 |     25 |    1 | 5.80e-09
   0.07060 | 707 |  184.13 |  0.104 |  1.100 |  1.030 |  1.000 |     2 |     25 |    1 | 3.80e-08
   0.07070 | 708 |  184.01 |  0.103 |  1.100 |  1.035 |  1.000 |     2 |     25 |    1 | 5.48e-09
   0.07080 | 709 |  183.90 |  0.103 |  1.100 |  1.040 |  1.000 |     2 |     25 |    1 | 4.49e-08
   0.07090 | 710 |  183.78 |  0.103 |  1.100 |  1.045 |  1.000 |     2 |     25 |    1 | 1.25e-08
   0.07100 | 711 |  183.67 |  0.103 |  1.101 |  1.050 |  1.000 |     2 |     25 |    1 | 7.49e-08

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.07110 | 712 |  183.55 |  0.103 |  1.101 |  1.055 |  1.000 |     2 |     25 |    1 | 2.86e-08
   0.07120 | 713 |  183.44 |  0.103 |  1.101 |  1.060 |  1.000 |     2 |     25 |    1 | 7.86e-09
   0.07130 | 714 |  183.34 |  0.103 |  1.101 |  1.065 |  1.000 |     2 |     25 |    1 | 6.78e-08
   0.07140 | 715 |  183.23 |  0.103 |  1.102 |  1.070 |  1.000 |     2 |     25 |    1 | 7.53e-09
   0.07150 | 716 |  183.13 |  0.103 |  1.102 |  1.075 |  1.000 |     2 |     25 |    1 | 9.86e-09
   0.07160 | 717 |  183.03 |  0.103 |  1.102 |  1.080 |  1.000 |     2 |     25 |    1 | 2.53e-08
   0.07170 | 718 |  182.93 |  0.103 |  1.103 |  1.085 |  1.000 |     2 |     25 |    1 | 1.05e-08
   0.07180 | 719 |  182.83 |  0.103 |  1.103 |  1.090 |  1.000 |     2 |     25 |    1 | 1.72e-08
   0.07190 | 720 |  182.74 |  0.102 |  1.103 |  1.095 |  1.000 |     2 |     25 |    1 | 2.48e-09
   0.07200 | 721 |  182.65 |  0.102 |  1.103 |  1.100 |  1.000 |     2 |     25 |    1 | 3.11e-08

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.07210 | 722 |  182.56 |  0.102 |  1.104 |  1.105 |  1.000 |     2 |     25 |    1 | 3.27e-08
   0.07220 | 723 |  182.48 |  0.102 |  1.104 |  1.110 |  1.000 |     2 |     25 |    1 | 2.50e-08
   0.07230 | 724 |  182.40 |  0.102 |  1.104 |  1.115 |  1.000 |     2 |     25 |    1 | 4.46e-07
   0.07240 | 725 |  182.32 |  0.102 |  1.105 |  1.120 |  1.000 |     2 |     25 |    1 | 1.39e-08
   0.07250 | 726 |  182.24 |  0.102 |  1.105 |  1.125 |  1.000 |     2 |     25 |    1 | 1.77e-08
   0.07260 | 727 |  182.16 |  0.102 |  1.105 |  1.130 |  1.000 |     2 |     25 |    1 | 1.05e-08
   0.07270 | 728 |  182.09 |  0.102 |  1.105 |  1.135 |  1.000 |     2 |     25 |    1 | 3.80e-08
   0.07280 | 729 |  182.02 |  0.102 |  1.106 |  1.140 |  1.000 |     2 |     25 |    1 | 1.15e-08
   0.07290 | 730 |  181.95 |  0.102 |  1.106 |  1.145 |  1.000 |     2 |     25 |    1 | 4.17e-08
   0.07300 | 731 |  181.89 |  0.102 |  1.106 |  1.150 |  1.000 |     2 |     25 |    1 | 1.24e-08

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.07310 | 732 |  181.82 |  0.102 |  1.107 |  1.155 |  1.000 |     2 |     25 |    1 | 4.57e-08
   0.07320 | 733 |  181.76 |  0.102 |  1.107 |  1.160 |  1.000 |     2 |     25 |    1 | 5.73e-09
   0.07330 | 734 |  181.70 |  0.102 |  1.107 |  1.165 |  1.000 |     2 |     25 |    1 | 1.68e-08
   0.07340 | 735 |  181.64 |  0.102 |  1.107 |  1.170 |  1.000 |     2 |     25 |    1 | 1.15e-08
   0.07350 | 736 |  181.58 |  0.102 |  1.108 |  1.175 |  1.000 |     2 |     25 |    1 | 5.05e-09
   0.07360 | 737 |  181.53 |  0.102 |  1.108 |  1.180 |  1.000 |     2 |     25 |    1 | 3.70e-07
   0.07370 | 738 |  181.47 |  0.102 |  1.108 |  1.185 |  1.000 |     2 |     25 |    1 | 2.01e-08
   0.07380 | 739 |  181.42 |  0.102 |  1.109 |  1.190 |  1.000 |     2 |     25 |    1 | 1.02e-08
   0.07390 | 740 |  181.37 |  0.102 |  1.109 |  1.195 |  1.000 |     2 |     25 |    1 | 6.49e-08
   0.07400 | 741 |  181.32 |  0.102 |  1.109 |  1.200 |  1.000 |     2 |     25 |    1 | 8.81e-08

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.07410 | 742 |  181.27 |  0.102 |  1.110 |  1.205 |  1.000 |     1 |     14 |    1 | 6.84e-05
   0.07420 | 743 |  181.22 |  0.102 |  1.110 |  1.210 |  1.000 |     2 |     25 |    1 | 2.14e-08
   0.07430 | 744 |  181.17 |  0.102 |  1.110 |  1.215 |  1.000 |     2 |     25 |    1 | 1.71e-08
   0.07440 | 745 |  181.13 |  0.102 |  1.110 |  1.220 |  1.000 |     2 |     25 |    1 | 1.90e-08
   0.07450 | 746 |  181.08 |  0.102 |  1.111 |  1.225 |  1.000 |     2 |     25 |    1 | 3.75e-08
   0.07460 | 747 |  181.04 |  0.102 |  1.111 |  1.230 |  1.000 |     2 |     25 |    1 | 2.59e-08
   0.07470 | 748 |  181.00 |  0.102 |  1.111 |  1.235 |  1.000 |     2 |     26 |    1 | 1.22e-08
   0.07480 | 749 |  180.96 |  0.102 |  1.111 |  1.240 |  1.000 |     2 |     27 |    1 | 1.76e-07
   0.07490 | 750 |  180.92 |  0.102 |  1.112 |  1.245 |  1.000 |     2 |     27 |    1 | 1.94e-07
   0.07500 | 751 |  180.88 |  0.102 |  1.112 |  1.250 |  1.000 |     2 |     28 |    1 | 4.21e-06

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.07510 | 752 |  180.84 |  0.102 |  1.112 |  1.255 |  1.000 |     2 |     26 |    1 | 3.36e-08
   0.07520 | 753 |  180.80 |  0.102 |  1.113 |  1.260 |  1.000 |     2 |     26 |    1 | 6.19e-09
   0.07530 | 754 |  180.77 |  0.102 |  1.113 |  1.265 |  1.000 |     2 |     26 |    1 | 1.28e-08
   0.07540 | 755 |  180.73 |  0.102 |  1.113 |  1.270 |  1.000 |     2 |     26 |    1 | 8.95e-09
   0.07550 | 756 |  180.70 |  0.102 |  1.113 |  1.275 |  1.000 |     2 |     26 |    1 | 6.35e-09
   0.07560 | 757 |  180.66 |  0.102 |  1.114 |  1.280 |  1.000 |     2 |     26 |    1 | 1.73e-09
   0.07570 | 758 |  180.63 |  0.102 |  1.114 |  1.285 |  1.000 |     2 |     26 |    1 | 1.23e-09
   0.07580 | 759 |  180.60 |  0.102 |  1.114 |  1.290 |  1.000 |     2 |     26 |    1 | 5.94e-09
   0.07590 | 760 |  180.57 |  0.102 |  1.114 |  1.295 |  1.000 |     2 |     26 |    1 | 3.15e-08
   0.07600 | 761 |  180.54 |  0.102 |  1.115 |  1.300 |  1.000 |     2 |     26 |    1 | 9.78e-09

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.07610 | 762 |  180.51 |  0.102 |  1.115 |  1.305 |  1.000 |     2 |     27 |    1 | 8.15e-08
   0.07620 | 763 |  180.48 |  0.102 |  1.115 |  1.310 |  1.000 |     2 |     27 |    1 | 4.23e-09
   0.07630 | 764 |  180.45 |  0.102 |  1.115 |  1.315 |  1.000 |     2 |     27 |    1 | 2.77e-08
   0.07640 | 765 |  180.43 |  0.102 |  1.115 |  1.320 |  1.000 |     2 |     27 |    1 | 3.91e-08
   0.07650 | 766 |  180.40 |  0.102 |  1.116 |  1.325 |  1.000 |     2 |     27 |    1 | 4.01e-08
   0.07660 | 767 |  180.38 |  0.102 |  1.116 |  1.330 |  1.000 |     2 |     27 |    1 | 1.07e-08
   0.07670 | 768 |  180.35 |  0.102 |  1.116 |  1.335 |  1.000 |     2 |     27 |    1 | 1.17e-08
   0.07680 | 769 |  180.33 |  0.102 |  1.116 |  1.340 |  1.000 |     2 |     27 |    1 | 6.64e-09
   0.07690 | 770 |  180.31 |  0.102 |  1.117 |  1.345 |  1.000 |     2 |     27 |    1 | 5.53e-09
   0.07700 | 771 |  180.29 |  0.102 |  1.117 |  1.350 |  1.000 |     2 |     27 |    1 | 4.43e-07

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.07710 | 772 |  180.26 |  0.102 |  1.117 |  1.355 |  1.000 |     2 |     27 |    1 | 1.73e-09
   0.07720 | 773 |  180.24 |  0.102 |  1.117 |  1.360 |  1.000 |     2 |     27 |    1 | 8.22e-09
   0.07730 | 774 |  180.22 |  0.102 |  1.117 |  1.365 |  1.000 |     2 |     27 |    1 | 5.33e-09
   0.07740 | 775 |  180.20 |  0.102 |  1.118 |  1.370 |  1.000 |     2 |     27 |    1 | 3.79e-09
   0.07750 | 776 |  180.18 |  0.102 |  1.118 |  1.375 |  1.000 |     2 |     27 |    1 | 4.32e-09
   0.07760 | 777 |  180.17 |  0.102 |  1.118 |  1.380 |  1.000 |     2 |     27 |    1 | 1.63e-09
   0.07770 | 778 |  180.15 |  0.102 |  1.118 |  1.385 |  1.000 |     2 |     27 |    1 | 1.05e-08
   0.07780 | 779 |  180.13 |  0.102 |  1.118 |  1.390 |  1.000 |     2 |     27 |    1 | 1.64e-09
   0.07790 | 780 |  180.11 |  0.102 |  1.119 |  1.395 |  1.000 |     2 |     27 |    1 | 5.56e-09
   0.07800 | 781 |  180.10 |  0.102 |  1.119 |  1.400 |  1.000 |     2 |     28 |    1 | 9.49e-10

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.07810 | 782 |  180.08 |  0.102 |  1.119 |  1.405 |  1.000 |     2 |     28 |    1 | 1.80e-09
   0.07820 | 783 |  180.07 |  0.102 |  1.119 |  1.410 |  1.000 |     2 |     28 |    1 | 1.64e-09
   0.07830 | 784 |  180.05 |  0.102 |  1.119 |  1.415 |  1.000 |     2 |     28 |    1 | 3.17e-09
   0.07840 | 785 |  180.04 |  0.102 |  1.120 |  1.420 |  1.000 |     2 |     28 |    1 | 4.40e-07
   0.07850 | 786 |  180.03 |  0.102 |  1.120 |  1.425 |  1.000 |     2 |     27 |    1 | 1.13e-08
   0.07860 | 787 |  180.01 |  0.102 |  1.120 |  1.430 |  1.000 |     2 |     27 |    1 | 4.11e-09
   0.07870 | 788 |  180.00 |  0.102 |  1.120 |  1.435 |  1.000 |     2 |     27 |    1 | 1.54e-09
   0.07880 | 789 |  179.99 |  0.102 |  1.120 |  1.440 |  1.000 |     2 |     27 |    1 | 6.03e-09
   0.07890 | 790 |  179.98 |  0.102 |  1.120 |  1.445 |  1.000 |     2 |     28 |    1 | 2.37e-09
   0.07900 | 791 |  179.97 |  0.102 |  1.121 |  1.450 |  1.000 |     2 |     28 |    1 | 1.13e-09

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.07910 | 792 |  179.95 |  0.102 |  1.121 |  1.455 |  1.000 |     2 |     29 |    1 | 1.18e-06
   0.07920 | 793 |  179.94 |  0.102 |  1.121 |  1.460 |  1.000 |     2 |     28 |    1 | 3.34e-09
   0.07930 | 794 |  179.94 |  0.102 |  1.121 |  1.465 |  1.000 |     2 |     28 |    1 | 1.21e-08
   0.07940 | 795 |  179.93 |  0.102 |  1.121 |  1.470 |  1.000 |     2 |     28 |    1 | 6.39e-09
   0.07950 | 796 |  179.92 |  0.102 |  1.121 |  1.475 |  1.000 |     2 |     28 |    1 | 1.39e-08
   0.07960 | 797 |  179.91 |  0.102 |  1.122 |  1.480 |  1.000 |     2 |     28 |    1 | 1.14e-08
   0.07970 | 798 |  179.91 |  0.102 |  1.122 |  1.485 |  1.000 |     2 |     28 |    1 | 1.08e-09
   0.07980 | 799 |  179.90 |  0.102 |  1.122 |  1.490 |  1.000 |     2 |     28 |    1 | 6.21e-09
   0.07990 | 800 |  179.89 |  0.102 |  1.122 |  1.495 |  1.000 |     2 |     28 |    1 | 4.62e-09
   0.08000 | 801 |  179.89 |  0.102 |  1.122 |  1.500 |  1.000 |     2 |     28 |    1 | 1.98e-08

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.08010 | 802 |  179.88 |  0.102 |  1.122 |  1.505 |  1.000 |     2 |     28 |    1 | 6.90e-09
   0.08020 | 803 |  179.88 |  0.102 |  1.122 |  1.510 |  1.000 |     2 |     28 |    1 | 6.49e-09
   0.08030 | 804 |  179.87 |  0.102 |  1.123 |  1.515 |  1.000 |     2 |     28 |    1 | 5.46e-09
   0.08040 | 805 |  179.87 |  0.102 |  1.123 |  1.520 |  1.000 |     2 |     28 |    1 | 4.58e-09
   0.08050 | 806 |  179.86 |  0.102 |  1.123 |  1.525 |  1.000 |     2 |     28 |    1 | 4.38e-09
   0.08060 | 807 |  179.86 |  0.102 |  1.123 |  1.530 |  1.000 |     2 |     28 |    1 | 3.01e-09
   0.08070 | 808 |  179.86 |  0.102 |  1.123 |  1.535 |  1.000 |     2 |     27 |    1 | 3.68e-09
   0.08080 | 809 |  179.85 |  0.102 |  1.123 |  1.540 |  1.000 |     2 |     28 |    1 | 1.62e-09
   0.08090 | 810 |  179.85 |  0.102 |  1.123 |  1.545 |  1.000 |     2 |     28 |    1 | 2.18e-09
   0.08100 | 811 |  179.84 |  0.102 |  1.123 |  1.550 |  1.000 |     2 |     29 |    1 | 3.51e-10

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.08110 | 812 |  179.84 |  0.102 |  1.124 |  1.555 |  1.000 |     2 |     28 |    1 | 4.64e-09
   0.08120 | 813 |  179.84 |  0.102 |  1.124 |  1.560 |  1.000 |     2 |     28 |    1 | 3.00e-09
   0.08130 | 814 |  179.83 |  0.102 |  1.124 |  1.565 |  1.000 |     2 |     28 |    1 | 4.44e-10
   0.08140 | 815 |  179.83 |  0.103 |  1.124 |  1.570 |  1.000 |     2 |     28 |    1 | 3.92e-09
   0.08150 | 816 |  179.83 |  0.103 |  1.124 |  1.575 |  1.000 |     2 |     29 |    1 | 7.60e-10
   0.08160 | 817 |  179.82 |  0.103 |  1.124 |  1.580 |  1.000 |     1 |     16 |    1 | 9.87e-05
   0.08170 | 818 |  179.82 |  0.103 |  1.124 |  1.585 |  1.000 |     1 |     16 |    1 | 8.35e-05
   0.08180 | 819 |  179.82 |  0.103 |  1.124 |  1.590 |  1.000 |     1 |     16 |    1 | 7.74e-05
   0.08190 | 820 |  179.81 |  0.103 |  1.124 |  1.595 |  1.000 |     2 |     29 |    1 | 1.51e-09
   0.08200 | 821 |  179.81 |  0.103 |  1.125 |  1.600 |  1.000 |     2 |     29 |    1 | 1.89e-09

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.08210 | 822 |  179.81 |  0.103 |  1.125 |  1.605 |  1.000 |     1 |     16 |    1 | 7.80e-05
   0.08220 | 823 |  179.81 |  0.103 |  1.125 |  1.610 |  1.000 |     1 |     16 |    1 | 9.98e-05
   0.08230 | 824 |  179.81 |  0.103 |  1.125 |  1.615 |  1.000 |     1 |     16 |    1 | 9.17e-05
   0.08240 | 825 |  179.81 |  0.103 |  1.125 |  1.620 |  1.000 |     1 |     16 |    1 | 9.67e-05
   0.08250 | 826 |  179.81 |  0.103 |  1.125 |  1.625 |  1.000 |     1 |     16 |    1 | 6.83e-05
   0.08260 | 827 |  179.81 |  0.103 |  1.125 |  1.630 |  1.000 |     1 |     16 |    1 | 8.73e-05
   0.08270 | 828 |  179.81 |  0.103 |  1.125 |  1.635 |  1.000 |     1 |     16 |    1 | 8.01e-05
   0.08280 | 829 |  179.81 |  0.103 |  1.125 |  1.640 |  1.000 |     1 |     16 |    1 | 7.82e-05
   0.08290 | 830 |  179.82 |  0.103 |  1.125 |  1.645 |  1.000 |     1 |     16 |    1 | 9.16e-05
   0.08300 | 831 |  179.82 |  0.103 |  1.126 |  1.650 |  1.000 |     1 |     16 |    1 | 7.35e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.08310 | 832 |  179.82 |  0.103 |  1.126 |  1.655 |  1.000 |     1 |     16 |    1 | 9.90e-05
   0.08320 | 833 |  179.83 |  0.103 |  1.126 |  1.660 |  1.000 |     1 |     16 |    1 | 8.87e-05
   0.08330 | 834 |  179.83 |  0.103 |  1.126 |  1.665 |  1.000 |     1 |     16 |    1 | 7.31e-05
   0.08340 | 835 |  179.83 |  0.103 |  1.126 |  1.670 |  1.000 |     1 |     17 |    1 | 7.89e-05
   0.08350 | 836 |  179.84 |  0.103 |  1.126 |  1.675 |  1.000 |     1 |     16 |    1 | 9.41e-05
   0.08360 | 837 |  179.85 |  0.103 |  1.126 |  1.680 |  1.000 |     2 |     29 |    1 | 9.48e-09
   0.08370 | 838 |  179.86 |  0.103 |  1.126 |  1.685 |  1.000 |     1 |     16 |    1 | 9.36e-05
   0.08380 | 839 |  179.87 |  0.103 |  1.126 |  1.690 |  1.000 |     2 |     29 |    1 | 1.63e-08
   0.08390 | 840 |  179.88 |  0.103 |  1.126 |  1.695 |  1.000 |     2 |     30 |    1 | 1.88e-08
   0.08400 | 841 |  179.89 |  0.103 |  1.126 |  1.700 |  1.000 |     2 |     30 |    1 | 1.01e-08

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.08410 | 842 |  179.90 |  0.103 |  1.127 |  1.705 |  1.000 |     1 |     17 |    1 | 9.58e-05
   0.08420 | 843 |  179.90 |  0.103 |  1.127 |  1.710 |  1.000 |     1 |     17 |    1 | 7.99e-05
   0.08430 | 844 |  179.91 |  0.103 |  1.127 |  1.715 |  1.000 |     1 |     17 |    1 | 9.74e-05
   0.08440 | 845 |  179.92 |  0.103 |  1.127 |  1.720 |  1.000 |     1 |     17 |    1 | 7.38e-05
   0.08450 | 846 |  179.93 |  0.103 |  1.127 |  1.725 |  1.000 |     1 |     17 |    1 | 8.80e-05
   0.08460 | 847 |  179.94 |  0.103 |  1.127 |  1.730 |  1.000 |     1 |     17 |    1 | 8.52e-05
   0.08470 | 848 |  179.94 |  0.103 |  1.127 |  1.735 |  1.000 |     1 |     17 |    1 | 8.25e-05
   0.08480 | 849 |  179.95 |  0.103 |  1.127 |  1.740 |  1.000 |     1 |     17 |    1 | 8.00e-05
   0.08490 | 850 |  179.96 |  0.103 |  1.127 |  1.745 |  1.000 |     1 |     17 |    1 | 6.42e-05
   0.08500 | 851 |  179.96 |  0.103 |  1.127 |  1.750 |  1.000 |     1 |     17 |    1 | 4.60e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.08510 | 852 |  179.97 |  0.103 |  1.127 |  1.755 |  1.000 |     1 |     17 |    1 | 5.24e-05
   0.08520 | 853 |  179.97 |  0.103 |  1.127 |  1.760 |  1.000 |     1 |     17 |    1 | 5.76e-05
   0.08530 | 854 |  179.98 |  0.103 |  1.127 |  1.765 |  1.000 |     1 |     17 |    1 | 4.20e-05
   0.08540 | 855 |  179.98 |  0.103 |  1.128 |  1.770 |  1.000 |     1 |     17 |    1 | 4.57e-05
   0.08550 | 856 |  179.98 |  0.103 |  1.128 |  1.775 |  1.000 |     1 |     17 |    1 | 4.88e-05
   0.08560 | 857 |  179.99 |  0.103 |  1.128 |  1.780 |  1.000 |     1 |     17 |    1 | 3.94e-05
   0.08570 | 858 |  179.99 |  0.103 |  1.128 |  1.785 |  1.000 |     1 |     17 |    1 | 4.16e-05
   0.08580 | 859 |  179.99 |  0.103 |  1.128 |  1.790 |  1.000 |     1 |     17 |    1 | 3.66e-05
   0.08590 | 860 |  180.00 |  0.103 |  1.128 |  1.795 |  1.000 |     1 |     17 |    1 | 3.72e-05
   0.08600 | 861 |  180.00 |  0.103 |  1.128 |  1.800 |  1.000 |     1 |     17 |    1 | 3.29e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.08610 | 862 |  180.00 |  0.103 |  1.128 |  1.805 |  1.000 |     1 |     17 |    1 | 2.58e-05
   0.08620 | 863 |  180.00 |  0.103 |  1.128 |  1.810 |  1.000 |     1 |     17 |    1 | 2.44e-05
   0.08630 | 864 |  180.01 |  0.103 |  1.128 |  1.815 |  1.000 |     1 |     17 |    1 | 2.33e-05
   0.08640 | 865 |  180.01 |  0.103 |  1.128 |  1.820 |  1.000 |     1 |     17 |    1 | 2.19e-05
   0.08650 | 866 |  180.01 |  0.103 |  1.128 |  1.825 |  1.000 |     1 |     17 |    1 | 1.32e-05
   0.08660 | 867 |  180.01 |  0.103 |  1.128 |  1.830 |  1.000 |     1 |     17 |    1 | 2.70e-05
   0.08670 | 868 |  180.01 |  0.103 |  1.128 |  1.835 |  1.000 |     1 |     17 |    1 | 1.27e-05
   0.08680 | 869 |  180.01 |  0.103 |  1.128 |  1.840 |  1.000 |     1 |     17 |    1 | 1.10e-05
   0.08690 | 870 |  180.01 |  0.103 |  1.128 |  1.845 |  1.000 |     1 |     17 |    1 | 1.06e-05
   0.08700 | 871 |  180.01 |  0.103 |  1.128 |  1.850 |  1.000 |     1 |     18 |    1 | 1.63e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.08710 | 872 |  180.01 |  0.103 |  1.128 |  1.855 |  1.000 |     1 |     18 |    1 | 1.38e-05
   0.08720 | 873 |  180.00 |  0.102 |  1.128 |  1.860 |  1.000 |     1 |     18 |    1 | 9.28e-06
   0.08730 | 874 |  180.00 |  0.102 |  1.128 |  1.865 |  1.000 |     1 |     17 |    1 | 2.03e-05
   0.08740 | 875 |  180.00 |  0.102 |  1.128 |  1.870 |  1.000 |     1 |     17 |    1 | 1.13e-05
   0.08750 | 876 |  179.99 |  0.102 |  1.128 |  1.875 |  1.000 |     1 |     17 |    1 | 1.67e-05
   0.08760 | 877 |  179.99 |  0.102 |  1.129 |  1.880 |  1.000 |     1 |     17 |    1 | 1.76e-05
   0.08770 | 878 |  179.99 |  0.102 |  1.129 |  1.885 |  1.000 |     1 |     17 |    1 | 1.44e-05
   0.08780 | 879 |  179.99 |  0.102 |  1.129 |  1.890 |  1.000 |     1 |     17 |    1 | 2.21e-05
   0.08790 | 880 |  179.98 |  0.102 |  1.129 |  1.895 |  1.000 |     1 |     17 |    1 | 9.91e-06
   0.08800 | 881 |  179.98 |  0.102 |  1.129 |  1.900 |  1.000 |     1 |     17 |    1 | 1.93e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.08810 | 882 |  179.98 |  0.102 |  1.129 |  1.905 |  1.000 |     1 |     18 |    1 | 2.16e-05
   0.08820 | 883 |  179.98 |  0.102 |  1.129 |  1.910 |  1.000 |     1 |     18 |    1 | 1.73e-05
   0.08830 | 884 |  179.97 |  0.102 |  1.129 |  1.915 |  1.000 |     1 |     18 |    1 | 1.52e-05
   0.08840 | 885 |  179.97 |  0.102 |  1.129 |  1.920 |  1.000 |     1 |     18 |    1 | 2.55e-05
   0.08850 | 886 |  179.96 |  0.102 |  1.129 |  1.925 |  1.000 |     2 |     33 |    1 | 1.29e-06
   0.08860 | 887 |  179.97 |  0.102 |  1.129 |  1.930 |  1.000 |     1 |     18 |    1 | 2.61e-05
   0.08870 | 888 |  179.97 |  0.102 |  1.129 |  1.935 |  1.000 |     1 |     18 |    1 | 5.46e-05
   0.08880 | 889 |  179.97 |  0.102 |  1.129 |  1.940 |  1.000 |     1 |     18 |    1 | 5.64e-05
   0.08890 | 890 |  179.97 |  0.102 |  1.129 |  1.945 |  1.000 |     1 |     18 |    1 | 3.50e-05
   0.08900 | 891 |  179.97 |  0.102 |  1.129 |  1.950 |  1.000 |     1 |     18 |    1 | 1.88e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.08910 | 892 |  179.97 |  0.102 |  1.129 |  1.955 |  1.000 |     1 |     18 |    1 | 1.95e-05
   0.08920 | 893 |  179.98 |  0.102 |  1.129 |  1.960 |  1.000 |     1 |     18 |    1 | 1.66e-05
   0.08930 | 894 |  179.99 |  0.102 |  1.129 |  1.965 |  1.000 |     1 |     18 |    1 | 2.99e-05
   0.08940 | 895 |  180.00 |  0.102 |  1.129 |  1.970 |  1.000 |     1 |     18 |    1 | 4.99e-05
   0.08950 | 896 |  180.00 |  0.102 |  1.129 |  1.975 |  1.000 |     1 |     18 |    1 | 3.03e-05
   0.08960 | 897 |  180.01 |  0.102 |  1.129 |  1.980 |  1.000 |     1 |     18 |    1 | 2.07e-05
   0.08970 | 898 |  180.02 |  0.102 |  1.129 |  1.985 |  1.000 |     1 |     18 |    1 | 6.50e-05
   0.08980 | 899 |  180.02 |  0.102 |  1.129 |  1.990 |  1.000 |     1 |     18 |    1 | 1.28e-05
   0.08990 | 900 |  180.03 |  0.102 |  1.129 |  1.995 |  1.000 |     1 |     18 |    1 | 1.78e-05
   0.09000 | 901 |  180.03 |  0.102 |  1.129 |  2.000 |  1.000 |     1 |     18 |    1 | 5.95e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.09010 | 902 |  179.96 |  0.102 |  1.129 |  2.000 |  1.000 |     3 |     51 |    1 | 4.22e-08
   0.09020 | 903 |  179.88 |  0.102 |  1.129 |  2.000 |  1.000 |     2 |     33 |    1 | 3.67e-08
   0.09030 | 904 |  179.80 |  0.102 |  1.128 |  2.000 |  1.000 |     2 |     33 |    1 | 1.84e-07
   0.09040 | 905 |  179.72 |  0.102 |  1.128 |  2.000 |  1.000 |     2 |     33 |    1 | 1.73e-07
   0.09050 | 906 |  179.64 |  0.102 |  1.128 |  2.000 |  1.000 |     2 |     34 |    1 | 1.58e-07
   0.09060 | 907 |  179.55 |  0.102 |  1.127 |  2.000 |  1.000 |     2 |     35 |    1 | 1.17e-06
   0.09070 | 908 |  179.47 |  0.102 |  1.127 |  2.000 |  1.000 |     2 |     34 |    1 | 5.44e-07
   0.09080 | 909 |  179.38 |  0.102 |  1.127 |  2.000 |  1.000 |     2 |     35 |    1 | 1.51e-07
   0.09090 | 910 |  179.31 |  0.102 |  1.126 |  2.000 |  1.000 |     2 |     36 |    1 | 1.48e-05
   0.09100 | 911 |  179.24 |  0.102 |  1.126 |  2.000 |  1.000 |     2 |     35 |    1 | 6.83e-08

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.09110 | 912 |  179.18 |  0.102 |  1.126 |  2.000 |  1.000 |     2 |     35 |    1 | 1.14e-07
   0.09120 | 913 |  179.11 |  0.102 |  1.125 |  2.000 |  1.000 |     2 |     35 |    1 | 3.63e-07
   0.09130 | 914 |  179.05 |  0.102 |  1.125 |  2.000 |  1.000 |     2 |     35 |    1 | 1.03e-07
   0.09140 | 915 |  179.00 |  0.102 |  1.125 |  2.000 |  1.000 |     2 |     35 |    1 | 3.73e-07
   0.09150 | 916 |  178.96 |  0.102 |  1.125 |  2.000 |  1.000 |     2 |     35 |    1 | 1.67e-07
   0.09160 | 917 |  178.92 |  0.102 |  1.124 |  2.000 |  1.000 |     2 |     35 |    1 | 9.63e-08
   0.09170 | 918 |  178.88 |  0.102 |  1.124 |  2.000 |  1.000 |     2 |     34 |    1 | 7.35e-08
   0.09180 | 919 |  178.84 |  0.102 |  1.124 |  2.000 |  1.000 |     2 |     35 |    1 | 3.32e-07
   0.09190 | 920 |  178.81 |  0.102 |  1.123 |  2.000 |  1.000 |     2 |     34 |    1 | 1.39e-07
   0.09200 | 921 |  178.79 |  0.102 |  1.123 |  2.000 |  1.000 |     2 |     35 |    1 | 1.13e-07

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.09210 | 922 |  178.77 |  0.102 |  1.123 |  2.000 |  1.000 |     2 |     34 |    1 | 1.98e-07
   0.09220 | 923 |  178.76 |  0.102 |  1.123 |  2.000 |  1.000 |     2 |     35 |    1 | 8.78e-07
   0.09230 | 924 |  178.75 |  0.101 |  1.123 |  2.000 |  1.000 |     2 |     33 |    1 | 1.44e-07
   0.09240 | 925 |  178.74 |  0.101 |  1.122 |  2.000 |  1.000 |     2 |     34 |    1 | 9.72e-07
   0.09250 | 926 |  178.74 |  0.101 |  1.122 |  2.000 |  1.000 |     2 |     33 |    1 | 7.15e-08
   0.09260 | 927 |  178.74 |  0.101 |  1.122 |  2.000 |  1.000 |     2 |     33 |    1 | 4.12e-08
   0.09270 | 928 |  178.75 |  0.101 |  1.122 |  2.000 |  1.000 |     2 |     33 |    1 | 1.65e-08
   0.09280 | 929 |  178.76 |  0.101 |  1.122 |  2.000 |  1.000 |     2 |     33 |    1 | 4.95e-07
   0.09290 | 930 |  178.78 |  0.101 |  1.121 |  2.000 |  1.000 |     2 |     34 |    1 | 2.33e-06
   0.09300 | 931 |  178.80 |  0.101 |  1.121 |  2.000 |  1.000 |     2 |     33 |    1 | 7.92e-09

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.09310 | 932 |  178.83 |  0.101 |  1.121 |  2.000 |  1.000 |     2 |     33 |    1 | 1.34e-08
   0.09320 | 933 |  178.86 |  0.101 |  1.121 |  2.000 |  1.000 |     2 |     33 |    1 | 3.45e-08
   0.09330 | 934 |  178.88 |  0.100 |  1.121 |  2.000 |  1.000 |     2 |     32 |    1 | 2.20e-08
   0.09340 | 935 |  178.91 |  0.100 |  1.121 |  2.000 |  1.000 |     2 |     33 |    1 | 7.50e-09
   0.09350 | 936 |  178.94 |  0.100 |  1.121 |  2.000 |  1.000 |     2 |     33 |    1 | 5.32e-09
   0.09360 | 937 |  178.98 |  0.100 |  1.121 |  2.000 |  1.000 |     1 |     18 |    1 | 5.20e-05
   0.09370 | 938 |  179.01 |  0.100 |  1.121 |  2.000 |  1.000 |     2 |     33 |    1 | 4.55e-09
   0.09380 | 939 |  179.04 |  0.100 |  1.121 |  2.000 |  1.000 |     1 |     18 |    1 | 5.38e-05
   0.09390 | 940 |  179.07 |  0.100 |  1.121 |  2.000 |  1.000 |     2 |     33 |    1 | 4.45e-09
   0.09400 | 941 |  179.11 |  0.100 |  1.121 |  2.000 |  1.000 |     1 |     18 |    1 | 6.24e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.09410 | 942 |  179.14 |  0.100 |  1.120 |  2.000 |  1.000 |     1 |     18 |    1 | 3.02e-05
   0.09420 | 943 |  179.17 |  0.100 |  1.120 |  2.000 |  1.000 |     1 |     18 |    1 | 4.14e-05
   0.09430 | 944 |  179.20 |  0.100 |  1.120 |  2.000 |  1.000 |     2 |     34 |    1 | 4.73e-07
   0.09440 | 945 |  179.24 |  0.100 |  1.120 |  2.000 |  1.000 |     1 |     18 |    1 | 5.78e-05
   0.09450 | 946 |  179.27 |  0.100 |  1.120 |  2.000 |  1.000 |     1 |     18 |    1 | 2.76e-05
   0.09460 | 947 |  179.30 |  0.100 |  1.120 |  2.000 |  1.000 |     1 |     18 |    1 | 2.63e-05
   0.09470 | 948 |  179.33 |  0.100 |  1.120 |  2.000 |  1.000 |     1 |     18 |    1 | 1.37e-05
   0.09480 | 949 |  179.36 |  0.100 |  1.120 |  2.000 |  1.000 |     1 |     18 |    1 | 5.92e-06
   0.09490 | 950 |  179.39 |  0.100 |  1.120 |  2.000 |  1.000 |     1 |     18 |    1 | 1.47e-05
   0.09500 | 951 |  179.42 |  0.100 |  1.120 |  2.000 |  1.000 |     1 |     18 |    1 | 4.28e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.09510 | 952 |  179.45 |  0.100 |  1.120 |  2.000 |  1.000 |     1 |     18 |    1 | 2.23e-05
   0.09520 | 953 |  179.48 |  0.100 |  1.120 |  2.000 |  1.000 |     1 |     18 |    1 | 5.12e-06
   0.09530 | 954 |  179.51 |  0.100 |  1.120 |  2.000 |  1.000 |     1 |     18 |    1 | 2.35e-05
   0.09540 | 955 |  179.53 |  0.100 |  1.120 |  2.000 |  1.000 |     1 |     18 |    1 | 3.93e-05
   0.09550 | 956 |  179.56 |  0.100 |  1.120 |  2.000 |  1.000 |     1 |     18 |    1 | 3.86e-05
   0.09560 | 957 |  179.58 |  0.100 |  1.119 |  2.000 |  1.000 |     1 |     18 |    1 | 4.89e-05
   0.09570 | 958 |  179.60 |  0.100 |  1.119 |  2.000 |  1.000 |     1 |     18 |    1 | 3.49e-05
   0.09580 | 959 |  179.63 |  0.100 |  1.119 |  2.000 |  1.000 |     1 |     18 |    1 | 1.99e-05
   0.09590 | 960 |  179.65 |  0.100 |  1.119 |  2.000 |  1.000 |     1 |     18 |    1 | 3.63e-05
   0.09600 | 961 |  179.67 |  0.100 |  1.119 |  2.000 |  1.000 |     1 |     18 |    1 | 6.36e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.09610 | 962 |  179.69 |  0.100 |  1.119 |  2.000 |  1.000 |     1 |     18 |    1 | 3.51e-05
   0.09620 | 963 |  179.71 |  0.100 |  1.119 |  2.000 |  1.000 |     1 |     18 |    1 | 5.75e-05
   0.09630 | 964 |  179.73 |  0.100 |  1.119 |  2.000 |  1.000 |     1 |     18 |    1 | 3.05e-05
   0.09640 | 965 |  179.75 |  0.100 |  1.119 |  2.000 |  1.000 |     1 |     18 |    1 | 5.26e-05
   0.09650 | 966 |  179.76 |  0.100 |  1.119 |  2.000 |  1.000 |     1 |     18 |    1 | 5.86e-05
   0.09660 | 967 |  179.78 |  0.100 |  1.119 |  2.000 |  1.000 |     1 |     18 |    1 | 6.66e-05
   0.09670 | 968 |  179.80 |  0.100 |  1.119 |  2.000 |  1.000 |     1 |     18 |    1 | 5.78e-05
   0.09680 | 969 |  179.81 |  0.101 |  1.118 |  2.000 |  1.000 |     1 |     18 |    1 | 6.01e-05
   0.09690 | 970 |  179.83 |  0.101 |  1.118 |  2.000 |  1.000 |     1 |     18 |    1 | 6.99e-05
   0.09700 | 971 |  179.84 |  0.101 |  1.118 |  2.000 |  1.000 |     1 |     18 |    1 | 6.85e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.09710 | 972 |  179.85 |  0.101 |  1.118 |  2.000 |  1.000 |     1 |     18 |    1 | 7.70e-05
   0.09720 | 973 |  179.87 |  0.101 |  1.118 |  2.000 |  1.000 |     1 |     18 |    1 | 7.44e-05
   0.09730 | 974 |  179.88 |  0.101 |  1.118 |  2.000 |  1.000 |     1 |     18 |    1 | 8.35e-05
   0.09740 | 975 |  179.89 |  0.101 |  1.118 |  2.000 |  1.000 |     1 |     18 |    1 | 8.59e-05
   0.09750 | 976 |  179.90 |  0.101 |  1.118 |  2.000 |  1.000 |     1 |     18 |    1 | 8.51e-05
   0.09760 | 977 |  179.91 |  0.101 |  1.118 |  2.000 |  1.000 |     1 |     18 |    1 | 8.19e-05
   0.09770 | 978 |  179.92 |  0.101 |  1.118 |  2.000 |  1.000 |     1 |     18 |    1 | 6.12e-05
   0.09780 | 979 |  179.93 |  0.101 |  1.118 |  2.000 |  1.000 |     1 |     18 |    1 | 8.82e-05
   0.09790 | 980 |  179.94 |  0.101 |  1.117 |  2.000 |  1.000 |     1 |     18 |    1 | 9.07e-05
   0.09800 | 981 |  179.95 |  0.101 |  1.117 |  2.000 |  1.000 |     1 |     18 |    1 | 7.58e-05

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.09810 | 982 |  179.96 |  0.101 |  1.117 |  2.000 |  1.000 |     1 |     18 |    1 | 9.57e-05
   0.09820 | 983 |  179.97 |  0.101 |  1.117 |  2.000 |  1.000 |     1 |     18 |    1 | 9.51e-05
   0.09830 | 984 |  179.97 |  0.101 |  1.117 |  2.000 |  1.000 |     1 |     18 |    1 | 9.89e-05
   0.09840 | 985 |  179.98 |  0.101 |  1.117 |  2.000 |  1.000 |     1 |     18 |    1 | 9.47e-05
   0.09850 | 986 |  179.99 |  0.101 |  1.117 |  2.000 |  1.000 |     2 |     33 |    1 | 5.12e-08
   0.09860 | 987 |  179.99 |  0.101 |  1.117 |  2.000 |  1.000 |     1 |     18 |    1 | 8.26e-05
   0.09870 | 988 |  180.00 |  0.101 |  1.117 |  2.000 |  1.000 |     2 |     34 |    1 | 1.01e-07
   0.09880 | 989 |  180.00 |  0.101 |  1.117 |  2.000 |  1.000 |     1 |     19 |    1 | 8.16e-05
   0.09890 | 990 |  180.01 |  0.101 |  1.116 |  2.000 |  1.000 |     2 |     34 |    1 | 1.36e-08
   0.09900 | 991 |  180.01 |  0.101 |  1.116 |  2.000 |  1.000 |     2 |     34 |    1 | 4.64e-08

         t |  it | Ip [kA] |  zA[m] |  rA[m] |     bp |     qA | niter | nfeval | conv |      res
  ===============================================================================================
   0.09910 | 992 |  180.02 |  0.101 |  1.116 |  2.000 |  1.000 |     2 |     34 |    1 | 1.16e-08
   0.09920 | 993 |  180.02 |  0.101 |  1.116 |  2.000 |  1.000 |     2 |     34 |    1 | 4.79e-08
   0.09930 | 994 |  180.03 |  0.101 |  1.116 |  2.000 |  1.000 |     1 |     19 |    1 | 9.54e-05
   0.09940 | 995 |  180.03 |  0.101 |  1.116 |  2.000 |  1.000 |     2 |     34 |    1 | 2.41e-08
   0.09950 | 996 |  180.03 |  0.101 |  1.116 |  2.000 |  1.000 |     1 |     19 |    1 | 9.87e-05
   0.09960 | 997 |  180.04 |  0.101 |  1.116 |  2.000 |  1.000 |     2 |     34 |    1 | 2.34e-08
   0.09970 | 998 |  180.04 |  0.101 |  1.116 |  2.000 |  1.000 |     1 |     19 |    1 | 9.34e-05
   0.09980 | 999 |  180.04 |  0.101 |  1.116 |  2.000 |  1.000 |     1 |     19 |    1 | 9.96e-05
   0.09990 | 1000 |  180.05 |  0.101 |  1.115 |  2.000 |  1.000 |     2 |     34 |    1 | 9.02e-09
   0.10000 | 1001 |  180.05 |  0.101 |  1.115 |  2.000 |  1.000 |     1 |     19 |    1 | 8.77e-05
Time elapsed (fget.m): 89.704[s]

Run RZP simulation

Lz.P.ctrlfct = @meqctrl;
Lz.P.ctrlpar = ctrlpar;
Lz.P.debug=0;
LYz = rzpt(Lz,LXz);
[..................................................] 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