fix(GwfGwtExchange): check node counts before comparison (MODFLOW-USG…

…S#1478)
geohackweek · Dec 5, 2023 · 7d0dc9c · 7d0dc9c
1 parent 4185096
commit 7d0dc9c
Showing 1 changed file with 8 additions and 4 deletions.
diff --git a/src/Exchange/GwfGwtExchange.f90 b/src/Exchange/GwfGwtExchange.f90
@@ -356,11 +356,15 @@ subroutine gwfconn2gwtconn(this, gwfModel, gwtModel)
           !
           if (associated(gwfEx%model1, gwfModel) .or. &
               associated(gwfEx%model2, gwfModel)) then
-            ! again, connecting the same nodes nrs will be
+
+            ! check exchanges have same node counts
+            areEqual = size(gwfEx%nodem1) == size(gwtConn%prim_exchange%nodem1)
+            ! then, connecting the same nodes nrs will be
             ! sufficient evidence of equality
-            areEqual = all(gwfEx%nodem1 == gwtConn%prim_exchange%nodem1)
-            areEqual = areEqual .and. &
-                       all(gwfEx%nodem2 == gwtConn%prim_exchange%nodem2)
+            if (areEqual) &
+              areEqual = all(gwfEx%nodem1 == gwtConn%prim_exchange%nodem1)
+            if (areEqual) &
+              areEqual = all(gwfEx%nodem2 == gwtConn%prim_exchange%nodem2)
             if (areEqual) then
               ! link exchange to connection
               write (iout, '(/6a)') 'Linking exchange ', &