integer i, j, kp1, two2
parameter(two2 = 2)
_RL plumekb(two2), SPdepth(two2), minusone
I two2,minusone,SPdepth,