Apollo lunar descent guidance demo

See also: ApolloLunarDescentGuidance

%--------------------------------------------------------------------------
%   Copyright (c) 2015 Princeton Satellite Systems, Inc.
%   All rights reserved.
%--------------------------------------------------------------------------

tEnd  = 300;
dT    = 1;
n     = ceil(tEnd/dT);
h     = 15;
r     = h+1738;
v     = sqrt(Constant('mu moon')/r)*1000;
x     = [1000*h;0;492000;0;0;-v];
xP    = zeros(9,n);
mode  = cell(n,1);
t     = 0;
dG    = ApolloLunarDescentGuidance;
fC    = 1/(r*1000);
gM    = fC*v^2;

for k = 1:n
  [aCG,dG]	= ApolloLunarDescentGuidance( x, t, dG );
  xP(:,k)   = [x;aCG];
  mode{k}   = dG.mode;
  x         = RK4(@RHS,x,dT,0,aCG,fC,gM);
  t         = t + dT;
end

[t,tL] = TimeLabl((0:(k-1))*dT);

yL = {'R_{G_x} (m)'     'R_{G_y} (m)',      'R_{G_z} (m)',...
      'V_{G_x} (m/s)'   'V_{G_y} (m/s)',    'V_{G_z} (m/s)',...
      'A_{G_x} (m/s^2)' 'A_{G_y} (m/s^2)',  'A_{G_z} (m/s^2)'};

h = xP(1,:);
z = xP(3,:);
NewFig('Trajectory')
plot(z,h);
grid on
text(z(1),1.01*h(1),mode{1});
text(z(1),0.99*h(1),sprintf('t = %4.1f min',t(1)));

for k = 2:n
  if( ~strcmp(mode{k},mode{k-1}) )
    text(z(k),1.01*h(k),mode{k});
    text(z(k),0.99*h(k),sprintf('t = %4.1f min',t(k)));
  end
end
text(z(n),h(n),sprintf('t = %4.1f min',t(n)));

Plot2D(t,xP(1:3,:),tL,yL(1:3),'R_G');
Plot2D(t,xP(4:6,:),tL,yL(4:6),'V_G');
Plot2D(t,xP(7:9,:),tL,yL(7:9),'A_G');

%--------------------------------------------------------------------------
% Right hand side
%--------------------------------------------------------------------------
function xDot = RHS(x,~,aCG,fC,gM)

a     = [fC*x(6)^2 - gM;0;0];
xDot	= [x(4:6);aCG+a];

end

%--------------------------------------
% $Date$
% $Id: e81e416d77672edbf7906fe2eee310623f3dc4e5 $
t =      0.0 mode: P63 rG = [15000.00     0.00 492000.00] aCG = [   -2.31     0.00    -2.26]
t =      1.0 mode: P63 rG = [14998.85     0.00 490326.51] aCG = [   -2.28     0.00    -2.20]
t =      2.0 mode: P63 rG = [14995.40     0.00 488650.78] aCG = [   -2.26     0.00    -2.15]
t =      3.0 mode: P63 rG = [14989.69     0.00 486972.88] aCG = [   -2.24     0.00    -2.10]
t =      4.0 mode: P63 rG = [14981.75     0.00 485292.85] aCG = [   -2.21     0.00    -2.04]
t =      5.0 mode: P63 rG = [14971.60     0.00 483610.75] aCG = [   -2.19     0.00    -1.99]
t =      6.0 mode: P63 rG = [14959.26     0.00 481926.64] aCG = [   -2.17     0.00    -1.94]
t =      7.0 mode: P63 rG = [14944.78     0.00 480240.56] aCG = [   -2.14     0.00    -1.88]
t =      8.0 mode: P63 rG = [14928.16     0.00 478552.58] aCG = [   -2.12     0.00    -1.83]
t =      9.0 mode: P63 rG = [14909.45     0.00 476862.73] aCG = [   -2.10     0.00    -1.78]
t =     10.0 mode: P63 rG = [14888.67     0.00 475171.09] aCG = [   -2.07     0.00    -1.72]
t =     11.0 mode: P63 rG = [14865.83     0.00 473477.69] aCG = [   -2.05     0.00    -1.67]
t =     12.0 mode: P63 rG = [14840.98     0.00 471782.60] aCG = [   -2.03     0.00    -1.62]
t =     13.0 mode: P63 rG = [14814.13     0.00 470085.86] aCG = [   -2.01     0.00    -1.57]
t =     14.0 mode: P63 rG = [14785.31     0.00 468387.53] aCG = [   -1.98     0.00    -1.52]
t =     15.0 mode: P63 rG = [14754.55     0.00 466687.65] aCG = [   -1.96     0.00    -1.46]
t =     16.0 mode: P63 rG = [14721.87     0.00 464986.29] aCG = [   -1.94     0.00    -1.41]
t =     17.0 mode: P63 rG = [14687.29     0.00 463283.49] aCG = [   -1.92     0.00    -1.36]
t =     18.0 mode: P63 rG = [14650.85     0.00 461579.30] aCG = [   -1.90     0.00    -1.31]
t =     19.0 mode: P63 rG = [14612.56     0.00 459873.78] aCG = [   -1.87     0.00    -1.26]
t =     20.0 mode: P63 rG = [14572.45     0.00 458166.97] aCG = [   -1.85     0.00    -1.21]
t =     21.0 mode: P63 rG = [14530.55     0.00 456458.93] aCG = [   -1.83     0.00    -1.16]
t =     22.0 mode: P63 rG = [14486.87     0.00 454749.70] aCG = [   -1.81     0.00    -1.11]
t =     23.0 mode: P63 rG = [14441.45     0.00 453039.34] aCG = [   -1.79     0.00    -1.06]
t =     24.0 mode: P63 rG = [14394.30     0.00 451327.90] aCG = [   -1.77     0.00    -1.01]
t =     25.0 mode: P63 rG = [14345.44     0.00 449615.43] aCG = [   -1.75     0.00    -0.96]
t =     26.0 mode: P63 rG = [14294.91     0.00 447901.98] aCG = [   -1.73     0.00    -0.91]
t =     27.0 mode: P63 rG = [14242.73     0.00 446187.59] aCG = [   -1.70     0.00    -0.86]
t =     28.0 mode: P63 rG = [14188.90     0.00 444472.32] aCG = [   -1.68     0.00    -0.81]
t =     29.0 mode: P63 rG = [14133.47     0.00 442756.21] aCG = [   -1.66     0.00    -0.76]
t =     30.0 mode: P63 rG = [14076.45     0.00 441039.32] aCG = [   -1.64     0.00    -0.71]
t =     31.0 mode: P63 rG = [14017.87     0.00 439321.69] aCG = [   -1.62     0.00    -0.66]
t =     32.0 mode: P63 rG = [13957.74     0.00 437603.38] aCG = [   -1.60     0.00    -0.61]
t =     33.0 mode: P63 rG = [13896.08     0.00 435884.43] aCG = [   -1.58     0.00    -0.57]
t =     34.0 mode: P63 rG = [13832.93     0.00 434164.88] aCG = [   -1.56     0.00    -0.52]
t =     35.0 mode: P63 rG = [13768.29     0.00 432444.80] aCG = [   -1.54     0.00    -0.47]
t =     36.0 mode: P63 rG = [13702.20     0.00 430724.22] aCG = [   -1.52     0.00    -0.42]
t =     37.0 mode: P63 rG = [13634.66     0.00 429003.19] aCG = [   -1.50     0.00    -0.38]
t =     38.0 mode: P63 rG = [13565.71     0.00 427281.76] aCG = [   -1.48     0.00    -0.33]
t =     39.0 mode: P63 rG = [13495.37     0.00 425559.98] aCG = [   -1.46     0.00    -0.28]
t =     40.0 mode: P63 rG = [13423.64     0.00 423837.90] aCG = [   -1.44     0.00    -0.23]
t =     41.0 mode: P63 rG = [13350.56     0.00 422115.55] aCG = [   -1.42     0.00    -0.19]
t =     42.0 mode: P63 rG = [13276.14     0.00 420393.00] aCG = [   -1.40     0.00    -0.14]
t =     43.0 mode: P63 rG = [13200.41     0.00 418670.28] aCG = [   -1.39     0.00    -0.10]
t =     44.0 mode: P63 rG = [13123.38     0.00 416947.44] aCG = [   -1.37     0.00    -0.05]
t =     45.0 mode: P63 rG = [13045.07     0.00 415224.54] aCG = [   -1.35     0.00    -0.00]
t =     46.0 mode: P63 rG = [12965.50     0.00 413501.60] aCG = [   -1.33     0.00     0.04]
t =     47.0 mode: P63 rG = [12884.69     0.00 411778.68] aCG = [   -1.31     0.00     0.09]
t =     48.0 mode: P63 rG = [12802.67     0.00 410055.83] aCG = [   -1.29     0.00     0.13]
t =     49.0 mode: P63 rG = [12719.44     0.00 408333.09] aCG = [   -1.27     0.00     0.18]
t =     50.0 mode: P63 rG = [12635.03     0.00 406610.51] aCG = [   -1.25     0.00     0.22]
t =     51.0 mode: P63 rG = [12549.45     0.00 404888.13] aCG = [   -1.23     0.00     0.27]
t =     52.0 mode: P63 rG = [12462.73     0.00 403165.99] aCG = [   -1.22     0.00     0.31]
t =     53.0 mode: P63 rG = [12374.87     0.00 401444.14] aCG = [   -1.20     0.00     0.36]
t =     54.0 mode: P63 rG = [12285.91     0.00 399722.63] aCG = [   -1.18     0.00     0.40]
t =     55.0 mode: P63 rG = [12195.86     0.00 398001.50] aCG = [   -1.16     0.00     0.44]
t =     56.0 mode: P63 rG = [12104.73     0.00 396280.79] aCG = [   -1.14     0.00     0.49]
t =     57.0 mode: P63 rG = [12012.54     0.00 394560.54] aCG = [   -1.12     0.00     0.53]
t =     58.0 mode: P63 rG = [11919.31     0.00 392840.81] aCG = [   -1.11     0.00     0.58]
t =     59.0 mode: P63 rG = [11825.06     0.00 391121.63] aCG = [   -1.09     0.00     0.62]
t =     60.0 mode: P63 rG = [11729.80     0.00 389403.05] aCG = [   -1.07     0.00     0.66]
t =     61.0 mode: P63 rG = [11633.55     0.00 387685.11] aCG = [   -1.05     0.00     0.70]
t =     62.0 mode: P63 rG = [11536.33     0.00 385967.85] aCG = [   -1.03     0.00     0.75]
t =     63.0 mode: P63 rG = [11438.15     0.00 384251.32] aCG = [   -1.02     0.00     0.79]
t =     64.0 mode: P63 rG = [11339.04     0.00 382535.55] aCG = [   -1.00     0.00     0.83]
t =     65.0 mode: P63 rG = [11238.99     0.00 380820.60] aCG = [   -0.98     0.00     0.87]
t =     66.0 mode: P63 rG = [11138.04     0.00 379106.50] aCG = [   -0.96     0.00     0.92]
t =     67.0 mode: P63 rG = [11036.20     0.00 377393.30] aCG = [   -0.95     0.00     0.96]
t =     68.0 mode: P63 rG = [10933.49     0.00 375681.03] aCG = [   -0.93     0.00     1.00]
t =     69.0 mode: P63 rG = [10829.91     0.00 373969.74] aCG = [   -0.91     0.00     1.04]
t =     70.0 mode: P63 rG = [10725.48     0.00 372259.46] aCG = [   -0.89     0.00     1.08]
t =     71.0 mode: P63 rG = [10620.23     0.00 370550.25] aCG = [   -0.88     0.00     1.12]
t =     72.0 mode: P63 rG = [10514.16     0.00 368842.14] aCG = [   -0.86     0.00     1.16]
t =     73.0 mode: P63 rG = [10407.30     0.00 367135.18] aCG = [   -0.84     0.00     1.20]
t =     74.0 mode: P63 rG = [10299.65     0.00 365429.39] aCG = [   -0.82     0.00     1.24]
t =     75.0 mode: P63 rG = [10191.23     0.00 363724.83] aCG = [   -0.81     0.00     1.28]
t =     76.0 mode: P63 rG = [10082.05     0.00 362021.53] aCG = [   -0.79     0.00     1.32]
t =     77.0 mode: P63 rG = [ 9972.14     0.00 360319.54] aCG = [   -0.77     0.00     1.36]
t =     78.0 mode: P63 rG = [ 9861.49     0.00 358618.89] aCG = [   -0.76     0.00     1.40]
t =     79.0 mode: P63 rG = [ 9750.14     0.00 356919.62] aCG = [   -0.74     0.00     1.44]
t =     80.0 mode: P63 rG = [ 9638.09     0.00 355221.77] aCG = [   -0.72     0.00     1.48]
t =     81.0 mode: P63 rG = [ 9525.36     0.00 353525.39] aCG = [   -0.71     0.00     1.52]
t =     82.0 mode: P63 rG = [ 9411.95     0.00 351830.51] aCG = [   -0.69     0.00     1.56]
t =     83.0 mode: P63 rG = [ 9297.89     0.00 350137.16] aCG = [   -0.67     0.00     1.60]
t =     84.0 mode: P63 rG = [ 9183.19     0.00 348445.39] aCG = [   -0.66     0.00     1.64]
t =     85.0 mode: P63 rG = [ 9067.87     0.00 346755.24] aCG = [   -0.64     0.00     1.67]
t =     86.0 mode: P63 rG = [ 8951.92     0.00 345066.75] aCG = [   -0.62     0.00     1.71]
t =     87.0 mode: P63 rG = [ 8835.38     0.00 343379.94] aCG = [   -0.61     0.00     1.75]
t =     88.0 mode: P63 rG = [ 8718.24     0.00 341694.87] aCG = [   -0.59     0.00     1.79]
t =     89.0 mode: P63 rG = [ 8600.54     0.00 340011.57] aCG = [   -0.57     0.00     1.82]
t =     90.0 mode: P63 rG = [ 8482.27     0.00 338330.07] aCG = [   -0.56     0.00     1.86]
t =     91.0 mode: P63 rG = [ 8363.45     0.00 336650.41] aCG = [   -0.54     0.00     1.90]
t =     92.0 mode: P63 rG = [ 8244.09     0.00 334972.64] aCG = [   -0.52     0.00     1.94]
t =     93.0 mode: P63 rG = [ 8124.21     0.00 333296.78] aCG = [   -0.51     0.00     1.97]
t =     94.0 mode: P63 rG = [ 8003.82     0.00 331622.87] aCG = [   -0.49     0.00     2.01]
t =     95.0 mode: P63 rG = [ 7882.93     0.00 329950.96] aCG = [   -0.48     0.00     2.04]
t =     96.0 mode: P63 rG = [ 7761.56     0.00 328281.07] aCG = [   -0.46     0.00     2.08]
t =     97.0 mode: P63 rG = [ 7639.71     0.00 326613.24] aCG = [   -0.44     0.00     2.12]
t =     98.0 mode: P63 rG = [ 7517.40     0.00 324947.52] aCG = [   -0.43     0.00     2.15]
t =     99.0 mode: P63 rG = [ 7394.64     0.00 323283.92] aCG = [   -0.41     0.00     2.19]
t =    100.0 mode: P63 rG = [ 7271.45     0.00 321622.50] aCG = [   -0.39     0.00     2.22]
t =    101.0 mode: P63 rG = [ 7147.83     0.00 319963.28] aCG = [   -0.38     0.00     2.26]
t =    102.0 mode: P63 rG = [ 7023.79     0.00 318306.31] aCG = [   -0.36     0.00     2.29]
t =    103.0 mode: P63 rG = [ 6899.35     0.00 316651.60] aCG = [   -0.35     0.00     2.33]
t =    104.0 mode: P63 rG = [ 6774.53     0.00 314999.21] aCG = [   -0.33     0.00     2.36]
t =    105.0 mode: P63 rG = [ 6649.32     0.00 313349.17] aCG = [   -0.31     0.00     2.40]
t =    106.0 mode: P63 rG = [ 6523.75     0.00 311701.50] aCG = [   -0.30     0.00     2.43]
t =    107.0 mode: P63 rG = [ 6397.82     0.00 310056.24] aCG = [   -0.28     0.00     2.46]
t =    108.0 mode: P63 rG = [ 6271.54     0.00 308413.44] aCG = [   -0.27     0.00     2.50]
t =    109.0 mode: P63 rG = [ 6144.94     0.00 306773.11] aCG = [   -0.25     0.00     2.53]
t =    110.0 mode: P63 rG = [ 6018.01     0.00 305135.30] aCG = [   -0.24     0.00     2.56]
t =    111.0 mode: P63 rG = [ 5890.77     0.00 303500.04] aCG = [   -0.22     0.00     2.60]
t =    112.0 mode: P63 rG = [ 5763.23     0.00 301867.35] aCG = [   -0.20     0.00     2.63]
t =    113.0 mode: P63 rG = [ 5635.40     0.00 300237.29] aCG = [   -0.19     0.00     2.66]
t =    114.0 mode: P63 rG = [ 5507.29     0.00 298609.87] aCG = [   -0.17     0.00     2.70]
t =    115.0 mode: P63 rG = [ 5378.92     0.00 296985.13] aCG = [   -0.16     0.00     2.73]
t =    116.0 mode: P63 rG = [ 5250.29     0.00 295363.10] aCG = [   -0.14     0.00     2.76]
t =    117.0 mode: P63 rG = [ 5121.41     0.00 293743.82] aCG = [   -0.13     0.00     2.79]
t =    118.0 mode: P63 rG = [ 4992.29     0.00 292127.31] aCG = [   -0.11     0.00     2.82]
t =    119.0 mode: P63 rG = [ 4862.95     0.00 290513.61] aCG = [   -0.09     0.00     2.86]
t =    120.0 mode: P63 rG = [ 4733.40     0.00 288902.76] aCG = [   -0.08     0.00     2.89]
t =    121.0 mode: P63 rG = [ 4603.64     0.00 287294.77] aCG = [   -0.06     0.00     2.92]
t =    122.0 mode: P63 rG = [ 4473.69     0.00 285689.69] aCG = [   -0.05     0.00     2.95]
t =    123.0 mode: P63 rG = [ 4343.55     0.00 284087.54] aCG = [   -0.03     0.00     2.98]
t =    124.0 mode: P63 rG = [ 4213.24     0.00 282488.36] aCG = [   -0.02     0.00     3.01]
t =    125.0 mode: P63 rG = [ 4082.77     0.00 280892.17] aCG = [   -0.00     0.00     3.04]
t =    126.0 mode: P63 rG = [ 3952.14     0.00 279299.01] aCG = [    0.01     0.00     3.07]
t =    127.0 mode: P63 rG = [ 3821.36     0.00 277708.91] aCG = [    0.03     0.00     3.10]
t =    128.0 mode: P63 rG = [ 3690.46     0.00 276121.89] aCG = [    0.04     0.00     3.13]
t =    129.0 mode: P63 rG = [ 3559.43     0.00 274537.99] aCG = [    0.06     0.00     3.16]
t =    130.0 mode: P63 rG = [ 3428.28     0.00 272957.24] aCG = [    0.08     0.00     3.19]
t =    131.0 mode: P63 rG = [ 3297.03     0.00 271379.66] aCG = [    0.09     0.00     3.22]
t =    132.0 mode: P63 rG = [ 3165.68     0.00 269805.29] aCG = [    0.11     0.00     3.25]
t =    133.0 mode: P63 rG = [ 3034.25     0.00 268234.15] aCG = [    0.12     0.00     3.28]
t =    134.0 mode: P63 rG = [ 2902.75     0.00 266666.28] aCG = [    0.14     0.00     3.31]
t =    135.0 mode: P63 rG = [ 2771.17     0.00 265101.70] aCG = [    0.15     0.00     3.34]
t =    136.0 mode: P63 rG = [ 2639.54     0.00 263540.44] aCG = [    0.17     0.00     3.36]
t =    137.0 mode: P63 rG = [ 2507.86     0.00 261982.53] aCG = [    0.18     0.00     3.39]
t =    138.0 mode: P63 rG = [ 2376.15     0.00 260427.99] aCG = [    0.20     0.00     3.42]
t =    139.0 mode: P63 rG = [ 2244.40     0.00 258876.87] aCG = [    0.21     0.00     3.45]
t =    140.0 mode: P63 rG = [ 2112.64     0.00 257329.17] aCG = [    0.23     0.00     3.47]
t =    141.0 mode: P63 rG = [ 1980.86     0.00 255784.94] aCG = [    0.24     0.00     3.50]
t =    142.0 mode: P63 rG = [ 1849.09     0.00 254244.20] aCG = [    0.26     0.00     3.53]
t =    143.0 mode: P63 rG = [ 1717.32     0.00 252706.97] aCG = [    0.28     0.00     3.56]
t =    144.0 mode: P63 rG = [ 1585.57     0.00 251173.28] aCG = [    0.29     0.00     3.58]
t =    145.0 mode: P63 rG = [ 1453.84     0.00 249643.17] aCG = [    0.31     0.00     3.61]
t =    146.0 mode: P63 rG = [ 1322.15     0.00 248116.65] aCG = [    0.32     0.00     3.64]
t =    147.0 mode: P63 rG = [ 1190.50     0.00 246593.76] aCG = [    0.34     0.00     3.66]
t =    148.0 mode: P63 rG = [ 1058.91     0.00 245074.51] aCG = [    0.35     0.00     3.69]
t =    149.0 mode: P63 rG = [  927.38     0.00 243558.94] aCG = [    0.37     0.00     3.71]
t =    150.0 mode: P63 rG = [  795.92     0.00 242047.07] aCG = [    0.38     0.00     3.74]
t =    151.0 mode: P63 rG = [  664.54     0.00 240538.93] aCG = [    0.40     0.00     3.77]
t =    152.0 mode: P63 rG = [  533.24     0.00 239034.54] aCG = [    0.41     0.00     3.79]
t =    153.0 mode: P63 rG = [  402.04     0.00 237533.92] aCG = [    0.43     0.00     3.82]
t =    154.0 mode: P63 rG = [  270.95     0.00 236037.11] aCG = [    0.44     0.00     3.84]
t =    155.0 mode: P63 rG = [  139.98     0.00 234544.13] aCG = [    0.46     0.00     3.86]
t =    156.0 mode: P63 rG = [    9.12     0.00 233055.00] aCG = [    0.47     0.00     3.89]
t =    157.0 mode: P63 rG = [ -121.60     0.00 231569.75] aCG = [    0.49     0.00     3.91]
t =    158.0 mode: P63 rG = [ -252.18     0.00 230088.40] aCG = [    0.50     0.00     3.94]
t =    159.0 mode: P63 rG = [ -382.62     0.00 228610.97] aCG = [    0.52     0.00     3.96]
t =    160.0 mode: P63 rG = [ -512.89     0.00 227137.50] aCG = [    0.53     0.00     3.99]
t =    161.0 mode: P63 rG = [ -643.01     0.00 225668.00] aCG = [    0.55     0.00     4.01]
t =    162.0 mode: P63 rG = [ -772.94     0.00 224202.49] aCG = [    0.56     0.00     4.03]
t =    163.0 mode: P63 rG = [ -902.70     0.00 222741.01] aCG = [    0.58     0.00     4.06]
t =    164.0 mode: P63 rG = [-1032.26     0.00 221283.57] aCG = [    0.59     0.00     4.08]
t =    165.0 mode: P63 rG = [-1161.62     0.00 219830.19] aCG = [    0.61     0.00     4.10]
t =    166.0 mode: P63 rG = [-1290.78     0.00 218380.91] aCG = [    0.62     0.00     4.12]
t =    167.0 mode: P63 rG = [-1419.72     0.00 216935.73] aCG = [    0.64     0.00     4.15]
t =    168.0 mode: P63 rG = [-1548.44     0.00 215494.69] aCG = [    0.65     0.00     4.17]
t =    169.0 mode: P63 rG = [-1676.92     0.00 214057.81] aCG = [    0.67     0.00     4.19]
t =    170.0 mode: P63 rG = [-1805.16     0.00 212625.11] aCG = [    0.68     0.00     4.21]
t =    171.0 mode: P63 rG = [-1933.16     0.00 211196.61] aCG = [    0.70     0.00     4.23]
t =    172.0 mode: P63 rG = [-2060.89     0.00 209772.32] aCG = [    0.71     0.00     4.25]
t =    173.0 mode: P63 rG = [-2188.36     0.00 208352.29] aCG = [    0.73     0.00     4.28]
t =    174.0 mode: P63 rG = [-2315.56     0.00 206936.52] aCG = [    0.74     0.00     4.30]
t =    175.0 mode: P63 rG = [-2442.47     0.00 205525.03] aCG = [    0.76     0.00     4.32]
t =    176.0 mode: P63 rG = [-2569.10     0.00 204117.85] aCG = [    0.77     0.00     4.34]
t =    177.0 mode: P63 rG = [-2695.43     0.00 202715.00] aCG = [    0.79     0.00     4.36]
t =    178.0 mode: P63 rG = [-2821.45     0.00 201316.50] aCG = [    0.81     0.00     4.38]
t =    179.0 mode: P63 rG = [-2947.15     0.00 199922.36] aCG = [    0.82     0.00     4.40]
t =    180.0 mode: P63 rG = [-3072.54     0.00 198532.61] aCG = [    0.84     0.00     4.42]
t =    181.0 mode: P63 rG = [-3197.59     0.00 197147.28] aCG = [    0.85     0.00     4.44]
t =    182.0 mode: P63 rG = [-3322.31     0.00 195766.37] aCG = [    0.87     0.00     4.46]
t =    183.0 mode: P63 rG = [-3446.68     0.00 194389.90] aCG = [    0.88     0.00     4.48]
t =    184.0 mode: P63 rG = [-3570.69     0.00 193017.91] aCG = [    0.90     0.00     4.50]
t =    185.0 mode: P63 rG = [-3694.34     0.00 191650.40] aCG = [    0.91     0.00     4.51]
t =    186.0 mode: P63 rG = [-3817.62     0.00 190287.39] aCG = [    0.93     0.00     4.53]
t =    187.0 mode: P63 rG = [-3940.53     0.00 188928.91] aCG = [    0.94     0.00     4.55]
t =    188.0 mode: P63 rG = [-4063.04     0.00 187574.98] aCG = [    0.95     0.00     4.57]
t =    189.0 mode: P63 rG = [-4185.16     0.00 186225.60] aCG = [    0.97     0.00     4.59]
t =    190.0 mode: P63 rG = [-4306.88     0.00 184880.80] aCG = [    0.98     0.00     4.61]
t =    191.0 mode: P63 rG = [-4428.19     0.00 183540.60] aCG = [    1.00     0.00     4.62]
t =    192.0 mode: P63 rG = [-4549.09     0.00 182205.02] aCG = [    1.01     0.00     4.64]
t =    193.0 mode: P63 rG = [-4669.55     0.00 180874.06] aCG = [    1.03     0.00     4.66]
t =    194.0 mode: P63 rG = [-4789.59     0.00 179547.76] aCG = [    1.04     0.00     4.68]
t =    195.0 mode: P63 rG = [-4909.18     0.00 178226.12] aCG = [    1.06     0.00     4.69]
t =    196.0 mode: P63 rG = [-5028.32     0.00 176909.17] aCG = [    1.07     0.00     4.71]
t =    197.0 mode: P63 rG = [-5147.00     0.00 175596.92] aCG = [    1.09     0.00     4.73]
t =    198.0 mode: P63 rG = [-5265.22     0.00 174289.39] aCG = [    1.10     0.00     4.74]
t =    199.0 mode: P63 rG = [-5382.97     0.00 172986.59] aCG = [    1.12     0.00     4.76]
t =    200.0 mode: P63 rG = [-5500.24     0.00 171688.54] aCG = [    1.13     0.00     4.77]
t =    201.0 mode: P63 rG = [-5617.02     0.00 170395.25] aCG = [    1.15     0.00     4.79]
t =    202.0 mode: P63 rG = [-5733.30     0.00 169106.75] aCG = [    1.16     0.00     4.81]
t =    203.0 mode: P63 rG = [-5849.08     0.00 167823.05] aCG = [    1.18     0.00     4.82]
t =    204.0 mode: P63 rG = [-5964.34     0.00 166544.16] aCG = [    1.19     0.00     4.84]
t =    205.0 mode: P63 rG = [-6079.09     0.00 165270.10] aCG = [    1.21     0.00     4.85]
t =    206.0 mode: P63 rG = [-6193.31     0.00 164000.88] aCG = [    1.22     0.00     4.87]
t =    207.0 mode: P63 rG = [-6307.00     0.00 162736.52] aCG = [    1.24     0.00     4.88]
t =    208.0 mode: P63 rG = [-6420.14     0.00 161477.04] aCG = [    1.25     0.00     4.90]
t =    209.0 mode: P63 rG = [-6532.73     0.00 160222.44] aCG = [    1.27     0.00     4.91]
t =    210.0 mode: P63 rG = [-6644.77     0.00 158972.75] aCG = [    1.28     0.00     4.92]
t =    211.0 mode: P63 rG = [-6756.24     0.00 157727.97] aCG = [    1.30     0.00     4.94]
t =    212.0 mode: P63 rG = [-6867.13     0.00 156488.12] aCG = [    1.31     0.00     4.95]
t =    213.0 mode: P63 rG = [-6977.45     0.00 155253.22] aCG = [    1.32     0.00     4.96]
t =    214.0 mode: P63 rG = [-7087.18     0.00 154023.27] aCG = [    1.34     0.00     4.98]
t =    215.0 mode: P63 rG = [-7196.31     0.00 152798.30] aCG = [    1.35     0.00     4.99]
t =    216.0 mode: P63 rG = [-7304.84     0.00 151578.31] aCG = [    1.37     0.00     5.00]
t =    217.0 mode: P63 rG = [-7412.75     0.00 150363.31] aCG = [    1.38     0.00     5.02]
t =    218.0 mode: P63 rG = [-7520.05     0.00 149153.33] aCG = [    1.40     0.00     5.03]
t =    219.0 mode: P63 rG = [-7626.72     0.00 147948.37] aCG = [    1.41     0.00     5.04]
t =    220.0 mode: P63 rG = [-7732.76     0.00 146748.44] aCG = [    1.43     0.00     5.05]
t =    221.0 mode: P63 rG = [-7838.16     0.00 145553.57] aCG = [    1.44     0.00     5.07]
t =    222.0 mode: P63 rG = [-7942.91     0.00 144363.75] aCG = [    1.45     0.00     5.08]
t =    223.0 mode: P63 rG = [-8047.00     0.00 143179.00] aCG = [    1.47     0.00     5.09]
t =    224.0 mode: P63 rG = [-8150.42     0.00 141999.33] aCG = [    1.48     0.00     5.10]
t =    225.0 mode: P63 rG = [-8253.18     0.00 140824.76] aCG = [    1.50     0.00     5.11]
t =    226.0 mode: P63 rG = [-8355.26     0.00 139655.30] aCG = [    1.51     0.00     5.12]
t =    227.0 mode: P63 rG = [-8456.65     0.00 138490.95] aCG = [    1.53     0.00     5.13]
t =    228.0 mode: P63 rG = [-8557.34     0.00 137331.73] aCG = [    1.54     0.00     5.14]
t =    229.0 mode: P63 rG = [-8657.34     0.00 136177.64] aCG = [    1.55     0.00     5.15]
t =    230.0 mode: P63 rG = [-8756.62     0.00 135028.71] aCG = [    1.57     0.00     5.16]
t =    231.0 mode: P63 rG = [-8855.19     0.00 133884.93] aCG = [    1.58     0.00     5.17]
t =    232.0 mode: P63 rG = [-8953.04     0.00 132746.32] aCG = [    1.60     0.00     5.18]
t =    233.0 mode: P63 rG = [-9050.16     0.00 131612.90] aCG = [    1.61     0.00     5.19]
t =    234.0 mode: P63 rG = [-9146.53     0.00 130484.65] aCG = [    1.63     0.00     5.20]
t =    235.0 mode: P63 rG = [-9242.17     0.00 129361.61] aCG = [    1.64     0.00     5.21]
t =    236.0 mode: P63 rG = [-9337.05     0.00 128243.78] aCG = [    1.65     0.00     5.22]
t =    237.0 mode: P63 rG = [-9431.16     0.00 127131.16] aCG = [    1.67     0.00     5.23]
t =    238.0 mode: P63 rG = [-9524.52     0.00 126023.77] aCG = [    1.68     0.00     5.24]
t =    239.0 mode: P63 rG = [-9617.09     0.00 124921.61] aCG = [    1.69     0.00     5.25]
t =    240.0 mode: P63 rG = [-9708.89     0.00 123824.69] aCG = [    1.71     0.00     5.26]
t =    241.0 mode: P63 rG = [-9799.89     0.00 122733.03] aCG = [    1.72     0.00     5.26]
t =    242.0 mode: P63 rG = [-9890.10     0.00 121646.62] aCG = [    1.74     0.00     5.27]
t =    243.0 mode: P63 rG = [-9979.51     0.00 120565.48] aCG = [    1.75     0.00     5.28]
t =    244.0 mode: P63 rG = [-10068.10     0.00 119489.62] aCG = [    1.76     0.00     5.29]
t =    245.0 mode: P63 rG = [-10155.88     0.00 118419.04] aCG = [    1.78     0.00     5.29]
t =    246.0 mode: P63 rG = [-10242.83     0.00 117353.75] aCG = [    1.79     0.00     5.30]
t =    247.0 mode: P63 rG = [-10328.95     0.00 116293.76] aCG = [    1.80     0.00     5.31]
t =    248.0 mode: P63 rG = [-10414.23     0.00 115239.07] aCG = [    1.82     0.00     5.31]
t =    249.0 mode: P63 rG = [-10498.67     0.00 114189.69] aCG = [    1.83     0.00     5.32]
t =    250.0 mode: P63 rG = [-10582.25     0.00 113145.63] aCG = [    1.84     0.00     5.33]
t =    251.0 mode: P63 rG = [-10664.97     0.00 112106.89] aCG = [    1.86     0.00     5.33]
t =    252.0 mode: P63 rG = [-10746.83     0.00 111073.49] aCG = [    1.87     0.00     5.34]
t =    253.0 mode: P63 rG = [-10827.81     0.00 110045.42] aCG = [    1.88     0.00     5.35]
t =    254.0 mode: P63 rG = [-10907.91     0.00 109022.69] aCG = [    1.90     0.00     5.35]
t =    255.0 mode: P63 rG = [-10987.12     0.00 108005.31] aCG = [    1.91     0.00     5.36]
t =    256.0 mode: P63 rG = [-11065.43     0.00 106993.29] aCG = [    1.92     0.00     5.36]
t =    257.0 mode: P63 rG = [-11142.85     0.00 105986.62] aCG = [    1.94     0.00     5.37]
t =    258.0 mode: P63 rG = [-11219.36     0.00 104985.32] aCG = [    1.95     0.00     5.37]
t =    259.0 mode: P63 rG = [-11294.95     0.00 103989.38] aCG = [    1.96     0.00     5.38]
t =    260.0 mode: P63 rG = [-11369.62     0.00 102998.82] aCG = [    1.97     0.00     5.38]
t =    261.0 mode: P63 rG = [-11443.36     0.00 102013.64] aCG = [    1.99     0.00     5.38]
t =    262.0 mode: P63 rG = [-11516.17     0.00 101033.84] aCG = [    2.00     0.00     5.39]
t =    263.0 mode: P63 rG = [-11588.03     0.00 100059.43] aCG = [    2.01     0.00     5.39]
t =    264.0 mode: P63 rG = [-11658.95     0.00 99090.41] aCG = [    2.02     0.00     5.40]
t =    265.0 mode: P63 rG = [-11728.91     0.00 98126.78] aCG = [    2.04     0.00     5.40]
t =    266.0 mode: P63 rG = [-11797.91     0.00 97168.54] aCG = [    2.05     0.00     5.40]
t =    267.0 mode: P63 rG = [-11865.95     0.00 96215.71] aCG = [    2.06     0.00     5.41]
t =    268.0 mode: P63 rG = [-11933.00     0.00 95268.29] aCG = [    2.07     0.00     5.41]
t =    269.0 mode: P63 rG = [-11999.08     0.00 94326.27] aCG = [    2.09     0.00     5.41]
t =    270.0 mode: P63 rG = [-12064.17     0.00 93389.66] aCG = [    2.10     0.00     5.41]
t =    271.0 mode: P63 rG = [-12128.27     0.00 92458.47] aCG = [    2.11     0.00     5.42]
t =    272.0 mode: P63 rG = [-12191.37     0.00 91532.68] aCG = [    2.12     0.00     5.42]
t =    273.0 mode: P63 rG = [-12253.46     0.00 90612.32] aCG = [    2.13     0.00     5.42]
t =    274.0 mode: P63 rG = [-12314.55     0.00 89697.38] aCG = [    2.14     0.00     5.42]
t =    275.0 mode: P63 rG = [-12374.61     0.00 88787.85] aCG = [    2.16     0.00     5.42]
t =    276.0 mode: P63 rG = [-12433.65     0.00 87883.75] aCG = [    2.17     0.00     5.42]
t =    277.0 mode: P63 rG = [-12491.66     0.00 86985.08] aCG = [    2.18     0.00     5.43]
t =    278.0 mode: P63 rG = [-12548.63     0.00 86091.82] aCG = [    2.19     0.00     5.43]
t =    279.0 mode: P63 rG = [-12604.57     0.00 85204.00] aCG = [    2.20     0.00     5.43]
t =    280.0 mode: P63 rG = [-12659.45     0.00 84321.60] aCG = [    2.21     0.00     5.43]
t =    281.0 mode: P63 rG = [-12713.28     0.00 83444.63] aCG = [    2.22     0.00     5.43]
t =    282.0 mode: P63 rG = [-12766.05     0.00 82573.09] aCG = [    2.24     0.00     5.43]
t =    283.0 mode: P63 rG = [-12817.76     0.00 81706.98] aCG = [    2.25     0.00     5.43]
t =    284.0 mode: P63 rG = [-12868.40     0.00 80846.30] aCG = [    2.26     0.00     5.43]
t =    285.0 mode: P63 rG = [-12917.96     0.00 79991.04] aCG = [    2.27     0.00     5.43]
t =    286.0 mode: P63 rG = [-12966.44     0.00 79141.22] aCG = [    2.28     0.00     5.43]
t =    287.0 mode: P63 rG = [-13013.83     0.00 78296.82] aCG = [    2.29     0.00     5.43]
t =    288.0 mode: P63 rG = [-13060.13     0.00 77457.85] aCG = [    2.30     0.00     5.43]
t =    289.0 mode: P63 rG = [-13105.33     0.00 76624.30] aCG = [    2.31     0.00     5.42]
t =    290.0 mode: P63 rG = [-13149.43     0.00 75796.19] aCG = [    2.32     0.00     5.42]
t =    291.0 mode: P63 rG = [-13192.42     0.00 74973.49] aCG = [    2.33     0.00     5.42]
t =    292.0 mode: P63 rG = [-13234.30     0.00 74156.22] aCG = [    2.34     0.00     5.42]
t =    293.0 mode: P63 rG = [-13275.06     0.00 73344.37] aCG = [    2.35     0.00     5.42]
t =    294.0 mode: P63 rG = [-13314.70     0.00 72537.94] aCG = [    2.36     0.00     5.42]
t =    295.0 mode: P63 rG = [-13353.21     0.00 71736.92] aCG = [    2.37     0.00     5.41]
t =    296.0 mode: P63 rG = [-13390.58     0.00 70941.32] aCG = [    2.38     0.00     5.41]
t =    297.0 mode: P63 rG = [-13426.82     0.00 70151.14] aCG = [    2.39     0.00     5.41]
t =    298.0 mode: P63 rG = [-13461.92     0.00 69366.36] aCG = [    2.40     0.00     5.41]
t =    299.0 mode: P63 rG = [-13495.86     0.00 68586.99] aCG = [    2.41     0.00     5.40]