& duconst(3,0:MAXRES),duxconst(3,0:MAXRES)
integer kstart,kend,lstart,lend,idummy
double precision delta /1.0d-7/
& duconst(3,0:MAXRES),duxconst(3,0:MAXRES)
integer kstart,kend,lstart,lend,idummy
double precision delta /1.0d-7/
write(iout,'(i5,3e15.5)') ii,(dxqwol(j,ii),j=1,3)
enddo
#endif
write(iout,'(i5,3e15.5)') ii,(dxqwol(j,ii),j=1,3)
enddo
#endif
-c Calculating numerical gradients of dU/dQi and dQi/dxi
-c call qwol_num(ifrag(1,i,iset),ifrag(2,i,iset),.true.
-c & ,idummy,idummy)
c do ii=1,nres
c write(iout,'(i5,3e15.5)') ii,(dxqwol(j,ii),j=1,3)
c enddo
c do ii=1,nres
c write(iout,'(i5,3e15.5)') ii,(dxqwol(j,ii),j=1,3)
c enddo
write(iout,'(i5,3e15.5)') ii,(duxconst(j,ii),j=1,3)
enddo
#endif
write(iout,'(i5,3e15.5)') ii,(duxconst(j,ii),j=1,3)
enddo
#endif