Merge pull request #614 from cgannonucm/dynamicalTimeUpdate

fix: floating point exception in nodePropertyExtractorSatelliteDynamicalTime
galacticusorg · May 19, 2024 · 382632c · 382632c
2 parents 1e63a6f + d75b746
commit 382632c
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                              &