Abort if matrixcn is .true. and nrepr is not 1

chrislxj · Jun 25, 2024 · 18e63e3 · 18e63e3
1 parent 47234a9
commit 18e63e3
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/biogeochem/CropReprPoolsMod.F90 b/src/biogeochem/CropReprPoolsMod.F90
@@ -98,6 +98,10 @@ subroutine crop_repr_pools_init()
     end if
 
     nrepr = nrepr_grain + nrepr_structure
+    ! matrixcn works with nrepr = 1 only
+    if (use_matrixcn .and. nrepr != 1) then
+       call endrun(msg="ERROR: nrepr should be 1 when use_matrixcn = .true."//errMsg(sourcefile, __LINE__))
+    end if
     allocate(repr_hist_fnames(nrepr))
     allocate(repr_rest_fnames(nrepr))
     allocate(repr_longnames(nrepr))