C
C Calculate the virtual-bond-angle energy.
C
- call ebend(ebe)
+ call ebend(ebe,ethetacnstr)
C
C Calculate the SC local energy.
C
energia(17)=estr
energia(19)=edihcnstr
energia(21)=esccor
+ energia(24)=ethetacnstr
c write (iout,*) "ETOTAL_SHORT before SUM_ENERGY"
call flush(iout)
call sum_energy(energia,.true.)