dumpErhaltungsgroeszen verträgt jetzt auch, wenn keine Teilchen da sind

1 files changed, 11 insertions, 9 deletions

diff --git a/Physikunit.pas b/Physikunit.pas
index 254c8d4..d546621 100644
--- a/Physikunit.pas
+++ b/Physikunit.pas
@@ -939,16 +939,18 @@ var
   i,j:  longint;
   dens: double;
 begin
-  dens:=0;
-  for i:=0 to length(massen)-1 do
-    dens:=dens+massen[i]/teilchen[i].eigenschaften[tsgMasse];
-  gitter.prot.schreibe('Gesamtdefizit: '+floattostr(gesamtDefizit)+' (Anteil '+floattostr(gesamtDefizit/dens*dX*dP)+')',true);
-  for i:=0 to length(massen)-1 do begin
+  if length(massen)>0 then begin
     dens:=0;
-    for j:=0 to aX-1 do
-      dens:=dens+impulsIntegral(j,i,msN);
-    dens:=dens*teilchen[i].eigenschaften[tsgMasse]*gitter.dX;
-    gitter.prot.schreibe('n['+inttostr(i+1)+'] = '+floattostr(dens)+' (relative Abweichung: '+floattostr(dens/massen[i]-1)+')',true);
+    for i:=0 to length(massen)-1 do
+      dens:=dens+massen[i]/teilchen[i].eigenschaften[tsgMasse];
+    gitter.prot.schreibe('Gesamtdefizit: '+floattostr(gesamtDefizit)+' (Anteil '+floattostr(gesamtDefizit/dens*dX*dP)+')',true);
+    for i:=0 to length(massen)-1 do begin
+      dens:=0;
+      for j:=0 to aX-1 do
+        dens:=dens+impulsIntegral(j,i,msN);
+      dens:=dens*teilchen[i].eigenschaften[tsgMasse]*gitter.dX;
+      gitter.prot.schreibe('n['+inttostr(i+1)+'] = '+floattostr(dens)+' (relative Abweichung: '+floattostr(dens/massen[i]-1)+')',true);
+    end;
   end;
 end;