if (me.eq.king) call cinfo
C Read force field parameters and job setup data
call readrtns
- call flush(iout)
C
if (me.eq.king .or. .not. out1file) then
write (iout,'(2a/)')
include 'COMMON.SETUP'
include 'COMMON.CONTROL'
include 'COMMON.IOUNITS'
- if (me.eq.king .or. .not. out1file)
- & write (iout,*) "Calling chainbuild"
+c if (me.eq.king .or. .not. out1file) then
+c write (iout,*) "Calling chainbuild"
+c call flush(iout)
+c endif
call chainbuild
+c if (me.eq.king .or. .not. out1file) then
+c write (iout,*) "Calling MD"
+c call flush(iout)
+c endif
call MD
return
end
time1=tcpu()
#endif
call minim_dc(etot,iretcode,nfun)
+ if(iretcode.eq.8) call check_ecartint
else
if (indpdb.ne.0) then
call bond_regular
include 'COMMON.SBRIDGE'
double precision varia(maxvar)
dimension ind(6)
- double precision energy(0:max_ene)
+ double precision energy(0:n_ene)
logical eof
eof=.false.
#ifdef MPI
include 'COMMON.MD'
include 'COMMON.SBRIDGE'
common /srutu/ icall
- double precision energy(0:max_ene)
+ double precision energy(0:n_ene)
c do i=2,nres
c vbld(i)=vbld(i)+ran_number(-0.1d0,0.1d0)
c if (itype(i).ne.10)
include 'DIMENSIONS'
include 'COMMON.IOUNITS'
include 'COMMON.CONTROL'
- double precision energy(0:max_ene)
+ double precision energy(0:n_ene)
call chainbuild
call etotal(energy(0))
call enerprint(energy(0))