From a2fe917e7cf1a6e9eacf92dccabe4114fd145332 Mon Sep 17 00:00:00 2001
From: Erich Eckner <git@eckner.net>
Date: Wed, 29 Jul 2015 16:13:41 +0200
Subject: Tagesendstand, alte (= inhaltliche) Bugs sind noch zu beseitigen

---
 input.plap | 39 ++++++++++++++++++---------------------
 1 file changed, 18 insertions(+), 21 deletions(-)

(limited to 'input.plap')

diff --git a/input.plap b/input.plap
index 34fcf17..5e98ccc 100644
--- a/input.plap
+++ b/input.plap
@@ -5,48 +5,45 @@ allgemein
   ortsschritt    10^-2 * λ
   zeitschritt    10^-2 * T
   zeit           20 * T
-  breite         5 * λ
+!setze $breite: (5 * λ)
+  breite         $breite
   mit Fortschrittsanzeige
 allgemeinEnde
 
 ausgaben
   prefix         /home_raid/erich/Dokumente/Prograemmchen/Plasmapropagation/Daten/
   suffix         _test.dat
-  felder         AX,AY,dAYDT,dPhiDX
-#  felder         AX,AY,dAYDT,N,dPhiDX,dPsiDX,Gamma
+  felder         AX,AY,dAYDT,dPhiDX,N1,dPsiDX1,Gamma1,N2,dPsiDX2,Gamma2
 ausgabenEnde
 
 !setze $tFwhm:   (2.5 * T)
 !setze $tMitte:  (1 * T)
 
-licht von links    0.5 * 2^(-2*((t-$tMitte)/$tFwhm)^2) * sin(ω*t)
-
-Dateiende
+licht von links      0.5 * 2^(-2*((t-$tMitte)/$tFwhm)^2) * sin(ω*t)
 
 teilchen1
-  ladung         -q
-  masse          me
-  maximaldichte  0
-  beweglich
-  !setze $profilbreite: 5 λ
-  !setze $randbreite: .1 λ
+  spezifische Ladung -q/me
+  maximaldichte      10
+!setze $profilbreite: (4 * λ)
+!setze $randbreite:   (0.1 * λ)
   verteilung stückweise
     0
-      (breite-$profilbreite)/2 - $randbreite
-    sin((x - (breite-$profilbreite)/2 - $randbreite)*π/2/$randbreite)^2
-      (breite-$profilbreite)/2
+      ($breite-$profilbreite)/2 - $randbreite
+    sin((x - ($breite-$profilbreite)/2 - $randbreite)*π/2/$randbreite)^2
+      ($breite-$profilbreite)/2
     1
-      (breite+$profilbreite)/2
-    sin((x - (breite+$profilbreite)/2 + $randbreite)*π/2/$randbreite)^2
-      (breite+$profilbreite)/2 + $randbreite
+      ($breite+$profilbreite)/2
+    sin((x - ($breite+$profilbreite)/2 + $randbreite)*π/2/$randbreite)^2
+      ($breite+$profilbreite)/2 + $randbreite
     0
   stückweiseEnde
 teilchen1Ende
 
 teilchen2
-  ladung         q
-  maximaldichte  nmax1
-  unbeweglich
+  spezifische Ladung (q/me)/2
+  maximaldichte      nmax1*2
   verteilung wie teilchen1
 teilchen2Ende
 
+Dateiende
+
-- 
cgit v1.2.3-54-g00ecf