c write(iout,*) 'b1=',b1(1,i-2)
c write (iout,*) 'theta=', theta(i-1)
enddo
+#else
+ b1(1,i-2)=b(3,iti)
+ b1(2,i-2)=b(5,iti)
+ b2(1,i-2)=b(2,iti)
+ b2(2,i-2)=b(4,iti)
+ b1tilde(1,i-2)=b1(1,i-2)
+ b1tilde(2,i-2)=-b1(2,i-2)
+ b2tilde(1,i-2)=b2(1,i-2)
+ b2tilde(2,i-2)=-b2(2,i-2)
+ EE(1,2,i-2)=eeold(1,2,iti)
+ EE(2,1,i-2)=eeold(2,1,iti)
+ EE(2,2,i-2)=eeold(2,2,iti)
+ EE(1,1,i-2)=eeold(1,1,iti)
+ enddo
+#endif
#ifdef PARMAT
do i=ivec_start+2,ivec_end+2
#else
do i=3,nres+1
#endif
-#endif
if (i .lt. nres+1) then
sin1=dsin(phi(i))
cos1=dcos(phi(i))
ghpbc(k,iii)=ghpbc(k,iii)-ggg(k)
enddo
endif
+ endif
enddo
ehpb=0.5D0*ehpb
return
gliptranc(3,i-1)=gliptranc(3,i-1)+ssgradlip*pepliptran/2.0d0
C gliptranc(3,i-2)=gliptranc(3,i)+ssgradlip*pepliptran
-C print *,"doing sccale for lower part"
+C print *,"doing sccale for lower part"
+C print *,i,sslip,fracinbuf,ssgradlip
elseif (positi.gt.bufliptop) then
fracinbuf=1.0d0-((bordliptop-positi)/lipbufthick)
sslip=sscalelip(fracinbuf)
gliptranc(3,i-1)=gliptranc(3,i-1)+ssgradlip*pepliptran/2.0d0
C gliptranc(3,i-2)=gliptranc(3,i)+ssgradlip*pepliptran
C print *, "doing sscalefor top part"
+C print *,i,sslip,fracinbuf,ssgradlip
else
eliptran=eliptran+pepliptran
C print *,"I am in true lipid"
C now multiply all by the peptide group transfer factor
C eliptran=eliptran*pepliptran
C now the same for side chains
-C do i=1,1
+CV do i=1,1
do i=ilip_start,ilip_end
if (itype(i).eq.ntyp1) cycle
positi=(mod(c(3,i+nres),boxzsize))
ssgradlip=-sscagradlip(fracinbuf)/lipbufthick
eliptran=eliptran+sslip*liptranene(itype(i))
gliptranx(3,i)=gliptranx(3,i)
- &+ssgradlip*liptranene(itype(i))/2.0d0
+ &+ssgradlip*liptranene(itype(i))
gliptranc(3,i-1)= gliptranc(3,i-1)
- &+ssgradlip*liptranene(itype(i))/2.0d0
+ &+ssgradlip*liptranene(itype(i))
C print *,"doing sccale for lower part"
elseif (positi.gt.bufliptop) then
fracinbuf=1.0d0-
ssgradlip=sscagradlip(fracinbuf)/lipbufthick
eliptran=eliptran+sslip*liptranene(itype(i))
gliptranx(3,i)=gliptranx(3,i)
- &+ssgradlip*liptranene(itype(i))/2.0d0
+ &+ssgradlip*liptranene(itype(i))
gliptranc(3,i-1)= gliptranc(3,i-1)
- &+ssgradlip*liptranene(itype(i))/2.0d0
+ &+ssgradlip*liptranene(itype(i))
C print *, "doing sscalefor top part",sslip,fracinbuf
else
eliptran=eliptran+liptranene(itype(i))