I two2,minusone,SPdepth,
integer two2
parameter(two2 = 2)
_RL plumekb(two2), SPdepth(two2)
_RL lon(two2), lat(two2)