Changeset 5931 for anuga_work/debug/will_power_example

Timestamp:

Nov 10, 2008, 4:44:39 PM (16 years ago)

Author:

duncan

Message:

example of adding segments and vertices to a mesh file, from a points file.

File:

• anuga_work/debug/will_power_example/make_mesh.py (modified) (4 diffs)

anuga_work/debug/will_power_example/make_mesh.py

@@ -1 @@
-""" This Script shows how a bad mesh can be generated causing
-excessive values for stage and momentum.
+""" This Script shows how exta points and segments can be added to a msh file.
 
+# WARNING #
+* this doesn't have tests assosiated with it.  For example that the
+resolutions stay as described in create_mesh_from_regions is not tested
+* Also, this code causes the mesh to be generated twice, so it is not
+an efficient way of doing this.
+
+If this is used a lot you may want to wrap the functionality into a function.
 """
 
@@ -36 +42 @@
 file.write("x,y\n\
 304200.0, 6185300.0\n\
-304250.0, 6185350.0\n\
-304300.0, 6185400.0\n")
+304300.0, 6185400.0\n\
+304400.0, 6185500.0\n")
 file.close()
 
@@ -67 +73 @@
 m = importMeshFromFile(meshname)
 
-m.add_vertices([(W +500,S +500),(W +1000,S +1000) ])
+#m.add_vertices([(W +500,S +500),(W +1000,S +1000) ])
 
 # load points from a points file
-m.add_vertices(Geospatial_data(ptfile))
+#m.add_vertices(Geospatial_data(ptfile))
 
 p1 = (W +500,N - 500)
 p2 = (W +1000,N - 1000)
 p3 = (W +1000,N - 1500) 
-m.add_points_and_segments([p1,p2,p3], [[0,1],[1,2]])
+#m.add_points_and_segments([p1,p2,p3], [[0,1],[1,2]])
 
 # add points and segments, where the segments are automatically added.
@@ -83 +89 @@
 p2 = (E -1000,N - 1500) 
 #m.add_points_and_segments([p1,p2])
+
+# load points from a points file
+m.add_points_and_segments(Geospatial_data(ptfile))
 
 m.generate_mesh()