X-Git-Url: https://hackdaworld.org/gitweb/?p=physik%2Fposic.git;a=blobdiff_plain;f=runmd;fp=runmd;h=609bcf316e8abdd5f3df8638582027df4f3014c8;hp=d50d9ae39867107745d1b1b24c5eb25928318c6f;hb=6e6d7126ea9a845f11637d8e1b8eb2b570ac4dc9;hpb=97dc63eb6a519b8e1f4fbfaa9760dd94539436b0

diff --git a/runmd b/runmd
index d50d9ae..609bcf3 100755
--- a/runmd
+++ b/runmd
@@ -8,13 +8,15 @@ if [ ! -f ./config ]; then
 	exit
 fi
 
-[ ! -d $1 ] && mkdir $1
+[ ! -d $1 ] && mkdir -p $1
 
 ./clean $1
 
 cp -v config $1/config
 
-time ./mdrun -c ./config -s $1
+mkdir -p logs
+logfile=logs/run_`basename $1`.log
+time ./mdrun -c ./config -s $1 | tee $logfile
 
 if [ "$?" == "0" ]; then
 	#./perms
@@ -22,9 +24,28 @@ if [ "$?" == "0" ]; then
 		# whole simulation cell
 		#./visualize -w 640 -h 480 -d $1
 
-		# center unit cell
+		# crt
+		mkdir -p $1/crt
+		check=`grep crt\ init $logfile`
+		if [ "$check" = "  crt init" ]; then
+			grep crt\ energy $logfile | \
+				awk -F: '{ print $2 }' | \
+				awk -F\ -\  '{ print $1 " " $2 }' > \
+				$1/crt/energy
+			for i in $1/s-crt_*.save; do
+				nr=`basename $i | awk -F_ '{ print $2 }' | \
+					awk -F. '{ print $1 }'`
+				./s2xyz $i > $1/crt/atomic_conf_${nr}.xyz
+			done
+			./visualize -w 640 -h 480 -d $1/crt \
+		       		    -nll -0.56 -0.56 -0.76 -fur 0.56 0.56 0.56 \
+				    -b -0.5 -0.5 -0.5 0.5 0.5 0.5 \
+				    -c -0.2 -2.0 0.6 -L 0 0 -0.1 \
+				    -r 0.6 -B 0.1
+		fi
+
 		./visualize -w 640 -h 480 -d $1 \
-		            -nll -0.56 -0.56 -0.56 -fur 0.56 0.56 0.56 \
+		            -nll -0.56 -0.56 -0.76 -fur 0.56 0.56 0.56 \
 			    -b -0.5 -0.5 -0.5 0.5 0.5 0.5 \
 			    -c -0.2 -2.0 0.6 -L 0 0 -0.1 \
 			    -r 0.6 -B 0.1
@@ -32,5 +53,7 @@ if [ "$?" == "0" ]; then
 		# old rasmol
 		#rasmol -32 -nodisplay < $1/visualize.scr > /dev/null 2>&1
 		./ppm2avi $1
+
 	fi
 fi
+