Plasmadichte aus Sergeys Trajektorien extrahieren

1 files changed, 13 insertions, 5 deletions

diff --git a/epostunit.pas b/epostunit.pas
index cf34038..1c446c5 100644
--- a/epostunit.pas
+++ b/epostunit.pas
@@ -1419,15 +1419,23 @@ begin
 
       sl:=tMyStringlist.create;
       sl.loadFromFile(dateien[i].Name+'../xrom.ini');
-      sl.grep('^\s*Pulse[12]\.a0\s*=.*[^0 ].*$');
-      if sl.count=0 then begin
+      t:=sl.grepFirst('^\s*Pulse[12]\.a0\s*=.*[^0 ].*$');
+      if t='' then begin
         gibAus(''''+s+''' enthält keinen Puls, den ich erkenne!',3);
         exit;
       end;
-      s:=sl[0];
+      erstesArgument(t,'Pulse');
+      dateien[i].groeszenFaktor:=strtoint(erstesArgument(t,'.'));
+
+      t:=sl.grepFirst('^\s*Plasma\.Density\s*=');
+      if t='' then begin
+        gibAus(''''+s+''' enthält keine Plasmadichte, die ich erkenne!',3);
+        exit;
+      end;
+      erstesArgument(t,'=');
+      dateien[i].params.np:=sqrt(strtofloat(t))/dateien[i].groeszenFaktor;
+
       sl.free;
-      erstesArgument(s,'Pulse');
-      dateien[i].groeszenFaktor:=strtoint(erstesArgument(s,'.'));
 
       assign(tf,dateien[i].Name+'dt');
       reset(tf);