#!/bin/sh

cd Re50
./run_mesh
touch Re50.foam

sed -i 's/simulation_type.*/simulation_type        laminar;/g' ./0.orig/include/initialConditions
sed -i 's/flowModel.*/flowModel        laminar;/g' ./0.orig/include/initialConditions
sed -i 's/flowVelocity.*/flowVelocity               (3.2979e-03 0.0 0.0);/g' ./0.orig/include/initialConditions
sed -i 's/flowUbar.*/flowUbar          (3.2979e-03 0.0 0.0);/g' ./0.orig/include/initialConditions
sed -i 's/turbulentKE.*/turbulentKE                 5.19e-09;/g' ./0.orig/include/initialConditions
sed -i 's/turbulentOmega.*/turbulentOmega           3.25e-06;/g' ./0.orig/include/initialConditions
sed -i 's/turbulentEpsilon.*/turbulentEpsilon       4.80e-15;/g' ./0.orig/include/initialConditions
sed -i 's/turbulentNut.*/turbulentNut               4.29e+01;/g' ./0.orig/include/initialConditions
sed -i 's/turbulentNuTilda.*/turbulentNuTilda       6.72e-05;/g' ./0.orig/include/initialConditions
cp -rf 0.orig 0
./run_solver
cd ..


# mkdir Re70
# cd Re70
# touch Re70.foam

# cp -rf ../Re50/0.orig 0.orig
# cp -rf ../Re50/constant constant
# cp -rf ../Re50/system system
# cp -rf ../Re50/run_mesh run_mesh
# cp -rf ../Re50/run_solver run_solver

# sed -i 's/simulation_type.*/simulation_type        laminar;/g' ./0.orig/include/initialConditions
# sed -i 's/flowModel.*/flowModel        laminar;/g' ./0.orig/include/initialConditions
# sed -i 's/flowVelocity.*/flowVelocity               (4.6171e-03 0.0 0.0);/g' ./0.orig/include/initialConditions
# sed -i 's/flowUbar.*/flowUbar          (4.6171e-03 0.0 0.0);/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentKE.*/turbulentKE                 6.77e-09;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentOmega.*/turbulentOmega           3.72e-06;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentEpsilon.*/turbulentEpsilon       7.17e-15;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentNut.*/turbulentNut               4.29e+01;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentNuTilda.*/turbulentNuTilda       6.72e-05;/g' ./0.orig/include/initialConditions
# cp -rf 0.orig 0
# ./run_solver
# cd ..

# mkdir Re100
# cd Re100
# touch Re100.foam

# cp -rf ../Re50/0.orig 0.orig
# cp -rf ../Re50/constant constant
# cp -rf ../Re50/system system
# cp -rf ../Re50/run_mesh run_mesh
# cp -rf ../Re50/run_solver run_solver

# sed -i 's/simulation_type.*/simulation_type        laminar;/g' ./0.orig/include/initialConditions
# sed -i 's/flowModel.*/flowModel        laminar;/g' ./0.orig/include/initialConditions
# sed -i 's/flowVelocity.*/flowVelocity               (6.5958e-03 0.0 0.0);/g' ./0.orig/include/initialConditions
# sed -i 's/flowUbar.*/flowUbar          (6.5958e-03 0.0 0.0);/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentKE.*/turbulentKE                 6.77e-09;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentOmega.*/turbulentOmega           3.72e-06;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentEpsilon.*/turbulentEpsilon       7.17e-15;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentNut.*/turbulentNut               4.29e+01;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentNuTilda.*/turbulentNuTilda       6.72e-05;/g' ./0.orig/include/initialConditions
# cp -rf 0.orig 0
# ./run_solver
# cd ..

# mkdir Re150
# cd Re150
# touch Re150.foam

# cp -rf ../Re50/0.orig 0.orig
# cp -rf ../Re50/constant constant
# cp -rf ../Re50/system system
# cp -rf ../Re50/run_mesh run_mesh
# cp -rf ../Re50/run_solver run_solver

# sed -i 's/simulation_type.*/simulation_type        laminar;/g' ./0.orig/include/initialConditions
# sed -i 's/flowModel.*/flowModel        laminar;/g' ./0.orig/include/initialConditions
# sed -i 's/flowVelocity.*/flowVelocity               (9.8937e-03 0.0 0.0);/g' ./0.orig/include/initialConditions
# sed -i 's/flowUbar.*/flowUbar          (9.8937e-03 0.0 0.0);/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentKE.*/turbulentKE                 2.38e-08;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentOmega.*/turbulentOmega           6.97e-06;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentEpsilon.*/turbulentEpsilon       4.72e-14;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentNut.*/turbulentNut               4.29e+01;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentNuTilda.*/turbulentNuTilda       6.72e-05;/g' ./0.orig/include/initialConditions
# cp -rf 0.orig 0
# ./run_solver
# cd ..

# mkdir Re200
# cd Re200
# touch Re200.foam

# cp -rf ../Re50/0.orig 0.orig
# cp -rf ../Re50/constant constant
# cp -rf ../Re50/system system
# cp -rf ../Re50/run_mesh run_mesh
# cp -rf ../Re50/run_solver run_solver

# sed -i 's/simulation_type.*/simulation_type        laminar;/g' ./0.orig/include/initialConditions
# sed -i 's/flowModel.*/flowModel        laminar;/g' ./0.orig/include/initialConditions
# sed -i 's/flowVelocity.*/flowVelocity               (1.3192e-02 0.0 0.0);/g' ./0.orig/include/initialConditions
# sed -i 's/flowUbar.*/flowUbar          (1.3192e-02 0.0 0.0);/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentKE.*/turbulentKE                 9.53e-08;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentOmega.*/turbulentOmega           1.39e-05;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentEpsilon.*/turbulentEpsilon       3.78e-13;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentNut.*/turbulentNut               4.29e+01;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentNuTilda.*/turbulentNuTilda       6.72e-05;/g' ./0.orig/include/initialConditions
# cp -rf 0.orig 0
# ./run_solver
# cd ..

# mkdir Re300
# cd Re300
# touch Re300.foam

# cp -rf ../Re50/0.orig 0.orig
# cp -rf ../Re50/constant constant
# cp -rf ../Re50/system system
# cp -rf ../Re50/run_mesh run_mesh
# cp -rf ../Re50/run_solver run_solver

# sed -i 's/simulation_type.*/simulation_type        laminar;/g' ./0.orig/include/initialConditions
# sed -i 's/flowModel.*/flowModel        laminar;/g' ./0.orig/include/initialConditions
# sed -i 's/flowVelocity.*/flowVelocity               (1.9787e-02 0.0 0.0);/g' ./0.orig/include/initialConditions
# sed -i 's/flowUbar.*/flowUbar          (1.9787e-02 0.0 0.0);/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentKE.*/turbulentKE                 9.53e-08;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentOmega.*/turbulentOmega           1.39e-05;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentEpsilon.*/turbulentEpsilon       3.78e-13;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentNut.*/turbulentNut               4.29e+01;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentNuTilda.*/turbulentNuTilda       6.72e-05;/g' ./0.orig/include/initialConditions
# cp -rf 0.orig 0
# ./run_solver
# cd ..

# mkdir Re400
# cd Re400
# touch Re400.foam

# cp -rf ../Re50/0.orig 0.orig
# cp -rf ../Re50/constant constant
# cp -rf ../Re50/system system
# cp -rf ../Re50/run_mesh run_mesh
# cp -rf ../Re50/run_solver run_solver

# sed -i 's/simulation_type.*/simulation_type        laminar;/g' ./0.orig/include/initialConditions
# sed -i 's/flowModel.*/flowModel        laminar;/g' ./0.orig/include/initialConditions
# sed -i 's/flowVelocity.*/flowVelocity               (2.6383e-02 0.0 0.0);/g' ./0.orig/include/initialConditions
# sed -i 's/flowUbar.*/flowUbar          (2.6383e-02 0.0 0.0);/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentKE.*/turbulentKE                 2.65e-07;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentOmega.*/turbulentOmega           2.32e-05;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentEpsilon.*/turbulentEpsilon       1.75e-12;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentNut.*/turbulentNut               4.29e+01;/g' ./0.orig/include/initialConditions
# sed -i 's/turbulentNuTilda.*/turbulentNuTilda       6.72e-05;/g' ./0.orig/include/initialConditions
# cp -rf 0.orig 0
# ./run_solver
# cd ..
