#!/bin/bash


qnumin=0
qnumax=1

v0qjmax=11
v1qjmax=7

qjstart=0

state=1
col=2
rm E0_2.dat
for (( qnu=$qnumin; qnu <= $qnumax; qnu++ ))
do
   #maximum J state selector -> qjmax is different per vibrational state
   if [ $qnu -eq 0 ]
   then
     qjmax=$v0qjmax
   fi
   if [ $qnu -eq 1 ]
   then
     qjmax=$v1qjmax
   fi
   let qj=$qjstart               #qj is rotational quantum number
   while [ $qj -le $qjmax ]
   do
     #satval=`tail -1 ALL.degeneracy.averaged | awk '{print $'$col'}'`


     NR=3; RP=0
     while (( $( echo "$RP <= ( 0.325 * 0.5 )" | bc -l ) ))
     do
       E=` awk ' NR == '$NR' {print $1}' ALL.degeneracy.averaged `
       RP=`awk ' NR == '$NR' {print $'$col'}' ALL.degeneracy.averaged `
       RP=`sed -E 's/([+-]?[0-9.]+)[eE]\+?(-?)([0-9]+)/(\1*10^\2\3)/g' <<<"$RP"`
       ((NR++))
       #echo $E $RP $NR $col
     done
     
     echo "$qj  $E  	$satval  $expE"
     echo "$qj  $E  	$satval  $expE" >> E0_2.dat
#     echo "$qj  $E" >> satval_${satval}
     ((col++))
     ((qj++))
     ((state++))
   done #qj
   echo 
  
done #qnu
echo
awk '{print $1, $2}' E0_2.dat