#! /bin/bash 

#path to the box model
box=/home/valorso/BOXMOD/


if [ $# -ne 2 ] 
then
  echo argument missing
  exit
fi

echo
echo --------------------------------
echo writing files for simulation ...
echo --------------------------------
echo ... cat fort.21 fort.17 '>' $1.mech
cat fort.21 fort.17 > $1.mech

echo ... cp fort.7 $1.dict
cp fort.7 $1.dict

echo ... cp pvap.dat $1.sat
cp pvap.dat $1.sat

echo ... cp fort.55 $1.55  AND  cp fort.8 $1.8
cp fort.55 $1.55
cp fort.8 $1.8

echo
echo --------------------------------
echo compteur ...
echo --------------------------------
cd COMPTEUR
rm fort.*
cp ../fort.7 ./
./count
cd ..

echo
echo --------------------------------
echo searching code name of the input formula
echo --------------------------------
name=$(./findname2)
echo code name found : $name

echo
echo --------------------------------
echo move files .....
echo --------------------------------
echo ... mv $1'.*' $box'CHEMDAT'
mv $1.* $box'CHEMDAT'

# echo copying files in BOXMOD/CHEMDAT/COMPTEUR/$1 ...

if test -d $box'CHEMDAT/COMPTEUR/'$2
  then
  echo
  echo '!' EXIT : le repertoire $2 existe déja dans CHEMDAT/COMPTEUR
  exit 
#  rm -r $box'CHEMDAT/COMPTEUR/'$2
fi

echo ... mkdir $box'CHEMDAT/COMPTEUR/'$2
mkdir $box'CHEMDAT/COMPTEUR/'$2

echo ... cp COMPTEUR/'X*' $box'CHEMDAT/COMPTEUR/'$2'/'
cp COMPTEUR/X* $box'CHEMDAT/COMPTEUR/'$2'/'

echo
echo ---------------------------------
echo interpreting $1.mech
echo ---------------------------------
cd $box'CHEMDAT'
./runintp.unix $1 > $1_interp.out

echo
echo ---------------------------------
echo creating $2 repertory for simulation
echo ---------------------------------

cd ..

if test -d $box$2
  then
  echo 
  echo '!' EXIT : le repertoire $2 existe déja dans $box
  exit
#  rm -r $box'CHEMDAT/COMPTEUR/'$2
fi

echo ... cp -r $box'XXXXX' $box$2
cp -r XXXXX $2

cd $2

echo ... modifying start.unix
sed -i s/xxxxx/$1/g start.unix
sed -i s/XXXXX/$2/g start.unix

exit

#echo ... modifying start.key
#sed -i "s/xxxxx/$name/g" start.key
#echo
#echo enter the initial concentration of the parent species in ppb
#read numpp
#let numpp2=$numpp*25
#sed -i "s/XX/$numpp2/g" start.key

#echo ... modifying plot.key
#sed -i "s/xxxxx/$name/g" plot.key
#echo
#echo enter the number of carbon of the parent species
#read numc
#let numc2=$numc*4
#sed -i "s/XX/$numc2/g" plot.key

echo ... done

echo
echo ----------------------------------
echo launch simulation :
echo ----------------------------------
echo enter yes for launching the simulation, anything else for no?
read answer
if [ $answer = yes ]
  then
  ulimit -s unlimited
  ./start.unix > simu_$1.out &
  echo ... simulation launched !
  else 
  exit
fi

echo
echo ---------------------------------
echo '!' WARNING
echo parent species name must be modified in the extracting programs
echo