Changeset 1837

Timestamp:

Sep 15, 2005, 5:37:47 PM (19 years ago)

Author:

ole

Message:

More karratha work

Location:

production/karratha_2005

Files:

• create_mesh.py (modified) (2 diffs)
• project.py (modified) (1 diff)
• run_karratha.py (modified) (2 diffs)

production/karratha_2005/create_mesh.py

@@ -12 +12 @@
 
 
-def create_mesh_boundary(polygon, tags, resolution, filename = None):
+def create_mesh(polygon, tags, resolution, filename = None):
     """Create mesh from bounding polygon, tags for all segments and resolution.
 
@@ -24 +24 @@
 
     FIXME - FUTURE:
-      Use class Point
+      Use class Point_set
       Take multiple polygons
 

production/karratha_2005/project.py

@@ -42 +42 @@
 #east = degminsec2decimal_degrees(117,0,0)
 
-south = degminsec2decimal_degrees(-20,45,0)
+south = degminsec2decimal_degrees(-20,50,0)
 north = degminsec2decimal_degrees(-20,15,0)
 west = degminsec2decimal_degrees(116,20,0)

production/karratha_2005/run_karratha.py

@@ -43 +43 @@
 
 #Mesh
-from create_mesh import create_mesh_boundary
+from create_mesh import create_mesh
 
-m = cache(create_mesh_boundary,
+m = cache(create_mesh,
           project.polygon,
           {'tags': {'back': [7, 8], 'side': [0, 6], 'ocean': [1, 2, 3, 4, 5]},
           'resolution': 52000,
+          #'resolution': 150000,           
           'filename': project.meshname + '.msh'},
           verbose = True)
-
-
 
 
@@ -92 +91 @@
 #Setup IC
 tide = 0
+domain.set_quantity('friction', 0.03)
 domain.set_quantity('stage', tide)
 domain.set_quantity('elevation',