# 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    5*10^-3 * λ
  zeitschritt    5*10^-3 * T
  maximalimpuls  10
  impulsschritt  10^-2
  zeit           7 * T
!setze $breite: (7 * λ)
  breite         $breite
  mit Fortschrittsanzeige
  Ausgabeverzeichnis aufräumen
allgemeinEnde

ausgaben
  prefix         /nfs/synologynfs/erich/Plasmapropagation/Daten/
  suffix         _test.dat
  felder         AX,AY,dAYDT,EX,BZ,VX1,PX1,PXSqr1,N1,N2,PXRipple1,Rho1,Rho2
  phasenraum     N1,N2
ausgabenEnde

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

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

!setze $IonenMassenFaktor: (1836.15267245 + 1838.68366158)

teilchen1
  ladung             -qe
  masse              me
  maximaldichte      10
  minimaldichte      0 # 10^-2
  maximales dLnN/dX  10
  maximales dLnN/dP  2
  impulsbreite       1/20
!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
  ladung             qe
  masse              $IonenMassenFaktor * me
  maximaldichte      nmax1
  minimaldichte      nmin1
  maximales dLnN/dX  10
  maximales dLnN/dP  2
  impulsbreite       1/20
  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