From 34f997cd33e08be1ac8b2db8e60c9097643836d2 Mon Sep 17 00:00:00 2001
From: rstiskalek <richard.stiskalek@protonmail.com>
Date: Sat, 5 Oct 2024 00:54:10 +0100
Subject: [PATCH] Add Quijote SPH

---
 scripts_independent/field_sph_gadget.sh | 41 +++++++++++++++++++------
 1 file changed, 32 insertions(+), 9 deletions(-)

diff --git a/scripts_independent/field_sph_gadget.sh b/scripts_independent/field_sph_gadget.sh
index ca14161..4839084 100755
--- a/scripts_independent/field_sph_gadget.sh
+++ b/scripts_independent/field_sph_gadget.sh
@@ -2,12 +2,12 @@
 nthreads=12
 memory=7
 on_login=${1}
-queue="berg"
+queue="cmb"
 env="/mnt/zfsusers/rstiskalek/csiborgtools/venv_csiborg/bin/python"
 file="field_sph_gadget.py"
 
 # Guilhem says higher resolution is better
-resolution=1024
+resolution=256
 SPH_executable="/mnt/users/rstiskalek/cosmotool/bld/sample/simple3DFilter"
 scratch_space="/mnt/extraspace/rstiskalek/dump/"
 
@@ -24,8 +24,8 @@ if [ -z "$on_login" ] || ! [[ "$on_login" =~ ^[0-1]$ ]]; then
     exit 1
 fi
 
-export OMP_NUM_THREADS={nthreads}
-export OMP_NESTED=true
+export OMP_NUM_THREADS=$nthreads
+export OMP_MAX_ACTIVE_LEVELS=4
 
 # pythoncm="$env $file --snapshot_path $snapshot_path --output_path $output_path --resolution $resolution --scratch_space $scratch_space --SPH_executable $SPH_executable --snapshot_kind $snapshot_kind"
 # if [ $on_login -eq 1 ]; then
@@ -40,11 +40,34 @@ export OMP_NESTED=true
 # fi
 
 
-# Manticore SWIFT submission loop
-snapshot_kind="swift"
-for k in {0..40}; do
-    snapshot_path="/mnt/extraspace/rstiskalek/MANTICORE/2MPP_N128_DES_V1/resimulations/R512/mcmc_$k/swift_monofonic/snap_0001/snap_0001.hdf5"
-    output_path="/mnt/extraspace/rstiskalek/MANTICORE/2MPP_N128_DES_V1/fields/R512/SPH_$k.hdf5"
+# ####### Manticore SWIFT submission loop #######
+# snapshot_kind="swift"
+# for k in {0..40}; do
+#     snapshot_path="/mnt/extraspace/rstiskalek/MANTICORE/2MPP_N128_DES_V1/resimulations/R512/mcmc_$k/swift_monofonic/snap_0001/snap_0001.hdf5"
+#     output_path="/mnt/extraspace/rstiskalek/MANTICORE/2MPP_N128_DES_V1/fields/R512/SPH_$k.hdf5"
+#
+#     pythoncm="$env $file --snapshot_path $snapshot_path --output_path $output_path --resolution $resolution --scratch_space $scratch_space --SPH_executable $SPH_executable --snapshot_kind $snapshot_kind"
+#     if [ $on_login -eq 1 ]; then
+#         echo $pythoncm
+#         $pythoncm
+#     else
+#         cm="addqueue -s -q $queue -n 1x$nthreads -m $memory $pythoncm"
+#         echo "Submitting:"
+#         echo $cm
+#         echo
+#         eval $cm
+#     fi
+#
+#     sleep 0.05
+# done
+
+####### Quijote submission loop #######
+snapshot_kind="gadget2"
+
+for k in {1..50}; do
+# for k in 0; do
+    snapshot_path="/mnt/extraspace/rstiskalek/quijote/Snapshots_fiducial/${k}/snapdir_004/snap_004"
+    output_path="/mnt/extraspace/rstiskalek/quijote/SPHField_fiducial/field_${k}_004.hdf5"
 
     pythoncm="$env $file --snapshot_path $snapshot_path --output_path $output_path --resolution $resolution --scratch_space $scratch_space --SPH_executable $SPH_executable --snapshot_kind $snapshot_kind"
     if [ $on_login -eq 1 ]; then