Commit 823ac083 authored by Matthieu Tristram's avatar Matthieu Tristram
Browse files

Force hpjul3

Force remove_cosmo_limits=true
Correct job name
parent dc80de71
......@@ -77,7 +77,8 @@ fi
file=$(readlink -e $parFile)
parbase=$(basename $file)
dirout=${parbase%".par"}_prof_"$var"
parbase=${parbase%".par"}
dirout=${parbase}_prof_"$var"
if [ -d "$dirout" ] ; then
echo "Output directory $dirout already exists: still want to run? (y/n)"
......@@ -131,26 +132,31 @@ type=\$(awk '(\$1 =="par" || \$1 =="fix") && \$2 =="${var/\*/\\*}" {print \$3}'
#if not found assume cosmo type
if [ -z "\$type" ] ; then
type="cosmo"
type="cosmo"
fi
echo "fix $var \$type $zeval" >> parfile_in
#force high precision settings
grep -v precisionFile parfile_in | grep -v remove_cosmo_limits > tmp
echo "precisionFile=class_pre/hpjul3.pre" >> tmp
echo "remove_cosmo_limits=true" >> tmp
mv tmp parfile_in
#special case of degenerate neutrinos (need to decrement Neff by 1 when mass=0)
if [ "$var" == "m_ncdm" ] ; then
grep -v ncdm parfile_in | grep -v N_ur > tmp
if [ "$zeval" == "0" ] ; then
echo "class N_ur 3.046" >> tmp
echo "fix N_ncdm cosmo 0" >> tmp
else
echo "class N_ncdm 1" >> tmp
echo "class N_ur 2.0328" >> tmp
echo "fix m_ncdm cosmo $zeval" >> tmp
fi
cp -f tmp parfile_in
grep -v ncdm parfile_in | grep -v N_ur > tmp
if [ "$zeval" == "0" ] ; then
echo "class N_ur 3.046" >> tmp
echo "fix N_ncdm cosmo 0" >> tmp
else
echo "class N_ncdm 1" >> tmp
echo "class N_ur 2.0328" >> tmp
echo "fix m_ncdm cosmo $zeval" >> tmp
fi
cp -f tmp parfile_in
fi
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment