include 'COMMON.SBRIDGE'
include 'COMMON.DISTFIT'
include 'COMMON.MD'
+ include 'COMMON.REMD'
include 'COMMON.SETUP'
integer itime
double precision energia(0:n_ene)
& rms,frac,frac_nn,kinetic_T,t_bath,gyrate(),
& potEcomp(23),me
format1="a133"
- print *,'A CHUJ',potEcomp(23)
- else
+ else
+C print *,'A CHUJ',potEcomp(23)
write (line1,'(i10,f15.2,7f12.3,i5,$)')
& itime,totT,EK,potE,totE,
& kinetic_T,t_bath,gyrate(),
& potEcomp(23),me
format1="a114"
endif
+ else if (selfguide.gt.0) then
+ distance=0.0
+ do j=1,3
+ distance=distance+(c(j,afmend)-c(j,afmbeg))**2
+ enddo
+ distance=dsqrt(distance)
+ if (refstr) then
+ call rms_nac_nnc(rms,frac,frac_nn,co,.false.)
+ write (line1,'(i10,f15.2,3f12.3,f7.2,2f6.3,f12.3,f10.1,2f8.2,
+ & f9.2,i5,$)')
+ & itime,totT,EK,potE,totE,
+ & rms,frac,frac_nn,kinetic_T,t_bath,gyrate(),
+ & distance,potEcomp(23),me
+ format1="a133"
+C print *,"CHUJOWO"
+ else
+C print *,'A CHUJ',potEcomp(23)
+ write (line1,'(i10,f15.2,8f12.3,i5,$)')
+ & itime,totT,EK,potE,totE,
+ & kinetic_T,t_bath,gyrate(),
+ & distance,potEcomp(23),me
+ format1="a114"
+ endif
else
if (refstr) then
call rms_nac_nnc(rms,frac,frac_nn,co,.false.)
line2=' '
endif
if (print_compon) then
+ if(itime.eq.0) then
+ write(format,'(a1,a4,a1,a4,a10)') "(",format1,",",format2,
+ & ",20a12)"
+ write (istat,format) "#","",
+ & (ename(print_order(i)),i=1,nprint_ene)
+ endif
write(format,'(a1,a4,a1,a4,a10)') "(",format1,",",format2,
& ",20f12.3)"
write (istat,format) line1,line2,