external rescode
logical iblnk
external iblnk
+ double precision cm(3)
ifreq=1
is=1
c(j,i+nres)=coord(j,ii+nres)
enddo
enddo
+c Calculate the CM
+ do j=1,3
+ cm(j)=0.0d0
+ do i=1,nres
+ cm(j)=cm(j)+c(j,i)
+ enddo
+ cm(j)=cm(j)/nres
+ enddo
+ do i=1,nres
+ do j=1,3
+ c(j,i)=c(j,i)-cm(j)
+ c(j,i+nres)=c(j,i+nres)-cm(j)
+ enddo
+ enddo
etot=potE
write (tytul,'(a,i6)') "Structure",kk
call pdbout(etot,tytul,9)