Changeset 2956
- Timestamp:
- May 24, 2006, 11:54:27 AM (19 years ago)
- Location:
- production/onslow_2006
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
production/onslow_2006/export_results.py
r2955 r2956 23 23 #which_var = int(raw_input('Stage = 0, Absolute Momentum = 1, Depth = 2, Speed = 3 ' )) 24 24 which_var = 2 25 sys.stderr.write(sys.stdout.data)25 #sys.stderr.write(sys.stdout.data) 26 26 if which_var == 0: # Stage 27 27 outname = name + '_stage' … … 47 47 48 48 print 'start sww2dem' 49 sys.stderr.write(sys.stdout.data)49 #sys.stderr.write(sys.stdout.data) 50 50 sww2dem(name, basename_out = outname, 51 51 quantity = quantityname, … … 60 60 format = 'asc') 61 61 62 sys.stderr.write(sys.stdout.data)62 #sys.stderr.write(sys.stdout.data) 63 63 64 64 #Check -
production/onslow_2006/run_onslow.py
r2902 r2956 25 25 from pyvolution.combine_pts import combine_rectangular_points_files 26 26 from pyvolution.pmesh2domain import pmesh_to_domain_instance 27 from geospatial_data import add_points_files27 #from geospatial_data import add_points_files 28 28 29 29 # Application specific imports … … 88 88 # print' most file', project.MOST_dir + project.boundary_basename 89 89 90 90 """ 91 91 # fine data (clipping the points file to smaller area) 92 92 # creates DEM from asc data … … 113 113 print'export G' 114 114 G.export_points_file(project.combined_dem_name + '.pts') 115 115 """ 116 116 117 117 #------------------------------------------------------------------------------- … … 122 122 123 123 from pmesh.mesh_interface import create_mesh_from_regions 124 ''' 125 # original 126 interior_res = 5000 127 high_res = 1500 128 interior_regions = [[project.poly_onslow, high_res], 129 [project.poly_thevenard, interior_res], 130 [project.poly_coast, interior_res]] 131 ''' 124 132 125 #new 133 126 region_res = 50000 … … 135 128 onslow_res = 500 136 129 interior_regions = [[project.poly_onslow, onslow_res], 137 [project.poly_coast, coast_res], 138 [project.poly_region, region_res]] 130 [project.poly_coast, coast_res]] 131 #, 132 # [project.poly_region, region_res]] 139 133 140 134 print 'number of interior regions', len(interior_regions) … … 147 141 'bottom': [4], 'bottomright': [5], 148 142 'topright':[6]}, 149 'maximum_triangle_area': 200000,143 'maximum_triangle_area': 1000000, 150 144 'filename': meshname, 151 145 'interior_regions': interior_regions},
Note: See TracChangeset
for help on using the changeset viewer.