Changeset 4310 for anuga_work/production/gold_coast_2007/project.py

Timestamp:

Mar 21, 2007, 6:25:37 PM (18 years ago)

Author:

sexton

Message:

update damage scripts for FESA areas

File:

• anuga_work/production/gold_coast_2007/project.py (modified) (6 diffs)

anuga_work/production/gold_coast_2007/project.py

@@ -28 +28 @@
 gc250 = 'dem_250'
 
+# better DEM
+gc5 = 'dem_5m_ascii'
+offshore = 'xy_AHD'
+coastline = 'coastline'
+
 #swollen/ all data output
 basename = 'source'
@@ -42 +47 @@
 outputtimedir = outputdir + local_time + sep
 polygondir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'polygons'+sep
+alex_datadir = home+sep+state+sep+scenario_dir_name+sep+'elevation_final'+sep+'points'+sep
 
-#gauge_filename = gaugedir + 'sydney_slide_gauges.csv'
 codedir = getcwd()+sep                            
 codedirname = codedir + 'project.py'
@@ -49 +54 @@
 
 # Necessary if using point datasets, rather than grid
-gc_dem_name = datadir + gc250
+gc_dem_name = alex_datadir + gc5
+offshore_name = datadir + offshore
+coast_name = datadir + coastline
 combined_dem_name   = datadir + 'gc_elevation'
 
+gauge_name = gaugedir + 'gc_gauges.csv'
+buildings_filename = gaugedir + 'gc_nexis.csv'
+buildings_filename_out = 'gc_nexis_modified.csv'
 ###############################
 # Domain definitions
@@ -57 +67 @@
 
 # bounding polygon for study area
-polyAll = read_polygon(polygondir+'extent.csv')
+polyAll = read_polygon(polygondir+'extentv2.csv')
 
 print 'Area of bounding polygon', polygon_area(polyAll)/1000000.0
@@ -70 +80 @@
 print 'Area of local polygon', polygon_area(poly_int)/1000000.0
 
+poly_bathy = read_polygon(polygondir+'mainland_only.csv')
+
 ###################################################################
 # Clipping regions for export to asc and regions for clipping data
@@ -75 +87 @@
 
 # exporting asc grid
-eastingmin = 538450
+eastingmin = 534300
 eastingmax = 546775
-northingmin = 6890540
-northingmax = 6902925
+northingmin = 6890000
+northingmax = 6912925