integer :: k,il,il1,i,j,nharp,ii,ierr,itime_master,irr,iex,&
i_set_temp,itmp,i_temp,i_mult,i_iset,i_mset,i_dir,i_temp1,&
i_mult1,i_iset1,i_mset1,ierror,itime,mnum
- integer,dimension(4,nres/3) :: iharp !(4,nres/3)(4,maxres/3)
+ integer,dimension(4,nres) :: iharp !(4,nres/3)(4,maxres/3)
!deb imin_itime_old=0
integer :: nres2 !el
WRITE(iout,*) "JUST AFTER CALL"
do i=nnt,nct-1
mnum=(molnum(i))
stdforcp(i)=stdfp(mnum)*dsqrt(gamp(mnum))
+! write(iout,*) "stdforcp=",stdforcp(i),itype(i,mnum),i
enddo
do i=nnt,nct
mnum=molnum(i)
if (itype(i,mnum).ne.ntyp1) stdforcsc(i)=stdfsc(iabs(itype(i,mnum)),mnum)&
*dsqrt(gamsc(iabs(itype(i,mnum)),mnum))
+! write(iout,*) "stdforcsc=",stdforcsc(i),itype(i,mnum),i
enddo
endif
call rescale_weights(t_bath)