anchor.diff:                  domain%lwup=stefan_boltzmann*EMISS*domain%skin_t**4
                   RAINBL=domain%rain
                   
  -                print*, "WARNING!! enforcing surface sensible heat flux greater than 0!!"
  -                where(domain%sensible_heat<0) domain%sensible_heat=0
  -                do j=1,ny
  -                    do i=1,nx
  -                        if (domain%landmask(i,j)==kLC_LAND) then
  -                            if (domain%sensible_heat(i,j)<0) domain%sensible_heat(i,j)=0
  -                        endif
  -                    end do
  -                end do
  +!                 print*, "WARNING!! enforcing surface sensible heat flux greater than 0!!"
  +!                 where(domain%sensible_heat<0) domain%sensible_heat=0
  +!                 do j=1,ny
  +!                     do i=1,nx
  +!                         if (domain%landmask(i,j)==kLC_LAND) then
  +!                             if (domain%sensible_heat(i,j)<0) domain%sensible_heat(i,j)=0
  +!                         endif
  +!                     end do
  +!                 end do
                   
               endif
               call surface_diagnostics(domain%sensible_heat, QFX, domain%skin_t, QSFC,  &
  
positive.code:                         cur_vegmonth=domain%current_month
                      endif
                  endif
                  call lsm_noah(domain%dz_inter,domain%qv,domain%p_inter,domain%th*domain%pii,domain%skin_t,  &
                              domain%sensible_heat,QFX,domain%latent_heat,domain%ground_heat, &
                              QGH,GSW,domain%swdown,domain%lwdown,SMSTAV,domain%soil_totalmoisture, &
  
negative.code:             endif
              cur_vegmonth=domain%current_month
              
              ! prevents init from failing when processing water poitns that may have "soil_t"=0
              where(domain%soil_t<200) domain%soil_t=200
              where(domain%soil_vwc<0.001) domain%soil_vwc=0.001
              call LSM_NOAH_INIT(VEGFRAC,SNOW,SNOWC,SNOWH,domain%canopy_water,domain%soil_t,    &
                                  domain%soil_vwc, SFCRUNOFF,UDRUNOFF,ACSNOW,  &
                                  ACSNOM,domain%veg_type,domain%soil_type,     &
  