summaryrefslogtreecommitdiff
path: root/input.plap
blob: 94a7c54712ad39ec2abb69c136004e1a595406fd (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
# Parameterdatei für Plasmapropagation

allgemein
#  runge-Kutta-14
#  runge-Kutta-12
#  runge-Kutta-10
  runge-Kutta-3/8
#  runge-Kutta-4
#  euler-Vorwärts
  ortsschritt    2*10^-3 * λ
  zeitschritt    2*10^-3 * T
  minZeitschritt 2*10^-3 * T
  zeit           40 * T
  diffusionsterm 2                # max. p/Lp
!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,N1,dPsiDX1,Gamma1,N2,dPsiDX2,Gamma2
ausgabenEnde

!setze $tFwhm:   (2.5 * T)
!setze $tMitte:  (1 * T)

# licht von links      0.1 * 2^(-2*((t-$tMitte)/$tFwhm)^2) * ω * cos(ω*t)    # Zeitableitung des A-Feldes

!setze $IonenMassenFaktor: (1836.15267245 + 1838.68366158)

teilchen1
  spezifische Ladung -q/me
  maximaldichte      10
  minimaldichte      0 # 10^-2
!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
    1
      ($breite+$profilbreite)/2
    sin((x - ($breite+$profilbreite)/2 + $randbreite)*π/2/$randbreite)^2
      ($breite+$profilbreite)/2 + $randbreite
    0
  stückweiseEnde
teilchen1Ende

teilchen2
  spezifische Ladung q/$IonenMassenFaktor/me
  maximaldichte      nmax1*$IonenMassenFaktor
  minimaldichte      nmin1*$IonenMassenFaktor
#  verteilung wie teilchen1
  verteilung stückweise
    0
      (0.1*λ + $breite-$profilbreite)/2 - $randbreite
    sin((x - (0.1*λ + $breite-$profilbreite)/2 - $randbreite)*π/2/$randbreite)^2
      (0.1*λ + $breite-$profilbreite)/2
    1
      (0.1*λ + $breite+$profilbreite)/2
    sin((x - (0.1*λ + $breite+$profilbreite)/2 + $randbreite)*π/2/$randbreite)^2
      (0.1*λ + $breite+$profilbreite)/2 + $randbreite
    0
  stückweiseEnde
teilchen2Ende

Dateiende