energy_dec=(index(controlcard,'ENERGY_DEC').gt.0)
unres_pdb = index(controlcard,'UNRES_PDB') .gt. 0
call readi(controlcard,'TORMODE',tor_mode,0)
energy_dec=(index(controlcard,'ENERGY_DEC').gt.0)
unres_pdb = index(controlcard,'UNRES_PDB') .gt. 0
call readi(controlcard,'TORMODE',tor_mode,0)