write giving a rodeo

Author: jginternational

kratos.gid/scripts/Mesh.tcl

@@ -72,11 +72,12 @@ proc Mesh::AddPointElementsIfNeeded {} {
         if {$cond eq ""} {
             continue
         }
+        
+        # if group element is point and has topology for points
         set topology [$cond getTopologyFeature "Point" 1]
         if {$topology eq ""} {
             continue
         }
-        W " Groups assigned to condition $condid will be meshed with Point elements."
 
         set group_nodes [dict get $condition_groups $condid]
         foreach node_tdom $group_nodes {
@@ -95,7 +96,4 @@ proc Mesh::AddPointElementsIfNeeded {} {
         }
     }
 
-    # if group element is point and has topology for points
-
-    # foreach nodeid [GiD_Mesh] GiD_Mesh create element append Point 1 [list $nodeid]
 }

kratos.gid/scripts/Writing/WriteGeometries.tcl

@@ -63,6 +63,33 @@ proc write::printGeometryConnectivities {group etype nnodes} {
         if {$etype == "Sphere" || $etype == "Circle"} {
             write::writeSphereRadiusOnGroup $group
         }
+    } else {
+        # Trick: GiD < 17.3.x return 0 if elements are of type Point
+        set elems [GiD_EntitiesGroups get $group elements -element_type point]
+        set num_elems [objarray length $elems]
+        if {$num_elems > 0} {
+            # Write header
+            set geometry_name ${etype}${nDim}
+            WriteString "${s}Begin Geometries $geometry_name // GUI group identifier: $group"
+            # increase indent (allows folding in text editor)
+            incr ::write::current_mdpa_indent_level
+            # Write the connectivities
+            set s1 [mdpaIndent]
+            objarray foreach elem $elems {
+                set node_id [GiD_Mesh get element $elem connectivities]
+                GiD_WriteCalculationFile puts "${s1}$elem $node_id"
+            }
+            # decrease indent
+            incr ::write::current_mdpa_indent_level -1
+            # Write footer
+            WriteString "${s}End Geometries"
+            WriteString ""
+
+            # Write the radius if it is a sphere or a circle
+            if {$etype == "Sphere" || $etype == "Circle"} {
+                write::writeSphereRadiusOnGroup $group
+            }
+        }
     }
     if {[GetConfigurationAttribute time_monitor]} {set endtime [clock seconds]; set ttime [expr {$endtime-$inittime}]; W "printGeometryConnectivities $geometry_name time: [Kratos::Duration $ttime]"}
 }