Changeset 3650 for anuga_work/production/hobart_2006/project.py

Timestamp:

Sep 22, 2006, 12:21:09 PM (18 years ago)

Author:

sexton

Message:

updates for Hobart and want to export elevation grid for Onslow

File:

• anuga_work/production/hobart_2006/project.py (modified) (4 diffs)

anuga_work/production/hobart_2006/project.py

@@ -13 +13 @@
 
 # data provided by Tas SES and checked by NM&I
-onshore_name = 'hob3_topo' 
+#onshore_name = 'hob3_topo' # original
+onshore_name = 'hob5_topo' # 12.5m grid and clipped to 100m elevation or 3000m from coast
+onshore_name_25 = 'hob5_topo_25' # 25m grid and clipped to 100m elevation or 3000m from coast
 offshore_name_tas1 = 'derwent_2m'
 offshore_name_tas2 = 'derwent_5m'
-offshore_name_tas3 = 'south_east_tas'
+offshore_name_tas3 = 'south_east_tas' #actually this is AHO
 offshore_name_tas4 = 'hobart_1m'
 
@@ -85 +87 @@
 
 onshore_dem_name = datadir + onshore_name
+onshore_dem_name_25 = datadir + onshore_name_25
+all_onshore_dem_name = datadir + 'combined_onshore'
 offshore_dem_name_local1 = datadir + offshore_name_tas1
 offshore_dem_name_local2 = datadir + offshore_name_tas2
@@ -110 +114 @@
 outputname = outputtimedir + basename  #Used by post processing
 
-# region to export
+# clipping 12.5m onshore data set
+eastingmin = 520000
+eastingmax = 536000
+northingmin = 5245000
+northingmax = 5260000
+
+# region to export for inundation map
 e_min_area = 500000
 e_max_area = 580000
@@ -132 +142 @@
 polyAll, zone = convert_points_from_latlon_to_utm([d0, d1, d2, d3])
 refzone = zone
+NWcorner = polyAll[0]
+SEcorner = polyAll[2]
+
+# region to export for Alex to make bathymetry map
+e_min_area = 480000#NWcorner[0]
+e_max_area = 600000#SEcorner[0]
+n_min_area = 5156000#NWcorner[1]
+n_max_area = 5285000#SEcorner[1]
 
 #Interior region - Hobart city area + Glenorchy, Kingston