dSPvolBelow2kLev = -dSPvolkLev2Above(i,j,kp1,bi,bj)
gTr_Below2kLev= dSPvolBelow2kLev * work(i,j,kp1)
& dSPvolBelow2kLev, gTr_totSurf2Below,