Fix seg-fault with DBEMT_Mod=-1

andrew-platt · Sep 18, 2024 · 040d6be · 040d6be
1 parent c67de45
commit 040d6be
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/modules/aerodyn/src/BEMT.f90 b/modules/aerodyn/src/BEMT.f90
@@ -929,7 +929,7 @@ subroutine BEMT_UpdateStates( t, n, u, utimes, p, x, xd, z, OtherState, AFInfo,
    !...............................................................................................................................
    !  update DBEMT states to step n+1
    !...............................................................................................................................
-   if (p%DBEMT_Mod /= DBEMT_none) then
+   if (p%DBEMT_Mod > DBEMT_none) then
 
       !........................
       ! update DBEMT states to t+dt
@@ -956,7 +956,7 @@ subroutine BEMT_UpdateStates( t, n, u, utimes, p, x, xd, z, OtherState, AFInfo,
             !............................................
             ! apply DBEMT correction to axInduction and tanInduction:
             !............................................
-            if (p%DBEMT_Mod /= DBEMT_none) then
+            if (p%DBEMT_Mod > DBEMT_none) then
                call calculate_Inductions_from_DBEMT_AllNodes(TimeIndex_t_plus_dt, uTimes(TimeIndex_t_plus_dt), u(TimeIndex_t_plus_dt), p, x, OtherState, m, m%axInduction, m%tanInduction)
             end if