Changeset 1086

Timestamp:

Mar 16, 2005, 2:05:39 PM (20 years ago)

Author:

duncan

Message:

4 regions

File:

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

inundation/ga/storm_surge/pmesh/create_gong.py

@@ -12 +12 @@
     bellambi_northing = 6195000.0
     north_beach_northing = 6183000.0
+    north_beach_easting = 308500
     border = 12000.0
     west = 301000.0 # From the map 
@@ -21 +22 @@
     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 = inner_center_easting,
@@ -33 +34 @@
     
     dict = {}
-    dict['pointlist'] = [[west, north_beach_northing - border],  #sw
+    dict['points'] = [[west, north_beach_northing - border],  #sw
                          [west, bellambi_northing + border],  #nw
                          [sea_boundary,
@@ -42 +43 @@
                           north_beach_northing - border]    #s_sea_boundary
                          ]
-    dict['segmentlist'] = [[0,1],[1,2],[2,3],
+    dict['segments'] = [[0,1],[1,2],[2,3],
                            [3,4],[4,5],[5,0], # the outer boarder
                            [2,5]   # the sea_boundary line
@@ -50 +51 @@
 
     dict = {}
-    dict['pointlist'] = [
+    dict['points'] = [
                          [inner_west_south,north_beach_northing], #sw
                          [inner_west_north,bellambi_northing],   #nw
@@ -56 +57 @@
                          [inner_east,north_beach_northing] #se
                          ]
-    dict['segmentlist'] = [[0,1],[1,2],[2,3],[3,0]] # the inner boarder
+    dict['segments'] = [[0,1],[1,2],[2,3],[3,0]] # the inner boarder
     m.addVertsSegs(dict)
 
+    #Even more inner points
+    flagstaff_x = 308000
+    flagstaff_y = 6189000
+    west_x = 305000
+    west_y = 6185000
+    
+    north_beach_northing = 6183000.0
+    
+    dict = {}
+    dict['points'] = [
+                         [flagstaff_x,flagstaff_y], #n
+                         [west_x,west_y],   #w
+                         [north_beach_easting,north_beach_northing], #s
+                         ]
+    dict['segments'] = [[0,1],[1,2],[2,0]] # the inner boarder
+    m.addVertsSegs(dict)
+    
     factor = 1000 #low res 10000, high res 1000
     low = m.addRegionEN(east-0.5, bellambi_northing + border-0.5)
@@ -68 +86 @@
     high = m.addRegionEN(inner_east-0.5,bellambi_northing-0.5)
     high.setMaxArea(5*factor) #doing 0.5, fac = 1000 gives degenerate error
+    
+    center = m.addRegionEN(flagstaff_x,flagstaff_y-0.5)
+    center.setMaxArea(0.05*factor) #doing 0.5, fac = 1000 gives degenerate error
     #m.generateMesh()
     m.generateMesh("pzq28.0za1000000a") 
     
-    m.export_mesh_file("wollongong_highres_a28.msh")
+    m.export_mesh_file("wollongong_4regions.msh")
     #m.export_mesh_file("wollongong_highres_a28.tsh")