R11910/run_model.sh
ceca456ea790cluster_updates

run_model.sh
View Options

			EXTRUDE=true
			PART=false
			GEOM_FILE="half_buchbinder_points.txt"
			THERMAL=true
			LAYERS=333
			START_LAYER=272
			LAYER_THICKNESS=0.03
			FREECAD_DIRECTORY="/home/ekinkubilay/anaconda3/envs/freecad_env/lib"

			source /home/ekinkubilay/anaconda3/bin/activate

			if $PART
			then
			echo "START PART GEOMETRY"
			conda activate freecad_env
			#go to Part_geometry directory
			cd Part_geometry/
			LOG_FILE="part_geometry.log"
			#extrude part using the python script if necessary
			if $EXTRUDE
			then
			python3 create_part_geometry.py $FREECAD_DIRECTORY $GEOM_FILE 5 ${GEOM_FILE%.txt}.step > $LOG_FILE
			cd ..
			conda deactivate
			echo "END PART GEOMETRY"
			fi
			fi



			if $THERMAL
			then
			for z in $(seq $START_LAYER $LAYERS)
			do
			echo "LAYER $z"
			cd Part_geometry/
			conda activate freecad_env
			LOG_FILE="part_geometry.log"
			#layer and mesh the part
			echo "START LAYER GEOMETRY"
			python3 layer_part.py $FREECAD_DIRECTORY ${GEOM_FILE%.txt}.step/${GEOM_FILE%.txt}.step $LAYER_THICKNESS $z layer mesh >> $LOG_FILE
			#convert .msh files into .xml files
			conda activate base
			sh convert_mesh.sh >> $LOG_FILE
			echo "END LAYER GEOMETRY"
			cd ..
			echo "START THERMAL MODEL"
			#go the Thermal_model directory
			cd Thermal_model/
			LOG_FILE="thermal_model.log"
			#adjust previous temperature field in series
			python adjust_previous_temperature.py $z >>$LOG_FILE
			#run thermal problem
			mpirun -np 4 python fem_layer_scan_final.py $z 8000 500 19.2 175 0.8 0.00015 >>$LOG_FILE
			echo "END THERMAL MODEL"
			cd ..

			done
			fi