esccor=0.0D0
do i=itau_start,itau_end
esccor_ii=0.0D0
-CC TU NIE MOZE BYC IABS DO ZMIANY JAK BEDA GOTOWE POTENCJALY AM1
- isccori=isccortyp(iabs(itype(i-2)))
- isccori1=isccortyp(iabs(itype(i-1)))
+ isccori=isccortyp((itype(i-2)))
+ isccori1=isccortyp((itype(i-1)))
phii=phi(i)
cccc Added 9 May 2012
cc Tauangle is torsional engle depending on the value of first digit