written moments should live at same time as fields and particles

[alecjohnson]

Currently, at the end of each cycle, we output the field and particle data at time tⁿ, but we output the moment data for time step t^n-1. Instead, at the end of each output cycle, all quantities saved should live at time step tⁿ.

The current algorithm cycle looks like this:

  for (int i = KCode.FirstCycle(); i < KCode.LastCycle(); i++) 
  {
      // computes moments M at t^(n-1) and implicit moments IM at t^n
      KCode.CalculateMoments();
      KCode.CalculateField(); // computes E at t^n
      KCode.ParticlesMover(); // computes pcls P at t^n
      KCode.CalculateB(); // computes B at t^n
      KCode.WriteOutput(i); // writes M, E, B, and P
  }   

To fix this problem, we can change the algorithm cycle to look like this:

  // computes moments M at t^0 and implicit moments IM at t^1
  KCode.CalculateMoments();
  for (int i = KCode.FirstCycle(); i < KCode.LastCycle(); i++) 
  {
      KCode.CalculateField(); // computes E at t^n
      KCode.ParticlesMover(); // computes pcls P at t^n
      KCode.CalculateB(); // computes B at t^n
      // computes moments M at t^n and implicit moments IM at t^(n+1)
      KCode.CalculateMoments();
      KCode.WriteOutput(i);  // writes M, E, B, and P
  }