_RL              cstar, cg

      common /kmixcbm/ cstar, cg

      cg      = cstar * vonk * (concs * vonk * epsilon)**p33

            ghat(i,ki) = (1.-stable(i)) * cg / max(phepsi,tempVar)

            cg(k) = land_dzF(k)*land_heatCs + mW(k)*land_CpWater

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

     &                          + 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

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

           eLoc = cg(1)*tg(1)

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

            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(2) = ( cg(2)*tg(2) + aLoc*tg(1) ) / (cg(2)+aLoc)

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

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

      _RL cg(land_nLev), mW(land_nLev)