program Plasmapropagation;

{$mode objfpc}{$H+}

uses
  {$IFDEF UNIX}{$IFDEF UseCThreads}
  cthreads,
  {$ENDIF}{$ENDIF}
  Classes, SysUtils, CustApp,
  { you can add units after this }
  math, Physikunit, crt, protokollunit;

var ErrorCode: longint;

type

  { TPlasmapropagation }

  TPlasmapropagation = class(TCustomApplication)
  protected
    procedure DoRun; override;
  public
    constructor Create(TheOwner: TComponent); override;
    destructor Destroy; override;
  end;

{ TPlasmapropagation }

procedure TPlasmapropagation.DoRun;
var
  Breite,i,j:                 longint;
  simulation:                 tSimulation;
  start,zeitDatei,zeitPhysik: extended;
  FI:                         TFeldInhalt;
  Abl:                        Boolean;
  Prot:                       TProtokollant;
  c:                          char;
begin
  prot:=tProtokollant.create('error');

  if paramcount<>1 then begin
    prot.schreibe('Bitte genau einen Parameter übergeben, nämlich die Parameterdatei!',true);
    halt(1);
  end;

  start:=now;


  zeitDatei:=0;
  zeitPhysik:=0;
  sT:=-Min(deltaT,sDT)/2;
  simulation:=tSimulation.create(paramstr(1));
//  Gitter.Al:=@InputFeld;
  while Gitter.t<Endzeit do begin
    if HasOption('F','Fortschrittsanzeige') then begin
      if errorCode<2 then
        s:=Gitter.gibErhaltungsgroessen;
      if (floor(100*Gitter.t/Endzeit) < floor(100*(Gitter.t+deltaT)/Endzeit)) or keypressed then begin
        if keypressed then c:=readkey
        else c:=#0;
        case c of
          #27,'q': begin
            errorCode:=3;
            break;
          end;
          ' ': begin
            writeln(' ... Pause (beliebige Taste drücken um fortzufahren) ...');
            readkey;
            writeln(' ... weiter geht''s ...');
          end;
          else begin
            Prot.schreibe(inttostr(round(100*Gitter.t/Endzeit))+'% (t='+floattostr(Gitter.t)+'T)',true);
            Prot.schreibe(timetostr(now-start)+' ('+floattostr(zeitPhysik/max(1e-11,zeitPhysik+zeitDatei))+')',true);
            Prot.schreibe('ETA: '+timetostr((now-start)*(Endzeit-Gitter.t)/max(Gitter.t,deltaT)),true);
            Prot.schreibe('aktueller Zeitschritt: '+floattostr(deltaT)+'T',true);
            Prot.schreibe(s);
          end;
        end{of case};
      end;
    end;

    zeitPhysik:=zeitPhysik-now;
    if errorCode<2 then
      Gitter.iteriereSchritt(deltaT);
    zeitPhysik:=zeitPhysik+now;
    zeitDatei:=zeitDatei-now;
    while Gitter.t>=sT do begin
      sT:=sT+sDT;
      for j:=0 to length(Ausgabedateien)-1 do
        Gitter.macheAusgabe(Ausgabedateien[j],sDX);
    end;
    zeitDatei:=zeitDatei+now;
  end;

  case errorCode of
    2: Prot.schreibe('Simulation wurde auf halbem Wege abgebrochen wegen Überlaufs.',true);
    3: Prot.schreibe('Simulation wurde auf halbem Wege abgebrochen wegen Benutzereingriffs.',true);
  end{of case};

  for i:=0 to length(Ausgabedateien)-1 do
    CloseFile(Ausgabedateien[i].Datei);

  Prot.schreibe('fertig!',true);
  s:=timetostr(now-start);
  t:=timetostr(zeitDatei);
  u:=timetostr(zeitPhysik);
  while length(s)<max(length(t),length(u)) do s:=' '+s;
  while length(t)<max(length(s),length(u)) do t:=' '+t;
  while length(u)<max(length(s),length(t)) do u:=' '+u;
  Prot.schreibe('Das hat '+s+' gedauert,',true);
  Prot.schreibe('  davon '+t+' für Dateizugriffe',true);
  Prot.schreibe('und nur '+u+' für die eigentliche Physik!',true);

  Gitter.Free;
  Prot.Free;
  // stop program loop
  Terminate;
  if errorCode=1 then errorCode:=0;
end;

constructor TPlasmapropagation.Create(TheOwner: TComponent);
begin
  inherited Create(TheOwner);
  errorCode:=1;
  StopOnException:=True;
end;

destructor TPlasmapropagation.Destroy;
begin
  inherited Destroy;
end;

var
  Application: TPlasmapropagation;
begin
  Application:=TPlasmapropagation.Create(nil);
  Application.Run;
  Application.Free;
  Halt(errorCode);
end.