indpdb=index(controlcard,'PDBSTART')
extconf=(index(controlcard,'EXTCONF').gt.0)
AFMlog=(index(controlcard,'AFM'))
- print *,'AFMlog',AFMlog
+ selfguide=(index(controlcard,'SELFGUIDE'))
+ print *,'AFMlog',AFMlog,selfguide,"KUPA"
call readi(controlcard,'IPRINT',iprint,0)
call readi(controlcard,'MAXGEN',maxgen,10000)
call readi(controlcard,'MAXOVERLAP',maxoverlap,1000)
bordliptop=(boxzsize+lipthick)/2.0
bordlipbot=bordliptop-lipthick
C endif
- if ((bordliptop.gt.boxzsize).or.(borlipbot.lt.0.0))
+ if ((bordliptop.gt.boxzsize).or.(bordlipbot.lt.0.0))
& write(iout,*) "WARNING WRONG SIZE OF LIPIDIC PHASE"
buflipbot=bordlipbot+lipbufthick
bufliptop=bordliptop-lipbufthick
call flush(iout)
if (constr_dist.gt.0) call read_dist_constr
write (iout,*) "After read_dist_constr nhpb",nhpb
- if (AFMlog.gt.0) call read_afminp
+ if ((AFMlog.gt.0).or.(selfguide.gt.0)) call read_afminp
call hpb_partition
if(me.eq.king.or..not.out1file)
& write (iout,*) 'Contact order:',co
include 'COMMON.MD'
open(irest2,file=rest2name,status='unknown')
read(irest2,*) totT,EK,potE,totE,t_bath
+ totTafm=totT
do i=1,2*nres
read(irest2,'(3e15.5)') (d_t(j,i),j=1,3)
enddo
call readi(afmcard,"BEG",afmbeg,0)
call readi(afmcard,"END",afmend,0)
call reada(afmcard,"FORCE",forceAFMconst,0.0d0)
+ call reada(afmcard,"VEL",velAFMconst,0.0d0)
print *,'FORCE=' ,forceAFMconst
CCCC NOW PROPERTIES FOR AFM
distafminit=0.0d0