wstrain=ww(15)
wbond=ww(18)
wsccor=ww(19)
+ wdfa_dist=ww(23)
+ wdfa_tor=ww(24)
+ wdfa_nei=ww(25)
+ wdfa_beta=ww(26)
+ write (iout,*) "wdfa_dist",wdfa_dist," wdfa_tor",wdfa_tor,
+ & " wdfa_nei",wdfa_nei," wdfa_beta",wdfa_beta
endif
cc write(iout,*) "KURWA", wstrain,akcm,akth,wsc,dyn_ss
if (lprint) then
write (iout,'(/a/)') 'Torsional constants:'
+ do l=1,maxinter
do i=1,nsccortyp
do j=1,nsccortyp
write (iout,*) 'ityp',i,' jtyp',j
enddo
enddo
enddo
+ enddo
endif
C