subroutine lipid_layer(xi,yi,zi,sslipi,ssgradlipi)
implicit none
include 'DIMENSIONS'
+ include 'COMMON.IOUNITS'
include 'COMMON.CHAIN'
double precision xi,yi,zi,sslipi,ssgradlipi
double precision fracinbuf
double precision sscalelip,sscagradlip
-
+#define DEBUG
+#ifdef DEBUG
+ write (iout,*) "bordlipbot",bordlipbot," bordliptop",bordliptop
+ write (iout,*) "buflipbot",buflipbot," lipbufthick",lipbufthick
+ write (iout,*) "xi yi zi",xi,yi,zi
+#endif
if ((zi.gt.bordlipbot).and.(zi.lt.bordliptop)) then
C the energy transfer exist
if (zi.lt.buflipbot) then
sslipi=0.0d0
ssgradlipi=0.0
endif
+#ifdef DEBUG
+ write (iout,*) "sslipi",sslipi," ssgradlipi",ssgradlipi
+#endif
+#undef DEBUG
return
end