make sure to define state_vector if not in land use mode

NGEET · Jul 15, 2024 · a2b71d0 · a2b71d0
1 parent fd7f343
commit a2b71d0
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/biogeochem/EDPatchDynamicsMod.F90 b/biogeochem/EDPatchDynamicsMod.F90
@@ -325,6 +325,8 @@ subroutine disturbance_rates( site_in, bc_in)
        call GetLUHStatedata(bc_in, state_vector)
        site_secondaryland_first_exceeding_min =  (state_vector(secondaryland) .gt. site_in%min_allowed_landuse_fraction) &
             .and. (.not. site_in%landuse_vector_gt_min(secondaryland))
+    else
+       state_vector = current_fates_landuse_state_vector
     end if
 
     currentPatch => site_in%oldest_patch