X-Git-Url: https://hackdaworld.org/gitweb/?a=blobdiff_plain;f=calc_delta_e;h=4ff46fb9f9084ba2a2fabeb3c72d53fd2a030299;hb=5680e09a2c7c8761d4c3eac199ee58ce20655b92;hp=bde20da40bdda2c4c9ea71947ece59c90d087561;hpb=0d2926a23ea7bbf7654da88770c827ad041357a5;p=physik%2Fposic.git

diff --git a/calc_delta_e b/calc_delta_e
index bde20da..4ff46fb 100755
--- a/calc_delta_e
+++ b/calc_delta_e
@@ -1,14 +1,9 @@
 #!/bin/bash
 
-# change these settings:
-atom_cnt=1000
-inj_time=10.00
-
+file=`ls $1/atomic_conf_* | tail -1`
+atom_cnt=`grep '# \[P\]' $file | awk '{ print $3 }'`
 e0=`grep ^0 $1/energy | awk '{ print $4 }'`
-e1=`grep ^${inj_time} $1/energy | awk '{ print $4 }'`
-e2=`tail -n 1 $1/energy | awk '{ print $4 }'`
+e1=`tail -n 1 $1/energy | awk '{ print $4 }'`
 
-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" }'
+echo "$e0 $e1 $atom_cnt" | \
+	awk '{ print ($2-$1)*$3" eV" }'