fix: floating point exception in nodePropertyExtractorSatelliteDynami…

…calTime Fixed floating point exception in nodePropertyExtractorSatelliteDynamicalTime. Check that both the tidal radius and tidal mass are non zero. If they are zero, return -1 to avoid divide by zero.
galacticusorg · May 18, 2024 · d75b746 · d75b746
1 parent 1e63a6f
commit d75b746
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/source/nodes.property_extractor.satellite.dynamical_time.F90 b/source/nodes.property_extractor.satellite.dynamical_time.F90
@@ -130,9 +130,11 @@ double precision function dynamicalTimeExtract(self,node,instance)
     type            (multiCounter                               ), intent(inout), optional :: instance
     double precision                                                                       :: radiusTidal, massTidal
     !$GLC attributes unused :: instance
-
+    dynamicalTimeExtract  =-1.0d0
     radiusTidal           =self%satelliteTidalStrippingRadius_%radius      (node            )
-    massTidal             =self%galacticStructure_            %massEnclosed(node,radiusTidal)
+    if (radiusTidal <= 0.0d0) return
+    massTidal             =self%galacticStructure_            %massEnclosed(node,radiusTidal) 
+    if (massTidal <= 0.0d0) return
     dynamicalTimeExtract  =+sqrt(                                    &  
          &                       +Pi                             **2 &
          &                       /4.0d0                              &