enddo
call flush(iout)
endif
- write (2,*) "Start reading THETA_PDB"
+cd write (2,*) "Start reading THETA_PDB"
do i=1,ntyp
read (ithep_pdb,*,err=111,end=111) a0thet(i),(athet(j,i),j=1,2),
& (bthet(j,i),j=1,2)
read (ithep_pdb,*,err=111,end=111) theta0(i),sig0(i),sigc0(i)
sigc0(i)=sigc0(i)**2
enddo
- write (2,*) "End reading THETA_PDB"
+cd write (2,*) "End reading THETA_PDB"
close (ithep_pdb)
#endif
close(ithep)