Changeset 1048

Timestamp:

Mar 9, 2005, 4:41:37 PM (20 years ago)

Author:

duncan

Message:

update

File:

• inundation/ga/storm_surge/pmesh/create_gong.py (modified) (2 diffs)

inundation/ga/storm_surge/pmesh/create_gong.py

@@ -9 +9 @@
 if __name__ == "__main__":
 
+
+    bellambi_northing = 6195000.0
+    north_beach_northing = 6183000.0
+    border = 12000.0
+    west = 301000.0 # From the map 
+    inner_east = 310000.0 # From the map 
+    inner_west_south = 303000.0 # From the map 
+    inner_west_north = 307000.0  # From the map 
+    east = 328421.0
+    sea_boundary = (east+inner_east)/2
+    inner_center_easting =  ( inner_east +inner_west_north)/2
+    inner_center_northing = (bellambi_northing + north_beach_northing)/2
+    
     # from 100m dem
-    geo = Geo_reference(xllcorner = 241191.35087962,
-                        yllcorner = 6111741.8729311,
+    geo = Geo_reference(xllcorner = inner_center_easting,
+                        yllcorner = inner_center_northing,
                         zone = 56)
+    print "***********************"
+    print "geo ref", geo
+    print "***********************"
+    
     m = Mesh(geo_reference=geo)
-    bellambi_northing = 6195000
-    north_beach_northing = 6183000
-    border = 12000
-    west = 301000 # From the map 
-    inner_east = 310000 # From the map 
-    inner_west_south = 303000 # From the map 
-    inner_west_north = 307000  # From the map 
-    east = 328421
-    center = (east+inner_east)/2
     
     dict = {}
     dict['pointlist'] = [[west, north_beach_northing - border],  #sw
                          [west, bellambi_northing + border],  #nw
-                         [center,
-                          bellambi_northing + border],  #n_center
+                         [sea_boundary,
+                          bellambi_northing + border],  #n_sea_boundary
                          [east, bellambi_northing + border], #ne
                          [east, north_beach_northing - border],    #se
-                         [center,
-                          north_beach_northing - border]    #s_center
+                         [sea_boundary,
+                          north_beach_northing - border]    #s_sea_boundary
                          ]
     dict['segmentlist'] = [[0,1],[1,2],[2,3],
                            [3,4],[4,5],[5,0], # the outer boarder
-                           [2,5]   # the center line
+                           [2,5]   # the sea_boundary line
                            ]
     
@@ -51 +59 @@
     m.addVertsSegs(dict)
 
-    factor = 10000 #low res 10000, high res 1000
+    factor = 1000 #low res 10000, high res 1000
     low = m.addRegionEN(east-0.5, bellambi_northing + border-0.5)
     low.setMaxArea(100*factor)
     
-    medium = m.addRegionEN(center-0.5, bellambi_northing + border-0.5)
+    medium = m.addRegionEN(sea_boundary-0.5, bellambi_northing + border-0.5)
     medium.setMaxArea(10*factor)
     
     high = m.addRegionEN(inner_east-0.5,bellambi_northing-0.5)
-    high.setMaxArea(1*factor)
+    high.setMaxArea(5*factor) #doing 0.5, fac = 1000 gives degenerate error
     m.generateMesh()
     
-    m.export_triangulation_file("wollongong_outline.msh")
+    m.export_triangulation_file("wollongong_highres.msh")