tg(k) = land_groundT(i,j,k,bi,bj)

     &         'LAND_IMPL_TEMP: 0 , ts,tg1&2,mSw=',0,tSurf,tg,mSnow

          tg(1) = ( cg(1)*tg(1) + fLoc*delT*beta

     &                          + cg(2)*tg(2)*aLoc/(cg(2)+aLoc)

          tg(2) = ( cg(2)*tg(2) + aLoc*tg(1) ) / (cg(2)+aLoc)

          tSurf = ( tg(1) + alpha*fLoc ) * beta

     &         'LAND_IMPL_TEMP: 1 , ts,tg1&2,mW1=',1,tSurf,tg,mW(1)

          IF ( tg(2)*land_groundT(i,j,2,bi,bj) .LE. 0. _d 0 

           tmp1 = tg(1)

           tmp2 = tg(2)

           tg(2) = 0.

           tg(1) = MIN( temp_bf, MAX(temp_af, 0. _d 0) )

           tSurf = ( tg(1) + alpha*fLoc ) * beta

             eg(1) = eLoc - (aLoc + bLoc*delT*beta)*tg(1)

             eg(2) = eg(2) + aLoc*tg(1) 

     &              tg,tSurf,mSnow

          IF ( tg(1)*land_groundT(i,j,1,bi,bj) .LE. 0. _d 0 

           tmp1 = tg(1)

           tg(1) = 0.

           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

     &             tg,tSurf,mSnow

           tg(1) = land_groundT(i,j,1,bi,bj)

           tg(2) = land_groundT(i,j,2,bi,bj)

           eLoc = cg(1)*tg(1)

     &          + delT*fLoc - land_Lfreez*mSnow + aLoc*tg(2)

             tg(1) = 0. _d 0

             tg(2) = cg(2)*tg(2) / (cg(2)+aLoc)

            tg(1) = ( cg(1)*tg(1) + cg(2)*tg(2)*aLoc/(cg(2)+aLoc) )

            tg(2) = ( cg(2)*tg(2) + aLoc*tg(1) ) / (cg(2)+aLoc)

            IF ( tg(2)*land_groundT(i,j,2,bi,bj).LE.0. _d 0 ) THEN

              tg(2) = 0.

              tg(1) = cg(1)*tg(1) / ( cg(1)+aLoc + delT/alpha )

            ELSEIF ( tg(1)*land_groundT(i,j,1,bi,bj).LE.0. _d 0 ) THEN

              tg(1) = 0.

              tg(2) = cg(2)*tg(2) / (cg(2)+aLoc)

            dMsn = ( fLoc+tg(1)/alpha )*delT*recip_Lfreez

     &          2, dMsn, fLoc,tg(1)/alpha

           eg(2) = eg(2) + aLoc*(tg(1)-tg(2))

           eg(1) = eg(1) - aLoc*(tg(1)-tg(2))

           eg(2) = eg(2) + aLoc*(tg(1)-tg(2))

           eg(1) = eg(1) - aLoc*(tg(1)-tg(2))

            land_groundT(i,j,k,bi,bj) = tg(k)

     &         tSurf, tg, dTsurf(i,j)

      _RL tSurf, tg(land_nLev), eg(land_nLev)

      tg     = 0.0

         tg  = tact

     &        'td, tg, tact, ifail = ', td, tg, tact, ifail

         test   = barr*(td-tg)

         left = tg+test

     &        'td, tg, tact, ifail = ', td, tg, tact, ifail

         if (td-tg .lt. tmin) go to 920

         if (tg .ne. 0.0) then

               xx(i) = xx(i) + tg*dd(i)

         xx(i) = xx(i) + tg*dd(i)

      double precision    tg, fg, td, ta