#!/bin/bash

Substatus="YES"

if [ $Substatus == "YES" ]; then

   for i in  "0.41" # "0.09"  "0.41" "0.81" "1.15" "1.86" 
   do
       rm -rf $i
       mkdir $i
       sed -i "8s/^.*/$i                # initial kinetic energy (eV)/" Input/input.txt 
       sed -i "3s/^.*/$(echo "451.98*$i+259.32" | bc -l)             # nozzle temperature (K)/" Input/input.txt 
       #python geninicon_1.py
       python gen_oriented_inicon_1.py
       cp  Output/ICON.txt Output/jv_QN_Inf.txt Output/veldist.txt $i/ 
       cat Input/input.txt
  done
fi

#for k in "vj00"   "vj01"  "vj10"  "vj20"  "vj30"
#do
#     if [ $Substatus == "YES" ]; then
#
#        jj=$(echo "${k:3:1}")
#        vv=$(echo "${k:2:1}")
#
#        sed -i "16s/^.*/$vv                  # v \(if state-resolved flag is on\)/" Input/input.txt
#        sed -i "17s/^.*/$jj                  # j \(if state-resolved flag is on\)/" Input/input.txt
#    
#        for i in  "1.50" "1.75"  "2.00" "2.25" "2.50" "2.75" "3.00" "3.25" 
#        do
#            if [ $(echo "${i:2:2}") == "00"  ]; then
#                j=$(echo $(echo "${i:0:1}-1" | bc)".95" )    
#            else    
#                j=$(echo "${i:0:2}"$(echo "${i:2:2}-5" | bc) )
#            fi
#            sed -i "8s/^.*/$j               # initial kinetic energy (eV)/" Input/input.txt 
#            python geninicon_1.py
#            cp  Output/ICON.txt Output/jv_QN_Inf.txt Output/veldist.txt $i/
#
#            cp  -r $i $k
#        done
#
#     fi
#done