#ifdef DEBUG
write (iout,*) "Opening file ",intinname(:ilen(intinname))
write (iout,*) "lenrec",lenrec_in
call flush(iout)
#endif
#ifdef DEBUG
write (iout,*) "Opening file ",intinname(:ilen(intinname))
write (iout,*) "lenrec",lenrec_in
call flush(iout)
#endif
#if (defined(AIX) && !defined(JUBL))
call xdrf3dfcoord_(ixdrf, csingle, itmp, prec, iret)
if (iret.eq.0) goto 101
call xdrfint_(ixdrf, nss, iret)
if (iret.eq.0) goto 101
do j=1,nss
#if (defined(AIX) && !defined(JUBL))
call xdrf3dfcoord_(ixdrf, csingle, itmp, prec, iret)
if (iret.eq.0) goto 101
call xdrfint_(ixdrf, nss, iret)
if (iret.eq.0) goto 101
do j=1,nss
call xdrfint_(ixdrf, ihpb(j), iret)
if (iret.eq.0) goto 101
call xdrfint_(ixdrf, jhpb(j), iret)
if (iret.eq.0) goto 101
call xdrfint_(ixdrf, ihpb(j), iret)
if (iret.eq.0) goto 101
call xdrfint_(ixdrf, jhpb(j), iret)
if (iret.eq.0) goto 101
call xdrfint(ixdrf, ihpb(k), iret)
if (iret.eq.0) goto 101
call xdrfint(ixdrf, jhpb(k), iret)
if (iret.eq.0) goto 101
call xdrfint(ixdrf, ihpb(k), iret)
if (iret.eq.0) goto 101
call xdrfint(ixdrf, jhpb(k), iret)
if (iret.eq.0) goto 101
#ifdef DEBUG
write (iout,'(5hREAD ,i5,3f15.4,i10)')
& jj+1,energy(jj+1),entfac(jj+1),
#ifdef DEBUG
write (iout,'(5hREAD ,i5,3f15.4,i10)')
& jj+1,energy(jj+1),entfac(jj+1),
write (iout,'(8f10.5)') ((c(j,i+nres),j=1,3),i=nnt,nct)
call flush(iout)
#endif
write (iout,'(8f10.5)') ((c(j,i+nres),j=1,3),i=nnt,nct)
call flush(iout)
#endif
call add_new_cconf(jjj,jj,jj_old,icount,Next)
enddo
101 continue
call add_new_cconf(jjj,jj,jj_old,icount,Next)
enddo
101 continue
if (j.gt.2) then
if (itel(j).ne.0 .and. itel(j-1).ne.0) then
write (iout,*) "Conformation",jjj,jj+1
if (j.gt.2) then
if (itel(j).ne.0 .and. itel(j-1).ne.0) then
write (iout,*) "Conformation",jjj,jj+1
write (iout,*) "Conformation",jjj,jj+1
write (iout,*) "Bad CA-SC bond length",j," ",vbld(nres+j)
write (iout,*) "The Cartesian geometry is:"
write (iout,*) "Conformation",jjj,jj+1
write (iout,*) "Bad CA-SC bond length",j," ",vbld(nres+j)
write (iout,*) "The Cartesian geometry is:"
#ifdef DEBUG
write (iout,*) "DAREAD_COORDS"
write (iout,*) "istart_conf",istart_conf," iend_conf",iend_conf
#ifdef DEBUG
write (iout,*) "DAREAD_COORDS"
write (iout,*) "istart_conf",istart_conf," iend_conf",iend_conf
write (iout,*) "Reading binary file, record",iii," ii",ii
call flush(iout)
#endif
write (iout,*) "Reading binary file, record",iii," ii",ii
call flush(iout)
#endif
+ if (dyn_ss) then
+ read(icbase,rec=iii) ((allcart(j,i,ij),j=1,3),i=1,nres),
+ & ((allcart(j,i,ij),j=1,3),i=nnt+nres,nct+nres),
+c & nss_all(ij),(ihpb_all(i,ij),jhpb_all(i,ij),i=1,nss),
+ & entfac(ii),rmstb(ii)
+ else
read(icbase,rec=iii) ((allcart(j,i,ij),j=1,3),i=1,nres),
& ((allcart(j,i,ij),j=1,3),i=nnt+nres,nct+nres),
& nss_all(ij),(ihpb_all(i,ij),jhpb_all(i,ij),i=1,nss),
& entfac(ii),rmstb(ii)
read(icbase,rec=iii) ((allcart(j,i,ij),j=1,3),i=1,nres),
& ((allcart(j,i,ij),j=1,3),i=nnt+nres,nct+nres),
& nss_all(ij),(ihpb_all(i,ij),jhpb_all(i,ij),i=1,nss),
& entfac(ii),rmstb(ii)
#ifdef DEBUG
write (iout,*) ii,iii,ij,entfac(ii)
write (iout,'(8f10.5)') ((allcart(j,i,ij),j=1,3),i=1,nres)
#ifdef DEBUG
write (iout,*) ii,iii,ij,entfac(ii)
write (iout,'(8f10.5)') ((allcart(j,i,ij),j=1,3),i=1,nres)
& jhpb_all(i,ij),i=1,nss)
call flush(iout)
#endif
& jhpb_all(i,ij),i=1,nss)
call flush(iout)
#endif
write (iout,*) "Writing binary file, record",iii," ii",ii
call flush(iout)
#endif
write (iout,*) "Writing binary file, record",iii," ii",ii
call flush(iout)
#endif
+ if (dyn_ss) then
+ write(unit_out,rec=iii) ((allcart(j,i,ij),j=1,3),i=1,nres),
+ & ((allcart(j,i,ij),j=1,3),i=nnt+nres,nct+nres),
+c & nss_all(ij),(ihpb_all(i,ij),jhpb_all(i,ij),i=1,nss_all(ij))
+ & entfac(ii),rmstb(ii)
+ else
write(unit_out,rec=iii) ((allcart(j,i,ij),j=1,3),i=1,nres),
& ((allcart(j,i,ij),j=1,3),i=nnt+nres,nct+nres),
& nss_all(ij),(ihpb_all(i,ij),jhpb_all(i,ij),i=1,nss_all(ij)),
& entfac(ii),rmstb(ii)
write(unit_out,rec=iii) ((allcart(j,i,ij),j=1,3),i=1,nres),
& ((allcart(j,i,ij),j=1,3),i=nnt+nres,nct+nres),
& nss_all(ij),(ihpb_all(i,ij),jhpb_all(i,ij),i=1,nss_all(ij)),
& entfac(ii),rmstb(ii)
#ifdef DEBUG
write (iout,'(8f10.5)') ((allcart(j,i,ij),j=1,3),i=1,nres)
write (iout,'(8f10.4)') ((allcart(j,i,ij),j=1,3),i=nnt+nres,
#ifdef DEBUG
write (iout,'(8f10.5)') ((allcart(j,i,ij),j=1,3),i=1,nres)
write (iout,'(8f10.4)') ((allcart(j,i,ij),j=1,3),i=nnt+nres,