C sscale=1.0d0
C else if(r.le.r_cut.and.r.ge.r_cut-rlamb) then
C gamm=(r-(r_cut-rlamb))/rlamb
- sscale=1.0d0+r*r*(2*r-3.0d0)
+ sscalelip=1.0d0+r*r*(2*r-3.0d0)
C else
C sscale=0d0
C endif
C sscagrad=0.0d0
C else if(r.le.r_cut.and.r.ge.r_cut-rlamb) then
C gamm=(r-(r_cut-rlamb))/rlamb
- sscagrad=r*(6*r-6.0d0)
+ sscagradlip=r*(6*r-6.0d0)
C else
C sscagrad=0.0d0
C endif