dd20 = (abs(SPDepth(i)))
IF ( dd20 .GT. zero ) THEN
S = (facz/dd20)
So = dd20*So
dd = dd20/SPovershoot
So=dd20-dd
IF( (facz.GE.dd).AND.(facz.LT.dd20) ) THEN
dd20 = one
So = 1. _d 0/dd20
plumek(i) = max(zero,two/dd20*facz-S/So)
& three/dd20*facz - three/(dd20*dd20)*facz*facz + S/So)
_RL facz, dd, dd20