From 280e85385b0f836749f9300bbaeb161e535eafef Mon Sep 17 00:00:00 2001
From: Shijie Shu <92333861+sshu88@users.noreply.github.com>
Date: Wed, 20 Nov 2024 11:03:59 -0800
Subject: [PATCH] Bug fix for issue #1289

Fix the issue of no wood product when forced by area-based harvest.
https://github.com/NGEET/fates/issues/1289
---
 biogeochem/EDLoggingMortalityMod.F90 | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/biogeochem/EDLoggingMortalityMod.F90 b/biogeochem/EDLoggingMortalityMod.F90
index 358c211e09..147eeab251 100644
--- a/biogeochem/EDLoggingMortalityMod.F90
+++ b/biogeochem/EDLoggingMortalityMod.F90
@@ -345,7 +345,7 @@ subroutine LoggingMortality_frac( currentSite, bc_in, pft_i, dbh, canopy_layer,
          ! transfer of area to secondary land is based on overall area affected, not just logged crown area
          ! l_degrad accounts for the affected area between logged crowns
          if(prt_params%woody(pft_i) == itrue)then ! only set logging rates for trees
-            if (cur_harvest_tag == 0) then
+            if (cur_harvest_tag == 0 .or. cur_harvest_tag == 2) then
                ! direct logging rates, based on dbh min and max criteria
                if (dbh >= logging_dbhmin .and. .not. &
                     ((logging_dbhmax < fates_check_param_set) .and. (dbh >= logging_dbhmax )) ) then