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)
do l=1,maxinter
do i=1,nsccortyp
do j=1,nsccortyp
- read (isccor,*,end=1113,err=1113) nterm_sccor(i,j),nlor_sccor(i,j)
+ read (isccor,*,end=1113,err=1113) nterm_sccor(i,j),
+ & nlor_sccor(i,j)
v0ijsccor=0.0d0
si=-1.0d0