do not write geometries if topology is empty

Author: jginternational

kratos.gid/apps/Fluid/write/write.tcl

@@ -79,7 +79,13 @@ proc ::Fluid::write::writeModelPartEvent { } {
             if {$condition ne "" && [$condition getGroupBy] eq "Condition"} {
                 dict lappend grouped_conditions $condition [$group @n]
             } else {
-                write::writeGroupSubModelPartAsGeometry [$group @n]
+                # Check topology in condition
+                set geoms 1
+                if {$condition ne "" && ![$condition hasTopologyFeatures]} {
+                    set geoms 0
+                }
+                write::writeGroupSubModelPartAsGeometry [$group @n] $geoms
+                
             }
         }
 
@@ -89,7 +95,11 @@ proc ::Fluid::write::writeModelPartEvent { } {
             set new_group_name "_HIDDEN_$condition_name"
             set groups [dict get $grouped_conditions $condition]
             set new_group [spdAux::MergeGroups $new_group_name $groups]
-            write::writeGroupSubModelPartAsGeometry $new_group_name 
+            set geoms 1
+            if {![$condition hasTopologyFeatures]} {
+                set geoms 0
+            }
+            write::writeGroupSubModelPartAsGeometry $new_group_name $geoms
             GiD_Groups delete $new_group_name
         }
 

kratos.gid/scripts/Writing/WriteSubModelPart.tcl

@@ -70,7 +70,7 @@ proc write::writeGroupSubModelPart { cid group {what "Elements"} {iniend ""} {ta
 }
 
 # 
-proc write::writeGroupSubModelPartAsGeometry { group } {
+proc write::writeGroupSubModelPartAsGeometry { group { write_geometries 1} } {
     variable submodelparts
     variable formats_dict
     variable geometry_cnd_name
@@ -106,10 +106,12 @@ proc write::writeGroupSubModelPartAsGeometry { group } {
         WriteString "${s1}Begin SubModelPartNodes"
         GiD_WriteCalculationFile nodes -sorted $gdict
         WriteString "${s1}End SubModelPartNodes"
-        # Print geometries
-        WriteString "${s1}Begin SubModelPartGeometries"
-        GiD_WriteCalculationFile elements -sorted $gdict
-        WriteString "${s1}End SubModelPartGeometries"
+        if {$write_geometries} {
+            # Print geometries only if requested, so conditions with no topology features are not printed
+            WriteString "${s1}Begin SubModelPartGeometries"
+            GiD_WriteCalculationFile elements -sorted $gdict
+            WriteString "${s1}End SubModelPartGeometries"
+        }
 
         WriteString "${s}End SubModelPart"
     }