O TG1, QG1, SE, VsurfSq, PSG, dpFac, kGround,
I kGround,bi,bj,myThid)
K = kGround(J)
INTEGER kGround(NGP)