RoVib="v1j7mj7.A.0.2eV-0.75eV.out"
starttable=`awk '/INFO FROM SETTING UP THE ASYMPTOTIC ROVIBRATIONAL STATES/ {print NR}' ../out/$RoVib` 
endtable=`awk '/INFO FROM SETTING UP THE DIFFRACTION BASIS/ {print NR}' ../out/$RoVib`
starttable=` echo "$starttable + 7" | bc `
endtable=` echo "$endtable -13" | bc `
#echo "$starttable   $endtable"

vmax=`awk ' NR=='$endtable' {print $1}' ../out/$RoVib`
jmax=`awk ' NR=='$endtable' {print $2}' ../out/$RoVib`
totstates=`echo "($vmax + 1) * $jmax " | bc `

echo "#RoVib levels from $RoVib"                                >> RoVib.levels
echo "$vmax          #nu_max"                                   >> RoVib.levels
echo "$jmax          #j_max"                                    >> RoVib.levels
echo "$totstates     #total states "                            >> RoVib.levels
while [ "$starttable" -le "$endtable" ];
do
  awk ' NR=='$starttable' {print $1, $2, $3}' ../out/$RoVib     >> RoVib.levels
  ((starttable++))
  ((starttable++))
done