Changeset 3951
- Timestamp:
- Nov 8, 2006, 6:41:19 PM (18 years ago)
- Location:
- anuga_work/production/broome_2006
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/broome_2006/project.py
r3944 r3951 56 56 coast_name = 'coastline' 57 57 58 # TIN model to fill in data gap 59 bathy_interp = 'interpolate' 60 58 61 boundary_basename = 'SU-AU' # Mw ? 59 62 … … 82 85 # Necessary if using point datasets, rather than grid 83 86 onshore_dem_name = datadir + onshore_name 87 offshore_interp_dem_name = datadir + bathy_interp 84 88 offshore_dem_name1 = datadir + offshore_name1 85 89 offshore_dem_name2 = datadir + offshore_name2 -
anuga_work/production/broome_2006/run_broome.py
r3897 r3951 63 63 # filenames 64 64 onshore_dem_name = project.onshore_dem_name 65 offshore_interp_dem_name = project.offshore_interp_dem_name 65 66 coast_points = project.coast_dem_name 66 67 meshname = project.meshname+'.msh' … … 71 72 #creates pts file for onshore DEM 72 73 dem2pts(onshore_dem_name, use_cache=True, verbose=True) 74 75 # creates DEM from asc data 76 convert_dem_from_ascii2netcdf(offshore_interp_dem_name, use_cache=True, verbose=True) 77 78 #creates pts file for offshore interpolated DEM 79 dem2pts(offshore_interp_dem_name, use_cache=True, verbose=True) 73 80 74 81 print 'create offshore' … … 94 101 Geospatial_data(file_name = project.offshore_dem_name20 + '.xya')+\ 95 102 Geospatial_data(file_name = project.offshore_dem_name21 + '.xya')+\ 96 Geospatial_data(file_name = project.offshore_dem_name22 + '.xya') 103 Geospatial_data(file_name = project.offshore_dem_name22 + '.xya')+\ 104 Geospatial_data(file_name = project.offshore_interp_dem_name + '.pts') 97 105 print 'create onshore' 98 106 G2 = Geospatial_data(file_name = project.onshore_dem_name + '.pts')
Note: See TracChangeset
for help on using the changeset viewer.