enddo
call flush(iout)
endif
- write (2,*) "Start reading THETA_PDB",ithep_pdb
+c write (2,*) "Start reading THETA_PDB",ithep_pdb
do i=1,ntyp
c write (2,*) 'i=',i
read (ithep_pdb,*,err=111,end=111)
gthet(j,i)=gthet(j,-i)
enddo
enddo
- write (2,*) "End reading THETA_PDB"
+c write (2,*) "End reading THETA_PDB"
close (ithep_pdb)
#endif
close(ithep)
C Read the parameters of the probability distribution/energy expression
C of the side chains.
C
- write (2,*) "Start reading ROTAM_PDB"
+c write (2,*) "Start reading ROTAM_PDB"
do i=1,ntyp
read (irotam_pdb,'(3x,i3,f8.3)',end=112,err=112) nlob(i),dsc(i)
if (i.eq.10) then
endif
enddo
close (irotam_pdb)
- write (2,*) "End reading ROTAM_PDB"
+c write (2,*) "End reading ROTAM_PDB"
#endif
close(irotam)
do i=-ntyp,-1
itortyp(i)=-itortyp(-i)
enddo
- write (iout,*) 'ntortyp',ntortyp
+c write (iout,*) 'ntortyp',ntortyp
do i=0,ntortyp-1
do j=-ntortyp+1,ntortyp-1
read (itorp,*,end=113,err=113) nterm(i,j,iblock),