Changeset 3650 for anuga_work/production/hobart_2006/run_hobart.py
- Timestamp:
- Sep 22, 2006, 12:21:09 PM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/hobart_2006/run_hobart.py
r3642 r3650 24 24 Dirichlet_boundary, Time_boundary, File_boundary 25 25 from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts 26 from anuga.abstract_2d_finite_volumes.combine_pts import combine_rectangular_points_files 26 27 from anuga.geospatial_data.geospatial_data import * 27 28 from anuga.abstract_2d_finite_volumes.util import Screen_Catcher … … 61 62 # filenames 62 63 onshore_dem_name = project.onshore_dem_name 64 onshore_dem_name_25 = project.onshore_dem_name_25 63 65 meshname = project.meshname+'.msh' 64 66 source_dir = project.boundarydir … … 69 71 convert_dem_from_ascii2netcdf(onshore_dem_name, use_cache=True, verbose=True) 70 72 71 #creates pts file for onshore DEM 72 dem2pts(onshore_dem_name, use_cache=True, verbose=True) 73 74 print 'local offshore data sets' 73 #creates pts file for onshore DEM - 12.5 74 dem2pts(onshore_dem_name, 75 easting_min=project.eastingmin, 76 easting_max=project.eastingmax, 77 northing_min=project.northingmin, 78 northing_max= project.northingmax, 79 use_cache=True, verbose=True) 80 81 #creates pts file for onshore DEM - 25 82 dem2pts(onshore_dem_name_25, use_cache=True, verbose=True) 83 84 combine_rectangular_points_files(project.onshore_dem_name + '.pts', 85 project.onshore_dem_name_25 + '.pts', 86 project.all_onshore_dem_name + '.pts') 87 #print 'local offshore data sets' 75 88 #G1 = Geospatial_data(file_name = project.offshore_dem_name_local1 + '.xya') 76 89 #G2 = Geospatial_data(file_name = project.offshore_dem_name_local2 + '.xya') … … 78 91 #G4 = Geospatial_data(file_name = project.offshore_dem_name_local4 + '.xya') 79 92 #G11 = G1 + G2 + G3 + G4 80 print 'aho offshore data sets'81 G5 = Geospatial_data(file_name = project.offshore_dem_name_aho1 + '.xya')82 G6 = Geospatial_data(file_name = project.offshore_dem_name_aho2 + '.xya')83 G7 = Geospatial_data(file_name = project.offshore_dem_name_aho3 + '.xya')84 G8 = Geospatial_data(file_name = project.offshore_dem_name_aho4 + '.xya')85 G9 = Geospatial_data(file_name = project.offshore_dem_name_aho5 + '.xya')86 G10 = Geospatial_data(file_name = project.offshore_dem_name_aho6 + '.xya')87 G11 = Geospatial_data(file_name = project.offshore_dem_name_aho7 + '.xya')88 G12 = Geospatial_data(file_name = project.offshore_dem_name_aho8 + '.xya')89 G13 = Geospatial_data(file_name = project.offshore_dem_name_aho9 + '.xya')90 G14 = Geospatial_data(file_name = project.offshore_dem_name_aho10 + '.xya')91 G15 = Geospatial_data(file_name = project.offshore_dem_name_aho11 + '.xya')92 G16 = Geospatial_data(file_name = project.offshore_dem_name_aho12 + '.xya')93 G17 = Geospatial_data(file_name = project.offshore_dem_name_aho13 + '.xya')94 G18 = Geospatial_data(file_name = project.offshore_dem_name_aho14 + '.xya')95 G19 = Geospatial_data(file_name = project.offshore_dem_name_aho15 + '.xya')96 G20 = Geospatial_data(file_name = project.offshore_dem_name_aho16 + '.xya')97 G12 = G5+G6+G7+G8+G9+G10+G11+G12+G13+G14+G15+G16+G17+G18+G19+G2098 print 'local onshore dem'99 G21 = Geospatial_data(file_name = project.onshore_dem_name + '.pts')93 #print 'aho offshore data sets' 94 #G5 = Geospatial_data(file_name = project.offshore_dem_name_aho1 + '.xya') 95 #G6 = Geospatial_data(file_name = project.offshore_dem_name_aho2 + '.xya') 96 #G7 = Geospatial_data(file_name = project.offshore_dem_name_aho3 + '.xya') 97 #G8 = Geospatial_data(file_name = project.offshore_dem_name_aho4 + '.xya') 98 #G9 = Geospatial_data(file_name = project.offshore_dem_name_aho5 + '.xya') 99 #G10 = Geospatial_data(file_name = project.offshore_dem_name_aho6 + '.xya') 100 #G11 = Geospatial_data(file_name = project.offshore_dem_name_aho7 + '.xya') 101 #G12 = Geospatial_data(file_name = project.offshore_dem_name_aho8 + '.xya') 102 #G13 = Geospatial_data(file_name = project.offshore_dem_name_aho9 + '.xya') 103 #G14 = Geospatial_data(file_name = project.offshore_dem_name_aho10 + '.xya') 104 #G15 = Geospatial_data(file_name = project.offshore_dem_name_aho11 + '.xya') 105 #G16 = Geospatial_data(file_name = project.offshore_dem_name_aho12 + '.xya') 106 #G17 = Geospatial_data(file_name = project.offshore_dem_name_aho13 + '.xya') 107 #G18 = Geospatial_data(file_name = project.offshore_dem_name_aho14 + '.xya') 108 #G19 = Geospatial_data(file_name = project.offshore_dem_name_aho15 + '.xya') 109 #G20 = Geospatial_data(file_name = project.offshore_dem_name_aho16 + '.xya') 110 #G12 = G5+G6+G7+G8+G9+G10+G11+G12+G13+G14+G15+G16+G17+G18+G19+G20 111 #print 'local onshore dem' 112 #G21 = Geospatial_data(file_name = project.onshore_dem_name + '.pts') 100 113 print 'coast' 101 G22 = Geospatial_data(file_name = project.coast_dem_name + '.xya')102 print 'adding data sets'114 #G22 = Geospatial_data(file_name = project.coast_dem_name + '.xya') 115 #print 'adding data sets' 103 116 #G = G11 + G12 + G21 + G22 104 G = G21 117 #G = G1+G2+G3+G4+G5+G6+G7+G8+G9+G10+G11+G12+G13+G14+G15+G16+G17+G18+G19+G20+G21+G22 118 #G = G5+G6+G7+G8+G9+G10+G11+G12+G13+G14+G15+G16+G17+G18+G19+G20+G21+G22 119 #G = G21+G22+G5+G6+G7+G8+G9+G10+G11+G12+G13+G14+G15+G16+G17+G18+G19+G20 120 #G = G21+G22+G5 121 #G = Geospatial_data(file_name = project.onshore_dem_name + '.pts') #+ \ 122 # Geospatial_data(file_name = project.coast_dem_name + '.xya') #+ \ 123 # Geospatial_data(file_name = project.offshore_dem_name_aho1 + '.xya') 124 G = Geospatial_data(file_name = project.offshore_dem_name_local1 + '.xya')+\ 125 Geospatial_data(file_name = project.offshore_dem_name_local2 + '.xya')+\ 126 Geospatial_data(file_name = project.offshore_dem_name_local3 + '.xya')+\ 127 Geospatial_data(file_name = project.offshore_dem_name_local4 + '.xya')+\ 128 Geospatial_data(file_name = project.offshore_dem_name_aho1 + '.xya')+\ 129 Geospatial_data(file_name = project.offshore_dem_name_aho2 + '.xya')+\ 130 Geospatial_data(file_name = project.offshore_dem_name_aho3 + '.xya')+\ 131 Geospatial_data(file_name = project.offshore_dem_name_aho4 + '.xya')+\ 132 Geospatial_data(file_name = project.offshore_dem_name_aho5 + '.xya')+\ 133 Geospatial_data(file_name = project.offshore_dem_name_aho6 + '.xya')+\ 134 Geospatial_data(file_name = project.offshore_dem_name_aho7 + '.xya')+\ 135 Geospatial_data(file_name = project.offshore_dem_name_aho8 + '.xya')+\ 136 Geospatial_data(file_name = project.offshore_dem_name_aho9 + '.xya')+\ 137 Geospatial_data(file_name = project.offshore_dem_name_aho10 + '.xya')+\ 138 Geospatial_data(file_name = project.offshore_dem_name_aho11 + '.xya')+\ 139 Geospatial_data(file_name = project.offshore_dem_name_aho12 + '.xya')+\ 140 Geospatial_data(file_name = project.offshore_dem_name_aho13 + '.xya')+\ 141 Geospatial_data(file_name = project.offshore_dem_name_aho14 + '.xya')+\ 142 Geospatial_data(file_name = project.offshore_dem_name_aho15 + '.xya')+\ 143 Geospatial_data(file_name = project.offshore_dem_name_aho16 + '.xya')+\ 144 Geospatial_data(file_name = project.all_onshore_dem_name + '.pts') 105 145 G.export_points_file(project.combined_dem_name + '.pts') 106 146
Note: See TracChangeset
for help on using the changeset viewer.