mirror of
https://asciireactor.com/otho/cloudy-agn.git
synced 2024-12-05 01:55:08 +00:00
116 lines
3.4 KiB
Bash
Executable File
116 lines
3.4 KiB
Bash
Executable File
#!/usr/bin/env bash
|
|
|
|
# This still isn't really generalized. Need to go to a
|
|
# specific directory of runs to run this. Feed it a list
|
|
# of column densities.
|
|
|
|
script_dir=$(
|
|
cd $(dirname $0) ;
|
|
pwd -P |sed 's@^\(.*\)/scripts.*@\1/scripts@'
|
|
)
|
|
base_dir=`pwd`
|
|
bin_dir="$script_dir/../bin"
|
|
|
|
griddir=$1
|
|
grid_id=$(echo ${griddir}|sed 's@^\.\/\(.*\)@\1@'|sed 's@.grids@@'|sed 's@\/@.@g')
|
|
|
|
if [[ -e "fortfiles_${grid_id}.tar.gz" ]]; then
|
|
echo "Package fortfiles_${grid_id}.tar.gz exists."
|
|
exit 5
|
|
else
|
|
echo "Setting up to compile fortfiles_${grid_id}.tar.gz."
|
|
fi
|
|
|
|
cd $griddir
|
|
|
|
if [[ -e "mpi_grid.out" ]]; then
|
|
echo -n ""
|
|
else
|
|
echo "mpi_grid.out doesn't exist"
|
|
cd $base_dir
|
|
exit 6
|
|
fi
|
|
|
|
mkdir -p fortfiles
|
|
cd fortfiles
|
|
# pwd
|
|
|
|
cp ../mpi_grid.in .
|
|
cp ../*.tab .
|
|
|
|
echo "Directory ready. Calling fort file creation."
|
|
$bin_dir/create_fort_files ../mpi_grid.out $script_dir/../reference/linelist.c17
|
|
|
|
# inwTfortnums=$(grep 'InwT' fort.*|sed 's@^fort\.\([0-9]\+\):InwT.*@\1@'|sort -nr)
|
|
# echo "Computing inward diffuse tables for $(echo "$inwTfortnums"|wc -l) coordinates."
|
|
# echo -n "Writing InwD tables to fortfiles and shifting files:"
|
|
# count=0
|
|
# for inwTnum in $(echo $inwTfortnums|sort -r)
|
|
# do
|
|
# inwCnum=$(($inwTnum + 1))
|
|
# inwDnum=$(($inwCnum + 1))
|
|
# allfortnums=$(ls fort.*|sed 's@fort\.\([0-9]\+\)@\1@g'|sort -n)
|
|
# fortnumstoproc=""
|
|
# for fortnum in $allfortnums
|
|
# do
|
|
# if (( $fortnum > $inwCnum ))
|
|
# then
|
|
# fortnumstoproc="$fortnum $fortnumstoproc"
|
|
# fi
|
|
# done
|
|
# # echo "Shifting $( echo "$fortnumstoproc"|wc -w ) files."
|
|
# for fortnum in $fortnumstoproc
|
|
# do
|
|
# newnum=$(( $fortnum + 1 ))
|
|
# mv fort.$fortnum fort.$newnum
|
|
# done
|
|
# #echo -n " $inwDnum"
|
|
# $bin_dir/subtract_fortfiles fort.$inwTnum fort.$inwCnum InwD > fort.$inwDnum
|
|
# echo -ne "\033[2K" ; printf "\r"
|
|
# echo -n "Writing InwD tables to fortfiles and shifting files: $(( ++count ))"
|
|
# done
|
|
# echo ""
|
|
#
|
|
# inwCfortnums=$(grep 'InwC' fort.*|sed 's@^fort\.\([0-9]\+\):InwC.*@\1@'|sort -nr)
|
|
# echo "Computing total diffuse tables for $(echo "$inwCfortnums"|wc -l) coordinates."
|
|
# echo -n "Writing nFnD tables to fortfiles and shifting files: 0"
|
|
# count=0
|
|
# for inwCnum in $(echo $inwCfortnums|sort -r)
|
|
# do
|
|
# nFnunum=$(($inwCnum - 3))
|
|
# nFnDnum=$(($inwCnum + 2))
|
|
# allfortnums=$(ls fort.*|sed 's@fort\.\([0-9]\+\)@\1@g'|sort -n)
|
|
# fortnumstoproc=""
|
|
# for fortnum in $allfortnums
|
|
# do
|
|
# if (( $fortnum >= $nFnDnum ))
|
|
# then
|
|
# fortnumstoproc="$fortnum $fortnumstoproc"
|
|
# fi
|
|
# done
|
|
# # echo "Shifting $( echo "$fortnumstoproc"|wc -w ) files."
|
|
# for fortnum in $fortnumstoproc
|
|
# do
|
|
# newnum=$(( $fortnum + 1 ))
|
|
# mv fort.$fortnum fort.$newnum
|
|
# done
|
|
#
|
|
# #echo -n " $nFnDnum"
|
|
# $bin_dir/subtract_fortfiles fort.$nFnunum fort.$inwCnum nFnD > fort.$nFnDnum
|
|
# echo -ne "\033[2K" ; printf "\r"
|
|
# echo -n "Writing nFnD tables to fortfiles and shifting files: $(( ++count ))"
|
|
# done
|
|
# echo ""
|
|
|
|
echo "Applying interpolative smoothing."
|
|
$script_dir/operations/bulk_interpolation_fix.sh fort.* > interpolation_report
|
|
|
|
|
|
grid_id=$(echo ${griddir}|sed 's@^\.\/\(.*\)@\1@'|sed 's@.grids@@'|sed 's@\/@.@g')
|
|
echo "Saving $grid_id grid"
|
|
tar cf fortfiles_${grid_id}.tar *
|
|
gzip fortfiles_${grid_id}.tar
|
|
mv fortfiles_${grid_id}.tar.gz $base_dir
|
|
|
|
cd $base_dir
|