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))