#!/bin/bash

poscar=0

while [ $poscar -lt 10 ]
do
	./MakeSlab-Temperature-032.py > Slab.out

	p=`grep 'Potential' Slab.out | awk '{print $9}'`
	k=`grep 'Kinetic'   Slab.out | awk '{print $10}'`

	P=`echo "scale=0; $p / 1" | bc`		# consider only the integer part
	K=`echo "scale=0; $k / 1" | bc`		#

	echo -n	"Potential: ${P} K, Kinetic: ${K} K	"

	if [[ $K -gt 530 && $K -lt 570 ]]
	then
		if [[ $P -gt 530 && $P -lt 570 ]]
		then
		let poscar++
		echo -n "... accepted"
		mkdir "$poscar"

		echo "Cu-111-, 550K surface equilibration"        > "$poscar"/POSCAR
		echo "1.000000000"                               >> "$poscar"/POSCAR
		echo "`grep -A4 Supercell Slab.out | tail -n 3`" >> "$poscar"/POSCAR
		echo "   Cu"					 >> "$poscar"/POSCAR
		echo "   48"                                     >> "$poscar"/POSCAR
		echo "Selective dynamics"                        >> "$poscar"/POSCAR
		echo "Cartesian"                                 >> "$poscar"/POSCAR
		for i in {1..36}
		do
			echo "`grep -A 49 'Coordinates with temperature' Slab.out | tail -n+3 | head -n$i | tail -n1` T T T" >> "$poscar"/POSCAR
		done
		for i in {37..48}
		do 
                        echo "`grep -A 49 'Coordinates with temperature' Slab.out | tail -n+3 | head -n$i | tail -n1` F F F" >> "$poscar"/POSCAR
		done
		
		echo ""                                                                 >> "$poscar"/POSCAR
		echo "`grep -A 49 'Velocities' Slab.out | tail -n+3`"                   >> "$poscar"/POSCAR

		echo "$poscar $K $P" >> check.dat
		cp INPUTS/* "$poscar"

		sed -i s/_JOBNAME_/surf_eq_${poscar}/g ${poscar}/Job

		cd ${poscar}
		qsub Job >JobID.txt
		echo "	`cat JobID.txt`"

		cd ../
		else
			echo -en "... rejected\r"
		fi
	else
		echo -en "... rejected\r"
	fi
done	
rm Slab.out