#!/bin/bash
#
# DIVYESH VARIYA
# Script for IITB


clear
rm -rf Solved*

for i in 1
do
	mkdir Solved_$i
	
	echo "Run Case with scalarTransportFoam"
	echo "----------------------------------------------"

	for k in upwind linearUpwind linear limitedLinear limitedLimitedLinear vanLeer limitedVanLeer cubic limitedCubic UMIST SFCD QUICK
	do

		if [ $i -eq 1 ]
		then
			echo "  - $k"
			if [ $i -eq 1 ] 
			then
				cp -r Case_Files Solved_$i/$k
			fi

				path="Solved_$i/$k"


				if [ "$k" == "linearUpwind" ]
				then
					sed "s/scheme/Gauss $k grad(T)/g" -i $path*/system/fvSchemes
				elif [ "$k" == "limitedLinear" ] || [ "$k" == "limitedCubic" ]
				then
					sed "s/scheme/Gauss $k 1/g" -i $path*/system/fvSchemes
				elif [ "$k" == "limitedLimitedLinear" ]
				then
					sed "s/scheme/Gauss $k 1 0 1/g" -i $path*/system/fvSchemes
				elif [ "$k" == "limitedVanLeer" ]
				then
					sed "s/scheme/Gauss $k 0 1/g" -i $path*/system/fvSchemes
				else
					sed "s/scheme/Gauss $k/g" -i $path*/system/fvSchemes
				fi

				blockMesh -case $path >> logFile
				
				setFields -case $path >> logFile

				scalarTransportFoam -case $path >> logFile

			 
		fi
done
	echo ""
done

