mirror of
https://asciireactor.com/otho/cloudy-agn.git
synced 2024-12-05 01:45:08 +00:00
25 lines
989 B
Bash
Executable File
25 lines
989 B
Bash
Executable File
#!/bin/bash
|
|
|
|
reference_file=`echo ${2}[0-9]*_${3}[0-9]*`
|
|
line_number=`grep -n "$1" ${reference_file}|tail -n1|sed 's#\([0-9]\): .*#\1#'`
|
|
avg_wavelength=`sed -n ${line_number}p ${reference_file}|sed 's#^......\(.......\).*$#\1#'`
|
|
num_potential_components=14
|
|
element=""
|
|
for line_num in `seq $num_potential_components`; do
|
|
line=`sed -n $((line_number + $line_num))p ${reference_file}`
|
|
# echo "$line"
|
|
if [[ `echo $line|grep Inwd` != "" ]]; then continue; fi
|
|
if [[ `echo $line|grep TOTL` != "" ]]; then break; fi
|
|
element_in_line=`echo "$line"|sed 's#^.\(....\).*$#\1#'`
|
|
if [[ ${element} == "" ]]; then element=${element_in_line}; fi
|
|
if [[ ${element} != ${element_in_line} ]]; then continue; fi
|
|
component_wavelength=`echo "$line"|sed 's#^......\(.......\).*$#\1#'`
|
|
# echo $avg_wavelength $component_wavelength
|
|
if [[ `bc <<< "scale=5; sqrt( ( ( ${avg_wavelength} - ${component_wavelength} ) / ${avg_wavelength} )^2 ) < .0162"` == 1 ]]; then
|
|
echo "$line"
|
|
else
|
|
break
|
|
fi
|
|
done
|
|
exit 0
|