X-Git-Url: https://hackdaworld.org/gitweb/?a=blobdiff_plain;f=calc_delta_e;fp=calc_delta_e;h=bde20da40bdda2c4c9ea71947ece59c90d087561;hb=e34be948a7efef2e183a68912046f529ff2e75cf;hp=9c5ecdafe8955c0dde7ab58345b59a57fb0de617;hpb=fd4ecbff76b2e4fc70ca88e64fb76c85fa953888;p=physik%2Fposic.git

diff --git a/calc_delta_e b/calc_delta_e
index 9c5ecda..bde20da 100755
--- a/calc_delta_e
+++ b/calc_delta_e
@@ -1,9 +1,14 @@
 #!/bin/bash
 
+# change these settings:
+atom_cnt=1000
+inj_time=10.00
+
 e0=`grep ^0 $1/energy | awk '{ print $4 }'`
-e1=`tail -n 1 $1/energy | awk '{ print $4 }'`
-ev=`grep ^19.000 $1/energy | awk '{ print $4 }'`
-en=`grep ^20.000 $1/energy | awk '{ print $4 }'`
+e1=`grep ^${inj_time} $1/energy | awk '{ print $4 }'`
+e2=`tail -n 1 $1/energy | awk '{ print $4 }'`
 
-echo "$e0 $e1 $2 $ev $en" | \
-	awk '{ print $5*($3+1)-$4*$3-($2*($3+1)-$1*$3) " eV" }'
+echo "$e0 $e1 $e2 $atom_cnt" | \
+	#awk '{ print ($3*($4+1))-($4*$1)" eV" }'
+	awk '{ print ($3*($4+1))-(($4+1)*$1)" eV" }'
+	#awk '{ print ($3-$1)*$4" eV" }'