Changeset 3660

Timestamp:

Sep 25, 2006, 8:41:52 AM (17 years ago)

Author:

sexton

Message:

Hobart updates

Location:

anuga_work/production/hobart_2006

Files:

• project.py (modified) (2 diffs)
• run_hobart.py (modified) (2 diffs)

anuga_work/production/hobart_2006/project.py

@@ -15 +15 @@
 #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
+onshore_name_25 = 'hob5_topo_25m' # 25m grid and clipped to 100m elevation or 3000m from coast
 offshore_name_tas1 = 'derwent_2m'
 offshore_name_tas2 = 'derwent_5m'
@@ -145 +145 @@
 SEcorner = polyAll[2]
 
-#Main Domain of Hobart: updated region from Alex JS 22/9/06
-j0 = [583819.54,5193158.37]
-j1 = [585025.05,5172148.01]
-j2 = [583991.76,5155443.05]
-j3 = [486000.80,5155443.05]
-j4 = [486173.02,5165603.80]
-j5 = [478939.94,5170770.28]
-j6 = [479628.80,5290787.67]
-j7 = [570731.12,5281488.00]
-j8 = [583991.76,5266160.77]
-j9 = [590363.75,5268399.58]
-j10 = [597252.40,5276838.17]
-#j11 = [602418.88,5276838.17]
-j11 = [600000,5276838.17]
-#j12 = [600868.93,5201924.17]
-j12 = [597600, 5201914]
-#j13 = [593119.21,5196241.04]
-j13 = [590000,5196241]
-#j14 = [583819.54,5193158.37]
-j14 = [580000, 5193160]ZZ
-polyAll = [j0,j1,j2,j3,j4,j5,j6,j7,j8,j9,j10,j11,j12,j13,j14]
-
 # region to export for Alex to make bathymetry map
 e_min_area = 480000#NWcorner[0]
 e_max_area = 600000#SEcorner[0]
-n_min_area = 5160000#NWcorner[1]
-n_max_area = 5280000#SEcorner[1]
+n_min_area = 5156000#NWcorner[1]
+n_max_area = 5285000#SEcorner[1]
 
 #Interior region - Hobart city area + Glenorchy, Kingston

anuga_work/production/hobart_2006/run_hobart.py

@@ -68 +68 @@
 copied_files = False
 
-# creates DEM from asc data
+# creates DEM from asc data - 12.5m
 convert_dem_from_ascii2netcdf(onshore_dem_name, use_cache=True, verbose=True)
 
@@ -78 +78 @@
         northing_max= project.northingmax,
         use_cache=True, verbose=True)
+
+# create DEM from asc data - 25m data
+convert_dem_from_ascii2netcdf(onshore_dem_name_25, use_cache=True, verbose=True)
 
 #creates pts file for onshore DEM - 25