aLoc = land_grdLambda*land_deltaT*land_rec_dzC(2)
& + cg(2)*tg(2)*aLoc/(cg(2)+aLoc)
& / ( cg(1) + aLoc + bLoc*delT*beta
& - aLoc*aLoc/(cg(2)+aLoc)
tg(2) = ( cg(2)*tg(2) + aLoc*tg(1) ) / (cg(2)+aLoc)
cLoc = cg(1) + aLoc + bLoc*delT*beta
eg(1) = eLoc - (aLoc + bLoc*delT*beta)*tg(1)
eg(2) = eg(2) + aLoc*tg(1)
tg(2) = cg(2)*tg(2) / (cg(2)+aLoc)
eg(2) = eg(2) - aLoc*tg(2)
eg(1) = eg(1) + aLoc*tg(2) + fLoc*delT*beta
& + delT*fLoc - land_Lfreez*mSnow + aLoc*tg(2)
tg(1) = ( cg(1)*tg(1) + cg(2)*tg(2)*aLoc/(cg(2)+aLoc) )
& / ( cg(1)+aLoc + delT/alpha - aLoc*aLoc/(cg(2)+aLoc) )
tg(1) = cg(1)*tg(1) / ( cg(1)+aLoc + delT/alpha )
eg(2) = eg(2) + aLoc*(tg(1)-tg(2))
eg(1) = eg(1) - aLoc*(tg(1)-tg(2))
_RL aLoc, bLoc, cLoc, eLoc, fLoc, alpha, beta