1 files changed, 475 insertions, 0 deletions
diff --git a/postpost/update_input b/postpost/update_input
new file mode 100755
index 0000000..6a3c414
--- /dev/null
+++ b/postpost/update_input
@@ -0,0 +1,475 @@
+#!/bin/sh
+
+# update_input input.lpi input.lpi.save
+
+path="$(dirname $0)"
+
+# neuer Wert:
+#  okey -> Name des Schlüssels
+#  obam -> Name des neuen Wertes
+#  olin -> einzufügende Zeile
+#  obnam -> Name des darüberliegenden Wertes
+
+okey[0]="lramp"
+onam[0]="form"
+olin[0]="form             = 0         # 0 = linear; 1 = exponential"
+obnam[0]="length"
+
+okey[1]="lramp"
+onam[1]="length"
+olin[1]="length           = 0         # cells in the linear ramp region (form=0) or scalelength"
+obnam[1]="------------------------------------------------------------------------------------------"
+
+okey[2]="lramp"
+onam[2]="length2"
+olin[2]="length2          = 100       # lefter scale length (form=2)"
+obnam[2]="form"
+
+okey[3]="lramp"
+onam[3]="cutoff"
+olin[3]="cutoff           = -1        # max ramplength (<0 <=> disabled)"
+obnam[3]="length2"
+
+okey[4]="lramp"
+onam[4]="gluepos"
+olin[4]="gluepos          = 250       # glue-point between the two exp-functions (form=2) ... length of the righter exp-ramp"
+obnam[4]="cutoff"
+
+okey[5]="rramp"
+onam[5]="length"
+olin[5]="length           = 0         # cells in the linear ramp region (form=0) or scalelength"
+obnam[5]="------------------------------------------------------------------------------------------"
+
+okey[6]="rramp"
+onam[6]="form"
+olin[6]="form             = 0         # 0 = linear; 1 = exponential"
+obnam[6]="length"
+
+okey[7]="rramp"
+onam[7]="length2"
+olin[7]="length2          = 100       # righter scale length (form=2)"
+obnam[7]="form"
+
+okey[8]="rramp"
+onam[8]="cutoff"
+olin[8]="cutoff          = -1        # max ramplength (<0 <=> disabled)"
+obnam[8]="length2"
+
+okey[9]="rramp"
+onam[9]="gluepos"
+olin[9]="gluepos         = 250       # glue-point between the two exp-functions (form=2) ... length of the lefter exp-ramp"
+obnam[9]="cutoff"
+
+okey[10]="box"
+onam[10]="ux0"
+olin[10]="ux0             = 0.0       # initial momentum in x-direction"
+obnam[10]="n_ion_over_nc"
+
+okey[11]="etx"
+onam[11]="Q"
+olin[11]="   Q            = 0         # thermal x-energy plots?"
+obnam[11]="------------------------------------------------------------------------------------------"
+
+okey[12]="etx"
+onam[12]="t_start"
+olin[12]="   t_start      = 0         # start time in periods"
+obnam[12]="Q"
+
+okey[13]="etx"
+onam[13]="t_stop"
+olin[13]="   t_stop       = 60        # stop time in periods"
+obnam[13]="t_start"
+
+okey[14]="etx"
+onam[14]="x_start"
+olin[14]="   x_start      = 0         # left boundary in cells"
+obnam[14]="t_stop"
+
+okey[15]="etx"
+onam[15]="x_stop"
+olin[15]="   x_stop       = 100000    # right boundary in cells"
+obnam[15]="x_start"
+
+okey[16]="ety"
+onam[16]="Q"
+olin[16]="   Q            = 0         # thermal y-energy plots?"
+obnam[16]="------------------------------------------------------------------------------------------"
+
+okey[17]="ety"
+onam[17]="t_start"
+olin[17]="   t_start      = 0         # start time in periods"
+obnam[17]="Q"
+
+okey[18]="ety"
+onam[18]="t_stop"
+olin[18]="   t_stop       = 60        # stop time in periods"
+obnam[18]="t_start"
+
+okey[19]="ety"
+onam[19]="x_start"
+olin[19]="   x_start      = 0         # left boundary in cells"
+obnam[19]="t_stop"
+
+okey[20]="ety"
+onam[20]="x_stop"
+olin[20]="   x_stop       = 100000    # right boundary in cells"
+obnam[20]="x_start"
+
+okey[21]="etz"
+onam[21]="Q"
+olin[21]="   Q            = 0         # thermal z-energy plots?"
+obnam[21]="------------------------------------------------------------------------------------------"
+
+okey[22]="etz"
+onam[22]="t_start"
+olin[22]="   t_start      = 0         # start time in periods"
+obnam[22]="Q"
+
+okey[23]="etz"
+onam[23]="t_stop"
+olin[23]="   t_stop       = 60        # stop time in periods"
+obnam[23]="t_start"
+
+okey[24]="etz"
+onam[24]="x_start"
+olin[24]="   x_start      = 0         # left boundary in cells"
+obnam[24]="t_stop"
+
+okey[25]="etz"
+onam[25]="x_stop"
+olin[25]="   x_stop       = 100000    # right boundary in cells"
+obnam[25]="x_start"
+
+okey[26]="lramp"
+onam[26]="middens"
+olin[26]="middens         = 0.5       # density of the constant portion in the middle (form=3) in fractions of ni"
+obnam[26]="gluepos"
+
+okey[27]="rramp"
+onam[27]="middens"
+olin[27]="middens         = 0.5       # density of the constant portion in the middle (form=3) in fractions of ni"
+obnam[27]="gluepos"
+
+for i in 0 1
+do
+ if [ ${i} -eq 0 ]
+ then
+  s="front"
+ else
+  s="rear"
+ fi
+ for j in 0 1 2 3
+ do
+  okey[28+36*${i}+9*${j}]="pulse_${s}_part_${j}"
+  onam[28+36*${i}+9*${j}]="amplitude"
+  olin[28+36*${i}+9*${j}]="amplitude         = 0                          # dimensionless laser field amplitude"
+  obnam[28+36*${i}+9*${j}]="Q"
+
+  okey[28+1+36*${i}+9*${j}]="pulse_${s}_part_${j}"
+  onam[28+1+36*${i}+9*${j}]="polarization"
+  olin[28+1+36*${i}+9*${j}]="polarization      = 2                          # s=1, p=2, c=3"
+  obnam[28+1+36*${i}+9*${j}]="amplitude"
+
+  okey[28+2+36*${i}+9*${j}]="pulse_${s}_part_${j}"
+  onam[28+2+36*${i}+9*${j}]="shape"
+  olin[28+2+36*${i}+9*${j}]="shape             = 3                          # linear=1, sin=2, sin^2=3"
+  obnam[28+2+36*${i}+9*${j}]="polarization"
+
+  okey[28+3+36*${i}+9*${j}]="pulse_${s}_part_${j}"
+  onam[28+3+36*${i}+9*${j}]="raise"
+  olin[28+3+36*${i}+9*${j}]="raise             = 15                         # pulse raise/fall time in periods"
+  obnam[28+3+36*${i}+9*${j}]="polarization"
+
+  okey[28+4+36*${i}+9*${j}]="pulse_${s}_part_${j}"
+  onam[28+4+36*${i}+9*${j}]="duration"
+  olin[28+4+36*${i}+9*${j}]="duration          = 30                         # pulse duration in fundamental periods"
+  obnam[28+4+36*${i}+9*${j}]="raise"
+
+  okey[28+5+36*${i}+9*${j}]="pulse_${s}_part_${j}"
+  onam[28+5+36*${i}+9*${j}]="delay"
+  olin[28+5+36*${i}+9*${j}]="delay             = 0                          # propagation delay in cycles"
+  obnam[28+5+36*${i}+9*${j}]="duration"
+
+  okey[28+6+36*${i}+9*${j}]="pulse_${s}_part_${j}"
+  onam[28+6+36*${i}+9*${j}]="phase"
+  olin[28+6+36*${i}+9*${j}]="phase             = 0                          # propagation delay in cycles"
+  obnam[28+6+36*${i}+9*${j}]="delay"
+
+  okey[28+7+36*${i}+9*${j}]="pulse_${s}_part_${j}"
+  onam[28+7+36*${i}+9*${j}]="chirp"
+  olin[28+7+36*${i}+9*${j}]="chirp             = 0                          # linear chirp (dw/dt) in 2 pi T^-2"
+  obnam[28+7+36*${i}+9*${j}]="phase"
+
+  okey[28+8+36*${i}+9*${j}]="pulse_${s}_part_${j}"
+  onam[28+8+36*${i}+9*${j}]="frequency"
+  olin[28+8+36*${i}+9*${j}]="frequency         = $[${j}+1]                          # frequency in fundamental frequencies"
+  obnam[28+8+36*${i}+9*${j}]="chirp"
+ done
+done
+
+# Neuer Schlüssel:
+#  nkey -> Name des neuen Schlüssels
+#  nakey -> Name des folgenden Schlüssels
+
+nkey[0]="&rramp"
+nakey[0]="&electrons"
+
+nkey[1]="&lramp"
+nakey[1]="&rramp"
+
+nkey[2]="&etz"
+nakey[2]="&edens"
+
+nkey[3]="&ety"
+nakey[3]="&etz"
+
+nkey[4]="&etx"
+nakey[4]="&ety"
+
+s="&pulse_rear"
+for i in 3 2 1 0
+do
+ nakey[$[8-${i}]]="$s"
+ s="&pulse_front_part_${i}"
+ nkey[$[8-${i}]]="$s"
+done
+
+s="&propagate"
+for i in 3 2 1 0
+do
+ nakey[$[12-${i}]]="$s"
+ s="&pulse_rear_part_${i}"
+ nkey[$[12-${i}]]="$s"
+done
+
+# umzuziehende Werte:
+#  toky -> alter Schlüssel
+#  tnky -> neuer Schlüssel
+#  tonm -> alter Name
+#  tnnm -> neuer Name
+#  tnonm -> Name des über der neuen Position liegenden Wertes
+
+toky[0]="box"
+tnky[0]="lramp"
+tonm[0]="cells_ramp"
+tnnm[0]="length    "
+tnonm[0]="------------------------------------------------------------------------------------------"
+
+toky[1]="box"
+tnky[1]="lramp"
+tonm[1]="cells_ramp_form"
+tnnm[1]="form           "
+tnonm[1]="length"
+twarn[1]="For form=0 the ramp now starts (consistently with form>=1) at cells_left with maximal density! You need to subtract &lramp -> length manually from &box -> cells_left"
+
+toky[2]="box"
+tnky[2]="lramp"
+tonm[2]="cells_ramp2"
+tnnm[2]="length2    "
+tnonm[2]="form"
+
+toky[3]="box"
+tnky[3]="lramp"
+tonm[3]="cells_rampcut"
+tnnm[3]="cutoff       "
+tnonm[3]="form"
+
+toky[4]="box"
+tnky[4]="lramp"
+tonm[4]="cells_gluepos"
+tnnm[4]="gluepos      "
+tnonm[4]="cutoff"
+
+for i in 0 1
+do
+ if [ $i -eq 0 ]
+ then
+  s="front"
+ else
+  s="rear"
+ fi
+ 
+ j=0
+ for t in "delay" "duration" "raise" "shape" "polarization" "amplitude"
+ do
+  toky[5+16*${i}+${j}]="pulse_${s}"
+  tnky[5+16*${i}+${j}]="pulse_${s}_part_0"
+  tonm[5+16*${i}+${j}]="${t}"
+  tnnm[5+16*${i}+${j}]="${t}"
+  tnonm[5+16*${i}+${j}]="------------------------------------------------------------------------------------------"
+  j=$[${j}+1]
+ done
+
+ toky[5+6+16*${i}]="pulse_${s}"
+ tnky[5+6+16*${i}]="pulse_${s}_part_3"
+ tonm[5+6+16*${i}]="seedfrequency"
+ tnnm[5+6+16*${i}]="frequency"
+ tnonm[5+6+16*${i}]="------------------------------------------------------------------------------------------"
+
+ k=1
+ l=0
+ for j in "2" "3" "n"
+ do
+  for t in "delay" "phase" "amplitude"
+  do
+   toky[5+7+16*${i}+${l}]="pulse_${s}"
+   tnky[5+7+16*${i}+${l}]="pulse_${s}_part_${k}"
+   tonm[5+7+16*${i}+${l}]="${t}${j}"
+   tnnm[5+7+16*${i}+${l}]="${t}"
+   tnonm[5+7+16*${i}+${l}]="------------------------------------------------------------------------------------------"
+   l=$[${l}+1]
+  done
+  k=$[${k}+1]
+ done
+done
+
+t1=0
+while [ -e "tmp$t1" ]
+do
+  t1=$[$t1+1]
+done
+t2=$[$t1+1]
+t1="tmp"$t1
+while [ -e "tmp$t2" ]
+do
+  t2=$[$t2+1]
+done
+t2="tmp"$t2
+
+cp $1 $t1
+mv $1 $2
+
+ # einfügen neuer Schlüssel
+for i in {0..12}
+do
+  if [ "$(grep -c "^${nkey[$i]}" $t1)" = "0" ]
+  then
+    alles=$(grep -c ".*" $t1)
+    stelle=$[$(grep -n "^${nakey[$i]}" $t1 | tr ":" " " | awk '{print $1}')-1]
+    head -n${stelle} $t1 > $t2
+    echo "${nkey[$i]}" >> $t2
+    echo "------------------------------------------------------------------------------------------" >> $t2
+    echo "" >> $t2
+    echo "" >> $t2
+    tail -n$[$alles-$stelle] $t1 >> $t2
+    rm $t1
+    mv $t2 $t1
+  fi
+  echo -ne "."
+done
+
+ # umziehen umzuziehender Werte
+for i in {0..36}
+do
+  if [ "$(${path}/finde_wert $t1 ${tnky[$i]} ${tnnm[$i]})" = "0" ] && [ "$(${path}/finde_wert $t1 ${toky[$i]} ${tonm[$i]})" != "0" ]
+  then
+    alles=$(grep -c ".*" $t1)
+    astelle=$(${path}/finde_wert $t1 ${toky[$i]} ${tonm[$i]})
+    nstelle=$(${path}/finde_wert $t1 ${tnky[$i]} ${tnonm[$i]})
+    zeile="$(head -n$astelle $t1 | tail -n1 | sed "s/${tonm[$i]}/${tnnm[$i]}/")"
+    if [ $astelle -lt $nstelle ]
+    then
+      head -n$[${astelle}-1] $t1 > $t2
+      head -n${nstelle} $t1 | tail -n$[${nstelle}-${astelle}] >> $t2
+      echo "${zeile}" >> $t2
+      tail -n$[$alles-$nstelle] $t1 >> $t2
+    else
+      head -n${nstelle} $t1 > $t2
+      echo "${zeile}" >> $t2
+      head -n$[${astelle}-1] $t1 | tail -n$[${astelle}-${nstelle}-1] >> $t2
+      tail -n$[$alles-$astelle] $t1 >> $t2
+    fi
+    rm $t1
+    mv $t2 $t1
+    if [ "${twarn[$i]}" != "" ]
+    then
+      echo
+      echo "Warning: ${twarn[$i]}"
+    fi
+  fi
+  echo -ne "."
+done
+
+ # einfügen neuer Werte
+for i in {0..99}
+do
+  if [ "$(${path}/finde_wert $t1 ${okey[$i]} ${onam[$i]})" == "0" ]
+  then
+    alles=$(grep -c ".*" $t1)
+    stelle=$(${path}/finde_wert $t1 ${okey[$i]} ${obnam[$i]})
+    head -n${stelle} $t1 > $t2
+    echo "${olin[$i]}" >> $t2
+    tail -n$[$alles-$stelle] $t1 >> $t2
+    rm $t1
+    mv $t2 $t1
+  fi
+  echo -ne "."
+done
+
+ # Ende auskommentieren
+alles=$(grep -c ".*" $t1)
+stelle=$(grep -n "^#\{0,1\}thermal velocity v/c vs. thermal energy" $t1 | sed "s/:.*$//" | tail -n1)
+head -n$[${stelle} - 1] $t1 > $t2
+tail -n$[${alles} - ${stelle} + 1] $t1 | sed "s/^/#/" | sed "s/^##/#/" | sed "s/^#\/\//\/\//" | sed "s/^#$//" >> $t2
+mv $t2 $t1
+
+# =s ausrichten
+minpos=$(
+for s in $(grep "^[^#]*=" $t1 | sed "s/ *=.*$//" | sort -u | tr " " "_")
+do
+  echo ${#s}
+done | sort -n | tail -n1)
+
+alles=$(grep -c ".*" $t1)
+i=1
+while [ ${i} -le ${alles} ]
+do
+  s="$(head -n${i} $t1 | tail -n1)"
+  if [ $(echo "$s" | grep -c "^[^=#]*=") -eq 1 ]
+  then
+    t="$(echo "$s" | sed "s/ *=.*$//")"
+    while [ ${#t} -lt ${minpos} ]
+    do
+      t="${t} "
+    done
+    echo "$t =$(echo "$s" | sed "s/^[^=]*=//")" >> $t2
+  else
+    echo "$s" >> $t2
+  fi
+  i=$[${i}+1]
+  echo -ne "."
+done
+mv $t2 $t1
+
+# #s ausrichten
+minpos=$(
+for s in $(grep "^[^#][^#]*#" $t1 | sed "s/ *#.*$//" | sort -u | tr " " "_")
+do
+  echo ${#s}
+done | sort -n | tail -n1)
+
+alles=$(grep -c ".*" $t1)
+i=1
+while [ ${i} -le ${alles} ]
+do
+  s="$(head -n${i} $t1 | tail -n1)"
+  if [ $(echo "$s" | grep -c "^[^#][^#]*#") -eq 1 ]
+  then
+    t="$(echo "$s" | sed "s/ *#.*$//")"
+    while [ ${#t} -lt ${minpos} ]
+    do
+      t="${t} "
+    done
+    echo "$t #$(echo "$s" | sed "s/^[^#]*#//")" >> $t2
+  else
+    echo "$s" >> $t2
+  fi
+  i=$[${i}+1]
+  echo -ne "."
+done
+echo
+
+mv $t2 $t1
+
+mv $t1 $1
+