Changeset 3683

Timestamp:

Oct 3, 2006, 5:28:47 PM (18 years ago)

Author:

sexton

Message:

updates (clipping 25m data into 50m grid, putting MOST boundary on relevant boundary segments, change saving unique vertices to False)

Location:

anuga_work/production/hobart_2006

Files:

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

anuga_work/production/hobart_2006/project.py

@@ -113 +113 @@
 coast_dem_name = datadir + coast_name
 # addition once total grid delivered
-#onshore_offshore_dem_name = datadir + '25m_se_tas' #25m grid
+onshore_offshore_dem_name_25 = datadir + '25m_se_tas' #25m grid
 onshore_offshore_dem_name = datadir + '50m_se_tas' #50m grid
-
+# output names
+bruny_dem_name_25 = datadir + 'bruny_25_dem'
+hobart_dem_name_25 = datadir + 'hobart_25_dem'
 combined_dem_name = datadir + 'hobart_combined_elevation'
-combined_dem_name2 = datadir + 'hobart_combined_elevation2' # for alpha checking
+combined_dem_name_2 = datadir + 'hobart_combined_elevation_2'
 
 outputname = outputtimedir + basename  #Used by post processing
@@ -167 +169 @@
 northingmin = 5245000
 northingmax = 5265000
-# clipping 25m onshore data set
-eastingmin25 = 520000
-eastingmax25 = 580000
-northingmin25 = 5170000
-northingmax25 = 5260000
+
+# clipping 25m data set - Hobart
+eastingmin25 = 524208.387
+eastingmax25 = 554867.24
+northingmin25 = 5229154.555
+northingmax25 = 5258511.857
+
+# clipping 25m data set - Bruny
+eastingmin25_2 = 523672.502
+eastingmax25_2 = 536020.057
+northingmin25_2 = 5203120.035
+northingmax25_2 = 5212052.309
+
+
 
 ###############################

anuga_work/production/hobart_2006/run_hobart.py

@@ -62 +62 @@
 # filenames
 onshore_offshore_dem_name = project.onshore_offshore_dem_name
+onshore_offshore_dem_name_25 = project.onshore_offshore_dem_name_25
 meshname = project.meshname+'.msh'
 source_dir = project.boundarydir
@@ -67 +68 @@
 copied_files = False
 
-# create DEM from asc data - 25m data
+# create DEM from 50m asc data
 convert_dem_from_ascii2netcdf(onshore_offshore_dem_name, use_cache=True, verbose=True)
 
-#creates pts file for combined DEM - 25
+# creates pts file for combined 50m DEM
 dem2pts(onshore_offshore_dem_name, use_cache=True, verbose=True)
 
+# 25m data (clipping the around the Hobart area)
+convert_dem_from_ascii2netcdf(onshore_offshore_dem_name_25, use_cache=True, verbose=True)
+# creates pts file for 25m data around Hobart
+dem2pts(onshore_offshore_dem_name_25, project.hobart_dem_name_25,
+        easting_min=project.eastingmin25,
+        easting_max=project.eastingmax25,
+        northing_min=project.northingmin25,
+        northing_max= project.northingmax25,
+        use_cache=True, 
+        verbose=True)
+
+# combining the 50m and Hobart 25m data
+combine_rectangular_points_files(project.hobart_dem_name_25 + '.pts',
+                                 project.onshore_offshore_dem_name + '.pts',
+                                 project.combined_dem_name + '.pts')
+
+# 25m data (clipping the around site 24 on Bruny Island)
+convert_dem_from_ascii2netcdf(onshore_offshore_dem_name_25, use_cache=True, verbose=True)
+# creates pts file for 25m data around site 24 at Bruny Island
+dem2pts(onshore_offshore_dem_name_25, project.bruny_dem_name_25,
+        easting_min=project.eastingmin25_2,
+        easting_max=project.eastingmax25_2,
+        northing_min=project.northingmin25_2,
+        northing_max= project.northingmax25_2,
+        use_cache=True, 
+        verbose=True)
+
+# combining the 50m and Hobart 25m data with Bruny Island 25m data
+combine_rectangular_points_files(project.bruny_dem_name_25 + '.pts',
+                                 project.combined_dem_name + '.pts',
+                                 project.combined_dem_name_2 + '.pts')
+
 # create geospatial data set and export 
-G = Geospatial_data(file_name = project.onshore_offshore_dem_name + '.pts')
-G.export_points_file(project.combined_dem_name + '.pts')
+#G = Geospatial_data(file_name = project.onshore_offshore_dem_name + '.pts')
+#G.export_points_file(project.combined_dem_name + '.pts')
 
 #----------------------------------------------------------------------------
@@ -104 +137 @@
                               'e12': [12], 'e13': [13], 'e14': [14],
                               'e15': [15]},
-           'maximum_triangle_area': 250000,
-           'filename': meshname,           
+           'maximum_triangle_area': 2000000,
+           'filename': meshname,
            'interior_regions': interior_regions},
           verbose = True, evaluate=False)
@@ -123 +156 @@
 domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
 domain.set_minimum_storable_height(0.01)
-domain.set_store_vertices_uniquely(True)  # for writting to sww
+domain.set_store_vertices_uniquely(False)  # for writting to sww
 
 #-------------------------------------------------------------------------------                                 
@@ -136 +169 @@
 domain.set_quantity('elevation', 
 #                    filename = project.onshore_dem_name + '.pts',
-                    filename = project.combined_dem_name + '.pts',
+                    filename = project.combined_dem_name_2 + '.pts',
 #                    filename = project.offshore_dem_name + '.pts',
                     use_cache = True,
@@ -195 +228 @@
 domain.set_boundary( {'e0': Bd, 'e1': Bd, 'e2': Bd, 'e3': Bd, 'e4': Bd,
                         'e5': Bd, 'e6': Bd, 'e7': Bd, 'e8': Bd, 'e9': Bd,
-                        'e10': Bd, 'e11': Bd, 'e12': Bd, 'e13': Bd, 'e14': Bf,
-                        'e15': Bd} )
+                        'e10': Bd, 'e11': Bd, 'e12': Bf, 'e13': Bf, 'e14': Bf,
+                        'e15': Bf} )
 
 #-------------------------------------------------------------------------------