Changeset 6975 for anuga_work/production/australia_ph2
- Timestamp:
- May 5, 2009, 6:51:54 PM (16 years ago)
- Location:
- anuga_work/production/australia_ph2/bamaga
- Files:
-
- 2 deleted
- 3 edited
- 1 copied
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/australia_ph2/bamaga/build_elevation.py
r6895 r6975 59 59 dem2pts(absolute_filename, use_cache=True, verbose=True) 60 60 61 G_grid= Geospatial_data(file_name=absolute_filename+'.pts',61 geospatial_data[filename] = Geospatial_data(file_name=absolute_filename+'.pts', 62 62 verbose=True) 63 print 'Clip geospatial object'64 geospatial_data[filename] = G_grid.clip(project.bounding_polygon)65 63 66 64 # Create Geospatial data from TXT files … … 68 66 for filename in project.point_filenames: 69 67 absolute_filename = join(project.topographies_folder, filename) 70 G_points= Geospatial_data(file_name=absolute_filename,68 geospatial_data[filename] = Geospatial_data(file_name=absolute_filename, 71 69 verbose=True) 72 print 'Clip geospatial object'73 geospatial_data[filename] = G_points.clip(project.bounding_polygon)74 70 75 71 #------------------------------------------------------------------------------- … … 86 82 G += geospatial_data[key] 87 83 88 #print 'Clip combined geospatial data'84 print 'Clip combined geospatial data' 89 85 ##G_clip = G.clip_outside(project.poly_aoi1) 90 86 ##G_all = G_clip + geospatial_data[project.point_filenames[4]] 91 87 #G_clipped = G_all.clip(project.poly_all) 92 #G_clip = G.clip(project.bounding_polygon)88 G_clip = G.clip(project.bounding_polygon) 93 89 94 90 95 91 print 'Export combined DEM file' 96 G .export_points_file(project.combined_elevation + '.pts')92 G_clip.export_points_file(project.combined_elevation + '.pts') 97 93 print 'Do txt version too' 98 94 # Use for comparision in ARC 99 G .export_points_file(project.combined_elevation + '.txt')95 G_clip.export_points_file(project.combined_elevation + '.txt') 100 96 -
anuga_work/production/australia_ph2/bamaga/export_results_max.py
r6939 r6975 22 22 directory = project.output_folder 23 23 24 ##time_dir1 = '20090408_152655_run_final_0_31975_2000_Tb__mhingee'25 ##time_dir2 = '20090409_123239_run_final_0_51390_2000_Tb__mhingee'26 ##time_dir3 = '20090410_110658_run_final_0_63735_2000_Tb__mhingee'24 time_dir1 = '20090408_152655_run_final_0_31975_2000_Tb__mhingee' 25 time_dir2 = '20090409_123239_run_final_0_51390_2000_Tb__mhingee' 26 time_dir3 = '20090410_110658_run_final_0_63735_2000_Tb__mhingee' 27 27 28 time_dir1 = '20090421_135412_run_final_0_31975_2000_Tb_internal_lfountai' 29 30 time_dirs = [time_dir1]#, time_dir2, time_dir3] 28 time_dirs = [time_dir1, time_dir2, time_dir3] 31 29 32 30 cellsize = 250 … … 84 82 85 83 name1 = join(directory, time_dir, project.scenario_name) 86 ## name2 = join(directory, time_dir, project.scenario_name)+'_time_37920_0' 87 name2 = join(directory, time_dir, project.scenario_name)+'_time_37860_0' 88 84 name2 = join(directory, time_dir, project.scenario_name)+'_time_37920_0' 85 89 86 names = [name1, name2] 90 87 -
anuga_work/production/australia_ph2/bamaga/run_model.py
r6911 r6975 87 87 # Number of boundary segments 88 88 num_ocean_segments = len(event_sts) - 1 89 print num_ocean_segments90 89 # Number of landward_boundary points 91 90 num_land_points = file_length(project.landward_boundary) 92 print num_land_points 91 93 92 # Boundary tags refer to project.landward_boundary 94 93 # 4 points equals 5 segments start at N 95 boundary_tags={'back': range(num_ocean_segments+ 2,96 num_ocean_segments+ 5),94 boundary_tags={'back': range(num_ocean_segments+1, 95 num_ocean_segments+num_land_points), 97 96 'side': [num_ocean_segments, 98 num_ocean_segments+1, 99 num_ocean_segments+5, 100 num_ocean_segments+6, 101 num_ocean_segments+7, 102 num_ocean_segments+8, 103 num_ocean_segments+9], 97 num_ocean_segments+num_land_points], 104 98 'ocean': range(num_ocean_segments)} 105 ##boundary_tags = {'back':[3],106 ## 'side':[2,4],107 ## 'ocean':[0,1,5]}108 109 print 'boundary tags', boundary_tags110 99 111 100 # Build mesh and domain
Note: See TracChangeset
for help on using the changeset viewer.