Changeset 5393
- Timestamp:
- Jun 4, 2008, 3:20:10 PM (17 years ago)
- Location:
- anuga_work/development/boxingday08
- Files:
-
- 5 added
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/development/boxingday08/combine_good_data.py
r5247 r5393 83 83 if not exists(dir+sep+filename+'_clipped'+'.csv'): 84 84 clip_csv(dir+sep+filename,project.bounding_polygon) 85 else: 86 print dir+sep+filename+'_clipped'+'.csv already exists ensure the .csv file has not changed' 85 87 86 88 print 'Creating geospatial data' … … 89 91 else: 90 92 G0 = Geospatial_data(file_name=dir+sep+filenames[0]+'.csv') 93 print 'read bay' 91 94 G1=Geospatial_data(file_name=dir+sep+filenames[1]+'_clipped.csv') 95 print 'read andaman' 92 96 G1.clip(project.bounding_polygon) 97 print 'clipped andaman' 93 98 G1.clip_outside(G0) 94 99 -
anuga_work/development/boxingday08/export_results.py
r5247 r5393 5 5 from os import sep 6 6 7 #name='good_simulation/good_simulation_polylinefriction20080529_115703' 8 #name='good_simulation/good_simulation_ursfriction20080529_160110' 9 name='good_simulation/good_simulationsourcefriction20080529_160144' 7 10 8 #time_dir = '20070322_035135' # birubi9 10 #directory = project_slide.outputdir11 #name = directory + time_dir + sep + project_slide.basename12 #name = 'test500000.035000'13 #name = 'test'14 name = 'good_simulationsourcefriction20080404_140016'15 is_parallel = False16 if is_parallel == True: nodes = 417 11 print 'output dir:', name 18 12 … … 38 32 quantityname = 'elevation' #Elevation 39 33 40 if is_parallel == True: 41 for i in range(0,nodes): 42 namei = name + '_P%d_%d' %(i,nodes) 43 outnamei = outname + '_P%d_%d' %(i,nodes) 44 print 'start sww2dem for sww file %d' %(i) 45 sww2dem(namei, basename_out = outnamei, 46 quantity = quantityname, 47 timestep = 1, 48 cellsize = 100, 49 #easting_min = project.e_min_area, 50 #easting_max = project.e_max_area, 51 #northing_min = project.n_min_area, 52 #northing_max = project.n_max_area, 53 reduction = max, 54 verbose = True, 55 format = 'asc') 56 else: 57 print 'start sww2dem' 58 sww2dem(name, basename_out = outname, 59 quantity = quantityname, 60 #timestep = 0, 61 cellsize = 50, 62 #easting_min = 416920, 63 #easting_max = 425310, 64 #northing_min = 877864, 65 #northing_max = 869570, 66 reduction = max, 67 verbose = True, 68 format = 'asc') 34 print 'start sww2dem' 35 sww2dem(name, basename_out = outname, 36 quantity = quantityname, 37 cellsize = 20, 38 easting_min = 417348, 39 easting_max = 425656, 40 northing_min = 870565, 41 northing_max = 877008, 42 reduction = max, 43 verbose = True, 44 format = 'asc') -
anuga_work/development/boxingday08/project.py
r5247 r5393 44 44 45 45 dir='mesh_polygons' 46 #extent = 'extent.csv'47 extent = 'boundary_new_pts.csv' 46 extent = 'extent.csv' 47 #extent = 'boundary_new_pts.csv'#old boundary 48 48 # bounding polygon for study area 49 49 bounding_polygon = read_polygon(dir+sep+extent) … … 75 75 #south = -5.0 76 76 77 #max_extent of ol simulation is why does this work and not below77 #max_extent of old simulation is why does this work and not below 78 78 #7.40911081272 8.71484358635 99.1683687224 97.513856322 79 79 80 80 #max extent of clipped MOST file 81 east = 100 82 west = 96.0 81 82 east = 97.7 83 west = 96.7 83 84 north = 9.0 84 south = 6.085 86 85 south = 5.9 86 87
Note: See TracChangeset
for help on using the changeset viewer.