diff --git a/mechanical_model/run_mechanical_model_helvetios.sh b/mechanical_model/run_mechanical_model_helvetios.sh
new file mode 100644
index 0000000..c966cdb
--- /dev/null
+++ b/mechanical_model/run_mechanical_model_helvetios.sh
@@ -0,0 +1,52 @@
+#!/bin/bash
+
+#SBATCH --partition debug
+#SBATCH --ntasks 1
+#SBATCH --cpus-per-task 36
+#SBATCH --mem-per-cpu 2G
+#SBATCH --time 00:59:00
+#SBATCH --error=slurm-%j.stderr
+#SBATCH --output=slurm-%j.stdout
+#SBATCH --job-name=mech_model
+#SBATCH --exclusive
+
+source /home/kubilay/anaconda3/bin/activate
+conda activate muspectre_env2
+export LD_LIBRARY_PATH=/home/kubilay/pnetcdf/lib/:$LD_LIBRARY_PATH
+
+#module load gcc mvapich2
+
+echo "*** STARTING JOB ***"
+
+LOG_FILE="mechanical_model.log"
+
+echo "*** STARTING JOB ***" > $LOG_FILE
+
+which python3 > $LOG_FILE
+which mpirun >> $LOG_FILE
+which srun >> $LOG_FILE
+ulimit -a >> $LOG_FILE
+START_LAYER=60
+LAYERS=61
+
+echo $OPENBLAS_NUM_THREAD >> $LOG_FILE 
+export OPENBLAS_NUM_THREADS=1
+echo $OPENBLAS_NUM_THREADS >> $LOG_FILE
+
+
+for z in $(seq $START_LAYER $LAYERS)
+do
+    echo $z
+    conda activate muspectre_env2
+	#valgrind --leak-check=full -s mpirun -np 16 python3 residual_stress_analysis.py $z 15 >> $LOG_FILE
+	mpirun -n 36 python3 residual_stress_analysis.py $z 35 >>$LOG_FILE
+	#mpirun -np 5 amplxe-cl -collect hotspots -r /scratch/kubilay/additive-manufacturing-work/mechanical_model/profiling -- python3 residual_stress_analysis.py 60 4 >> $LOG_FILE
+	#python3 netcdf2h5.py $z >> $LOG_FILE
+    conda deactivate
+done
+
+conda deactivate
+
+echo "*** JOB FINISHED ***"
+
+echo "*** JOB FINISHED ***" >> $LOG_FILE