diff --git a/biogeochem/EDPatchDynamicsMod.F90 b/biogeochem/EDPatchDynamicsMod.F90
index 9622669409..5c3db014b9 100644
--- a/biogeochem/EDPatchDynamicsMod.F90
+++ b/biogeochem/EDPatchDynamicsMod.F90
@@ -570,7 +570,6 @@ subroutine spawn_patches( currentSite, bc_in, bc_out)
     ! zero the diagnostic disturbance rate fields
     currentSite%disturbance_rates(:,:,:) = 0._r8
     bc_out%fire_closs_to_atm_si = 0._r8
-    bc_out%grazing_closs_to_atm_si = 0._r8
     ! get rules for vegetation clearing during land use change
     call GetLanduseChangeRules(clearing_matrix)
     
diff --git a/main/EDMainMod.F90 b/main/EDMainMod.F90
index afb3c64a54..c089c3fd9d 100644
--- a/main/EDMainMod.F90
+++ b/main/EDMainMod.F90
@@ -420,7 +420,8 @@ subroutine ed_integrate_state_variables(currentSite, bc_in, bc_out )
     ! Clear site GPP and AR passing to HLM
     bc_out%gpp_site = 0._r8
     bc_out%ar_site = 0._r8
-
+    bc_out%grazing_closs_to_atm_si = 0._r8
+    
     ! Patch level biomass are required for C-based harvest
     call get_harvestable_carbon(currentSite, bc_in%site_area, bc_in%hlm_harvest_catnames, harvestable_forest_c)