C & rad2deg*theta_constr0(i), rad2deg*theta_drange(i),
C & rad2deg*difi,0.25d0*for_thet_constr(i)*difi**4,
C & gloc(itheta+nphi-2,icg)
- endif
+C endif
enddo
C Ufff.... We've done all this!!!
return
logical seq_comp
integer i,j,k,nres_pdb,iaux
double precision ddsc,dist
+ integer nnt_old,nct_old
integer ilen,kkk
external ilen
C
nres0=nres
+ nnt_old=nnt
+ nct_old=nct
write (iout,*) "pdbref",pdbref
if (pdbref) then
read(inp,'(a)') pdbfile
& 'Number of residues to be superposed:',nsup,
& ' (from residue',nstart_sup,' to residue',
& nend_sup,').'
+ nnt=nnt_old
+ nct=nct_old
return
end