hs(i) = 0.0

           grid(i,lats) = hs(i)

           hs(i) = grid(i,lats)

        hs(i) = hs(i)*scale

         hs  = hSnow(i,j)

         albedo = albsno + (albice-albsno)*EXP(-hs/hAlbSnow)

          sAlbNIR(i,j) = albNIR_ice * ( 1. _d 0 - hs/(hs + 0.02 _d 0) )

     &         * hs/(hs + 0.02 _d 0)

      _RL  hs                  ! snow height

      COMMON /exfl_turb_r/ hs, hl

      _RL hs        (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)

           hs(i,j,bi,bj) = atmcp*tau(i,j)*tstar(i,j)

           hs     (i,j,bi,bj) = 0. _d 0

         CALL DIAGNOSTICS_FILL(hs    ,     'EXFhs   ',0,1,0,1,1,myThid)

     &              - hs(i,j,bi,bj) 

          hs(i,j,bi,bj) = 0.