From df6c5e7c9c6ecf1a998c199a72e031175ae856ac 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 73ea7231e7..a091e4f3d6 100644
--- a/biogeochem/EDLoggingMortalityMod.F90
+++ b/biogeochem/EDLoggingMortalityMod.F90
@@ -291,7 +291,7 @@ subroutine LoggingMortality_frac( pft_i, dbh, canopy_layer, lmort_direct, &
          ! 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