Merge pull request #1223 from glemieux/landuse-v2-nanfix

Make sure to define `state_vector` if not in land use mode
NGEET · Jul 17, 2024 · 1982b00 · 1982b00
2 parents fd7f343 + a2b71d0
commit 1982b00
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