Changeset 4063 for anuga_work/production/newcastle_2006/project_slide.py

Timestamp:

Dec 6, 2006, 5:35:00 PM (17 years ago)

Author:

sexton

Message:

updates for nsw slide modelling and cairns demo

File:

• anuga_work/production/newcastle_2006/project_slide.py (modified) (4 diffs)

anuga_work/production/newcastle_2006/project_slide.py

@@ -62 +62 @@
 
 # Necessary if using point datasets, rather than grid
-on_offshore_dem_name = datadir + on_offshore_name
+nsw_dem_name = datadir + nsw100
 offshore_dem_name2 = datadir + offshore_name2
 offshore_dem_name3 = datadir + offshore_name3
@@ -78 +78 @@
 
 # bounding polygon for study area
-polyAll = read_polygon(polygondir+'extentC.csv')
+polyAll = read_polygon(polygondir+'extentC2.csv')
 
 print 'Area of bounding polygon', polygon_area(polyAll)/1000000.0
@@ -87 +87 @@
 
 # newcastle digitized polygons
-poly_newcastle1 = read_polygon(polygondir+'.csv')
-poly_newcastle2 = read_polygon(polygondir+'.csv')
-poly_newcastle3 = read_polygon(polygondir+'.csv')
+poly_local = read_polygon(polygondir+'local2.csv')
+poly_newcastle = read_polygon(polygondir+'newcastle.csv')
 
-print 'Area of local polygon', polygon_area(poly_newcastle1)/1000000.0
-print 'Area of close polygon', polygon_area(poly_newcastle2)/1000000.0
-print 'Area of coastal polygon', polygon_area(poly_newcastle3)/1000000.0
+print 'Area of local polygon', polygon_area(poly_local)/1000000.0
+print 'Area of CBD polygon', polygon_area(poly_newcastle)/1000000.0
 
+#plot_polygons([polyAll,poly_local,poly_newcastle],'fig',verbose=True)
 ###################################################################
 # Clipping regions for export to asc and regions for clipping data
@@ -100 +99 @@
 
 # clipping for nsw grid for conversion
-eastingmin_nsw = 546000
-eastingmax_nsw = 350390
+eastingmin_nsw = 350390
+eastingmax_nsw = 546000
 northingmin_nsw = 6253970
 northingmax_nsw = 6424600
 
+poly_surveyclip = read_polygon(polygondir+'surveyclip.csv')
+
 # exporting asc grid
-eastingmin = 
-eastingmax = 
-northingmin = 
-northingmax = 
+#eastingmin = 
+#eastingmax = 
+#northingmin = 
+#northingmax = 
 
 ###################################################################