Changeset 3683 for anuga_work/production/hobart_2006/run_hobart.py
- Timestamp:
- Oct 3, 2006, 5:28:47 PM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/hobart_2006/run_hobart.py
r3679 r3683 62 62 # filenames 63 63 onshore_offshore_dem_name = project.onshore_offshore_dem_name 64 onshore_offshore_dem_name_25 = project.onshore_offshore_dem_name_25 64 65 meshname = project.meshname+'.msh' 65 66 source_dir = project.boundarydir … … 67 68 copied_files = False 68 69 69 # create DEM from asc data - 25mdata70 # create DEM from 50m asc data 70 71 convert_dem_from_ascii2netcdf(onshore_offshore_dem_name, use_cache=True, verbose=True) 71 72 72 # creates pts file for combined DEM - 2573 # creates pts file for combined 50m DEM 73 74 dem2pts(onshore_offshore_dem_name, use_cache=True, verbose=True) 74 75 76 # 25m data (clipping the around the Hobart area) 77 convert_dem_from_ascii2netcdf(onshore_offshore_dem_name_25, use_cache=True, verbose=True) 78 # creates pts file for 25m data around Hobart 79 dem2pts(onshore_offshore_dem_name_25, project.hobart_dem_name_25, 80 easting_min=project.eastingmin25, 81 easting_max=project.eastingmax25, 82 northing_min=project.northingmin25, 83 northing_max= project.northingmax25, 84 use_cache=True, 85 verbose=True) 86 87 # combining the 50m and Hobart 25m data 88 combine_rectangular_points_files(project.hobart_dem_name_25 + '.pts', 89 project.onshore_offshore_dem_name + '.pts', 90 project.combined_dem_name + '.pts') 91 92 # 25m data (clipping the around site 24 on Bruny Island) 93 convert_dem_from_ascii2netcdf(onshore_offshore_dem_name_25, use_cache=True, verbose=True) 94 # creates pts file for 25m data around site 24 at Bruny Island 95 dem2pts(onshore_offshore_dem_name_25, project.bruny_dem_name_25, 96 easting_min=project.eastingmin25_2, 97 easting_max=project.eastingmax25_2, 98 northing_min=project.northingmin25_2, 99 northing_max= project.northingmax25_2, 100 use_cache=True, 101 verbose=True) 102 103 # combining the 50m and Hobart 25m data with Bruny Island 25m data 104 combine_rectangular_points_files(project.bruny_dem_name_25 + '.pts', 105 project.combined_dem_name + '.pts', 106 project.combined_dem_name_2 + '.pts') 107 75 108 # create geospatial data set and export 76 G = Geospatial_data(file_name = project.onshore_offshore_dem_name + '.pts')77 G.export_points_file(project.combined_dem_name + '.pts')109 #G = Geospatial_data(file_name = project.onshore_offshore_dem_name + '.pts') 110 #G.export_points_file(project.combined_dem_name + '.pts') 78 111 79 112 #---------------------------------------------------------------------------- … … 104 137 'e12': [12], 'e13': [13], 'e14': [14], 105 138 'e15': [15]}, 106 'maximum_triangle_area': 2 50000,107 'filename': meshname, 139 'maximum_triangle_area': 2000000, 140 'filename': meshname, 108 141 'interior_regions': interior_regions}, 109 142 verbose = True, evaluate=False) … … 123 156 domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum']) 124 157 domain.set_minimum_storable_height(0.01) 125 domain.set_store_vertices_uniquely( True) # for writting to sww158 domain.set_store_vertices_uniquely(False) # for writting to sww 126 159 127 160 #------------------------------------------------------------------------------- … … 136 169 domain.set_quantity('elevation', 137 170 # filename = project.onshore_dem_name + '.pts', 138 filename = project.combined_dem_name + '.pts',171 filename = project.combined_dem_name_2 + '.pts', 139 172 # filename = project.offshore_dem_name + '.pts', 140 173 use_cache = True, … … 195 228 domain.set_boundary( {'e0': Bd, 'e1': Bd, 'e2': Bd, 'e3': Bd, 'e4': Bd, 196 229 'e5': Bd, 'e6': Bd, 'e7': Bd, 'e8': Bd, 'e9': Bd, 197 'e10': Bd, 'e11': Bd, 'e12': B d, 'e13': Bd, 'e14': Bf,198 'e15': B d} )230 'e10': Bd, 'e11': Bd, 'e12': Bf, 'e13': Bf, 'e14': Bf, 231 'e15': Bf} ) 199 232 200 233 #-------------------------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.