gTr_Below2kLev= dSPvolBelow2kLev * work(i,j,kp1)
& ( gTr_Surf2kLev + gTr_Below2kLev + gTr_kLev2Above )
& ( gTr_Surf2kLev + gTr_Below2kLev)
_RL gTr_Surf2kLev, gTr_Below2kLev, gTr_kLev2Above,