do not store sub-maps

recursive_gmconvert.sh

@@ -50,6 +50,11 @@ then
 	exit 1
 fi
 
+if [[ -z "$n_gaussians0" ]]
+then
+	n_gaussians0=$n_gaussians
+fi
+
 MaxArraySize=100000
 get_max_jobs(){
 	echo ${MaxArraySize}
@@ -70,6 +75,7 @@ cd $map_root_dir
 mkdir -p converged
 for ((i=$i0 ; i<=N ; i++))
 do
+	echo "i=$i"
 	n=$((i))
 	jobname=${map_name}_${n}
 
@@ -91,8 +97,8 @@ do
 			n_jobs=$((n_jobs + $(awk '/NGAUSS/{s+=$3}END{print s}' ${gmm_file}) ))
 		done
 		echo echo ${gmm_name} '$((SLURM_ARRAY_TASK_ID-1+offset))' \> 'tmp_$((SLURM_ARRAY_TASK_ID+offset)).txt'  
-		echo ${profile_cmd} ${gmconvert} -imap ${map_name}  -gmml 'tmp_$((SLURM_ARRAY_TASK_ID+offset)).txt' -oimap $n/'sub_$((SLURM_ARRAY_TASK_ID+offset-1)).mrc' -ogmm $n/${n}_'$((SLURM_ARRAY_TASK_ID-1+offset))'.gmm -ng ${n_gaussians} -zth ${threshold}
-		echo ${gmconvert}  -igmm $n/${n}_'$((SLURM_ARRAY_TASK_ID-1+offset)).gmm' -imap $n/'sub_$((SLURM_ARRAY_TASK_ID+offset-1)).mrc' -omap /dev/null
+		echo ${profile_cmd} ${gmconvert} -imap ${map_name}  -gmml 'tmp_$((SLURM_ARRAY_TASK_ID+offset)).txt' -ogmm $n/${n}_'$((SLURM_ARRAY_TASK_ID-1+offset))'.gmm -ng ${n_gaussians} -zth ${threshold}
+		echo ${gmconvert}  -igmm $n/${n}_'$((SLURM_ARRAY_TASK_ID-1+offset)).gmm' -imap ${map_name}  -gmml 'tmp_$((SLURM_ARRAY_TASK_ID+offset)).txt' -omap /dev/null
 	fi > commands_$i.sh
 	mkdir -p $n
 	max_jobs=$(get_max_jobs)